specify bound for minimization

This commit is contained in:
Valentin Boettcher 2022-06-21 13:34:13 +02:00
parent aec44f2248
commit f14ff4bd77

View file

@ -49,9 +49,11 @@ def find_integral_boundary(integrand, tol, ref_val, max_val, x0):
x0 = float(x0)
while integrand(ref_val) < 1e-2:
ref_val += 0.1
ref_val += 0.01
ref_val = minimize(lambda x: -abs(integrand(x)), x0=ref_val).x[0]
ref_val = minimize(
lambda x: -abs(integrand(x)), x0=ref_val, bounds=[(0, np.inf)]
).x[0]
I_ref = integrand(ref_val)
log.debug("find_integral_boundary")