Changeset 9808

Timestamp:

Oct 31, 2006, 9:38:04 AM (20 years ago)

Author:

magnier

Message:

if-def out some test code, comments describing SXY fitting issues

File:

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

trunk/psModules/src/objects/pmPSF.c

@@ -6 +6 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.13 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-10-30 03:04:03 $
+ *  @version $Revision: 1.14 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-10-31 19:38:04 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -219 +219 @@
 
         // XXX Test output
-        psPolynomial2D *poly = psf->params_NEW->data[i];
+        # if (0)
+
+            psPolynomial2D *poly = psf->params_NEW->data[i];
         fprintf (stderr, "stats: %f +/- %f\n", stats->robustMedian, stats->robustStdev);
         fprintf (stderr, "PO: %g %g %g\n", poly->coeff[0][0], poly->coeff[1][0], poly->coeff[0][1]);
         fprintf (stderr, "PO: %g %g %g\n", poly->coeff[0][2], poly->coeff[1][1], poly->coeff[2][0]);
+        fprintf (stderr, "PO: %g %g %g\n", poly->coeff[0][0], poly->coeff[1][0], poly->coeff[0][1]);
+        fprintf (stderr, "PO: %g %g %g\n", poly->coeff[0][2], poly->coeff[1][1], poly->coeff[2][0]);
+        # endif
+
     }
 
@@ -298 +304 @@
 // these are the values of the model->params.  the psf->params term for sxy is actually fitted
 // to sxy/(sxx^-2 + syy^-2)^2
+
+// XXX this is only an approximate solution.  A better solution would be to fit the second moment, Mxy. the
+// problem here is that converting from Mxy,SXX,SYY -> SXY is a third order problem:
+// Mxy = SXY * (SXX^-4 + SYY^-4 - 2 SXY ^2)
 
 // input: model->param, output: psf->param[PM_PAR_SXY]