From b06a2ab4e54f5d1c9087667e8bece541e4924820 Mon Sep 17 00:00:00 2001 From: ftong Date: Mon, 25 Aug 2025 13:32:59 +0200 Subject: [PATCH] Update src/seismic_hazard_forecasting.py --- src/seismic_hazard_forecasting.py | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py index dd6d5ea..6b5052b 100644 --- a/src/seismic_hazard_forecasting.py +++ b/src/seismic_hazard_forecasting.py @@ -185,9 +185,10 @@ def compute_IMT_exceedance(rx_lat, rx_lon, r, fr, p, lambdas, D, percentages_D, start = timer() try: - method='bisect' + method='brentq' logger.debug("Now trying Scipy " + method) - output = root_scalar(exceedance_root_function, bracket=[IMT_min, IMT_max], rtol=rtol, method=method) + #output = root_scalar(exceedance_root_function, bracket=[IMT_min, IMT_max], rtol=rtol, method=method) + output = root_scalar(exceedance_root_function, method=method) gm_est = output.root except Exception as error: logger.error(f"An exception occurred: {error}")