Trisurf Monte Carlo simulator
Samo Penic
2018-10-11 b4e2e6dbd5049ae279f0f74838cd21b625148c72
src/vesicle.c
@@ -38,6 +38,7 @@
    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;
@@ -78,3 +79,15 @@
    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;
}