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

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

@@ -2 +2 @@
 
 /* local private functions */
-float fmaxwellOD (float, float *, int, float *);
+opihi_flt fmaxwellOD (opihi_flt, opihi_flt *, int, opihi_flt *);
 
 # define GET_VAR(V,A) \
@@ -16 +16 @@
 
   int i, N, Npts, Npar, Quiet;
-  float par[5], *v1, *v2, *dy, **covar;
-  float chisq, ochisq, dchisq;
+  opihi_flt par[5], *v1, *v2, *dy, **covar;
+  opihi_flt chisq, ochisq, dchisq;
   Vector *xvec, *yvec, *svec, *ovec;
   char *c, name[16];
@@ -43 +43 @@
 
   Npts = xvec[0].Nelements;
-  ALLOCATE (dy, float, Npts);
-  REALLOCATE (ovec[0].elements, float, Npts);
+  ALLOCATE (dy, opihi_flt, Npts);
+  REALLOCATE (ovec[0].elements.Flt, opihi_flt, Npts);
 
   GET_VAR (par[0], "C0");
@@ -54 +54 @@
   /* careful of variable renomalization */
 
-  v1 = svec[0].elements;
+  v1 = svec[0].elements.Flt;
   v2 = dy;
   for (i = 0; i < Npts; i++, v1++, v2++) *v2 = 1.0 / (*v1 * *v1);
   
-  ochisq = mrqinit (xvec[0].elements, yvec[0].elements, dy, Npts, par, Npar, fmaxwellOD, !Quiet);
+  ochisq = mrqinit (xvec[0].elements.Flt, yvec[0].elements.Flt, dy, Npts, par, Npar, fmaxwellOD, !Quiet);
   dchisq = ochisq + 2*Npts;
 
   for (i = 0; (i < 20) && ((dchisq > 0.1*(Npts - Npar)) || (dchisq <= 0.0)); i++) {
-    chisq = mrqmin (xvec[0].elements, yvec[0].elements, dy, Npts, par, Npar, fmaxwellOD, !Quiet);
+    chisq = mrqmin (xvec[0].elements.Flt, yvec[0].elements.Flt, dy, Npts, par, Npar, fmaxwellOD, !Quiet);
     dchisq = ochisq - chisq;
     ochisq = chisq;
@@ -70 +70 @@
 
   for (i = 0; i < Npts; i++) {
-    ovec[0].elements[i] = fmaxwellOD (xvec[0].elements[i], par, Npar, dy);
+    ovec[0].elements.Flt[i] = fmaxwellOD (xvec[0].elements.Flt[i], par, Npar, dy);
   }
   ovec[0].Nelements = Npts;
@@ -90 +90 @@
 
 /* pars: x_o, -0.5/sigma^2, I, back, ref */
-float fmaxwellOD (float x, float *par, int Npar, float *dpar) {
+opihi_flt fmaxwellOD (opihi_flt x, opihi_flt *par, int Npar, opihi_flt *dpar) {
 
-  float z, r, f;
+  opihi_flt z, r, f;
 
   z = (x - par[0])/par[1];