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: MultisegmentWellEquations::residual()

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

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

@ -110,6 +110,12 @@ public:
const int seg_pressure_var_ind, const int seg_pressure_var_ind,
const WellState& well_state) const; const WellState& well_state) const;
//! \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_;

10
opm/simulators/wells/MultisegmentWellEval.cpp
View File

@ -113,7 +113,7 @@ getWellConvergence(const WellState& well_state,
std::vector<std::vector<double>> abs_residual(this->numberOfSegments(), std::vector<double>(numWellEq, 0.0)); std::vector<std::vector<double>> abs_residual(this->numberOfSegments(), std::vector<double>(numWellEq, 0.0));
for (int seg = 0; seg < this->numberOfSegments(); ++seg) { for (int seg = 0; seg < this->numberOfSegments(); ++seg) {
for (int eq_idx = 0; eq_idx < numWellEq; ++eq_idx) { for (int eq_idx = 0; eq_idx < numWellEq; ++eq_idx) {
abs_residual[seg][eq_idx] = std::abs(linSys_.resWell_[seg][eq_idx]); abs_residual[seg][eq_idx] = std::abs(linSys_.residual()[seg][eq_idx]);
} }
} }
@ -178,7 +178,7 @@ getWellConvergence(const WellState& well_state,
tolerance_wells, tolerance_wells,
tolerance_wells, tolerance_wells,
max_residual_allowed}, max_residual_allowed},
std::abs(linSys_.resWell_[0][SPres]), std::abs(linSys_.residual()[0][SPres]),
report, report,
deferred_logger); deferred_logger);
@ -1484,10 +1484,10 @@ getFiniteWellResiduals(const std::vector<Scalar>& B_avg,
for (int eq_idx = 0; eq_idx < numWellEq; ++eq_idx) { for (int eq_idx = 0; eq_idx < numWellEq; ++eq_idx) {
double residual = 0.; double residual = 0.;
if (eq_idx < baseif_.numComponents()) { if (eq_idx < baseif_.numComponents()) {
residual = std::abs(linSys_.resWell_[seg][eq_idx]) * B_avg[eq_idx]; residual = std::abs(linSys_.residual()[seg][eq_idx]) * B_avg[eq_idx];
} else { } else {
if (seg > 0) { if (seg > 0) {
residual = std::abs(linSys_.resWell_[seg][eq_idx]); residual = std::abs(linSys_.residual()[seg][eq_idx]);
} }
} }
if (std::isnan(residual) || std::isinf(residual)) { if (std::isnan(residual) || std::isinf(residual)) {
@ -1504,7 +1504,7 @@ getFiniteWellResiduals(const std::vector<Scalar>& B_avg,
// handling the control equation residual // handling the control equation residual
{ {
const double control_residual = std::abs(linSys_.resWell_[0][numWellEq - 1]); const double control_residual = std::abs(linSys_.residual()[0][numWellEq - 1]);
if (std::isnan(control_residual) || std::isinf(control_residual)) { if (std::isnan(control_residual) || std::isinf(control_residual)) {
deferred_logger.debug("nan or inf value for control residal get for well " + baseif_.name()); deferred_logger.debug("nan or inf value for control residal get for well " + baseif_.name());
return {false, residuals}; return {false, residuals};