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::solve()

Browse Source 
this applies the inverted D matrix to the residual vector.
use the new method in the well implementation.
This commit is contained in:
Arne Morten Kvarving 2022-11-11 21:41:24 +01:00
parent 83c41e3f9c
commit daf0f90fe5
3 changed files with 10 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -160,6 +160,12 @@ void StandardWellEquations<Scalar,numEq>::invert()
}
}
template<class Scalar, int numEq>
void StandardWellEquations<Scalar,numEq>::solve(BVectorWell& dx_well) const
{
invDuneD_.mv(resWell_, dx_well);
}
#define INSTANCE(N) \
template class StandardWellEquations<double,N>;

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

@ -79,6 +79,9 @@ public:
//! \brief Apply linear operator to vector.
void apply(BVector& r) const;
//! \brief Apply inverted D matrix to residual and store in vector.
void solve(BVectorWell& dx_well) const;
//! \brief Invert D matrix.
void invert();

2
opm/simulators/wells/StandardWell_impl.hpp
View File

@ -1679,7 +1679,7 @@ namespace Opm
// which is why we do not put the assembleWellEq here.
BVectorWell dx_well(1);
dx_well[0].resize(this->numWellEq_);
this->linSys_.invDuneD_.mv(this->linSys_.resWell_, dx_well);
this->linSys_.solve( dx_well);
updateWellState(dx_well, well_state, deferred_logger);
}