From 2b054ce2a2748577e598fb1499b4838f7b5dc72a Mon Sep 17 00:00:00 2001
From: Kai Bao <kai.bao@sintef.no>
Date: Mon, 1 May 2023 16:00:05 +0200
Subject: [PATCH] updating the interval when doing bisectBracket()

to avoid getting stuck and iterate for no purpose. It might not affect the
result much, while from the code, it looks like it should be updated
iteratively.
---
 opm/simulators/wells/WellBhpThpCalculator.cpp | 1 +
 1 file changed, 1 insertion(+)

diff --git a/opm/simulators/wells/WellBhpThpCalculator.cpp b/opm/simulators/wells/WellBhpThpCalculator.cpp
index 03b0e1c18..b3f1a056d 100644
--- a/opm/simulators/wells/WellBhpThpCalculator.cpp
+++ b/opm/simulators/wells/WellBhpThpCalculator.cpp
@@ -695,6 +695,7 @@ bisectBracket(const std::function<double(const double)>& eq,
                 eq_high = eq(high);
                 abs_high = std::fabs(eq_high);
             }
+            interval = high - low;
             ++bracket_attempts;
         }