OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-22 23:43:28 -06:00
Code Issues Packages Projects Releases Wiki Activity

Add comments at the spots where we multiply with D^(-1) - we actually can do this on all processes

Browse Source
This commit is contained in:
Lisa Julia Nebel 2024-10-08 15:02:27 +02:00
parent 627b9c98ba
commit 4801a17703
1 changed files with 15 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
opm/simulators/wells/MultisegmentWellEquations.cpp
View File

@ -157,6 +157,9 @@ apply(const BVector& x, BVector& Ax) const
this->pw_info_.communication().sum(Bx.data(), Bx.size()); this->pw_info_.communication().sum(Bx.data(), Bx.size());
} }
// It is ok to do this on each process instead of only on one,
// because the other processes would remain idle while waiting for
// the single process to complete the computation.
// invDBx = duneD^-1 * Bx_ // invDBx = duneD^-1 * Bx_
const BVectorWell invDBx = mswellhelpers::applyUMFPack(*duneDSolver_, Bx); const BVectorWell invDBx = mswellhelpers::applyUMFPack(*duneDSolver_, Bx);
@ -168,6 +171,9 @@ template<class Scalar, int numWellEq, int numEq>
void MultisegmentWellEquations<Scalar,numWellEq,numEq>:: void MultisegmentWellEquations<Scalar,numWellEq,numEq>::
apply(BVector& r) const apply(BVector& r) const
{ {
// It is ok to do this on each process instead of only on one,
// because the other processes would remain idle while waiting for
// the single process to complete the computation.
// invDrw_ = duneD^-1 * resWell_ // invDrw_ = duneD^-1 * resWell_
const BVectorWell invDrw = mswellhelpers::applyUMFPack(*duneDSolver_, resWell_); const BVectorWell invDrw = mswellhelpers::applyUMFPack(*duneDSolver_, resWell_);
// r = r - duneC_^T * invDrw // r = r - duneC_^T * invDrw
@ -198,6 +204,9 @@ template<class Scalar, int numWellEq, int numEq>
typename MultisegmentWellEquations<Scalar,numWellEq,numEq>::BVectorWell typename MultisegmentWellEquations<Scalar,numWellEq,numEq>::BVectorWell
MultisegmentWellEquations<Scalar,numWellEq,numEq>::solve() const MultisegmentWellEquations<Scalar,numWellEq,numEq>::solve() const
{ {
// It is ok to do this on each process instead of only on one,
// because the other processes would remain idle while waiting for
// the single process to complete the computation.
return mswellhelpers::applyUMFPack(*duneDSolver_, resWell_); return mswellhelpers::applyUMFPack(*duneDSolver_, resWell_);
} }
@ -205,6 +214,9 @@ template<class Scalar, int numWellEq, int numEq>
typename MultisegmentWellEquations<Scalar,numWellEq,numEq>::BVectorWell typename MultisegmentWellEquations<Scalar,numWellEq,numEq>::BVectorWell
MultisegmentWellEquations<Scalar,numWellEq,numEq>::solve(const BVectorWell& rhs) const MultisegmentWellEquations<Scalar,numWellEq,numEq>::solve(const BVectorWell& rhs) const
{ {
// It is ok to do this on each process instead of only on one,
// because the other processes would remain idle while waiting for
// the single process to complete the computation.
return mswellhelpers::applyUMFPack(*duneDSolver_, rhs); return mswellhelpers::applyUMFPack(*duneDSolver_, rhs);
} }
@ -227,6 +239,9 @@ recoverSolutionWell(const BVector& x, BVectorWell& xw) const
} }
// xw = D^-1 * resWell // xw = D^-1 * resWell
// It is ok to do this on each process instead of only on one,
// because the other processes would remain idle while waiting for
// the single process to complete the computation.
xw = mswellhelpers::applyUMFPack(*duneDSolver_, resWell); xw = mswellhelpers::applyUMFPack(*duneDSolver_, resWell);
} }