Changeset 25980 for trunk/psModules/src/objects/pmSourcePhotometry.c

Timestamp:

Nov 1, 2009, 3:53:37 PM (17 years ago)

Author:

eugene

Message:

add covarianceFactor to pmSourceChisq, pmSourceModelWeight, pmSourceModelDotModel, pmSourceDataDotModel

File:

• trunk/psModules/src/objects/pmSourcePhotometry.c (modified) (8 diffs)

trunk/psModules/src/objects/pmSourcePhotometry.c

@@ -545 +545 @@
 
 // determine chisq, etc for linear normalization-only fit
-bool pmSourceChisq (pmModel *model, psImage *image, psImage *mask, psImage *variance,
-                    psImageMaskType maskVal)
+bool pmSourceChisq (pmModel *model, psImage *image, psImage *mask, psImage *variance, psImageMaskType maskVal, const float covarFactor)
 {
     PS_ASSERT_PTR_NON_NULL(model, false);
@@ -561 +560 @@
             if (variance->data.F32[j][i] <= 0)
                 continue;
-            dC += PS_SQR (image->data.F32[j][i]) / variance->data.F32[j][i];
+            dC += PS_SQR (image->data.F32[j][i]) / (covarFactor * variance->data.F32[j][i]);
             Npix ++;
         }
@@ -573 +572 @@
 
 
-double pmSourceModelWeight(const pmSource *Mi,
-                      int term,
-                      const bool unweighted_sum) // should the cross product be weighted?
+double pmSourceModelWeight(const pmSource *Mi, int term, const bool unweighted_sum, const float covarFactor)
 {
     PS_ASSERT_PTR_NON_NULL(Mi, NAN);
@@ -594 +591 @@
                 continue;
             if (!unweighted_sum) {
-                wt = Wi->data.F32[yi][xi];
+                wt = covarFactor * Wi->data.F32[yi][xi];
                 if (wt == 0)
                     continue;
@@ -623 +620 @@
 }
 
-double pmSourceModelDotModel (const pmSource *Mi,
-                              const pmSource *Mj,
-                              const bool unweighted_sum) // should the cross product be weighted?
+double pmSourceModelDotModel (const pmSource *Mi, const pmSource *Mj, const bool unweighted_sum, const float covarFactor)
 {
     PS_ASSERT_PTR_NON_NULL(Mi, NAN);
@@ -677 +672 @@
                 flux += (Pi->data.F32[yi][xi] * Pj->data.F32[yj][xj]);
             } else {
-                wt = Wi->data.F32[yi][xi];
+                wt = covarFactor * Wi->data.F32[yi][xi];
                 if (wt > 0) {
                     flux += (Pi->data.F32[yi][xi] * Pj->data.F32[yj][xj]) / wt;
@@ -687 +682 @@
 }
 
-double pmSourceDataDotModel (const pmSource *Mi,
-                             const pmSource *Mj,
-                             const bool unweighted_sum) // should the cross product be weighted?
+double pmSourceDataDotModel (const pmSource *Mi, const pmSource *Mj, const bool unweighted_sum, const float covarFactor)
 {
     PS_ASSERT_PTR_NON_NULL(Mi, NAN);
@@ -741 +734 @@
                 flux += (Pi->data.F32[yi][xi] * Pj->data.F32[yj][xj]);
             } else {
-                wt = Wi->data.F32[yi][xi];
+                wt = covarFactor * Wi->data.F32[yi][xi];
                 if (wt > 0) {
                     flux += (Pi->data.F32[yi][xi] * Pj->data.F32[yj][xj]) / wt;