TEST: testing degridding function in development

miguelcarcamov · Nov 12, 2024 · 9992b09 · 9992b09
1 parent f125906
commit 9992b09
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 9 deletions.
diff --git a/src/functions.cu b/src/functions.cu
@@ -1993,13 +1993,13 @@ __host__ void degridding(std::vector<Field>& fields,
               blockSizeV);
         }
 
-        hermitianSymmetry<<<
+        /*hermitianSymmetry<<<
             fields[f].device_visibilities[i][s].numBlocksUV,
             fields[f].device_visibilities[i][s].threadsPerBlockUV>>>(
             fields[f].device_visibilities[i][s].uvw,
             fields[f].device_visibilities[i][s].Vo, fields[f].nu[i],
             fields[f].numVisibilitiesPerFreqPerStoke[i][s]);
-        checkCudaErrors(cudaDeviceSynchronize());
+        checkCudaErrors(cudaDeviceSynchronize());*/
 
         // Interpolation / Degridding
         vis_mod2<<<fields[f].device_visibilities[i][s].numBlocksUV,

diff --git a/src/mfs.cu b/src/mfs.cu
@@ -1104,17 +1104,16 @@ void MFS::writeResiduals() {
         "Measurement Set File\n");
     // In the de-gridding procedure weights are also restored to the original
     for (int d = 0; d < nMeasurementSets; d++) {
-      // degridding(datasets[d].fields, datasets[d].data, deltau, deltav,
-      // num_gpus,
-      //            firstgpu, variables.blockSizeV, M, N, this->ckernel);
-      getOriginalVisibilitiesBack(datasets[d].fields, datasets[d].data,
-                                  num_gpus, firstgpu, variables.blockSizeV);
+      degridding(datasets[d].fields, datasets[d].data, deltau, deltav, num_gpus,
+                 firstgpu, variables.blockSizeV, M, N, this->ckernel);
+      /*getOriginalVisibilitiesBack(datasets[d].fields, datasets[d].data,
+                                  num_gpus, firstgpu, variables.blockSizeV);*/
     }
-    Fi* chi2 = optimizer->getObjectiveFunction()->getFiByName("Chi2");
+    /*Fi* chi2 = optimizer->getObjectiveFunction()->getFiByName("Chi2");
     float res = chi2->calcFi(image->getImage());
     printf(
         "Non-gridded chi2 after de-gridding using bilinear interpolation %f\n",
-        res);
+        res);*/
   }
 
   for (int d = 0; d < nMeasurementSets; d++)