Changeset 7917 for trunk/Ohana/src/opihi/cmd.data/fit2d.c

Timestamp:

Jul 16, 2006, 10:58:49 PM (20 years ago)

Author:

eugene

Message:

major updates to use gprint for buffered / threaded printing

File:

• trunk/Ohana/src/opihi/cmd.data/fit2d.c (modified) (7 diffs)

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

@@ -48 +48 @@
 
   if (argc != 5) {
-    fprintf (stderr, "USAGE: fit x y z order [-dz wt]\n");
+    gprint (GP_ERR, "USAGE: fit x y z order [-dz wt]\n");
     return (FALSE);
   }
@@ -57 +57 @@
 
   if (xvec[0].Nelements != yvec[0].Nelements) {
-    fprintf (stderr, "vectors must have same length\n");
+    gprint (GP_ERR, "vectors must have same length\n");
     return (FALSE);
   }
   if (xvec[0].Nelements != zvec[0].Nelements) {
-    fprintf (stderr, "vectors must have same length\n");
+    gprint (GP_ERR, "vectors must have same length\n");
     return (FALSE);
   }
   if (Weight && xvec[0].Nelements != dzvec[0].Nelements) {
-    fprintf (stderr, "vectors must have same length\n");
+    gprint (GP_ERR, "vectors must have same length\n");
     return (FALSE);
   }
@@ -160 +160 @@
     }
 
-    fprintf (stderr, "matrix:\n");
+    gprint (GP_ERR, "matrix:\n");
     for (k = 0; k < wterm; k++) {
         for (K = 0; K < wterm; K++) {
-            fprintf (stderr, "%f ", c[k][K]);
-        }
-        fprintf (stderr, "\n");
+            gprint (GP_ERR, "%f ", c[k][K]);
+        }
+        gprint (GP_ERR, "\n");
     }
 
     gaussj (c, wterm, b, 1);
 
-    fprintf (stderr, "inverse:\n");
+    gprint (GP_ERR, "inverse:\n");
     for (k = 0; k < wterm; k++) {
         for (K = 0; K < wterm; K++) {
-            fprintf (stderr, "%f ", c[k][K]);
-        }
-        fprintf (stderr, "\n");
+            gprint (GP_ERR, "%f ", c[k][K]);
+        }
+        gprint (GP_ERR, "\n");
     }
 
@@ -182 +182 @@
       for (i = ny = 0; ny < nterm; ny++) {
         for (nx = 0; nx < nterm - ny; nx++, i++) {
-          fprintf (stderr, "x^%d y^%d: %g\n", nx, ny, b[i][0]);
+          gprint (GP_ERR, "x^%d y^%d: %g\n", nx, ny, b[i][0]);
         }
       }
@@ -228 +228 @@
     maxsigma = ClipNSigma * sigma;
 
-    if (VERBOSE) fprintf (stderr, "mean: %g, sigma: %g, maxsigma: %g\n", mean, sigma, maxsigma);
+    if (VERBOSE) gprint (GP_ERR, "mean: %g, sigma: %g, maxsigma: %g\n", mean, sigma, maxsigma);
 
     /* mask outlier points */
@@ -245 +245 @@
       } 
     }
-    if (VERBOSE) fprintf (stderr, "pass: %d, Nmask: %d\n", N, Nmask);
-  }
-
-  if (!Quiet) fprintf (stderr, "z = ");
+    if (VERBOSE) gprint (GP_ERR, "pass: %d, Nmask: %d\n", N, Nmask);
+  }
+
+  if (!Quiet) gprint (GP_ERR, "z = ");
   for (N = ny = 0; ny < nterm; ny++) {
     for (nx = 0; nx < nterm - ny; nx++, N++) {
       sprintf (name, "CX%dY%d", nx, ny);
       set_variable (name, b[N][0]);
-      if (!Quiet) fprintf (stderr, "%f x^%d y^%d  ", b[N][0], nx, ny);
+      if (!Quiet) gprint (GP_ERR, "%f x^%d y^%d  ", b[N][0], nx, ny);
     }
   }
@@ -259 +259 @@
   set_variable (name, (double) order);
   
-  if (!Quiet) fprintf (stderr, "\n");
+  if (!Quiet) gprint (GP_ERR, "\n");
 
   /* free internal data */