Trisurf Monte Carlo simulator
Samo Penic
2018-06-21 0a1e99fea2e3d92e073231707361932b5af9cd81
A fix in bondflip energy calculation
1 files modified
2 ■■■ changed files
src/bondflip.c 2 ●●● patch | view | raw | blame | history
src/bondflip.c
@@ -199,7 +199,7 @@
        delta_energy+=lm->energy+lp->energy;
    */
        temp_area=temp_area+lm->area+lp->area;
        delta_energy+=stretchenergy2(temp_area,vesicle->tlist->a0*vesicle->tlist->n)-stretchenergy2(vesicle->area,vesicle->tlist->a0*vesicle->tlist->n);
        delta_energy+=0.5*vesicle->tape->xkA0*(stretchenergy2(temp_area,vesicle->tlist->a0*vesicle->tlist->n)-stretchenergy2(vesicle->area,vesicle->tlist->a0*vesicle->tlist->n));
    }
    delta_energy-=oldenergy;