Comment faire symétrique inverse fft avec CUFFT

Question

Je dois convertir cette ligne (Matlab) à CUDA:

picTimeFiltered = ifft((picFFT_filt), size(I,3), 3 ,'symmetric'); 

Mon implémentation actuelle est de cette ligne (sans drapeau 'symétrique'):

picTimeFiltered = ifft((picFFT_filt), size(I,3), 3); 

Ceci est ma mise en œuvre de CUDA:

void ifftDouble_many(cuDoubleComplex*& input, cuDoubleComplex*& outputMatrixAfterIFFT, const int width, const int height, const int depth) 
{ 
    const int NX = depth; 
    const int NY = width * height; 

    // Allocate and set the CUDA Input 
    cuDoubleComplex *d_input; 
    cudaMalloc(&d_input, NX*NY * sizeof(cuDoubleComplex)); 
    cudaMemcpy(d_input, input, NX * NY * sizeof(cuDoubleComplex), cudaMemcpyHostToDevice); 

    // Allocate the CUDA output 
    cufftDoubleComplex* d_output = nullptr; 
    cudaMalloc((void**)&d_output, sizeof(cuDoubleComplex)*NX*NY); 

    // CUDA FFT 
    cufftHandle plan; 
    int n[1] = { NX }; 
    int inembed[] = { NY, NX }; 
    int onembed[] = { NY, NX }; 
    cufftPlanMany(&plan, 1, n, inembed, 1, NX, onembed, 1, NX, CUFFT_Z2Z, NY); 
    cufftExecZ2Z(plan, d_input, d_output, CUFFT_INVERSE); 

    // Devide the results by depth 
    devideCufftDoubleComplexArrayByScalar_CUDA(d_output, NX * NY, depth); 

    cudaMemcpy(outputMatrixAfterIFFT, d_output, NY*NX * sizeof(cuDoubleComplex), cudaMemcpyDeviceToHost); 

    /* Destroy the CUFFT plan. */ 
    cufftDestroy(plan); 
    cudaFree(d_input); 
    cudaFree(d_output); 
} 

S'il vous plaît conseiller - comment faire symétrique inverse fft via CUDA?

Answer 1

Utilisez cufftExecC2R() ou cufftExecZ2D() pour calculer l'inverse symétrique fft pour la précision simple/double.

http://docs.nvidia.com/cuda/cufft/index.html#axzz4l0egR62U