hiro/stocproc
1
0
Fork 0
mirror of https://github.com/vale981/stocproc synced 2025-03-04 17:21:42 -05:00
Code Issues Projects Releases Packages Wiki Activity Actions

optimize the reference value better

Browse source
This commit is contained in:
Valentin Boettcher 2022-06-17 14:41:24 +02:00
parent 95b8c8e629
commit aec44f2248
1 changed files with 3 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
stocproc/method_ft.py
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 += 1
ref_val += 0.1
ref_val = minimize(lambda x: -abs(integrand(x)), x0=ref_val).x[0]
I_ref = integrand(ref_val)
log.debug("find_integral_boundary")
log.debug("I_ref: {:.2e} at ref_val: {:.2e}".format(I_ref, ref_val))
if I_ref < tol: