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

fix ft_ref check of shifted bcf

Browse source
This commit is contained in:
Valentin Boettcher 2022-07-19 16:24:19 +02:00
parent 0035c0185a
commit 3a9bd00a45
1 changed files with 6 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

7
stocproc/method_ft.py
View file

@ -502,9 +502,14 @@ def get_N_a_b_for_accurate_fourier_integral(
""" """
if opt_b_only:
I0 = quad(integrand, 0, np.inf)[0]
# sometimes the integrand is zero for a while...
x0 = 0
while integrand(x0 + 0.1) == 0:
x0 += 0.1
I0 = quad(integrand, x0, np.inf)[0]
else:
I0 = quad(integrand, -np.inf, np.inf)[0]
ft_ref_0 = ft_ref(0)
rd = np.abs(ft_ref_0 - I0) / np.abs(ft_ref_0)
log.debug("ft_ref check yields rd {:.3e}".format(rd))