Changeset 20839 for branches/eam_branch_20081124/Ohana/src/opihi/cmd.data/svd.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.data/svd.c (modified) (2 diffs)

branches/eam_branch_20081124/Ohana/src/opihi/cmd.data/svd.c

@@ -4 +4 @@
   
   int i, Nx, Ny, status;
-  float *in, *out, *A, *U, *W, *V;
+  float *in, *out, *A, *U, *V;
+  Buffer *Ma, *Mu, *Mv;
   Vector *Vw;
-  Buffer *Ma, *Mu, *Mv;
+  opihi_flt *W;
 
   if (argc != 6) goto usage;
@@ -43 +44 @@
   /* w is Nx */
   Vw[0].Nelements = Nx;
-  REALLOCATE (Vw[0].elements, float, Nx);
+  REALLOCATE (Vw[0].elements.Flt, opihi_flt, Nx);
 
   /* pointers to the various arrays */
   A = (float *) Ma[0].matrix.buffer;
   U = (float *) Mu[0].matrix.buffer;
-  W = (float *) Vw[0].elements;
   V = (float *) Mv[0].matrix.buffer;
+  W = Vw[0].elements.Flt;
 
   /* copy A to U (svdcmp replaces A with U) */