dangerous bug in Metropolis
When summing over the neighbours, averagevector
must not in general be a StateVector. Why? In models where StateVector uses an int
data type, this line can still produce doubles
averagevector = averagevector + callbonds<Grid>(grid, a, rsite, i, nbr);
This should be fixed as soon as possible.