Changeset 36680 for branches/eam_branches/ps2-tc3-20130727/Ohana/src/libohana/src

Timestamp:

Apr 21, 2014, 5:42:34 AM (12 years ago)

Author:

eugene

Message:

merge changes from trunk

Location:

branches/eam_branches/ps2-tc3-20130727

Files:

• . (modified) (1 prop)
• Ohana (modified) (1 prop)
• Ohana/src/libohana/src (modified) (1 prop)
• Ohana/src/libohana/src/gaussj.c (modified) (3 diffs)

branches/eam_branches/ps2-tc3-20130727/Ohana/src/libohana/src/gaussj.c

@@ -1 +1 @@
 # 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 +91 @@
     /* 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 +118 @@
     }
   }
+
+# 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