From 7d326cfdd6f204eef7099186babe808ee4e5487a Mon Sep 17 00:00:00 2001 From: fjosw Date: Mon, 7 Mar 2022 11:28:15 +0000 Subject: [PATCH] Documentation updated --- docs/pyerrors/obs.html | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/docs/pyerrors/obs.html b/docs/pyerrors/obs.html index cc4f0571..97b5bbf2 100644 --- a/docs/pyerrors/obs.html +++ b/docs/pyerrors/obs.html @@ -571,6 +571,7 @@ gamma_div = np.zeros(w_max) for r_name in e_content[e_name]: gamma_div += self._calc_gamma(np.ones((self.shape[r_name])), self.idl[r_name], self.shape[r_name], w_max, fft) + gamma_div[gamma_div < 1] = 1.0 e_gamma[e_name] /= gamma_div[:w_max] if np.abs(e_gamma[e_name][0]) < 10 * np.finfo(float).tiny: # Prevent division by zero @@ -1767,7 +1768,7 @@ continue gamma_div += calc_gamma(np.ones(obs1.shape[r_name]), np.ones(obs2.shape[r_name]), obs1.idl[r_name], obs2.idl[r_name], idl_d[r_name]) e_N += len(idl_d[r_name]) - gamma /= gamma_div + gamma /= max(gamma_div, 1.0) # Bias correction hep-lat/0306017 eq. (49) dvalue += (1 + 1 / e_N) * gamma / e_N @@ -2139,6 +2140,7 @@ gamma_div = np.zeros(w_max) for r_name in e_content[e_name]: gamma_div += self._calc_gamma(np.ones((self.shape[r_name])), self.idl[r_name], self.shape[r_name], w_max, fft) + gamma_div[gamma_div < 1] = 1.0 e_gamma[e_name] /= gamma_div[:w_max] if np.abs(e_gamma[e_name][0]) < 10 * np.finfo(float).tiny: # Prevent division by zero @@ -3162,6 +3164,7 @@ list of ranges or lists on which the samples are defined gamma_div = np.zeros(w_max) for r_name in e_content[e_name]: gamma_div += self._calc_gamma(np.ones((self.shape[r_name])), self.idl[r_name], self.shape[r_name], w_max, fft) + gamma_div[gamma_div < 1] = 1.0 e_gamma[e_name] /= gamma_div[:w_max] if np.abs(e_gamma[e_name][0]) < 10 * np.finfo(float).tiny: # Prevent division by zero