Complex floating point numbers in ccsv5

#include <stdio.h>
#include <math.h>
#include <complex.h>
#define SIZE 12

void CircularShift(double complex p[SIZE]);
void Print(double complex p[SIZE]);



/*--------------------------------------------------------------
Difference between float and double: The maximum value of float is only about 3e38, but double is
about 1.7e308, so using float can hit Infinity much easier than
double for something simple e.g. computing 60!.

----------------------------------------------------------------*/

int main()
{

    double tau_rms=991e-9;
    double delta_f = 15e3;                       //pilot spacing
    double Tsym=71.7e-6;
    double fc=2.4e9;
    double c=3e8;
    double v = 10.0;                              //kmp
    double fd=(v/c)*fc;
    double pi = 3.142857143;
    int n,i,j;
    double Nzc=11.0;
    double complex primary_seq[SIZE];
    double complex p[SIZE];
    double complex rs_0[SIZE],rs_1[SIZE],rs_2[SIZE];
    double complex pilot_f[3*SIZE];


printf("\n-----------------------------------------------------------\n");
printf("\n Primary sequence is:");

    for (n =1; n <=SIZE; n++)
    {
       if (n<=SIZE-1)
       {
            primary_seq[n] = cos(pi*n*(n-1)/Nzc) - I*sin(pi*n*(n-1)/Nzc);
            printf("\n primary_seq[%d]: %.15g %+.15gi\n",n,creal(primary_seq[n]),cimag(primary_seq[n]));
       }

       else
       {
            primary_seq[SIZE] = primary_seq[1];
            printf("\n primary_seq[%d]: %.15g %+.15gi\n",n,creal(primary_seq[n]),cimag(primary_seq[n]));
       }


    }

    //Print(primary_seq); - Cannot use this because array starts from 1
printf("\n-----------------------------------------------------------\n");
printf("\n p array is:");


    for ( n=0; n<=SIZE-1; n++)
    {
        p[n] = primary_seq[n+1];
    }

    Print(p);



//Saving the Derived sequence in a temp array

printf("\n-----------------------------------------------------------\n");
printf("\nGeneration of reference signals (frequency domain)\n");
printf("\n rs_0 array:");
for ( n=0; n<=SIZE-1; n++)
    {
        rs_0[n] = p[n];
    }


    for (j=0;j<=1;j++)
    {
        CircularShift(rs_0);
    }

    Print(rs_0);

printf("\n-----------------------------------------------------------\n");
printf("\n rs_1 array:");
for ( n=0; n<=SIZE-1; n++)
    {
        rs_1[n] = p[n];

    }


    for (j=0;j<=2;j++)
    {
        CircularShift(rs_1);
    }

    Print(rs_1);

    printf("\n-----------------------------------------------------------\n");
    printf("\n rs_2 array:");
    for ( n=0; n<=SIZE-1; n++)
        {
            rs_2[n] = p[n];
        }


    for (j=0;j<=3;j++)
        {
            CircularShift(rs_2);
        }

    Print(rs_2);


    getchar();

    return 0;

}

void CircularShift(double complex p[SIZE])
{
    double complex temp;
    temp=p[SIZE-1];//put the last element in the temp variable
    int i;
    for(i=SIZE-1;i>0;i--)//shift each value to the next value so that a hole can be created in the beginning
    p[i]=p[i-1];//shift the values,index 0 is not changed because faida nahi hai
    p[0]=temp;
}

void Print(double complex p[SIZE])
{
    int i;
    for(i=0;i<SIZE;i++)
    printf("\n p[%d]: %.15g %+.15gi\n",i,creal(p[i]),cimag(p[i]));
}