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

Fix reduction calculation

Browse Source
This commit is contained in:
Tong Dong Qiu 2022-11-09 10:27:10 +01:00
parent 1dd6af15dc
commit b51a4a863e
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
opm/simulators/linalg/bda/rocalutionSolverBackend.cpp
View File

@ -175,6 +175,8 @@ SolverStatus rocalutionSolverBackend<block_size>::solve_system(std::shared_ptr<B
// so it just calls ILU::Build() everytime
roc_solver->ReBuildNumeric();
double norm_0 = roc_rhs.Norm(); // since the initial guess is a vector with 0s, initial error is norm(b)
// actually solve
Dune::Timer t_solve;
roc_solver->Solve(roc_rhs, &roc_x);
@ -188,7 +190,7 @@ SolverStatus rocalutionSolverBackend<block_size>::solve_system(std::shared_ptr<B
res.elapsed = t_solve.stop();
res.iterations = roc_solver->GetIterationCount();
res.reduction = roc_solver->GetCurrentResidual();
res.reduction = roc_solver->GetCurrentResidual() / norm_0;
res.conv_rate = static_cast<double>(pow(res.reduction, 1.0 / res.iterations));
res.converged = (roc_solver->GetSolverStatus() == 2);