Commit 733cebcf authored by Claudio Gheller's avatar Claudio Gheller
Browse files

few bug fixed

parent b854572e

inverse-imaging.c

+17 −9
Original line number Diff line number Diff line
@@ -84,7 +84,7 @@ int main(int argc, char * argv[]) 
        // Image related files

	char imagepath[900] = "./";

        #ifdef FITSIO

	char imagename[FILENAMELENGTH] = "image_name.fits";

	char imagename[FILENAMELENGTH] = "test.fits";

        #else

	char imagename[FILENAMELENGTH] = "image_name.bin";

        #endif
@@ -113,8 +113,8 @@ int main(int argc, char * argv[]) 
	double resolution;



        // Mesh related parameters: global size

	int grid_size_x = 2048;

	int grid_size_y = 2048;

	int grid_size_x = 1024;

	int grid_size_y = 1024;

	// Split Mesh size (auto-calculated)

	int local_grid_size_x;

	int local_grid_size_y;
@@ -215,7 +215,7 @@ if(rank == 0){ 


	// INPUT FILES (only the first ndatasets entries are used)

	int ndatasets = 1;

        strcpy(datapath_multi[0],"/m100_scratch/userexternal/cgheller/gridding/Lofar/L798046_SB244_uv.uncorr_130B27932t_146MHz.pre-cal.binMS/");

        strcpy(datapath_multi[0],"/m100_scratch/userexternal/cgheller/inverse_imaging/L798046_SB244_uv.uncorr_130B27932t_146MHz.pre-cal.binMS/");



	strcpy(datapath,datapath_multi[0]);

	// Read metadata
@@ -363,8 +363,10 @@ if(rank == 0){ 
	int naxis;

	long * naxes;

	fits_get_img_dim(fptr, &naxis,  &status);

	printf("number of axis: %d\n",naxis);

	naxes = (long *)malloc(sizeof(long)*naxis);

	fits_get_img_size(fptr, naxis, naxes, &status);

	printf("image size (from FITS file): %d\n",naxes[0]);

	if (naxes[0] != xaxis){

		printf("Wrong Image Size : %d vs %d\n",naxes[0],xaxis);

		exit(2);
@@ -378,11 +380,11 @@ if(rank == 0){ 
        fpixel[0] = 1;

	lpixel[0] = xaxis;

        inc[0] = 1;

	fpixel[1] = 1;

	lpixel[1] = rank*yaxis+1;

	fpixel[1] = rank*yaxis+1;

	lpixel[1] = (rank+1)*yaxis+1;

        inc[1] = 1;

	

	fits_read_subset(fptr, TFLOAT, fpixel, lpixel, inc, NULL, image_real, &anynul, &status);

	fits_read_subset(fptr, TDOUBLE, fpixel, lpixel, inc, NULL, image_real, &anynul, &status);

	//fits_read_img(fptr, TDOUBLE, 1, xaxis*yaxis, NULL, image_real, &anynul, &status);

	fits_close_file(fptr, &status);



	#else
@@ -394,7 +396,12 @@ if(rank == 0){ 
        fclose(pFilereal);

        #endif

	// image read

	// We read binary images, however we can easily extend to fits files

	

        #ifdef WRITE_DATA

        pFilereal = fopen ("revtest.bin","wb");

	fwrite(image_real, xaxis*yaxis, sizeof(double), pFilereal);

	fclose(pFilereal);

        #endif

	

        if(rank == 0)printf("FFT TRANSFORMING (from Real to Complex Fourier space)\n");
@@ -521,6 +528,7 @@ if(rank == 0){ 
        #endif



#endif //WRITE_DATA

	exit(3);



        if(rank == 0)printf("CREATING LINKED LISTS\n");