CUDA Replaces FFT Part in Library

 unsigned character imageIN [number of pixels] [number of pixels]

for(i=0;i<width;i++){
  for(j=0;j<width;j++){
    data[i][j]=(float)imageIN[i][j];
    jdata[i][j] = 0.0;
  }
}

int main(intargc, char**argv){
  inti,k,d;
  int n = 512; // Size of one side of the image
  float*v = NULL;
  cufftComplex*dv;
  cufftHandle plan;

  cudaMalloc(void**)&dv,n*n*sizeof(float2));// Securing device memory
  cudaMallocHost(void**)&v,n*n*sizeof(float2)); // Host memory availability
  for(i=0;i<n*n*2;i++)v[i]=(float)(i%(n-1)); // Initialize data

  load_image_data(); // Load original image
  make_original_data(); // Create original data

  // Do you want to put the image data in v here?
  cudaMemcpy(dv,v,n*n*sizeof(double2),cudaMemcpyHostToDevice); // Host to device transfer
  Preparing for 2D FFT for cufftPlan2d(&plan,n,n,CUFFT_R2C);//n×n
  cufftExecC2C (plan, dv, dv, CUFFT_FORWARD); // Forward Conversion

  cudaMemcpy(data,dv,n*n*sizeof(double2),cudaMemcpyDeviceToHost);
                                  // Transfer from device to host
  // cuftDestroy(plan); // Free up resources
  // (I want blur to be done by CPU)
  Blur(); // Blur filtering in frequency domain

  /*
  cudaMemcpy(dv,v,n*n*sizeof(double2),cudaMemcpyHostToDevice); // Host to device transfer
  Preparing for 2D FFT for cufftPlan2d(&plan,n,n,CUFFT_C2R);//n×n
  cufftExecC2C(plan, dv, dv, CUFFT_INVERSE); // Reverse Conversion
  cuffTestroy(plan);// Free up resources
  cudaMemcpy(v,dv,n*n*sizeof(double2),cudaMemcpyDeviceToHost);
                          // Transfer from device to host
  */

  cudaFree(v);
  save_image_data(); // Save output image
  return 0;
}