| | |
| | | /* vim: set ts=4 sts=4 sw=4 noet : */ |
| | | #include<general.h> |
| | | #include "vesicle.h" |
| | | #include "vertex.h" |
| | |
| | | triangle_list_free(vesicle->tlist); |
| | | cell_list_free(vesicle->clist); |
| | | poly_list_free(vesicle->poly_list); |
| | | poly_list_free(vesicle->filament_list); |
| | | complex_sph_free(vesicle->sphHarmonics); |
| | | free(vesicle); |
| | | return TS_SUCCESS; |
| | |
| | | vesicle->area=area; |
| | | return TS_SUCCESS; |
| | | } |
| | | |
| | | ts_double vesicle_meancurvature(ts_vesicle *vesicle){ |
| | | // Integrates (H dA) over vesicle area A, where H=(C1+C2)/2. |
| | | // (To be devided by A outside of function) |
| | | ts_double mc; |
| | | ts_uint i; |
| | | mc=0; |
| | | for(i=0;i<vesicle->vlist->n;i++){ |
| | | mc=mc+vesicle->vlist->vtx[i]->curvature; |
| | | } |
| | | return mc/2.0; |
| | | } |