Changeset 30761

Timestamp:

Feb 25, 2011, 10:30:31 AM (15 years ago)

Author:

bills

Message:

when sanity checking dRA take declination into account.

File:

• trunk/Ohana/src/dvomerge/src/merge_catalogs_old.c (modified) (1 diff)

trunk/Ohana/src/dvomerge/src/merge_catalogs_old.c

@@ -219 +219 @@
       }
       if (fabs(output[0].measure[Nmeas].dR) > 10*RADIUS) {
-        fprintf (stderr, "error: %10.6f,%10.6f vs %10.6f,%10.6f (%f,%f vs %f,%f)\n", 
-                 output[0].average[n].R, output[0].average[n].D, Rin, Din,
-                 X1[i], X2[Jmin], Y1[i], Y2[Jmin]);
-        // XXX abort on this? -- this is a bad failure...
+        // ok take declination into account and check again.
+        double cosD = cos(RAD_DEG*Din);
+        if (fabs(output[0].measure[Nmeas].dR*cosD) > 10*RADIUS) {
+
+            fprintf (stderr, "error: %10.6f,%10.6f vs %10.6f,%10.6f (%f,%f vs %f,%f)\n", 
+             output[0].average[n].R, output[0].average[n].D, Rin, Din,
+             X1[i], X2[Jmin], Y1[i], Y2[Jmin]);
+
+            // XXX abort on this? -- this is a bad failure...
+        }
       }
       input[0].found[N] = Nmeas;