diff --git a/opm/core/simulator/EquilibrationHelpers.hpp b/opm/core/simulator/EquilibrationHelpers.hpp index 2ea49e59e..8619cf640 100644 --- a/opm/core/simulator/EquilibrationHelpers.hpp +++ b/opm/core/simulator/EquilibrationHelpers.hpp @@ -350,7 +350,7 @@ namespace Opm const double temp, const double sat_oil = 0.0 ) const { - if (sat_oil > 0.0) { + if (std::abs(sat_oil) > 1e-16) { return satRv(press, temp); } else { return std::min(satRv(press, temp), linearInterpolation(depth_, rv_, depth));