Skip to main content

Author: Dirk Hesse <herr.dirk.hesse@gmail.com>

We implement the method to estimate autocorrelation times of Monte Carlo data presented in

U. Wolff [ALPHA Collaboration], Monte Carlo errors with less errors, Comput. Phys. Commun. 156, 143 (2004) [hep-lat/0306017].

PUBLICATIONS MAKING USE OF THIS CODE MUST CITE THE PAPER.

The main objective is the following: Data coming from a Monte Carlo simulation usually suffers from autocorrelation. It is not straight-forward to estimate this autocorrelation, which is required to give robust estimates for errors. This program implements a method proposed by Wolff to estimate autocorrelations in a safe way.

Quick start

This package contains code to generate correlated data, so we can conveniently demonstrate the basic functionality of the code in a short example:

>>> from puwr import tauint, correlated_data
>>> correlated_data(2, 10)
[[array([ 1.02833043,  1.08615234,  1.16421776,  1.15975754,
          1.23046603,  1.13941114,  1.1485227 ,  1.13464388,
          1.12461557,  1.15413354])]]
>>> mean, delta, tint, d_tint = tauint(correlated_data(10, 200), 0)
>>> print "mean = {0} +/- {1}".format(mean, delta)
mean = 1.42726267057 +/- 0.03013853
>>> print "tau_int = {0} +/- {1}".format(tint, d_tint)
tau_int = 9.89344869217 +/- 4.10466090332

The data is expected to be in the format data[observable][replicum][measurement]. See the documentation that comes with this code for more information.

License

See LICENSE file.

Release files for py-uwerr 0.1

For a detailed explanation of source distributions (sdists) and built distributions (wheels), please see the package formats documentation.

Built distribution (wheel)

Table of built distributions (wheels) for py-uwerr 0.1
File Interpreter ABI Platform
py_uwerr-0.1-py3-none-any.whl Python 3 none any Details

Release files / py_uwerr-0.1-py3-none-any.whl

Download URL py_uwerr-0.1-py3-none-any.whl
Size 6.9 kB
Tags Python 3
SHA-256 checksum
How to use checksums
18d56950ad4ffeaf750ec15a256a3a261b42a20327ea7ecd7aef742b7d476483
BLAKE2b-256 checksum
How to use checksums
4d321fe59cf9cdd3512b1dbe326aff8cf93e466258e828af7b842d2da98ae548
Upload date
Uploaded using Trusted Publishing?
What is trusted publishing?
No
Uploaded via python-requests/2.26.0

Release history Release notifications | RSS feed

This release

0.1 This release

1 release file

Anthropic, PBC Visionary sponsor Bloomberg Visionary sponsor Hudson River Trading Visionary sponsor Meta Visionary sponsor NVIDIA Visionary sponsor Microsoft Sustainability sponsor Depot Continuous Integration AWS Cloud computing and Security Sponsor Datadog Monitoring Fastly CDN Google Download Analytics Sentry Error logging StatusPage Status page