Index: /branches/eam_branches/ipp-20130711/Ohana/src/libohana/src/gaussj.c
===================================================================
--- /branches/eam_branches/ipp-20130711/Ohana/src/libohana/src/gaussj.c	(revision 36064)
+++ /branches/eam_branches/ipp-20130711/Ohana/src/libohana/src/gaussj.c	(revision 36065)
@@ -1,5 +1,5 @@
 # include <ohana.h>
 # define GROWTHTEST 0
-# define MAX_RANGE 1.0e7
+# define MAX_RANGE 1.0e9
 
 // Gauss-Jordan elimination using full pivots based on Press et al's description.  Substantially
@@ -91,4 +91,5 @@
     /* rescale by pivot reciprocal */
     double tmpval = 1.0 / A[maxcol][maxcol];
+    // XXX why is this here (don't I double count this element?) A[maxcol][maxcol] = 1.0;
     A[maxcol][maxcol] = 1.0;
     for (col = 0; col < N; col++) A[maxcol][col] *= tmpval;
@@ -117,4 +118,15 @@
     }
   }
+
+# if (GROWTHTEST)
+    fprintf (stderr, "final A matrix:\n");
+    for (row = 0; row < N; row++) {
+	for (col = 0; col < N; col++) {
+	    fprintf (stderr, "%10.3e ", A[row][col]);
+	}
+	fprintf (stderr, "\n");
+    }
+    fprintf (stderr, "\n");
+# endif
 
   // swap back the inverse matrix based on the row swaps above