Changeset 16058

Timestamp:

Jan 14, 2008, 1:48:06 PM (18 years ago)

Author:

eugene

Message:

convert to new gauss-jordan code (non-Press)

Location:

trunk/Ohana/src

Files:

• gastro/src/gfit.c (modified) (1 diff)
• gastro/src/gfitpoly.c (modified) (1 diff)
• gastro2/src/lumfunc.c (modified) (1 diff)
• gastro2/src/polyfit.c (modified) (2 diffs)
• mosastro/src/GetGradients.c (modified) (1 diff)
• mosastro/src/fitpoly.c (modified) (1 diff)

trunk/Ohana/src/gastro/src/gfit.c

@@ -307 +307 @@
       vector[2][1] = Yy;
 
-      dgaussj (matrix, NR, vector, NC); 
+      dgaussjordan (matrix, vector, NR, NC); 
       
       /* 

trunk/Ohana/src/gastro/src/gfitpoly.c

@@ -212 +212 @@
       }
       /* svd (matrix, NPARS, vector, 2);  */
-      dgaussj (matrix, NPARS, vector, 2); 
+      dgaussjordan (matrix, vector, NPARS, 2); 
       i = 0;
       for (m = 0; m < NPOWR; m++) {

trunk/Ohana/src/gastro2/src/lumfunc.c

@@ -107 +107 @@
     b[1][0] += y[i]*x[i];
   }
-  dgaussj (c, 2, b, 1);
+  dgaussjordan (c, b, 2, 1);
   *C0 = b[0][0];
   *C1 = b[1][0];

trunk/Ohana/src/gastro2/src/polyfit.c

@@ -101 +101 @@
     vector[i][1] /= max;
   }
-  dgaussj (matrix, NPARS, vector, 2); 
+  dgaussjordan (matrix, vector, NPARS, 2); 
   i = 0;
   for (m = 0; m < NPOWR; m++) {
@@ -271 +271 @@
       B[0][0] = -Fx;
       B[1][0] = -Fy;
-      dgaussj (A, 2, B, 1);
+      dgaussjordan (A, B, 2, 1);
       Xo += B[0][0]; 
       Yo += B[1][0];

trunk/Ohana/src/mosastro/src/GetGradients.c

@@ -90 +90 @@
 
         if (Npts < 5) continue;
-        if (!dgaussj (a, 3, b, 2)) continue;
+        if (!dgaussjordan (a, b, 3, 2)) continue;
 
         /* we only care about the slopes, not the offsets */

trunk/Ohana/src/mosastro/src/fitpoly.c

@@ -109 +109 @@
 # endif
 
-  dgaussj (matrix, NPARS, vector, 2); 
+  dgaussjordan (matrix, vector, NPARS, 2); 
 
 # if (0)