Index: /trunk/Ohana/src/gastro/src/gfit.c
===================================================================
--- /trunk/Ohana/src/gastro/src/gfit.c	(revision 16057)
+++ /trunk/Ohana/src/gastro/src/gfit.c	(revision 16058)
@@ -307,5 +307,5 @@
       vector[2][1] = Yy;
 
-      dgaussj (matrix, NR, vector, NC); 
+      dgaussjordan (matrix, vector, NR, NC); 
       
       /* 
Index: /trunk/Ohana/src/gastro/src/gfitpoly.c
===================================================================
--- /trunk/Ohana/src/gastro/src/gfitpoly.c	(revision 16057)
+++ /trunk/Ohana/src/gastro/src/gfitpoly.c	(revision 16058)
@@ -212,5 +212,5 @@
       }
       /* svd (matrix, NPARS, vector, 2);  */
-      dgaussj (matrix, NPARS, vector, 2); 
+      dgaussjordan (matrix, vector, NPARS, 2); 
       i = 0;
       for (m = 0; m < NPOWR; m++) {
Index: /trunk/Ohana/src/gastro2/src/lumfunc.c
===================================================================
--- /trunk/Ohana/src/gastro2/src/lumfunc.c	(revision 16057)
+++ /trunk/Ohana/src/gastro2/src/lumfunc.c	(revision 16058)
@@ -107,5 +107,5 @@
     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];
Index: /trunk/Ohana/src/gastro2/src/polyfit.c
===================================================================
--- /trunk/Ohana/src/gastro2/src/polyfit.c	(revision 16057)
+++ /trunk/Ohana/src/gastro2/src/polyfit.c	(revision 16058)
@@ -101,5 +101,5 @@
     vector[i][1] /= max;
   }
-  dgaussj (matrix, NPARS, vector, 2); 
+  dgaussjordan (matrix, vector, NPARS, 2); 
   i = 0;
   for (m = 0; m < NPOWR; m++) {
@@ -271,5 +271,5 @@
       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];
Index: /trunk/Ohana/src/mosastro/src/GetGradients.c
===================================================================
--- /trunk/Ohana/src/mosastro/src/GetGradients.c	(revision 16057)
+++ /trunk/Ohana/src/mosastro/src/GetGradients.c	(revision 16058)
@@ -90,5 +90,5 @@
 
 	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 */
Index: /trunk/Ohana/src/mosastro/src/fitpoly.c
===================================================================
--- /trunk/Ohana/src/mosastro/src/fitpoly.c	(revision 16057)
+++ /trunk/Ohana/src/mosastro/src/fitpoly.c	(revision 16058)
@@ -109,5 +109,5 @@
 # endif
 
-  dgaussj (matrix, NPARS, vector, 2); 
+  dgaussjordan (matrix, vector, NPARS, 2); 
 
 # if (0)