Index: /trunk/Ohana/src/gastro/Makefile
===================================================================
--- /trunk/Ohana/src/gastro/Makefile	(revision 8298)
+++ /trunk/Ohana/src/gastro/Makefile	(revision 8299)
@@ -37,5 +37,4 @@
 $(SRC)/sort.$(ARCH).o \
 $(SRC)/gfitpoly.$(ARCH).o \
-$(SRC)/gaussj.$(ARCH).o \
 $(SRC)/get_region_coords.$(ARCH).o \
 $(SRC)/greference.$(ARCH).o \
Index: unk/Ohana/src/gastro/src/gaussj.c
===================================================================
--- /trunk/Ohana/src/gastro/src/gaussj.c	(revision 8298)
+++ 	(revision )
@@ -1,77 +1,0 @@
-# include <ohana.h>
-
-int gaussj (double **a, int n, double **b, int m) {
-
-  int *indxc,*indxr,*ipiv;
-  int i, icol, irow, j, k, l, ll;
-  double big,dum,pivinv;
-  
-  irow = icol = 0;
-
-  ALLOCATE (indxc, int, n);
-  ALLOCATE (indxr, int, n);
-  ALLOCATE (ipiv, int, n);
-  for (j = 0; j < n; j++) 
-    ipiv[j] = 0;
-
-  for (i = 0; i < n; i++) {
-    big = 0.0;
-    for (j = 0; j < n; j++) {
-      if (ipiv[j] != 1) {
-	for (k = 0; k < n; k++) {
-	  if (ipiv[k] == 0) {
-	    if (fabs (a[j][k]) >= big) {
-	      big  = fabs (a[j][k]);
-	      irow = j;
-	      icol = k;
-	    }
-	  } 
-	  else 
-	    if (ipiv[k] > 1) {
-	      fprintf (stderr, "GAUSSJ: Singular Matrix! (1)\n");
-	      return (0);
-	    }
-	}
-      }
-    }
-    ipiv[icol]++;
-    if (irow != icol) {
-      for (l = 0; l < n; l++) 
-	SWAP (a[irow][l], a[icol][l]);
-      for (l = 0; l < m; l++) 
-	SWAP (b[irow][l], b[icol][l]);
-    }
-    indxr[i] = irow;
-    indxc[i] = icol;
-    if (a[icol][icol] == 0.0) {
-      fprintf (stderr, "GAUSSJ: Singular Matrix! (2)\n");
-      return (0);
-    }
-    pivinv = 1.0 / a[icol][icol];
-    a[icol][icol] = 1.0;
-    for (l = 0; l < n; l++) 
-      a[icol][l] *= pivinv;
-    for (l = 0; l < m; l++) 
-      b[icol][l] *= pivinv;
-    for (ll = 0; ll < n; ll++) {
-      if (ll != icol) {
-	dum = a[ll][icol];
-	a[ll][icol] = 0.0;
-	for (l = 0; l < n; l++) 
-	  a[ll][l] -= a[icol][l]*dum;
-	for (l = 0; l < m; l++) 
-	  b[ll][l] -= b[icol][l]*dum;
-      }
-    }
-  }
-
-  for (l = n - 1; l >= 0; l--) {
-    if (indxr[l] != indxc[l])
-      for (k = 0; k < n; k++)
-	SWAP (a[k][indxr[l]], a[k][indxc[l]]);
-  }
-  free (ipiv);
-  free (indxr);
-  free (indxc);
-  return (1);
-}
Index: /trunk/Ohana/src/gastro/src/gfit.c
===================================================================
--- /trunk/Ohana/src/gastro/src/gfit.c	(revision 8298)
+++ /trunk/Ohana/src/gastro/src/gfit.c	(revision 8299)
@@ -307,5 +307,5 @@
       vector[2][1] = Yy;
 
-      gaussj (matrix, NR, vector, NC); 
+      dgaussj (matrix, NR, vector, NC); 
       
       /* 
Index: /trunk/Ohana/src/gastro/src/gfitpoly.c
===================================================================
--- /trunk/Ohana/src/gastro/src/gfitpoly.c	(revision 8298)
+++ /trunk/Ohana/src/gastro/src/gfitpoly.c	(revision 8299)
@@ -212,5 +212,5 @@
       }
       /* svd (matrix, NPARS, vector, 2);  */
-      gaussj (matrix, NPARS, vector, 2); 
+      dgaussj (matrix, NPARS, vector, 2); 
       i = 0;
       for (m = 0; m < NPOWR; m++) {