Changeset 39168 for trunk/Ohana/src/opihi/cmd.astro/fitpm_irls.c

Timestamp:

Nov 19, 2015, 5:43:11 PM (11 years ago)

Author:

watersc1

Message:

Fixed transposition that was critically broken for fitplx_irls, but seemed to have dodged through in fitpm_irls.

File:

• trunk/Ohana/src/opihi/cmd.astro/fitpm_irls.c (modified) (2 diffs)

trunk/Ohana/src/opihi/cmd.astro/fitpm_irls.c

@@ -231 +231 @@
   sigma_ols = 0.0;
   for (i = 0; i < Npts; i++) {
-    rx[i] = X[i] - (T[i] * B[0][0] + B[1][0]);
-    ry[i] = Y[i] - (T[i] * B[2][0] + B[3][0]);
+    rx[i] = X[i] - (T[i] * B[1][0] + B[0][0]);
+    ry[i] = Y[i] - (T[i] * B[3][0] + B[2][0]);
     //    u[i] = r[i] /
     sigma_ols += SQ(rx[i]) + SQ(ry[i]);
@@ -280 +280 @@
     sigma_hat = 0.0;
     for (i = 0; i < Npts; i++) {
-      rx[i] = X[i] - (T[i] * B[0][0] + B[1][0]);
-      ry[i] = Y[i] - (T[i] * B[2][0] + B[3][0]);
+      rx[i] = X[i] - (T[i] * B[1][0] + B[0][0]);
+      ry[i] = Y[i] - (T[i] * B[3][0] + B[2][0]);
       u[i] = sqrt(SQ(rx[i] / dX[i]) + SQ(ry[i] / dY[i]));
     }