Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/modules/Solverx/SolverxSeq.cpp
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/modules/Solverx/SolverxSeq.cpp	(revision 13319)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/modules/Solverx/SolverxSeq.cpp	(revision 13320)
@@ -74,6 +74,10 @@
 int EDF_fos_forward_for_solverx(int n, IssmPDouble *inVal, IssmPDouble *inDeriv, int m, IssmPDouble *outVal, IssmPDouble *outDeriv) { /*{{{*/
 #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);
   /*Initialize gsl matrices and vectors: */
@@ -188,6 +192,10 @@
 	gsl_vector      *x = NULL;
 	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);