Changeset 19894

Timestamp:

Oct 3, 2008, 4:29:19 PM (18 years ago)

Author:

eugene

Message:

clip negative fluxes in moments analysis

File:

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

trunk/psModules/src/objects/pmSourceMoments.c

@@ -6 +6 @@
  *  @author EAM, IfA: significant modifications.
  *
- *  @version $Revision: 1.2 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2008-10-03 21:55:42 $
+ *  @version $Revision: 1.3 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2008-10-04 02:29:19 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -211 +211 @@
             // XXX EAM : should this limit be user-defined?
             if (PS_SQR(pDiff) < wDiff) continue;
+            if (pDiff < 0) continue;
 
             Sum += pDiff;
@@ -269 +270 @@
              source->moments->Mxxxx, source->moments->Mxxxy, source->moments->Mxxyy, source->moments->Mxyyy, source->moments->Myyyy);
 
+    if (source->moments->Mxx < 0) {
+        fprintf (stderr, "error: neg second moment??\n");
+    }
+    if (source->moments->Myy < 0) {
+        fprintf (stderr, "error: neg second moment??\n");
+    }
+
     // XXX TEST:
-    pmSourceMoments_Old (source, radius);
+    // pmSourceMoments_Old (source, radius);
     return(true);
 }