correct an assertion error in getWellConvergence in StandardWell

and adding comments why it is the case for polymer case.
This commit is contained in:
Kai Bao
2017-07-31 14:20:02 +02:00
parent 5ed9f4d497
commit d470edfb0b

View File

@@ -1705,7 +1705,9 @@ namespace Opm
const int np = numberOfPhases();
const int numComp = numComponents();
assert(int(B_avg.size()) == numComp);
// the following implementation assume that the polymer is always after the w-o-g phases
// For the polymer case, there is one more mass balance equations of reservoir than wells
assert((int(B_avg.size()) == numComp) || has_polymer);
const double tol_wells = param.tolerance_wells_;
const double maxResidualAllowed = param.max_residual_allowed_;