Changeset 40838 for branches/eam_branches/ohana.20190329/src/opihi/mana/deimos_mkobj.c

Timestamp:

Jul 19, 2019, 5:33:23 PM (7 years ago)

Author:

eugene

Message:

avoid NANs; testing fitalt

File:

• branches/eam_branches/ohana.20190329/src/opihi/mana/deimos_mkobj.c (modified) (2 diffs)

branches/eam_branches/ohana.20190329/src/opihi/mana/deimos_mkobj.c

@@ -223 +223 @@
       dxr = 1 - dxf;
       Npof += dxi;
+      // fprintf (stderr, "%d: %f : %d\n", iy, dx, Npof);
       if (VERBOSE && (iy % 200 == 0)) { gprint (GP_ERR, "%d : %f : %d : %f %f\n", iy, dx, dxi, dxf, dxr); }
     }
@@ -241 +242 @@
           if (doInterp) {
             if ((n > 0) && (n < NpsfFull)) {
-              value += objVy*psfV[n]*dxr + objVy*psfV[n-1]*dxf;
+              if (isfinite(psfV[n]) && isfinite(psfV[n-1])) {
+                value += objVy*psfV[n]*dxr + objVy*psfV[n-1]*dxf;
+              }
             }
             if (n == 0) {
-              value += objVy*psfV[n]*dxr;
+              if (isfinite(psfV[n])) {
+                value += objVy*psfV[n]*dxr;
+              }
             }
           } else {
-            value += objVy*psfV[n];
+            if (isfinite(psfV[n])) {
+              value += objVy*psfV[n];
+            }
           }
         }