107 cusp::array1d<ValueType, cusp::host_memory> resid(
numRHS);
108 std::cout <<
"Successfully converged after " <<
iteration_count() <<
" iterations to tolerance " <<
tolerance(0) << std::endl;
109 std::cout <<
"with max residual norm ";
111 for (
int i = 0; i <
numRHS; i++) {
112 resid[i] = cusp::blas::nrm2(r.column(i));
113 if (resid[i] > norm_max) norm_max = resid[i];
115 std::cout << norm_max << std::endl;
125 cusp::array1d<ValueType, cusp::host_memory> resid(
numRHS);
126 std::cout <<
"Failed to converge after " <<
iteration_count() <<
" iterations." << std::endl;
127 std::cout <<
"with max residual norm ";
129 for (
int i = 0; i <
numRHS; i++) {
130 resid[i] = cusp::blas::nrm2(r.column(i));
131 if (resid[i] > norm_max) norm_max = resid[i];
133 std::cout << norm_max << std::endl;