maybe bugs for code "equil"

Alloy Theoretic Automated Toolkit (ATAT)
1

In the file `equil’, the following code-snippet does not seem to do any averaging…
for (int i=bin1+1; i<=cur_bin; i++) {
avg_corr+=corr;
}
avg_corr/=(Real)(cur_bin-bin1);
Am I overlooking something basic? is it intentional?

2

I’m very sorry for the very delayed answer! (It is proposal writing season in the US and I have a high teaching load this semester!)

Indeed

avg_corr+=corr;

should have been

avg_corr+=buf_corr[i];

Fortunately, this would have had a minor effect on the results.
The code only used the correlation evaluated from the latest time step buffer instead of all available ones. This correlation should not change over time, so there is no systematic error introduced.

Sorry for the bug and thanks for pointing it out!!