Mark the spots where we multiply with D^(-1) - do this only on one pr…

…ocess later

opm/simulators/wells/MultisegmentWellEquations.cpp

@@ -156,6 +156,7 @@ apply(const BVector& x, BVector& Ax) const
     this->pw_info_.communication().sum(Bx.data(), Bx.size());
 
     // invDBx = duneD^-1 * Bx_
+    // TODO: do this only on one process and distribute the solution
     const BVectorWell invDBx = mswellhelpers::applyUMFPack(*duneDSolver_, Bx);
 
     // Ax = Ax - duneC_^T * invDBx
@@ -167,6 +168,7 @@ void MultisegmentWellEquations<Scalar,numWellEq,numEq>::
 apply(BVector& r) const
 {
     // invDrw_ = duneD^-1 * resWell_
+    // TODO: do this only on one process and distribute the solution
     const BVectorWell invDrw = mswellhelpers::applyUMFPack(*duneDSolver_, resWell_);
     // r = r - duneC_^T * invDrw
     duneC_.mmtv(invDrw, r);
@@ -196,13 +198,15 @@ template<class Scalar, int numWellEq, int numEq>
 typename MultisegmentWellEquations<Scalar,numWellEq,numEq>::BVectorWell
 MultisegmentWellEquations<Scalar,numWellEq,numEq>::solve() const
 {
+    // TODO: do this only on one process and distribute the solution
     return mswellhelpers::applyUMFPack(*duneDSolver_, resWell_);
 }
 
 template<class Scalar, int numWellEq, int numEq>
 typename MultisegmentWellEquations<Scalar,numWellEq,numEq>::BVectorWell
 MultisegmentWellEquations<Scalar,numWellEq,numEq>::solve(const BVectorWell& rhs) const
 {
+    // TODO: do this only on one process and distribute the solution
     return mswellhelpers::applyUMFPack(*duneDSolver_, rhs);
 }
 
@@ -220,6 +224,7 @@ recoverSolutionWell(const BVector& x, BVectorWell& xw) const
     resWell -= Bx;
 
     // xw = D^-1 * resWell
+    // TODO: do this only on one process and distribute the solution
     xw = mswellhelpers::applyUMFPack(*duneDSolver_, resWell);
 }