lattice/pyerrors
1
0
Fork 0
mirror of https://github.com/fjosw/pyerrors.git synced 2025-05-14 19:43:41 +02:00
Code Issues Projects Releases Packages Wiki Activity

Merge branch 'develop'

Browse source
This commit is contained in:
Fabian Joswig 2020-11-03 15:12:50 +01:00
commit 839a690725
3 changed files with 10 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

5
CHANGELOG.md
View file

@ -2,6 +2,11 @@
All notable changes to this project will be documented in this file. All notable changes to this project will be documented in this file.
## [1.0.1] - 2020-11-03
### Fixed
- Bug in pyerrors.covariance fixed that appeared when working with several
replica of different length.
## [1.0.0] - 2020-10-13 ## [1.0.0] - 2020-10-13
### Added ### Added
- Compatibility with the BDIO Native format outlined [here](https://ific.uv.es/~alramos/docs/ADerrors/tutorial/). Read and write function added to input.bdio - Compatibility with the BDIO Native format outlined [here](https://ific.uv.es/~alramos/docs/ADerrors/tutorial/). Read and write function added to input.bdio

7
pyerrors/pyerrors.py
View file

@ -852,12 +852,12 @@ def covariance(obs1, obs2, correlation=False, **kwargs):
r_length.append(len(obs1.deltas[r_name])) r_length.append(len(obs1.deltas[r_name]))
gamma += np.sum(obs1.deltas[r_name] * obs2.deltas[r_name]) / len(obs1.deltas[r_name]) gamma += np.sum(obs1.deltas[r_name] * obs2.deltas[r_name])
e_N = np.sum(r_length) e_N = np.sum(r_length)
tau_combined = (obs1.e_tauint[e_name] + obs2.e_tauint[e_name]) / 2 tau_combined = (obs1.e_tauint[e_name] + obs2.e_tauint[e_name]) / 2
dvalue += gamma * (1 + 1 / e_N) / e_N * 2 * tau_combined dvalue += gamma / e_N * (1 + 1 / e_N) / e_N * 2 * tau_combined
if np.abs(dvalue / obs1.dvalue / obs2.dvalue) > 1.0: if np.abs(dvalue / obs1.dvalue / obs2.dvalue) > 1.0:
dvalue = np.sign(dvalue) * obs1.dvalue * obs2.dvalue dvalue = np.sign(dvalue) * obs1.dvalue * obs2.dvalue
@ -1137,7 +1137,8 @@ def plot_corrs(observables, **kwargs):
y[i] = observables[j][i].value y[i] = observables[j][i].value
y_err[i] = observables[j][i].dvalue y_err[i] = observables[j][i].dvalue
plt.errorbar(x, y, yerr=y_err, ls='none', fmt='o', capsize=3, markersize=5, label=label[j]) plt.errorbar(x, y, yerr=y_err, ls='none', fmt='o', capsize=3,
markersize=5, lw=1, label=label[j])
if kwargs.get('logscale'): if kwargs.get('logscale'):
plt.yscale('log') plt.yscale('log')

2
setup.py
View file

@ -3,7 +3,7 @@
from setuptools import setup, find_packages from setuptools import setup, find_packages
setup(name='pyerrors', setup(name='pyerrors',
version='1.0.0', version='1.0.1',
description='Error analysis for lattice QCD', description='Error analysis for lattice QCD',
author='Fabian Joswig', author='Fabian Joswig',
author_email='fabian.joswig@wwu.de', author_email='fabian.joswig@wwu.de',