From d397fd847ba765bc8b6dea51d2e6959947f9a20b Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Atgeirr=20Fl=C3=B8=20Rasmussen?= <atgeirr@sintef.no>
Date: Mon, 2 Apr 2012 13:24:57 +0200
Subject: [PATCH] Changed single-cell solver call.  - Using [0,1] interval
 instead of [smin, smax] interval to handle compressible case.  - Using new
 version of Regula Falsi function which exploits initial guess.

---
 opm/core/transport/reorder/TransportModelTwophase.cpp | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/opm/core/transport/reorder/TransportModelTwophase.cpp b/opm/core/transport/reorder/TransportModelTwophase.cpp
index dc9cd5316..33a91905a 100644
--- a/opm/core/transport/reorder/TransportModelTwophase.cpp
+++ b/opm/core/transport/reorder/TransportModelTwophase.cpp
@@ -167,7 +167,8 @@ namespace Opm
 	//     return;
 	// }
 	int iters_used;
-	saturation_[cell] = modifiedRegulaFalsi(res, smin_[2*cell], smax_[2*cell], maxit_, tol_, iters_used);
+	// saturation_[cell] = modifiedRegulaFalsi(res, smin_[2*cell], smax_[2*cell], maxit_, tol_, iters_used);
+	saturation_[cell] = modifiedRegulaFalsi(res, saturation_[cell], 0.0, 1.0, maxit_, tol_, iters_used);
 	fractionalflow_[cell] = fracFlow(saturation_[cell], cell);
     }