Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Unsure about update information in HDF5 output #10

Open
fsohn opened this issue Oct 30, 2017 · 1 comment
Open

Unsure about update information in HDF5 output #10

fsohn opened this issue Oct 30, 2017 · 1 comment

Comments

@fsohn
Copy link
Collaborator

fsohn commented Oct 30, 2017

Hi,
recently I had some trouble understanding the information about the updates in the output HDF5 file, e.g. the number of attempts and acceptances of an update. For example, in

/simulation/results/1-pair_insertion_remover_accepted_scalar Group
/simulation/results/1-pair_insertion_remover_accepted_scalar/count Dataset {SCALAR}
/simulation/results/1-pair_insertion_remover_accepted_scalar/mean Group
/simulation/results/1-pair_insertion_remover_accepted_scalar/mean/error Dataset {SCALAR}
/simulation/results/1-pair_insertion_remover_accepted_scalar/mean/value Dataset {SCALAR}

I found an integer number in the 'count' dataset, which I expected to be the number of accepted updates of this type. However, the same number showed up in the 'count' dataset for attempted updates. So, I am not sure what it actually means because not every attempted update will be accepted. Furthermore, how exactly is the 'mean/value' (and the 'mean/error') obtained?

I assume, that a certain number of updates of different types is performed, where the number of attempts of a certain update is determined by an update rate, before the information is passed on to the ALPSCore library and the 'count' is increased by 1. The 'mean/value' then is the mean value of the 'count' times-accumulated data. Is that correct?

Could you implement a counter for the number of attempted, valid and accepted updates for each type of update?

@shinaoka
Copy link
Collaborator

shinaoka commented Nov 1, 2017

Hi, all data below /simulation/results are not intended to be public to users.
These data are sort of raw data before postparocess (basis rotation, normalization etc.)
You can see how high the acceptance rates are in output messages instead.

Internally, we measure 0 for a rejected update and 1 for an accepted update by means of the accumulator of the ALPSCore libraries.
"count" is the number of these measurements.
"mean/value" is the mean value of this observable, which should be the actual acceptance.
"mean/error" is the statistical error bar of this observable.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

2 participants