Changeset 20010 for trunk/psModules/src/objects/pmSource.c

Timestamp:

Oct 8, 2008, 2:58:39 PM (18 years ago)

Author:

Paul Price

Message:

Adding check for NAN moments.

File:

• trunk/psModules/src/objects/pmSource.c (modified) (3 diffs)

trunk/psModules/src/objects/pmSource.c

@@ -6 +6 @@
  *  @author EAM, IfA: significant modifications.
  *
- *  @version $Revision: 1.61 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2008-10-08 04:23:31 $
+ *  @version $Revision: 1.62 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2008-10-09 00:58:39 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -347 +347 @@
             float ar = Mxx / Myy;       // Radius
 
+            if (!isfinite(Mxx) || !isfinite(Myy)) {
+                psTrace("psModules.objects", 10,
+                        "Rejecting source from clump because of non-finite moments (%f,%f)\n",
+                        Mxx, Myy);
+                continue;
+            }
+
             // Sx,Sy are limited at 0.  a peak at 0,0 is artificial
             if (fabs(Mxx) < 0.05 || fabs(Myy < 0.05)) {
@@ -612 +619 @@
 
     if (starsn_moments->n) {
-        psStats *stats = NULL;
-        stats = psStatsAlloc (PS_STAT_MIN | PS_STAT_MAX);
-
-        if (!psVectorStats (stats, starsn_moments, NULL, NULL, 0)) {
-            // Don't care about this error
-            psErrorClear();
-        }
-        psLogMsg ("pmObjects", 3, "SN range (moments): %f - %f\n", stats->min, stats->max);
-        psFree (stats);
+        psStats *stats = NULL;
+        stats = psStatsAlloc (PS_STAT_MIN | PS_STAT_MAX);
+
+        if (!psVectorStats (stats, starsn_moments, NULL, NULL, 0)) {
+            // Don't care about this error
+            psErrorClear();
+        }
+        psLogMsg ("pmObjects", 3, "SN range (moments): %f - %f\n", stats->min, stats->max);
+        psFree (stats);
     }
     psFree (starsn_moments);
 
     if (starsn_peaks->n) {
-        psStats *stats = NULL;
-        stats = psStatsAlloc (PS_STAT_MIN | PS_STAT_MAX);
-        if (!psVectorStats (stats, starsn_peaks, NULL, NULL, 0)) {
-            // Don't care about this error
-            psErrorClear();
-        }
-        psLogMsg ("pmObjects", 3, "SN range (peaks)  : %f - %f\n", stats->min, stats->max);
-        psFree (stats);
+        psStats *stats = NULL;
+        stats = psStatsAlloc (PS_STAT_MIN | PS_STAT_MAX);
+        if (!psVectorStats (stats, starsn_peaks, NULL, NULL, 0)) {
+            // Don't care about this error
+            psErrorClear();
+        }
+        psLogMsg ("pmObjects", 3, "SN range (peaks)  : %f - %f\n", stats->min, stats->max);
+        psFree (stats);
     }
     psFree (starsn_peaks);