@@ -66,7 +66,7 @@ NewtonSolve(const T & x,
6666
6767 using std::isnan;
6868 if (verbose)
69- std::cout << " Target value for 1D Newton inversion:\n " << y << std::endl;
69+ Moose::out << " Target value for 1D Newton inversion:\n " << y << std::endl;
7070
7171 do
7272 {
@@ -92,10 +92,10 @@ NewtonSolve(const T & x,
9292
9393 if (verbose)
9494 {
95- std::cout << " Iteration " << iteration << std::endl;
96- std::cout << " Current solution vector: " << z << std::endl;
97- std::cout << " Current (minus) residual: " << -R << std::endl;
98- std::cout << " Current Jacobian: " << dy_dz << std::endl;
95+ Moose::out << " Iteration " << iteration << std::endl;
96+ Moose::out << " Current solution vector: " << z << std::endl;
97+ Moose::out << " Current (minus) residual: " << -R << std::endl;
98+ Moose::out << " Current Jacobian: " << dy_dz << std::endl;
9999 }
100100
101101 // Check for NaNs
@@ -198,7 +198,7 @@ NewtonSolve2D(const T & f,
198198 return ss;
199199 };
200200 if (debug)
201- std::cout << " Target values for 2D Newton inversion:\n " << targets << std::endl;
201+ Moose::out << " Target values for 2D Newton inversion:\n " << targets << std::endl;
202202
203203 using std::isnan, std::max, std::abs;
204204
@@ -226,10 +226,10 @@ NewtonSolve2D(const T & f,
226226
227227 if (debug)
228228 {
229- std::cout << " Iteration " << iteration << std::endl;
230- std::cout << " Current solution vector:\n " << u << std::endl;
231- std::cout << " Current (minus) residual:\n " << minus_R << std::endl;
232- std::cout << " Current Jacobian:\n " << J << std::endl;
229+ Moose::out << " Iteration " << iteration << std::endl;
230+ Moose::out << " Current solution vector:\n " << u << std::endl;
231+ Moose::out << " Current (minus) residual:\n " << minus_R << std::endl;
232+ Moose::out << " Current Jacobian:\n " << J << std::endl;
233233 }
234234
235235 // Do some Jacobi (rowmax) preconditioning and check for an empty row
@@ -264,7 +264,7 @@ NewtonSolve2D(const T & f,
264264 raw_J.svd (svs);
265265 raw_J2.evd (evs_real, evs_imag);
266266 if (debug)
267- std::cout << " Jacobian singular values:\n " << svs << std::endl;
267+ Moose::out << " Jacobian singular values:\n " << svs << std::endl;
268268#endif
269269
270270 if (degenerate_row == -1 )
0 commit comments