Changeset 34521 for branches/eam_branches/ipp-20120905/Ohana/src/opihi

Timestamp:

Oct 7, 2012, 1:12:18 PM (14 years ago)

Author:

eugene

Message:

do not modify input vectors when rationalizing ra range

File:

• branches/eam_branches/ipp-20120905/Ohana/src/opihi/cmd.astro/cplot.c (modified) (2 diffs)

branches/eam_branches/ipp-20120905/Ohana/src/opihi/cmd.astro/cplot.c

@@ -3 +3 @@
 int cplot (int argc, char **argv) {
   
+  double ra_prev;
   int i, kapa, Npts, status, leftside;
   opihi_flt *x, *y, *r, *d, Rmin, Rmax, Rmid;
@@ -41 +42 @@
   Npts = 0;
   for (i = 0; i < Xvec.Nelements; i++, r++, d++) {
-    *r = ohana_normalize_angle (*r);
-    while (*r < Rmin) *r += 360.0;
-    while (*r > Rmax) *r -= 360.0;
+    double ra = ohana_normalize_angle (*r);
+    while (ra < Rmin) ra += 360.0;
+    while (ra > Rmax) ra -= 360.0;
 
     // for pair-by-pair connections, check on second point if we straddle the back midline
-    if ((graphmode.ptype == 100) && (i % 2)) {
-      leftside = (r[-1] < Rmid); // if first of the pair is left, second must be as well
-      if ( leftside && (r[0] > Rmid + 90)) { r[0] -= 360.0; }
-      if (!leftside && (r[0] < Rmid - 90)) { r[0] += 360.0; }
+    if (graphmode.ptype == 100) {
+      if (i % 2) {
+        leftside = (ra_prev < Rmid); // if first of the pair is left, second must be as well
+        if ( leftside && (ra > Rmid + 90)) { ra -= 360.0; }
+        if (!leftside && (ra < Rmid - 90)) { ra += 360.0; }
+      } else {
+        ra_prev = ra;
+      }
     }
-    status = RD_to_XY (x, y, *r, *d, &graphmode.coords);
+    status = RD_to_XY (x, y, ra, *d, &graphmode.coords);
 
     // if we fail on one of the points, drop the corresponding pair