Merge pull request #4861 from hnil/fixed_msw_addwellcontricutions

fixed add-well-contributions for multisegement wells
2025-02-25 18:55:30 -06:00 · 2023-09-13 11:40:40 +02:00
parent 3673c05242 29e968c232
commit 708401ca14
1 changed files with 2 additions and 2 deletions
--- a/opm/simulators/wells/MultisegmentWellEquations.cpp
+++ b/opm/simulators/wells/MultisegmentWellEquations.cpp
@@ -266,8 +266,8 @@ template<class SparseMatrixAdapter>
 void MultisegmentWellEquations<Scalar,numWellEq,numEq>::
 extract(SparseMatrixAdapter& jacobian) const
 {
-    const auto invDuneD = mswellhelpers::invertWithUMFPack<BVectorWell>(numWellEq,
-                                                                        numEq,
+    const auto invDuneD = mswellhelpers::invertWithUMFPack<BVectorWell>(duneD_.M(),
+                                                                        numWellEq,
                                                                        *duneDSolver_);

    // We need to change matrix A as follows