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

added: StandardWellEquations::residual()

Browse Source 
this returns a const reference to the residual vector.
use this in StandardWellEval::getWellConvergence
This commit is contained in:
Arne Morten Kvarving 2022-11-18 12:59:16 +01:00
parent 5bab68cb1f
commit 5eac7cae4e
2 changed files with 8 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
opm/simulators/wells/StandardWellEquations.hpp
View File

@ -117,6 +117,12 @@ public:
//! \brief Sum with off-process contribution. //! \brief Sum with off-process contribution.
void sumDistributed(Parallel::Communication comm); void sumDistributed(Parallel::Communication comm);
//! \brief Returns a const reference to the residual.
const BVectorWell& residual() const
{
return resWell_;
}
// two off-diagonal matrices // two off-diagonal matrices
OffDiagMatWell duneB_; OffDiagMatWell duneB_;
OffDiagMatWell duneC_; OffDiagMatWell duneC_;

4
opm/simulators/wells/StandardWellEval.cpp
View File

@ -679,7 +679,7 @@ getWellConvergence(const WellState& well_state,
res.resize(numWellEq_); res.resize(numWellEq_);
for (int eq_idx = 0; eq_idx < numWellEq_; ++eq_idx) { for (int eq_idx = 0; eq_idx < numWellEq_; ++eq_idx) {
// magnitude of the residual matters // magnitude of the residual matters
res[eq_idx] = std::abs(this->linSys_.resWell_[0][eq_idx]); res[eq_idx] = std::abs(this->linSys_.residual()[0][eq_idx]);
} }
std::vector<double> well_flux_residual(baseif_.numComponents()); std::vector<double> well_flux_residual(baseif_.numComponents());
@ -716,7 +716,7 @@ getWellConvergence(const WellState& well_state,
WellConvergence(baseif_). WellConvergence(baseif_).
checkConvergenceControlEq(well_state, checkConvergenceControlEq(well_state,
{1.e3, 1.e4, 1.e-4, 1.e-6, maxResidualAllowed}, {1.e3, 1.e4, 1.e-4, 1.e-6, maxResidualAllowed},
std::abs(this->linSys_.resWell_[0][Bhp]), std::abs(this->linSys_.residual()[0][Bhp]),
report, report,
deferred_logger); deferred_logger);