Changeset 20839 for branches/eam_branch_20081124/Ohana/src/opihi/cmd.astro/profile.c

Timestamp:

Nov 25, 2008, 4:51:39 PM (18 years ago)

Author:

eugene

Message:

convert opihi vectors to double type, add option for int type as well

File:

• branches/eam_branch_20081124/Ohana/src/opihi/cmd.astro/profile.c (modified) (2 diffs)

branches/eam_branch_20081124/Ohana/src/opihi/cmd.astro/profile.c

@@ -29 +29 @@
 
   xvec[0].Nelements = yvec[0].Nelements = (int)SQ(2*N+1);
-  REALLOCATE (xvec[0].elements, float, 5*(int)SQ(2*N+1));
-  REALLOCATE (yvec[0].elements, float, 5*(int)SQ(2*N+1));
-  bzero (yvec[0].elements, (int)SQ(2*N+1)*sizeof(float)+1);
+  REALLOCATE (xvec[0].elements.Flt, opihi_flt, 5*(int)SQ(2*N+1));
+  REALLOCATE (yvec[0].elements.Flt, opihi_flt, 5*(int)SQ(2*N+1));
+  bzero (yvec[0].elements.Flt, (int)SQ(2*N+1)*sizeof(opihi_flt)+1);
   V = (float *)(buf[0].matrix.buffer); 
   Npt = 0;
@@ -37 +37 @@
   for (i = sx - N; i < sx + N; i++) {
     for (j = sy - N; j < sy + N; j++, Npt++) {
-      yvec[0].elements[Npt] = V[i + j*Nx];
-      xvec[0].elements[Npt] = hypot (i - sx, j - sy);
+      yvec[0].elements.Flt[Npt] = V[i + j*Nx];
+      xvec[0].elements.Flt[Npt] = hypot (i - sx, j - sy);
     }
   }
 
-  fsortpair (xvec[0].elements, yvec[0].elements, xvec[0].Nelements);
+  dsortpair (xvec[0].elements.Flt, yvec[0].elements.Flt, xvec[0].Nelements);
 
   return (TRUE);