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

checking magnitude of the residual for convergence

Browse Source 
in StandardWell_impl.
This commit is contained in:
Kai Bao 2017-07-19 16:10:50 +02:00
parent ddb21e1f2c
commit 168205a993
2 changed files with 3 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
opm/autodiff/StandardWell_impl.hpp
View File

@ -1722,7 +1722,8 @@ namespace Opm
std::vector<Scalar> res(numComp);
for (int comp = 0; comp < numWellEq; ++comp) {
res[comp] = resWell_[0][comp];
// magnitude of the residual matters
res[comp] = std::abs(resWell_[0][comp]);
}
Vector well_flux_residual(numComp);

4
opm/autodiff/StandardWellsDense_impl.hpp
View File

@ -237,8 +237,6 @@ namespace Opm {
if (param_.solve_welleq_initially_ && iterationIdx == 0) {
// solve the well equations as a pre-processing step
report = solveWellEq(ebosSimulator, dt, well_state);
std::cout << " solvWellEq is done ! " << std::endl;
std::cin.ignore();
}
assembleWellEq(ebosSimulator, dt, well_state, false);
@ -989,8 +987,8 @@ namespace Opm {
for (const auto& well : well_container_) {
if ( !well->getWellConvergence(ebosSimulator, B_avg, param_) ) {
converged_well = false;
// break; // TODO: no need to check other wells?
}
break;
}
// TODO: to think about the output here.