Changeset 5813 for trunk/psLib/src/astro/psCoord.c

Timestamp:

Dec 19, 2005, 1:58:47 PM (20 years ago)

Author:

gusciora

Message:

Implemented the chi-squared algorithm for calculated Chebyshev poly
coefficients.

File:

• trunk/psLib/src/astro/psCoord.c (modified) (4 diffs)

trunk/psLib/src/astro/psCoord.c

@@ -10 +10 @@
 *  @author GLG, MHPCC
 *
-*  @version $Revision: 1.96 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2005-12-05 21:33:29 $
+*  @version $Revision: 1.97 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2005-12-19 23:58:47 $
 *
 *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -61 +61 @@
 should rename this.
  
-XXX: Use the ADD version which is based on determinants.
- 
 XXX: This code no longer makes sense.  The merge must be reviewed.
  *****************************************************************************/
@@ -85 +83 @@
 
     if (0) {
-        //This is sample code from IfA.  It didn't work initially, and I did not
-        //spend any time debugging it.
-        psF64 a = transform->x->coeff[1][0];
-        psF64 b = transform->x->coeff[0][1];
-        psF64 c = transform->y->coeff[1][0];
-        psF64 d = transform->y->coeff[0][1];
-        psF64 e = transform->x->coeff[0][0];
-        psF64 f = transform->y->coeff[0][0];
-
-        psF64 invDet = 1.0 / (a * d - b * c); // Inverse of the determinant
-
-        // Not entirely sure why this works, but it appears to do so
-        out->x->coeff[1][0] = invDet * a;
-        out->x->coeff[0][1] = - invDet * b;
-        out->y->coeff[1][0] = - invDet * c;
-        out->y->coeff[0][1] = invDet * d;
-
-        out->x->coeff[0][0] = - invDet * (d * e + c * f);
-        out->y->coeff[0][0] = - invDet * (b * e + a * f);
-    }
-
-    if (1) {
-        // XXX: There is no reason to execute this code and the following code.
+        // XXX: Get rid of this code.
         psF64 A = transform->x->coeff[1][0];
         psF64 B = transform->x->coeff[0][1];
@@ -129 +105 @@
     }
 
-
-    // XXX: There is no reason to execute this code and the previous codes.
-    // unless the cross terms are available, set these matrix elements to 0
     psF64 r12 = 0.0;
     if (transform->x->nY == 1) {