Why does the mkl fft routine return unreasonably large number?

Why does the mkl fft routine return unreasonably large number?

The code is like this

Status = DftiCreateDescriptor(FFT_Handle32, DFTI_DOUBLE, DFTI_REAL, 1, k32+1)

Status = DftiCommitDescriptor(FFT_Handle32)

Do

Status = DftiComputeBackward(FFT_Handle32, temptv32)

End do

I have initialized the FFT_Handle32, and run the loop for many times for different arrays temptv32. One of the arrays has exactly zero, but the backward FFT routine returns numbers like 6.187317420301585E+195. Why is that?

Thank you in advance.

3 posts / 0 new
Last post
For more complete information about compiler optimizations, see our Optimization Notice.

Hi,

Please send us small testcase to analyze your problem together with information: MKL-version, architecture, OS, etc.

Thanks,
-- Victor

Hellohadesmajesty,how about the test case for reproducing the problem? The existing cases don't show the problem on our side.--Gennady

Leave a Comment

Please sign in to add a comment. Not a member? Join today