Index: /branches/eam_branches/ipp-20130509/psModules/src/astrom/pmAstrometryUtils.c
===================================================================
--- /branches/eam_branches/ipp-20130509/psModules/src/astrom/pmAstrometryUtils.c	(revision 35647)
+++ /branches/eam_branches/ipp-20130509/psModules/src/astrom/pmAstrometryUtils.c	(revision 35648)
@@ -78,7 +78,10 @@
             Alpha->data.F32[1][1] = psPolynomial2DEval (YdY, Xo, Yo);
 
+	    // Beta = (-L,-M) for the current guess
             Beta->data.F32[0] = -1.0 * psPolynomial2DEval (trans->x, Xo, Yo);
             Beta->data.F32[1] = -1.0 * psPolynomial2DEval (trans->y, Xo, Yo);
-
+            // fprintf (stderr, " Beta: %f %f\n", Beta->data.F32[0], Beta->data.F32[1]);
+
+	    // since we want (L,M) = (0,0), Beta is also the offset
             if (!psMatrixGJSolve (Alpha, Beta)) {
 		psError(PS_ERR_UNKNOWN, false, "Unable to solve for center.");