Changeset 3896 for trunk/Ohana/src/opihi/cmd.data

Timestamp:

May 11, 2005, 9:39:17 PM (21 years ago)

Author:

eugene

Message:

fixed nterm vs norder error

File:

• trunk/Ohana/src/opihi/cmd.data/fit.c (modified) (4 diffs)

trunk/Ohana/src/opihi/cmd.data/fit.c

@@ -56 +56 @@
   }
   
+  /* nterm is number of polynomial terms, starting at x^0 */
   order = atof (argv[3]);
   nterm = order + 1;
-  mterm = 2*order + 1;
+  mterm = 2*nterm;
 
   ALLOCATE (yfit, float, xvec[0].Nelements);
@@ -64 +65 @@
   memset (mask, 0, xvec[0].Nelements);
 
-  ALLOCATE (s, double, 2*mterm);
+  ALLOCATE (s, double, mterm);
   ALLOCATE (b, double *, nterm);
   ALLOCATE (c, double *, nterm);
@@ -75 +76 @@
 
     /* init registers for current pass */
-    memset (s, 0, 2*mterm*sizeof(double));
+    memset (s, 0, mterm*sizeof(double));
     for (i = 0; i < nterm; i++) {
       memset (c[i], 0, nterm*sizeof(double));
@@ -85 +86 @@
     y = yvec[0].elements;
     if (Weight) dy = dyvec[0].elements;
+
     for (i = 0; i < xvec[0].Nelements; i++, x++, y++) {
       if (mask[i]) continue;