for (x=0L; x < (long) fourier_info->center; x++)
{
magnitude[i]=cabs(fourier[i]);
phase[i]=carg(fourier[i]);
phase[i]=atan2(cimag(fourier[i]),creal(fourier[i]));
i++;
}
else