| | |
| | | ts_bool centermass(ts_vesicle *vesicle){ |
| | | ts_uint i,j, n=vesicle->vlist->n; |
| | | ts_vertex **vtx=vesicle->vlist->vtx; |
| | | ts_double temp_z_cm=0; |
| | | vesicle->cm[0]=0; |
| | | vesicle->cm[1]=0; |
| | | vesicle->cm[2]=0; |
| | |
| | | vesicle->cm[0]/=(ts_float)n; |
| | | vesicle->cm[1]/=(ts_float)n; |
| | | vesicle->cm[2]/=(ts_float)n; |
| | | |
| | | //center mass for z component does not change if we confine the vesicle with plates |
| | | if(vesicle->tape->plane_confinement_switch){ |
| | | temp_z_cm=vesicle->cm[2]; |
| | | vesicle->cm[2]=0; |
| | | } |
| | | |
| | | for(i=0;i<n;i++){ |
| | | vtx[i]->x-=vesicle->cm[0]; |
| | |
| | | |
| | | vesicle->cm[0]=0.0; |
| | | vesicle->cm[1]=0.0; |
| | | vesicle->cm[2]=0.0; |
| | | if(vesicle->tape->plane_confinement_switch){ |
| | | vesicle->cm[2]=temp_z_cm; |
| | | } |
| | | |
| | | for(i=0;i<vesicle->tlist->n;i++){ |
| | | triangle_normal_vector(vesicle->tlist->tria[i]); |