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

Add auxiliary module in BlackoilModelEbos constructor.

Browse Source 
Previously, we did this during assemble. Unfortunately, this resulted
in the sparsity pattern being recomputed every time instead of just
once for a non-linear solve.
This commit is contained in:
Markus Blatt 2017-05-15 20:33:36 +02:00
parent 799cbb4b62
commit 07f61c2a2b
1 changed files with 9 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
opm/autodiff/BlackoilModelEbos.hpp
View File

@ -156,6 +156,15 @@ namespace Opm {
{
OPM_THROW(std::logic_error,"solver down cast to ISTLSolver failed");
}
if ( param_.matrix_add_well_contributions_ )
{
// This might be dangerous?!
ebosSimulator_.model().clearAuxiliaryModules();
auto auxMod = std::make_shared<WellConnectionAuxiliaryModule<TypeTag> >( wellModel().wells() );
ebosSimulator_.model().addAuxiliaryModule(auxMod);
}
}
bool isParallel() const
@ -354,15 +363,6 @@ namespace Opm {
SimulatorReport assemble(const SimulatorTimerInterface& timer,
const int iterationIdx)
{
if ( param_.matrix_add_well_contributions_ )
{
// This might be dangerous?!
ebosSimulator_.model().clearAuxiliaryModules();
auto auxMod = std::make_shared<WellConnectionAuxiliaryModule<TypeTag> >( wellModel().wells() );
ebosSimulator_.model().addAuxiliaryModule(auxMod);
}
// -------- Mass balance equations --------
ebosSimulator_.model().newtonMethod().setIterationIndex(iterationIdx);
ebosSimulator_.problem().beginIteration();