Trisurf Monte Carlo simulator
Samo Penic
2022-05-04 06afc729f9061c1cfa14c78728d61d518324c2f0
src/energy.c
@@ -213,7 +213,7 @@
}
ts_double direct_force_energy(ts_vesicle *vesicle, ts_vertex *vtx, ts_vertex *vtx_old){
   if(fabs(vtx->c)<1e-15) return 0.0;
   if(fabs(vtx->direct_interaction_force)<1e-15) return 0.0;
//   printf("was here");
   if(fabs(vesicle->tape->F)<1e-15) return 0.0;
@@ -240,5 +240,5 @@
}
void stretchenergy(ts_vesicle *vesicle, ts_triangle *triangle){
   triangle->energy=vesicle->tape->xkA0/vesicle->tlist->n/2*pow((triangle->area/vesicle->tlist->a0-1),2);
   triangle->energy=vesicle->tape->xkA0/2.0*pow((triangle->area/vesicle->tlist->a0-1.0),2);
}