Changeset 13361
- Timestamp:
- 09/13/12 13:37:43 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/modules/Solverx/SolverxSeq.cpp
r13334 r13361 73 73 #ifdef _HAVE_GSL_ 74 74 // the matrix will be modified by LU decomposition. Use gsl_A copy 75 for (int i=0; i<m*m; ++i)76 std::cout << "EDF_fos_forward_for_solverx A["<< i << "]=" << inVal[i] << std::endl;77 for (int i=0; i<m; ++i)78 std::cout << "EDF_fos_forward_for_solverx b["<< i << "]=" << inVal[i+m*m] << std::endl;79 75 double* Acopy = xNew<double>(m*m); 80 76 xMemCpy(Acopy,inVal,m*m); … … 187 183 gsl_permutation *p = NULL; 188 184 /*A will be modified by LU decomposition. Use copy*/ 189 for (int i=0; i<n*n; ++i)190 std::cout << "SolverxSeq A["<< i << "]=" << A[i] << std::endl;191 for (int i=0; i<n; ++i)192 std::cout << "SolverxSeq b["<< i << "]=" << B[i] << std::endl;193 185 double* Acopy = xNew<double>(n*n); 194 186 xMemCpy(Acopy,A,n*n);
Note:
See TracChangeset
for help on using the changeset viewer.