From 79048ff443472f339a95f070940a6df4b14a9035 Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@gmail.com>
Date: Thu, 07 Jun 2012 21:29:41 +0000
Subject: [PATCH] Ylmi and Ulm seems to be working... needs to be further debugged. Triangle area is missing, thus ulm are always 0.
---
src/vesicle.c | 18 ++++++++++++++++++
1 files changed, 18 insertions(+), 0 deletions(-)
diff --git a/src/vesicle.c b/src/vesicle.c
index d6c29e1..e1cc21d 100644
--- a/src/vesicle.c
+++ b/src/vesicle.c
@@ -36,3 +36,21 @@
free(vesicle);
return TS_SUCCESS;
}
+
+ts_bool vesicle_volume(ts_vesicle *vesicle){
+ ts_double volume;
+ ts_double vol;
+ ts_uint i;
+ ts_triangle **tria=vesicle->tlist->tria;
+ volume=0;
+ for(i=0; i<vesicle->tlist->n;i++){
+ vol=(tria[i]->vertex[0]->x+ tria[i]->vertex[1]->x + tria[i]->vertex[2]->x) * tria[i]->xnorm +
+ (tria[i]->vertex[0]->y+ tria[i]->vertex[1]->y + tria[i]->vertex[2]->y) * tria[i]->ynorm +
+ (tria[i]->vertex[0]->z+ tria[i]->vertex[1]->z + tria[i]->vertex[2]->z) *
+tria[i]->znorm;
+ volume=volume-vol/18.0;
+ }
+
+ vesicle->volume=volume;
+ return TS_SUCCESS;
+}
--
Gitblit v1.9.3