DftiComputeForward vs dzfft1dc

[melbinmathew1]

Hi,

I have the old code of MKL which was using the dzfft1dc for forward transform and zdfft1dc for backward transform. Now i want to use the new MKL having the DftiComputeForward and DftiComputeBackward as the previous MKL library was not supporting 64-bit. I tried doing it but the ouput is having a huge varition from the older version.

This is the usage in the old code for forward transform:
double *x, *wsave;

const long n=1<<11;
x = new double[n+2];
wsave = new double[2*n+4];
ZeroMemory(x, (n+2)*sizeof(double));
CopyMemory(x, vIn,vIn.size()*sizeof(double));
dzfft1dc(x, n, 0, wsave);
dzfft1dc(x, n, 1, wsave);

This is the usage in the old code forBackward transform:

zdfft1dc(x, n, 0, wsave);
zdfft1dc(x, n, 1, wsave);

With the new code i tried like this for Forward transform:
----------------------------------------------------------------

DFTI_DESCRIPTOR_HANDLE handle;

long ret = DftiCreateDescriptor(&handle,

(DFTI_CONFIG_VALUE)DFTI_DOUBLE , (DFTI_CONFIG_VALUE)DFTI_REAL , 1, length);

long setval1 = DftiSetValue(handle, (DFTI_CONFIG_PARAM)DFTI_PLACEMENT, (DFTI_CONFIG_VALUE)DFTI_INPLACE);

long setval2 = DftiSetValue(handle, (DFTI_CONFIG_PARAM)DFTI_PACKED_FORMAT, (DFTI_CONFIG_VALUE)DFTI_PACK_FORMAT);

long comval = DftiCommitDescriptor(handle);

long retfor = DftiComputeForward(handle,vIn);

long retfree = DftiFreeDescriptor(&handle);

With the new code i tried like this for backward transform:
-------------------------------------------------------------------

DFTI_DESCRIPTOR_HANDLE handlebak;

long retbak = DftiCreateDescriptor(&handlebak,

(DFTI_CONFIG_VALUE)DFTI_DOUBLE , (DFTI_CONFIG_VALUE)DFTI_REAL , 1, length);

long transform_size = length;

double scale_factor = 1.0 / transform_size;

int setval = DftiSetValue(handlebak, (DFTI_CONFIG_PARAM)DFTI_BACKWARD_SCALE, scale_factor);

long setvalbak1 = DftiSetValue(handlebak, (DFTI_CONFIG_PARAM)DFTI_PLACEMENT, (DFTI_CONFIG_VALUE)DFTI_NOT_INPLACE);

long setvalbak2 = DftiSetValue(handlebak, (DFTI_CONFIG_PARAM)DFTI_PACKED_FORMAT, (DFTI_CONFIG_VALUE)DFTI_PACK_FORMAT);

long comvalbak = DftiCommitDescriptor(handlebak);

long retbakr = DftiComputeBackward(handlebak,vIn,vOut);

long retfreebak = DftiFreeDescriptor(&handlebak);


What i wanted to know is the setting which i doin the old code and the new code are same or not?
Like the settings which i do using the DftiSetValue is proper or not when we compare the old code and the new code.

Thanks,
Melbin

[Evgueni_P_Intel]

Hi melbinmathew1,

The issue that you describe may be caused by an inconsistent link line -- 32-bit integers at compile time but 64-bit integers at run time.
Did you check Intel MKL Link Line Advisor?

[melbinmathew1]

Hi,

Thanks for the reply.

I tried compliling and running both versions n 32-bit itself and compared the values just to see if what i am doing is correct or not.

What i wanted to know is in the older verion(dzfft1dc) we had to give the array size as n+2. Is it the same in the new version also?

Is there any site where it gives the actual comparison between the old and the new ones?

Thanks & Regards,
Melbin

[Evgueni_P_Intel]

Hi melbinmathew1,
Yes,in-place real-to-complex FFT of length n=2^11 requires a realarray of length n+2.
No, we don't provide any documentthat compares DFTI and FFT API in MKL 8.0 and beforeexcept the fft2dfti wrappersin MKL 9.0.

BTW, dzfft1dc and DftiComputeForward for real-to-complex in-place FFT store outputs differently.
dzfft1dc stores real parts in elements 0, 1, ..., n/2 and imaginary partsin elements n/2+1, n/2+2, ..., 2*(n/2)+1.
DftiComputeForward stores real parts in elements 0, 2, 4, ..., 2*(n/2) and imaginary parts in elements 1, 3, 5, ..., 2*(n/2)+1.

Evgueni