From 15def8846646ee1dd4d2f96b6f65efe510412d07 Mon Sep 17 00:00:00 2001
From: nandhanas <s.nandhana@gmail.com>
Date: Thu, 3 Mar 2022 14:35:23 +0100
Subject: [PATCH] Changes in parameters

---
 fourier_transform.c |  4 ++--
 gridding.c          | 23 ++++++++++++-----------
 2 files changed, 14 insertions(+), 13 deletions(-)

diff --git a/fourier_transform.c b/fourier_transform.c
index ba8c652..7648e78 100644
--- a/fourier_transform.c
+++ b/fourier_transform.c
@@ -104,7 +104,7 @@ void write_fftw_data(){
                 for (int iv=0; iv<yaxis; iv++)
                 for (int iu=0; iu<xaxis; iu++)
                 {
-                          long global_index = (iu + (iv+isector*yaxis)*xaxis + iw*grid_size_x*grid_size_y)*sizeof(double);
+                          long global_index = (iu + (iv+isector*yaxis)*xaxis + iw*param.grid_size_x*param.grid_size_y)*sizeof(double);
                           long index = iu + iv*xaxis + iw*xaxis*yaxis;
                           fseek(file.pFilereal, global_index, SEEK_SET);
                           fwrite(&gridss_real[index], 1, sizeof(double), file.pFilereal);
@@ -113,7 +113,7 @@ void write_fftw_data(){
                 for (int iv=0; iv<yaxis; iv++)
                 for (int iu=0; iu<xaxis; iu++)
                 {
-                          long global_index = (iu + (iv+isector*yaxis)*xaxis + iw*grid_size_x*grid_size_y)*sizeof(double);
+                          long global_index = (iu + (iv+isector*yaxis)*xaxis + iw*param.grid_size_x*param.grid_size_y)*sizeof(double);
                           long index = iu + iv*xaxis + iw*xaxis*yaxis;
                           fseek(file.pFileimg, global_index, SEEK_SET);
                           fwrite(&gridss_img[index], 1, sizeof(double), file.pFileimg);
diff --git a/gridding.c b/gridding.c
index 52f7811..f9e7ba7 100644
--- a/gridding.c
+++ b/gridding.c
@@ -333,22 +333,22 @@ void write_grided_data()
                       gridss_real[i] = gridss[2*i];
                       gridss_img[i] = gridss[2*i+1];
               }
-              if (num_w_planes > 1)
+              if (param.num_w_planes > 1)
               {
-                      for (int iw=0; iw<num_w_planes; iw++)
+                      for (int iw=0; iw<param.num_w_planes; iw++)
                         for (int iv=0; iv<yaxis; iv++)
                           for (int iu=0; iu<xaxis; iu++)
                           {
-                               long global_index = (iu + (iv+isector*yaxis)*xaxis + iw*grid_size_x*grid_size_y)*sizeof(double);
+                               long global_index = (iu + (iv+isector*yaxis)*xaxis + iw*param.grid_size_x*param.grid_size_y)*sizeof(double);
                                long index = iu + iv*xaxis + iw*xaxis*yaxis;
                                fseek(file.pFilereal, global_index, SEEK_SET);
                                fwrite(&gridss_real[index], 1, sizeof(double), file.pFilereal);
                           }
-                      for (int iw=0; iw<num_w_planes; iw++)
+                      for (int iw=0; iw<param.num_w_planes; iw++)
                         for (int iv=0; iv<yaxis; iv++)
                           for (int iu=0; iu<xaxis; iu++)
                           {
-                               long global_index = (iu + (iv+isector*yaxis)*xaxis + iw*grid_size_x*grid_size_y)*sizeof(double);
+                               long global_index = (iu + (iv+isector*yaxis)*xaxis + iw*param.grid_size_x*param.grid_size_y)*sizeof(double);
                                long index = iu + iv*xaxis + iw*xaxis*yaxis;
                                fseek(file.pFileimg, global_index, SEEK_SET);
                                fwrite(&gridss_img[index], 1, sizeof(double), file.pFileimg);
@@ -359,9 +359,9 @@ void write_grided_data()
               }
               else
               {
-                      for (int iw=0; iw<num_w_planes; iw++)
+                      for (int iw=0; iw<param.num_w_planes; iw++)
                       {
-                          long global_index = (xaxis*isector*yaxis + iw*grid_size_x*grid_size_y)*sizeof(double);
+                          long global_index = (xaxis*isector*yaxis + iw*param.grid_size_x*param.grid_size_y)*sizeof(double);
                           long index = iw*xaxis*yaxis;
                           fseek(file.pFilereal, global_index, SEEK_SET);
                           fwrite(&gridss_real[index], xaxis*yaxis, sizeof(double), file.pFilereal);
@@ -371,11 +371,11 @@ void write_grided_data()
               }
           }
        #else
-          for (int iw=0; iw<num_w_planes; iw++)
-             for (int iv=0; iv<grid_size_y; iv++)
-               for (int iu=0; iu<grid_size_x; iu++)
+          for (int iw=0; iw<param.num_w_planes; iw++)
+             for (int iv=0; iv<param.grid_size_y; iv++)
+               for (int iu=0; iu<param.grid_size_x; iu++)
                 {
-                      long index = 2*(iu + iv*grid_size_x + iw*grid_size_x*grid_size_y);
+                      long index = 2*(iu + iv*param.grid_size_x + iw*param.grid_size_x*param.grid_size_y);
                       fwrite(&gridtot[index], 1, sizeof(double), file.pFilereal);
                       fwrite(&gridtot[index+1], 1, sizeof(double), file.pFileimg);
                       //double v_norm = sqrt(gridtot[index]*gridtot[index]+gridtot[index+1]*gridtot[index+1]);
@@ -389,6 +389,7 @@ void write_grided_data()
      #ifdef USE_MPI
         MPI_Win_fence(0,slabwin);
      #endif
+
    #endif //WRITE_DATA 
 
 }
-- 
GitLab