Changeset 9561

Timestamp:

Oct 13, 2006, 2:54:37 PM (20 years ago)

Author:

magnier

Message:

test for NAN chisq, chisqNorm values

File:

• trunk/psModules/src/objects/pmSourceIO_RAW.c (modified) (2 diffs)

trunk/psModules/src/objects/pmSourceIO_RAW.c

@@ -3 +3 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.6 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-10-12 21:41:22 $
+ *  @version $Revision: 1.7 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-10-14 00:54:37 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -97 +97 @@
             fprintf (f, "%9.6f ", dPAR[j]);
         }
+
+        float logChi = ((model[0].chisq == 0.0) || (model[0].nDOF == 0)) ? NAN : log10(model[0].chisq/model[0].nDOF);
+        float logChiNorm = ((model[0].chisqNorm == 0.0) || (model[0].nDOF == 0)) ? NAN : log10(model[0].chisqNorm/model[0].nDOF);
+
         fprintf (f, ": %8.4f %2d %#5x %7.3f %7.3f  %7.1f %7.2f %4.2f %4d %2d\n",
                  source[0].apMag, source[0].type, source[0].mode,
-                 log10(model[0].chisq/model[0].nDOF),
-                 log10(model[0].chisqNorm/model[0].nDOF),
+                 logChi, logChiNorm,
                  source[0].moments->SN,
                  model[0].radiusTMP,