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-22 15:09:01 +02:00
parent 5f0b397793
commit 49cc8a5c7d
1 changed files with 3 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
src/seismic_hazard_forecasting.py
View File

@@ -185,11 +185,10 @@ def compute_IMT_exceedance(rx_lat, rx_lon, r, fr, p, lambdas, D, percentages_D,
start = timer()
try:
method='brenth'
method='brentq'
logger.debug("Now trying Scipy " + method)
#output = root_scalar(exceedance_root_function, bracket=[IMT_min, IMT_max], rtol=rtol, method=method)
#gm_est = output.root
gm_est = np.nan
output = root_scalar(exceedance_root_function, bracket=[IMT_min, IMT_max], rtol=rtol, method=method)
gm_est = output.root
except Exception as error:
logger.error(f"An exception occurred: {error}")
logger.info("Set ground motion value to nan")