Changeset 16059 for trunk/Ohana/src/opihi/lib.data/mrqmin.c

Timestamp:

Jan 14, 2008, 1:50:11 PM (19 years ago)

Author:

eugene

Message:

convert to new gauss-jordan code (non-Press); fixing Wall,Werror warnings; fixed bug in shell command; fixed bug in matrix element lookup

File:

• trunk/Ohana/src/opihi/lib.data/mrqmin.c (modified) (2 diffs)

trunk/Ohana/src/opihi/lib.data/mrqmin.c

@@ -60 +60 @@
   }
 
-  fgaussj (talpha, Npar, tbeta, 1);
+  fgaussjordan (talpha, tbeta, Npar, 1);
 
   for (j = 0; j < Npar; j++) partry[j] = par[j] - tbeta[j][0];
@@ -137 +137 @@
 /* don't invoke this in the middle of a run, only near the end */ 
 float **mrqcovar (int Npar) {
-  fgaussj (alpha, Npar, beta, 1);
+  fgaussjordan (alpha, beta, Npar, 1);
   return (alpha);
 }