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

Removes rest of the commented out code introduced with applyUMFPack.

Browse Source
This commit is contained in:
Markus Blatt
2020-08-20 10:13:53 +02:00
parent cca38ddc36
commit 49ed3d6468
1 changed files with 0 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@@ -640,7 +640,6 @@ namespace Opm
apply(BVector& r) const apply(BVector& r) const
{ {
// invDrw_ = duneD^-1 * resWell_ // invDrw_ = duneD^-1 * resWell_
// const BVectorWell invDrw = mswellhelpers::invDXDirect(duneD_, resWell_);
const BVectorWell invDrw = mswellhelpers::applyUMFPack(duneD_, duneDSolver_, resWell_); const BVectorWell invDrw = mswellhelpers::applyUMFPack(duneD_, duneDSolver_, resWell_);
// r = r - duneC_^T * invDrw // r = r - duneC_^T * invDrw
duneC_.mmtv(invDrw, r); duneC_.mmtv(invDrw, r);
@@ -1000,7 +999,6 @@ namespace Opm
// resWell = resWell - B * x // resWell = resWell - B * x
duneB_.mmv(x, resWell); duneB_.mmv(x, resWell);
// xw = D^-1 * resWell // xw = D^-1 * resWell
//xw = mswellhelpers::invDXDirect(duneD_, resWell);
xw = mswellhelpers::applyUMFPack(duneD_, duneDSolver_, resWell); xw = mswellhelpers::applyUMFPack(duneD_, duneDSolver_, resWell);
} }
@@ -2398,7 +2396,6 @@ namespace Opm
assembleWellEqWithoutIteration(ebosSimulator, dt, inj_controls, prod_controls, well_state, deferred_logger); assembleWellEqWithoutIteration(ebosSimulator, dt, inj_controls, prod_controls, well_state, deferred_logger);
// const BVectorWell dx_well = mswellhelpers::invDXDirect(duneD_, resWell_);
const BVectorWell dx_well = mswellhelpers::applyUMFPack(duneD_, duneDSolver_, resWell_); const BVectorWell dx_well = mswellhelpers::applyUMFPack(duneD_, duneDSolver_, resWell_);
if (it > param_.strict_inner_iter_ms_wells_) if (it > param_.strict_inner_iter_ms_wells_)