Trisurf Monte Carlo simulator
Samo Penic
2014-04-28 9166cbcd0e28d61a69646911af35bb7895ff9203
src/vertexmove.c
@@ -9,6 +9,7 @@
#include "timestep.h"
#include "cell.h"
//#include "io.h"
#include "io.h"
#include<stdio.h>
#include "vertexmove.h"
#include <string.h>
@@ -90,7 +91,7 @@
   memcpy((void *)&backupvtx[i+1],(void *)vtx->neigh[i],sizeof(ts_vertex));
   }
   if(vesicle->pswitch == 1){
   if(vesicle->pswitch == 1 || vesicle->tape->constvolswitch == 1){
      for(i=0;i<vtx->tristar_no;i++) dvol-=vtx->tristar[i]->volume;
   };
@@ -107,9 +108,9 @@
        delta_energy+=vtx->neigh[i]->xk*(vtx->neigh[i]->energy-oenergy);
    }
   if(vesicle->pswitch == 1){
   if(vesicle->pswitch == 1 || vesicle->tape->constvolswitch == 1){
      for(i=0;i<vtx->tristar_no;i++) dvol+=vtx->tristar[i]->volume;
      delta_energy-=vesicle->pressure*dvol;
        if(vesicle->pswitch == 1) delta_energy-=vesicle->pressure*dvol;
   };
/* No poly-bond energy for now!