From f977e8fe32b6a1a2ae85662ab35fd7c9f647cc3e Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@gmail.com>
Date: Mon, 02 Jul 2018 10:50:58 +0000
Subject: [PATCH] Added a funciton that calculates everything needed for spherical harmonics

---
 src/vesicle.c |   12 ++++++++++++
 1 files changed, 12 insertions(+), 0 deletions(-)

diff --git a/src/vesicle.c b/src/vesicle.c
index bc69e62..fe4c4b1 100644
--- a/src/vesicle.c
+++ b/src/vesicle.c
@@ -79,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;
+}

--
Gitblit v1.9.3