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; @@ -238,3 +238,7 @@ return vesicle->tape->F*ddp; } void stretchenergy(ts_vesicle *vesicle, ts_triangle *triangle){ triangle->energy=vesicle->tape->xkA0/2.0*pow((triangle->area/vesicle->tlist->a0-1.0),2); }