fixed add-well-contributions for multisegement wells

This commit is contained in:
hnil
2023-09-11 20:01:29 +02:00
parent 87b9dd3ec2
commit 29e968c232

View File

@@ -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