From 142a67fe82b830e5c7816914afa62445959c87ca Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@fe.uni-lj.si>
Date: Tue, 05 Nov 2013 14:04:21 +0000
Subject: [PATCH] changes in bondflip call. No need to bondflip all the bonds, but only as many bonds as there are vertices. Also, rnvec seems to be not needed for bondflip, so it is commented out

---
 src/general.h |   50 +++++++++++++++++++++++++++++++++++++++++---------
 1 files changed, 41 insertions(+), 9 deletions(-)

diff --git a/src/general.h b/src/general.h
index 6d2ab06..f3328f6 100644
--- a/src/general.h
+++ b/src/general.h
@@ -110,12 +110,29 @@
 
 /* STRUCTURES */
 
+
+/** @brief Data structure for keeping the coordinates in selected coordinate
+ * system
+ */
+#define TS_COORD_CARTESIAN 0
+#define TS_COORD_SPHERICAL 1
+#define TS_COORD_CYLINDRICAL 2
+
+typedef struct {
+    ts_double e1;
+    ts_double e2;
+    ts_double e3;
+    ts_uint coord_type;
+} ts_coord;
+
+
+
 /** @brief Data structure of all data connected to a vertex
  *
- *  ts_vertex_data holds the data for one single point (bead, vertex). To understand how to use it
+ *  ts_vertex holds the data for one single point (bead, vertex). To understand how to use it
  *  here is a detailed description of the fields in the data structure. */
-struct ts_vertex_data {
-        ts_uint idx; /**< Represents index of the vertex point. Should become obsolete, since it is also present in ts_vertex structure. */        
+struct ts_vertex {
+        ts_uint idx;
         ts_double x; /**< The x coordinate of vertex. */
         ts_double y; /**< The y coordinate of vertex. */
         ts_double z; /**< The z coordinate of vertex. */
@@ -134,12 +151,9 @@
         ts_double xk;
         ts_double c;
         ts_uint id;
-};
-typedef struct ts_vertex_data ts_vertex_data;
-
-struct ts_vertex {
-        ts_uint idx;
-        ts_vertex_data *data;
+        ts_double projArea;
+        ts_double relR;
+        ts_double solAngle;
 };
 typedef struct ts_vertex ts_vertex;
 
@@ -155,6 +169,7 @@
 	ts_vertex *vtx2;
     ts_double bond_length;
     ts_double bond_length_dual;
+	ts_bool tainted;
 };
 typedef struct ts_bond ts_bond;
 
@@ -172,6 +187,8 @@
 	ts_double xnorm;
 	ts_double ynorm;
 	ts_double znorm;
+    ts_double area; // firstly needed for sh.c
+    ts_double volume; // firstly needed for sh.c
 };
 typedef struct ts_triangle ts_triangle;
 
@@ -199,6 +216,17 @@
 
 
 typedef struct {
+    ts_uint l;
+    ts_double **ulm;
+    ts_double **sumUlm2;
+    ts_uint N;
+    ts_double **co;
+    ts_double ***Ylmi;
+} ts_spharm;
+
+
+
+typedef struct {
 	ts_vertex_list *vlist;
 	ts_bond_list *blist;
 	ts_triangle_list *tlist;
@@ -208,9 +236,13 @@
     ts_double dmax;
     ts_double stepsize;
     ts_double cm[3];
+    ts_double volume;
+    ts_spharm *sphHarmonics;
 } ts_vesicle;
 
 
+
+
 /* GLOBAL VARIABLES */
 
 int quiet;

--
Gitblit v1.9.3