Changeset 20936 for trunk/Ohana/src/opihi/cmd.astro/precess.c

Timestamp:

Dec 7, 2008, 3:31:01 PM (18 years ago)

Author:

eugene

Message:

big update from eam_branch_20081124 with updates to Opihi math

File:

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

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

@@ -90 +90 @@
   }
   
-  for (i = 0; i < xvec[0].Nelements; i++) {
-    A = xvec[0].elements[i];
-    D = yvec[0].elements[i];
+  // cast xvec or yvec to FLT if it is type INT
+  CastVector (xvec, OPIHI_FLT);
+  CastVector (yvec, OPIHI_FLT);
+
+  opihi_flt *Av = xvec[0].elements.Flt;
+  opihi_flt *Dv = yvec[0].elements.Flt;
+
+  for (i = 0; i < xvec[0].Nelements; i++, Av, Dv) {
+    A = *Av;
+    D = *Dv;
     SD =  cos(RAD_DEG*A + zeta)*sin(theta)*cos(RAD_DEG*D) + cos(theta)*sin(RAD_DEG*D);
     CD = sqrt (1 - SD*SD);
@@ -101 +108 @@
     RA  = DEG_RAD*atan2(SA, CA) + z;
 
-    if (RA < 0)
-      RA += 360;
+    if (RA < 0) RA += 360;
     
-    xvec[0].elements[i] = RA;
-    yvec[0].elements[i] = DEC; 
+    *Av = RA;
+    *Dv = DEC; 
   }