OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Improve solver output in Newton iterations.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2012-05-21 14:08:38 +02:00
parent bccb7381fd
commit 182c4718b0
1 changed files with 6 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
opm/core/pressure/CompressibleTpfa.cpp
View File

@ -122,10 +122,9 @@ namespace Opm
double inc_norm = 0.0; double inc_norm = 0.0;
int iter = 0; int iter = 0;
double res_norm = residualNorm(); double res_norm = residualNorm();
const int width = 18; std::cout << "\nIteration Residual Change in p\n"
std::cout << "\nIteration Residual Change" << std::setw(9) << iter
<< std::setw(width) << iter << std::setw(18) << res_norm << std::endl;
<< std::setw(width) << res_norm << std::endl;
for (; (iter < maxiter_) && (res_norm > residual_tol_); ++iter) { for (; (iter < maxiter_) && (res_norm > residual_tol_); ++iter) {
// Solve for increment in Newton method: // Solve for increment in Newton method:
// incr = x_{n+1} - x_{n} = -J^{-1}F // incr = x_{n+1} - x_{n} = -J^{-1}F
@ -151,10 +150,9 @@ namespace Opm
// Update residual norm. // Update residual norm.
res_norm = residualNorm(); res_norm = residualNorm();
std::cout << "\nIteration Residual Change" std::cout << std::setw(9) << iter + 1
<< std::setw(width) << iter + 1 << std::setw(18) << res_norm
<< std::setw(width) << res_norm << std::setw(18) << inc_norm << std::endl;
<< std::setw(width) << inc_norm << std::endl;
} }
if ((iter == maxiter_) && (res_norm > residual_tol_) && (inc_norm > change_tol_)) { if ((iter == maxiter_) && (res_norm > residual_tol_) && (inc_norm > change_tol_)) {