Quantcast
Channel: Intel® Software - Intel® oneAPI Math Kernel Library & Intel® Math Kernel Library
Viewing all articles
Browse latest Browse all 3005

FFT not working for N=2048

$
0
0

Hello,

I have a problem with mkl FFT. It stops working without a hint for N=2048.

Here is my implementation:

	MKL_LONG status;
	int nThreads = omp_get_num_threads();
	int th=0;

	DFTI_DESCRIPTOR_HANDLE handle = 0;
	status=DftiCreateDescriptor(&handle, DFTI_SINGLE, DFTI_COMPLEX,1,(MKL_LONG) nChannels);
	status=DftiSetValue(handle, DFTI_PLACEMENT, DFTI_INPLACE);
	status=DftiSetValue (handle, DFTI_NUMBER_OF_USER_THREADS, nThreads);
	status=DftiCommitDescriptor(handle);
	#pragma omp parallel for shared(nChannels,nThreads, input) private(th)
	for (th = 0; th < nSpectra; th++) {
		DftiComputeForward(handle, &input[th*nChannels]);
	}
	DftiFreeDescriptor(&handle);

where N=nChannels. It works fine for N=1024 as well as for N=1536. All statuses are zero until  DftiComputeForward which stops and will not end.

Thanks.

 


Viewing all articles
Browse latest Browse all 3005

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>