From 6c3bb91fea92c5013507bdf4f5ac49e889e213ff Mon Sep 17 00:00:00 2001 From: Samo Penic <samo.penic@fe.uni-lj.si> Date: Tue, 06 Oct 2015 12:07:39 +0000 Subject: [PATCH] Unkonown merge difference between different directories are merged together. --- src/vesicle.c | 59 ++++++++++++++++++++++++++++++++++++++++++++++++++--------- 1 files changed, 50 insertions(+), 9 deletions(-) diff --git a/src/vesicle.c b/src/vesicle.c index c87ead1..f63eb12 100644 --- a/src/vesicle.c +++ b/src/vesicle.c @@ -1,21 +1,62 @@ #include<general.h> #include "vesicle.h" +#include "vertex.h" +#include "triangle.h" +#include "bond.h" +#include "cell.h" +#include "stdlib.h" +#include "poly.h" +#include "sh.h" +#include "shcomplex.h" + +ts_vesicle *init_vesicle(ts_uint N, ts_uint ncmax1, ts_uint ncmax2, ts_uint +ncmax3, ts_double stepsize){ + ts_vesicle *vesicle=(ts_vesicle *)calloc(1,sizeof(ts_vesicle)); + vesicle->vlist=init_vertex_list(N); + vesicle->blist=init_bond_list(); + vesicle->tlist=init_triangle_list(); + vesicle->clist=init_cell_list(ncmax1, ncmax2, ncmax3, stepsize); + return vesicle; +} + ts_bool vesicle_translate(ts_vesicle *vesicle,ts_double x, ts_double y, ts_double z){ ts_uint i; - ts_vertex *vtx=vesicle->vlist.vertex; - ts_uint nn=vesicle->vlist.n; + ts_vertex **vtx=vesicle->vlist->vtx; + ts_uint nn=vesicle->vlist->n; for(i=0;i<nn;i++){ - vtx[i].x+=x; - vtx[i].y+=y; - vtx[i].z+=z; + vtx[i]->x+=x; + vtx[i]->y+=y; + vtx[i]->z+=z; } return TS_SUCCESS; } ts_bool vesicle_free(ts_vesicle *vesicle){ - vertex_list_free(&vesicle->vlist); - bond_list_free(&vesicle->blist); - triangle_list_free(&vesicle->tlist); - cell_list_free(&vesicle->clist); + vtx_list_free(vesicle->vlist); + bond_list_free(vesicle->blist); + triangle_list_free(vesicle->tlist); + cell_list_free(vesicle->clist); + poly_list_free(vesicle->poly_list); + complex_sph_free(vesicle->sphHarmonics); + free(vesicle); + return TS_SUCCESS; +} + +/* @brief Function makes a sum of partial volumes of each triangle. Volumes of + * + * Partial volumes are calculated when we calculate normals of triangles. It is + * relatively easy to calculate the volume of vesicle if we take into account + * that the volume of the whole vertex is simply sum of all partial volumes of + * all the triangles. + */ +ts_bool vesicle_volume(ts_vesicle *vesicle){ + ts_double volume; + ts_uint i; + ts_triangle **tria=vesicle->tlist->tria; + volume=0; + for(i=0; i<vesicle->tlist->n;i++){ + volume=volume+tria[i]->volume; + } + vesicle->volume=volume; return TS_SUCCESS; } -- Gitblit v1.9.3