Sorry.... for error in posting the code
Code is :
Code:
#define re(c) ((c)[0]) ///< real part of a complex number as defined by fftw_complex in fftw
#define im(c) ((c)[1])
#include <fftw3.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
unsigned int fs = 16367600;
double fif = 4.1304E6;
double range = fif;
double span = 10E3;
double fd_step = 500;
double Tint = 1;
unsigned int l,C1,D1;
fftw_complex *SUM_SS;
int main()
{
C1 = ceil((double) fs * 1E-3);
D1 = ceil(span / fd_step) + 1;
SUM_SS = (fftw_complex *) fftw_malloc(sizeof(fftw_complex) * C1 * D1);
if (SUM_SS == NULL) {
fprintf(stderr, "Error allocationg memory fftw operations \n");
}
for (l = 0; l < C1 *D1; l++) {
re(SUM_SS[l]) = 0.0;
im(SUM_SS[l]) = 0.0;
}
}