montecarlo- все статьи тега

Алгоритм Метрополиса модели Изинга: решетка не уравновешивается

У меня есть некоторый код для модели Изинга в python (2d), и решетка не достигнет равновесия. Здесь код выводит число вращений, которые переворачиваются для каждой развертки Монте-Карло, и то же самое число переворачивается для каждой развертки. Если я прав, то число, которое переворачивается, должно уменьшаться с каждой разверткой, поскольку решетка достигает равновесия. Может ли кто-нибудь увидеть какие-либо ошибки в коде? import numpy as np from numpy import random as rn N=20 a=[1,-1] b=[N,N ...