Gennady,

Thnax for th equick response.
setting SSE4.2 worked,

Now I could run more tests, and now the next example crashes for DFTI_COMPLEX_COMPLEX (not for DFTI_COMPLEX_REAL
(crash happens typically at nrOfTransforms 3, nrOfSamples 2658):

for (unsigned nrOfTransforms = 1; nrOfTransforms <= 5; ++nrOfTransforms)
{
for (unsigned nrOfSamples = 1; nrOfSamples <= 10000; ++nrOfSamples)
{
std::cout << "Test 3c, Forward FFT Real-2-complex out-of-place nrOfTransforms " << nrOfTransforms << ", nrOfSamples " << nrOfSamples << std::endl;

MKL_LONG status;
DFTI_DESCRIPTOR_HANDLE _fft;

// allocate buffer (make buffer too big, just to be sure that inplace FFT does not go beyond allocate memory
float *x_in = new float [nrOfSamples*nrOfTransforms*10];
std::complex *x_out = new std::complex[nrOfSamples*nrOfTransforms*10];

status = DftiCreateDescriptor( &_fft, DFTI_SINGLE, DFTI_REAL, 1, nrOfSamples);
checkStatus(status);

status = DftiSetValue(_fft, DFTI_PLACEMENT, DFTI_NOT_INPLACE);
checkStatus(status);

// Specify the number of transforms
status = DftiSetValue(_fft, DFTI_NUMBER_OF_TRANSFORMS, nrOfTransforms);
checkStatus(status);

//status = DftiSetValue(_fft, DFTI_CONJUGATE_EVEN_STORAGE, DFTI_COMPLEX_REAL);
status = DftiSetValue(_fft, DFTI_CONJUGATE_EVEN_STORAGE, DFTI_COMPLEX_COMPLEX);
checkStatus(status);

// The FFT is now fully specified
status = DftiCommitDescriptor( _fft );

// Calculate forward FFT
status = DftiComputeForward(_fft, x_in, x_out);
checkStatus(status);

// cleanup
delete[] x_in, x_out;
status = DftiFreeDescriptor(&_fft);
checkStatus(status);
}
}

Dima,

you are correct that one should specify the input/output distance,

non-the-less the example code still crashes at the same position...

Dirk-Jan

Dirk-Jan