Commit d63dce1e authored by Claudio Gheller's avatar Claudio Gheller
Browse files

...

parent 6c831c1c

phase_correction.cu

+8 −8
Original line number Diff line number Diff line
@@ -109,18 +109,18 @@ void phase_correction(double* gridss, double* image_real, double* image_imag, in 
	double * gridss_g;



        mmm=cudaMalloc(&gridss_g, 2*num_w_planes*xaxis*yaxis*sizeof(double));

	printf("CUDA ERROR 1 %s\n",cudaGetErrorString(mmm));

	//printf("CUDA ERROR 1 %s\n",cudaGetErrorString(mmm));

	mmm=cudaMalloc(&image_real_g, xaxis*yaxis*sizeof(double));

	printf("CUDA ERROR 2 %s\n",cudaGetErrorString(mmm));

	//printf("CUDA ERROR 2 %s\n",cudaGetErrorString(mmm));

	mmm=cudaMalloc(&image_imag_g, xaxis*yaxis*sizeof(double));

	printf("CUDA ERROR 3 %s\n",cudaGetErrorString(mmm));

	//printf("CUDA ERROR 3 %s\n",cudaGetErrorString(mmm));



	mmm=cudaMemcpy(gridss_g, gridss, 2*num_w_planes*xaxis*yaxis*sizeof(double), cudaMemcpyHostToDevice);

	printf("CUDA ERROR 4 %s\n",cudaGetErrorString(mmm));

	//printf("CUDA ERROR 4 %s\n",cudaGetErrorString(mmm));

	mmm=cudaMemset(image_real_g, 0.0, xaxis*yaxis*sizeof(double));

	printf("CUDA ERROR 5 %s\n",cudaGetErrorString(mmm));

	//printf("CUDA ERROR 5 %s\n",cudaGetErrorString(mmm));

	mmm=cudaMemset(image_imag_g, 0.0, xaxis*yaxis*sizeof(double));

	printf("CUDA ERROR 6 %s\n",cudaGetErrorString(mmm));

	//printf("CUDA ERROR 6 %s\n",cudaGetErrorString(mmm));



	// call the phase correction kernel

	phase_g <<<Nbl,Nth>>> (xaxis,
@@ -138,9 +138,9 @@ void phase_correction(double* gridss, double* image_real, double* image_imag, in 
			       nbucket);



	mmm = cudaMemcpy(image_real, image_real_g, xaxis*yaxis*sizeof(double), cudaMemcpyDeviceToHost);

	printf("CUDA ERROR 7 %s\n",cudaGetErrorString(mmm));

	//printf("CUDA ERROR 7 %s\n",cudaGetErrorString(mmm));

	mmm = cudaMemcpy(image_imag, image_imag_g, xaxis*yaxis*sizeof(double), cudaMemcpyDeviceToHost);

	printf("CUDA ERROR 8 %s\n",cudaGetErrorString(mmm));

	//printf("CUDA ERROR 8 %s\n",cudaGetErrorString(mmm));



#else