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

Timestamp:

Oct 7, 2008, 6:21:47 PM (18 years ago)

Author:

eugene

Message:

skip stats for empty vectors

File:

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

trunk/psModules/src/objects/pmSource.c

@@ -6 +6 @@
  *  @author EAM, IfA: significant modifications.
  *
- *  @version $Revision: 1.59 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2008-10-07 20:12:27 $
+ *  @version $Revision: 1.60 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2008-10-08 04:21:47 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -611 +611 @@
     }
 
-    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[0].min, stats[0].max);
-    psFree (starsn_moments);
-    psFree (stats);
-
-    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[0].min, stats[0].max);
-    psFree (starsn_peaks);
-    psFree (stats);
+    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 (starsn_moments);
+        psFree (stats);
+    }
+
+    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 (starsn_peaks);
+        psFree (stats);
+    }
 
     psTrace ("psModules.objects", 2, "Nstar:    %3d\n", Nstar);