Index: trunk/psModules/src/detrend/pmPattern.c
===================================================================
--- trunk/psModules/src/detrend/pmPattern.c	(revision 33324)
+++ trunk/psModules/src/detrend/pmPattern.c	(revision 33326)
@@ -803,6 +803,15 @@
     
     // 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; */