Trisurf Monte Carlo simulator
Samo Penic
2020-04-13 2a1e3d528265c3d366dab9ad4b0c4d3d65e1c4ed
src/bondflip.c
@@ -32,7 +32,7 @@
    ts_vertex *k=bond->vtx2;
    ts_uint nei,neip,neim;
    ts_uint i,j;
    ts_double oldenergy, delta_energy, dvol=0.0, darea=0.0;
    ts_double oldenergy=0.0, delta_energy=0.0, dvol=0.0, darea=0.0;
    ts_triangle *lm=NULL,*lp=NULL, *lp1=NULL, *lm2=NULL;
    ts_vertex *kp,*km;
@@ -163,7 +163,7 @@
/* end backup vertex */
/* Save old energy */
  oldenergy=0;
  oldenergy=0.0;
  oldenergy+=k->xk* k->energy;
  oldenergy+=kp->xk* kp->energy;
  oldenergy+=km->xk* km->energy;
@@ -181,7 +181,7 @@
/* Calculating the new energy */
  delta_energy=0;
  delta_energy=0.0;
  delta_energy+=k->xk* k->energy;
  delta_energy+=kp->xk* kp->energy;
  delta_energy+=km->xk* km->energy;