Changeset 10076 for trunk/psphot/src/psphotEnsemblePSF.c

Timestamp:

Nov 17, 2006, 9:35:22 PM (20 years ago)

Author:

eugene

Message:

working on new linear fit process

File:

• trunk/psphot/src/psphotEnsemblePSF.c (modified) (6 diffs)

trunk/psphot/src/psphotEnsemblePSF.c

@@ -10 +10 @@
 // the analysis is performed wrt the simulated pixel values
 
-static bool SetBorderT (psSparseBorder *border, pmReadout *readout, psArray *refSources, bool fitSlope);
+static bool SetBorderT (psSparseBorder *border, pmReadout *readout, psArray *refSources, psArray *fitSources, psVector *index, bool fitSlope);
 
 bool psphotEnsemblePSF (pmReadout *readout, psArray *refSources, psMetadata *recipe, pmPSF *psf, bool final) {
@@ -162 +162 @@
     // set the sky, sky_x, sky_y components of border matrix
     // XXX ignore sky slope terms for now
-    SetBorderT (border, readout, refSources, false);
+    SetBorderT (border, readout, refSources, fitSources, index, false);
 
     psSparseConstraint constraint;
@@ -170 +170 @@
 
     // solve for normalization terms (need include local sky?)
+    # if (0)
     psVector *norm = NULL;
     psVector *skyfit = NULL;
     psSparseBorderSolve (&norm, &skyfit, constraint, border, 3);
+    # else
+    psVector *norm = psSparseSolve (NULL, constraint, sparse, 3);
+    # endif
 
     // adjust fitSources, set refSources and subtract
@@ -215 +219 @@
     psFree (weight);
     psFree (errors);
+    psFree (border);
 
     psLogMsg (__func__, 3, "measure ensemble of PSFs: %f sec\n", psTimerMark ("psphot"));
@@ -221 +226 @@
 
 // calculate the weight terms for the sky fit component of the matrix
-static bool SetBorderT (psSparseBorder *border, pmReadout *readout, psArray *refSources, bool fitSlope) {
+static bool SetBorderT (psSparseBorder *border, pmReadout *readout, psArray *refSources, psArray *fitSources, psVector *index, bool fitSlope) {
 
     // generate the image-wide weight terms
@@ -231 +236 @@
     // turn off MARK for all object pixels
     // we must use the refSources here since their pixels point at the science image
-    for (int i = 0; i < refSources->n; i++) {
-        pmSource *refSource = refSources->data[i];
+    for (int i = 0; i < fitSources->n; i++) {
+        int N = index->data.U32[i];
+        pmSource *refSource = refSources->data[N];
         pmModel *model = pmSourceGetModel (NULL, refSource);
+        if (model == NULL) continue;
         float x = model->params->data.F32[PM_PAR_XPOS];
         float y = model->params->data.F32[PM_PAR_YPOS];