Index: trunk/Ohana/src/opihi/cmd.astro/multifit.c
===================================================================
--- trunk/Ohana/src/opihi/cmd.astro/multifit.c	(revision 16059)
+++ trunk/Ohana/src/opihi/cmd.astro/multifit.c	(revision 20936)
@@ -46,18 +46,25 @@
     sprintf (name, "c%d", i);
     if ((Nc[i]  = SelectVector (name, ANYVECTOR, TRUE)) == NULL) goto escape;
+    ResetVector (Nc[i], OPIHI_FLT, 1); 
     sprintf (name, "Mb%d", i);
     if ((NMb[i] = SelectVector (name, OLDVECTOR, TRUE)) == NULL) goto escape;
+    REQUIRE_VECTOR_FLT (NMb[i], FALSE); 
     sprintf (name, "mb%d", i);
     if ((Nmb[i] = SelectVector (name, OLDVECTOR, TRUE)) == NULL) goto escape;
+    REQUIRE_VECTOR_FLT (Nmb[i], FALSE); 
     sprintf (name, "wb%d", i);
     if ((Nwb[i] = SelectVector (name, OLDVECTOR, TRUE)) == NULL) goto escape;
+    REQUIRE_VECTOR_FLT (Nwb[i], FALSE); 
   }
   for (i = 0; i < Norder - 1; i++) {
     sprintf (name, "ml%d", i);
     if ((Nml[i] = SelectVector (name, OLDVECTOR, TRUE)) == NULL) goto escape;
+    REQUIRE_VECTOR_FLT (Nml[i], FALSE); 
     sprintf (name, "mh%d", i);
     if ((Nmh[i] = SelectVector (name, OLDVECTOR, TRUE)) == NULL) goto escape;
+    REQUIRE_VECTOR_FLT (Nmh[i], FALSE); 
     sprintf (name, "wo%d", i);
     if ((Nwo[i] = SelectVector (name, OLDVECTOR, TRUE)) == NULL) goto escape;
+    REQUIRE_VECTOR_FLT (Nwo[i], FALSE); 
   }
 
@@ -72,5 +79,5 @@
 	    v = 0;
 	    for (n = 0; n < Nwo[i-1][0].Nelements; n++) {
-	      v += -pow (Nwo[i-1][0].elements[n], (double)(j+J));
+	      v += -pow (Nwo[i-1][0].elements.Flt[n], (double)(j+J));
 	    }
 	    a[Ny][Nx] = v;
@@ -83,5 +90,5 @@
 	    v = 0;
 	    for (n = 0; n < Nwo[i][0].Nelements; n++) {
-	      v += -pow (Nwo[i][0].elements[n], (double)(j+J));
+	      v += -pow (Nwo[i][0].elements.Flt[n], (double)(j+J));
 	    }
 	    a[Ny][Nx] = v;
@@ -94,14 +101,14 @@
 	    v = 0;
 	    for (n = 0; n < Nwb[i][0].Nelements; n++) {
-	      v += pow (Nwb[i][0].elements[n], (double)(j+J));
+	      v += pow (Nwb[i][0].elements.Flt[n], (double)(j+J));
 	    }
 	    if (i > 0) {
 	      for (n = 0; n < Nwo[i-1][0].Nelements; n++) {
-		v += pow (Nwo[i-1][0].elements[n], (double)(j+J));
+		v += pow (Nwo[i-1][0].elements.Flt[n], (double)(j+J));
 	      }
 	    }
 	    if (i < Norder - 1) {
 	      for (n = 0; n < Nwo[i][0].Nelements; n++) {
-		v += pow (Nwo[i][0].elements[n], (double)(j+J));
+		v += pow (Nwo[i][0].elements.Flt[n], (double)(j+J));
 	      }
 	    }
@@ -123,17 +130,17 @@
       v = 0;
       for (n = 0; n < Nwb[i][0].Nelements; n++) {
-	v += NMb[i][0].elements[n]*pow (Nwb[i][0].elements[n], (double)j);
-	v -= Nmb[i][0].elements[n]*pow (Nwb[i][0].elements[n], (double)j);
+	v += NMb[i][0].elements.Flt[n]*pow (Nwb[i][0].elements.Flt[n], (double)j);
+	v -= Nmb[i][0].elements.Flt[n]*pow (Nwb[i][0].elements.Flt[n], (double)j);
       }
       if (i > 0) {
 	for (n = 0; n < Nwo[i-1][0].Nelements; n++) {
-	  v += Nmh[i-1][0].elements[n] * pow (Nwo[i-1][0].elements[n], (double)j);
-	  v -= Nml[i-1][0].elements[n] * pow (Nwo[i-1][0].elements[n], (double)j);
+	  v += Nmh[i-1][0].elements.Flt[n] * pow (Nwo[i-1][0].elements.Flt[n], (double)j);
+	  v -= Nml[i-1][0].elements.Flt[n] * pow (Nwo[i-1][0].elements.Flt[n], (double)j);
 	}
       }
       if (i < Norder - 1) {
 	for (n = 0; n < Nwo[i][0].Nelements; n++) {
-	  v += Nml[i][0].elements[n] * pow (Nwo[i][0].elements[n], (double)j);
-	  v -= Nmh[i][0].elements[n] * pow (Nwo[i][0].elements[n], (double)j);
+	  v += Nml[i][0].elements.Flt[n] * pow (Nwo[i][0].elements.Flt[n], (double)j);
+	  v -= Nmh[i][0].elements.Flt[n] * pow (Nwo[i][0].elements.Flt[n], (double)j);
 	}
       }
@@ -146,7 +153,7 @@
   for (i = 0; i < Norder; i++) {
     Nc[i][0].Nelements = nterm[i];
-    REALLOCATE (Nc[i][0].elements, float, nterm[i]);
+    REALLOCATE (Nc[i][0].elements.Flt, opihi_flt, nterm[i]);
     for (j = 0; j < nterm[i]; j++, Ny++) {
-      Nc[i][0].elements[j] = b[Ny][0];
+      Nc[i][0].elements.Flt[j] = b[Ny][0];
     }
   }