assembleWellEq

This commit is contained in:
steink 2023-03-20 13:15:59 +01:00
parent af30df5dff
commit c997b438dd

View File

@ -243,7 +243,7 @@ assemblePressureEq(const int seg,
// contribution from the outlet segment
eqns.residual()[seg][SPres] -= outlet_pressure.value();
for (int pv_idx = 0; pv_idx < numWellEq; ++pv_idx) {
eqns.D()[seg][outlet_segment_index][SPres][pv_idx] = -outlet_pressure.derivative(pv_idx + Indices::numEq);
eqns.D()[seg][outlet_segment_index][SPres][pv_idx] -= outlet_pressure.derivative(pv_idx + Indices::numEq);
}
}