official-apps/SeismicHazardForecasting
10
0
Fork 1
Code Issues Pull Requests Activity

Update src/seismic_hazard_forecasting.py

Browse Source
This commit is contained in:
ftong
2025-08-25 13:32:59 +02:00
parent 2002114425
commit b06a2ab4e5
1 changed files with 3 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
src/seismic_hazard_forecasting.py
View File

@@ -185,9 +185,10 @@ def compute_IMT_exceedance(rx_lat, rx_lon, r, fr, p, lambdas, D, percentages_D,
start = timer() start = timer()
try: try:
method='bisect' method='brentq'
logger.debug("Now trying Scipy " + method) 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 gm_est = output.root
except Exception as error: except Exception as error:
logger.error(f"An exception occurred: {error}") logger.error(f"An exception occurred: {error}")