Merge pull request #84 from fjosw/refactor/covariance_e_N_determination

Simplification and speed up of the computation of e_N in _covariance_element
This commit is contained in:
Fabian Joswig 2022-03-03 11:55:47 +00:00 committed by GitHub
commit 329e4abfc3
No known key found for this signature in database
GPG key ID: 4AEE18F83AFDEB23

View file

@ -1446,7 +1446,7 @@ def _covariance_element(obs1, obs2):
if r_name not in obs2.e_content[e_name]:
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 += np.sum(np.ones_like(idl_d[r_name]))
e_N += len(idl_d[r_name])
gamma /= gamma_div
# Bias correction hep-lat/0306017 eq. (49)