From 5e8b577cd02f94cd178a7692c3e9e6b8d793099f Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Atgeirr=20Fl=C3=B8=20Rasmussen?= Date: Tue, 5 Nov 2019 11:10:43 +0100 Subject: [PATCH] Use the RegulaFalsiBisection method instead of RegulaFalsi. --- opm/simulators/wells/StandardWell_impl.hpp | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/opm/simulators/wells/StandardWell_impl.hpp b/opm/simulators/wells/StandardWell_impl.hpp index a8c47fae7..b5dbb9667 100644 --- a/opm/simulators/wells/StandardWell_impl.hpp +++ b/opm/simulators/wells/StandardWell_impl.hpp @@ -3749,7 +3749,7 @@ namespace Opm const double flo_tolerance = 1e-6 * std::fabs(flo_samples.back()); int iteration = 0; try { - const double solved_bhp = RegulaFalsi<>:: + const double solved_bhp = RegulaFalsiBisection<>:: solve(eq, low, high, max_iteration, flo_tolerance, iteration); bhp_samples.push_back(solved_bhp); } @@ -3821,7 +3821,7 @@ namespace Opm return boost::optional(); } try { - const double solved_bhp = RegulaFalsi<>:: + const double solved_bhp = RegulaFalsiBisection<>:: solve(eq, low, high, max_iteration, bhp_tolerance, iteration); #ifdef EXTRA_THP_DEBUGGING OpmLog::debug("***** " + name() + " solved_bhp = " + std::to_string(solved_bhp) @@ -3949,7 +3949,7 @@ namespace Opm const double flo_tolerance = 1e-6 * std::fabs(flo_samples.back()); int iteration = 0; try { - const double solved_bhp = RegulaFalsi<>:: + const double solved_bhp = RegulaFalsiBisection<>:: solve(eq, low, high, max_iteration, flo_tolerance, iteration); bhp_samples.push_back(solved_bhp); } @@ -4021,7 +4021,7 @@ namespace Opm return boost::optional(); } try { - const double solved_bhp = RegulaFalsi<>:: + const double solved_bhp = RegulaFalsiBisection<>:: solve(eq, low, high, max_iteration, bhp_tolerance, iteration); #ifdef EXTRA_THP_DEBUGGING OpmLog::debug("***** " + name() + " solved_bhp = " + std::to_string(solved_bhp)