Index: /trunk/psModules/src/objects/pmPSF.c
===================================================================
--- /trunk/psModules/src/objects/pmPSF.c	(revision 9807)
+++ /trunk/psModules/src/objects/pmPSF.c	(revision 9808)
@@ -6,6 +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,8 +219,14 @@
 
         // 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,4 +304,8 @@
 // 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]