Index: /issm/trunk-jpl/src/c/modules/Solverx/SolverxSeq.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/Solverx/SolverxSeq.cpp	(revision 13319)
+++ /issm/trunk-jpl/src/c/modules/Solverx/SolverxSeq.cpp	(revision 13320)
@@ -75,4 +75,8 @@
 #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,4 +193,8 @@
 	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);