Expand concentration interval in which we look for solution.

This is an ad-hoc fix. We should investigate if (or how much)
the concentration may be above from the injected value.
This commit is contained in:
Atgeirr Flø Rasmussen 2012-04-10 14:08:32 +02:00
parent dce3068343
commit d71742871b

View File

@ -570,7 +570,7 @@ namespace Opm
{ {
ResidualC res(*this, cell); ResidualC res(*this, cell);
const double a = 0.0; const double a = 0.0;
const double b = polyprops_.cMax(); const double b = polyprops_.cMax()*1.1; // Add 10% to account for possible non-monotonicity of hyperbolic system.
int iters_used; int iters_used;
// Check if current state is an acceptable solution. // Check if current state is an acceptable solution.