Changeset 33327

Timestamp:

Feb 21, 2012, 4:52:51 PM (14 years ago)

Author:

watersc1

Message:

Switching from LUSolve to GJSolve. See commit log on the trunk for details.

Location:

tags/ipp-20120216/psModules/src/detrend

Files:

• . (modified) (1 prop)
• pmPattern.c (modified) (1 diff)

tags/ipp-20120216/psModules/src/detrend/pmPattern.c

@@ -803 +803 @@
     
     // Solve the Ax=b equation
-    psMatrixLUSolve(XX,solution);
-
+    //    psMatrixLUSolve(XX,solution);
+    psMatrixGJSolve(XX,solution);
+#if (0)
+    for (int i = 0; i < numCells; i++) { // print vector b
+      psTrace("psModules.detrend.cont",5,"b: %d %f",
+              i,
+              solution->data.F64[i]
+              );
+    }
+#endif
+    
     /* old code to remove the minimum solution value from the set, to give a "minimal set of offsets." Mathematically unnecessary. */
 /*     double min = 99e99; */