fixing a typo bug in addWellControlEq

verified with a case with a simple horizontal multi-segment well.
This commit is contained in:
Kai Bao 2015-10-14 13:14:29 +02:00
parent 50938effc4
commit 92c5927f30

View File

@ -1196,7 +1196,7 @@ namespace Opm {
for (int w = 0; w < nw; ++w) {
WellMultiSegmentConstPtr well = wellsMultiSegment()[w];
const int nseg = well->numberOfSegments();
ADB segp = subset(state.segqs, Span(nseg, 1, start_segment));
ADB segp = subset(state.segp, Span(nseg, 1, start_segment));
ADB well_residual = segp - well->wellOps().s2s_outlet * segp;
// - subset(well_segment_pressures_delta_, Span(nseg, 1, start_segment));
ADB others_well_residual = subset(well_residual, Span(nseg - 1, 1, 1));