src/energy.c
@@ -240,5 +240,10 @@ } void stretchenergy(ts_vesicle *vesicle, ts_triangle *triangle){ triangle->energy=vesicle->tape->xkA0/vesicle->tlist->n*pow((triangle->area/vesicle->tlist->a0-1),2); triangle->energy=vesicle->tape->xkA0/2.0*pow((triangle->area/vesicle->tlist->a0-1.0),2); } ts_double stretchenergy2(ts_double current_area,ts_double tensionless_area){ return pow((current_area-tensionless_area),2)/tensionless_area; }