-

diff --git a/ising.c b/ising.c
index 0b342332c8cf0dde84c4ef50679e791f48753e2a..d3690247f9a94136609a8495d27ad637656ba65a 100644 (file)
--- a/ising.c
+++ b/ising.c
@@ -200,7 +200,7 @@ int main(int argc, char **argv)
    if((*(atom+((max_x+x_c-1)%max_x)+y_c*max_x))&1) ++count_p;
    if(((*(atom+x_c+y_c*max_x))&1)==0) count_p=4-count_p;
    delta_e=(2*count_p-4)*s;
    if((*(atom+((max_x+x_c-1)%max_x)+y_c*max_x))&1) ++count_p;
    if(((*(atom+x_c+y_c*max_x))&1)==0) count_p=4-count_p;
    delta_e=(2*count_p-4)*s;

-   if(delta_e<0) *(atom+x_c+y_c*max_x)=(*(atom+x_c+y_c*max_x)+1)&1;
+   if(delta_e<=0) *(atom+x_c+y_c*max_x)=(*(atom+x_c+y_c*max_x)+1)&1;

    else
    {
     if(1.0*rand()/RAND_MAX<exp(-1.0*delta_e*beta))
    else
    {
     if(1.0*rand()/RAND_MAX<exp(-1.0*delta_e*beta))