Changeset 10031

Timestamp:

Nov 16, 2006, 4:39:34 PM (20 years ago)

Author:

magnier

Message:

moved pmSourceChisq into here

Location:

trunk/psModules/src/objects

Files:

• pmSourcePhotometry.c (modified) (2 diffs)
• pmSourcePhotometry.h (modified) (2 diffs)

trunk/psModules/src/objects/pmSourcePhotometry.c

@@ -3 +3 @@
  *  @author EAM, IfA; GLG, MHPCC
  *
- *  @version $Revision: 1.17 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-11-09 20:17:10 $
+ *  @version $Revision: 1.18 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-11-17 02:39:34 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -487 +487 @@
 }
 
+bool pmSourceChisq (pmModel *model, psImage *image, psImage *mask, psImage *weight)
+{
+
+    double dC = 0.0;
+    int Npix = 0;
+    for (int j = 0; j < image->numRows; j++) {
+        for (int i = 0; i < image->numCols; i++) {
+            if (mask->data.U8[j][i])
+                continue;
+            if (weight->data.F32[j][i] <= 0)
+                continue;
+            dC += PS_SQR (image->data.F32[j][i]) / weight->data.F32[j][i];
+            Npix ++;
+        }
+    }
+    model->nDOF = Npix - 1;
+    model->chisq = dC;
+
+    return (true);
+}

trunk/psModules/src/objects/pmSourcePhotometry.h

@@ -3 +3 @@
  *  @author EAM, IfA; GLG, MHPCC
  *
- *  @version $Revision: 1.5 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-11-07 09:11:38 $
+ *  @version $Revision: 1.6 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-11-17 02:39:34 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -51 +51 @@
 double pmSourceCrossWeight(const pmSource *Mi, const pmSource *Mj, const bool unweighted_sum);
 bool pmSourcePixelWeight (float *pixWeight, pmModel *model, psImage *image, psImage *mask);
-
-psImage *tmp_psImageShift (psImage *in, float dx, float dy);
+bool pmSourceChisq (pmModel *model, psImage *image, psImage *mask, psImage *weight);
 
 # endif /* PM_SOURCE_PHOTOMETRY_H */