Changeset 41159 for trunk/Ohana/src/opihi/mana/deimos_arclines.c

Timestamp:

Nov 27, 2019, 11:04:51 AM (7 years ago)

Author:

eugene

Message:

improvments to deimos extraction tools

File:

• trunk/Ohana/src/opihi/mana/deimos_arclines.c (modified) (2 diffs)

trunk/Ohana/src/opihi/mana/deimos_arclines.c

@@ -53 +53 @@
   deimos_make_LSF_kernel (LSF, stilt, Nx);
 
+  // FILE *Fout = NULL;
+
   // loop over the rows
   for (int iy = 0; iy < Ny; iy++) {
     float Fsum = 0;
+    float Ksum = 0;
+
+    // if (iy == 5802) Fout = fopen ("test.5802.dat", "w");
+    // if (iy == 5815) Fout = fopen ("test.5815.dat", "w");
+
     // cross-correlation of buffer values and kernel values
     for (int ix = 0; ix < Nx; ix++) {
@@ -62 +69 @@
         if ((iy + ky) < 0) continue;
         if ((iy + ky) >= Ny) continue;
+        if (!isfinite(Fin[ix + (iy + ky)*Nx])) continue;
+        if (!isfinite(kernel[ix + ko*Nx])) continue;
         Fsum += Fin[ix + (iy + ky)*Nx] * kernel[ix + ko*Nx];
+        Ksum += kernel[ix + ko*Nx];
+        // if ((iy == 5815) || (iy == 5802)) {
+        //   fprintf (Fout, "%d %d %d : %f %f\n", iy, ix, ky, Fin[ix + (iy + ky)*Nx], kernel[ix + ko*Nx]);
+        // }
       }
     }
+
+    // if ((iy == 5815) || (iy == 5802)) fclose (Fout);
+
     coord->elements.Flt[iy] = iy;
     flux->elements.Flt[iy] = Fsum;