Changeset 27488 for branches/eam_branches/relastro.20100326/src/fitpoly.c

Timestamp:

Mar 26, 2010, 4:43:50 PM (16 years ago)

Author:

eugene

Message:

updates to handle image coordinate failures

File:

• branches/eam_branches/relastro.20100326/src/fitpoly.c (modified) (6 diffs)

branches/eam_branches/relastro.20100326/src/fitpoly.c

@@ -103 +103 @@
 
 /* convert the xsum,ysum,sum terms into vector,matrix and solve */
-void fit_eval (CoordFit *fit) {
+int fit_eval (CoordFit *fit) {
 
   int i, j, ix, iy, jx, jy;
@@ -147 +147 @@
   }     
 
-  dgaussjordan (matrix, vector, fit[0].Nelems, 2); 
+  status = dgaussjordan (matrix, vector, fit[0].Nelems, 2); 
+  if (!status) {
+    return (FALSE);
+  }
 
   for (i = 0; i < fit[0].Nelems; i++) {
@@ -166 +169 @@
   array_free (matrix, fit[0].Nelems);
   array_free (vector, fit[0].Nelems);
+  return (TRUE);
 }
 
@@ -271 +275 @@
 /* this should only apply to the polynomial, not the projection terms */
 /* compare with psastro supporting code */
-CoordFit *fit_apply_coords (CoordFit *fit, Coords *coords) {
+int fit_apply_coords (CoordFit *fit, Coords *coords) {
 
   double Xo, Yo, R1, R2;
@@ -281 +285 @@
   // L = pc1_1*cd1*(x - cp1) + pc1_2*cd2*(y - cp2) + ...
 
-  CoordsGetCenter (fit, 0.001, &Xo, &Yo);
+  if (!CoordsGetCenter (fit, 0.001, &Xo, &Yo)) {
+    fprintf (stderr, "failed to modify model\n");
+    return (FALSE);
+  }
   coords[0].crpix1 = Xo;
   coords[0].crpix2 = Yo;
@@ -329 +336 @@
   /* keep the order and type from initial values */
   
-  // XXX if desired in the future, return modfit (and free above)
-  return (NULL);
-}
-
+  return (TRUE);
+}
+