Index: trunk/Ohana/src/opihi/lib.data/mrqmin.c
===================================================================
--- trunk/Ohana/src/opihi/lib.data/mrqmin.c	(revision 16059)
+++ trunk/Ohana/src/opihi/lib.data/mrqmin.c	(revision 20936)
@@ -7,15 +7,15 @@
 */
 
-static float **alpha, **talpha;
-static float **beta, **tbeta;
-static float *partry, *dyda;
-static float ochisq, lambda;
+static opihi_flt **alpha, **talpha;
+static opihi_flt **beta, **tbeta;
+static opihi_flt *partry, *dyda;
+static opihi_flt ochisq, lambda;
 
-float mrqcof (float *x, float *y, float *dy, int Npts, 
-	      float *par, int Npar, float **ta, float **tb, 
-	      float (funcs)(float, float *, int, float *)) {
+opihi_flt mrqcof (opihi_flt *x, opihi_flt *y, opihi_flt *dy, int Npts, 
+	      opihi_flt *par, int Npar, opihi_flt **ta, opihi_flt **tb, 
+	      opihi_flt (funcs)(opihi_flt, opihi_flt *, int, opihi_flt *)) {
 
   int k, j, i;
-  float ydiff, wt, chisq;
+  opihi_flt ydiff, wt, chisq;
 
   for (j = 0; j < Npar; j++) {
@@ -45,11 +45,11 @@
 }
 
-float mrqmin (float *x, float *y, float *dy, int Npts, 
-	      float *par, int Npar, 
-	      float (funcs)(float, float *, int, float *), int VERBOSE) {
+opihi_flt mrqmin (opihi_flt *x, opihi_flt *y, opihi_flt *dy, int Npts, 
+	      opihi_flt *par, int Npar, 
+	      opihi_flt (funcs)(opihi_flt, opihi_flt *, int, opihi_flt *), int VERBOSE) {
 
   int j, k;
-  float chisq;
-  float rho, dX, dL;
+  opihi_flt chisq;
+  opihi_flt rho, dX, dL;
 
   /* set up test matrixes for this run */
@@ -60,5 +60,5 @@
   }
 
-  fgaussjordan (talpha, tbeta, Npar, 1);
+  dgaussjordan (talpha, tbeta, Npar, 1);
 
   for (j = 0; j < Npar; j++) partry[j] = par[j] - tbeta[j][0];
@@ -102,21 +102,21 @@
 }
 
-float mrqinit (float *x, float *y, float *dy, int Npts, 
-	      float *par, int Npar, 
-	      float (funcs)(float, float *, int, float *), int VERBOSE) {
+opihi_flt mrqinit (opihi_flt *x, opihi_flt *y, opihi_flt *dy, int Npts, 
+	      opihi_flt *par, int Npar, 
+	      opihi_flt (funcs)(opihi_flt, opihi_flt *, int, opihi_flt *), int VERBOSE) {
 
   int i;
 
-  ALLOCATE (dyda, float, Npar);
-  ALLOCATE (partry, float, Npar);
-  ALLOCATE (alpha, float *, Npar);
-  ALLOCATE (beta, float *, Npar);
-  ALLOCATE (talpha, float *, Npar);
-  ALLOCATE (tbeta, float *, Npar);
+  ALLOCATE (dyda, opihi_flt, Npar);
+  ALLOCATE (partry, opihi_flt, Npar);
+  ALLOCATE (alpha, opihi_flt *, Npar);
+  ALLOCATE (beta, opihi_flt *, Npar);
+  ALLOCATE (talpha, opihi_flt *, Npar);
+  ALLOCATE (tbeta, opihi_flt *, Npar);
   for (i = 0; i < Npar; i++) {
-    ALLOCATE (alpha[i], float, Npar);
-    ALLOCATE (beta[i], float, Npar);
-    ALLOCATE (talpha[i], float, Npar);
-    ALLOCATE (tbeta[i], float, Npar);
+    ALLOCATE (alpha[i], opihi_flt, Npar);
+    ALLOCATE (beta[i], opihi_flt, Npar);
+    ALLOCATE (talpha[i], opihi_flt, Npar);
+    ALLOCATE (tbeta[i], opihi_flt, Npar);
   }
   
@@ -136,6 +136,6 @@
 
 /* don't invoke this in the middle of a run, only near the end */ 
-float **mrqcovar (int Npar) {
-  fgaussjordan (alpha, beta, Npar, 1);
+opihi_flt **mrqcovar (int Npar) {
+  dgaussjordan (alpha, beta, Npar, 1);
   return (alpha);
 }