Changeset 36065

Timestamp:

Aug 30, 2013, 4:37:29 PM (13 years ago)

Author:

eugene

Message:

some verbosity on dgaussjordan for testing

File:

• branches/eam_branches/ipp-20130711/Ohana/src/libohana/src/gaussj.c (modified) (3 diffs)

branches/eam_branches/ipp-20130711/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