Trisurf Monte Carlo simulator
Samo Penic
2017-01-04 96b0a0ac783d35fbf4e764fe8e610ddc2a79af27
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;
}