| | |
| | | } |
| | | |
| | | 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; |
| | | |
| | |
| | | } |
| | | |
| | | 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); |
| | | } |