Changeset 13320

Timestamp:

09/11/12 11:30:54 (13 years ago)

Author:

utke

Message:

CHG add in temporary debug output

File:

• issm/trunk-jpl/src/c/modules/Solverx/SolverxSeq.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/modules/Solverx/SolverxSeq.cpp

@@ -75 +75 @@
 #ifdef _HAVE_GSL_
   // the matrix will be modified by LU decomposition. Use gsl_A copy
+  for (int i=0; i<m*m; ++i)
+    std::cout << "EDF_fos_forward_for_solverx A["<< i << "]=" << inVal[i] << std::endl;
+  for (int i=0; i<m; ++i)
+    std::cout << "EDF_fos_forward_for_solverx b["<< i << "]=" << inVal[i+m*m] << std::endl;
   double* Acopy = xNew<double>(m*m);
   xMemCpy(Acopy,inVal,m*m);
@@ -189 +193 @@
         gsl_permutation *p = NULL;
         /*A will be modified by LU decomposition. Use copy*/
+        for (int i=0; i<n*n; ++i)
+          std::cout << "SolverxSeq A["<< i << "]=" << A[i] << std::endl;
+        for (int i=0; i<n; ++i)
+          std::cout << "SolverxSeq b["<< i << "]=" << B[i] << std::endl;
         double* Acopy = xNew<double>(n*n);
         xMemCpy(Acopy,A,n*n);