Skip to content

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.

Edited by Manuel Schrauth