From 3197854e89eb9ffe5ff3137e274d16806c2028a5 Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@gmail.com>
Date: Tue, 27 Jul 2021 12:51:09 +0000
Subject: [PATCH] Found one bug in the code when calculating mprod. Fixed.

---
 src/general.h |   20 +++++++++++++++++---
 1 files changed, 17 insertions(+), 3 deletions(-)

diff --git a/src/general.h b/src/general.h
index 9a8feab..4388e2c 100644
--- a/src/general.h
+++ b/src/general.h
@@ -158,6 +158,10 @@
         ts_double solAngle;
 		struct ts_poly *grafted_poly;
 		struct ts_cluster *cluster;
+	ts_double edge_x[14];
+	ts_double edge_y[14];
+	ts_double edge_z[14];
+
 };
 typedef struct ts_vertex ts_vertex;
 
@@ -201,6 +205,7 @@
 
 struct ts_triangle_list{
     ts_uint n;
+	ts_double a0;
     ts_triangle **tria;
 };
 typedef struct ts_triangle_list ts_triangle_list;
@@ -218,7 +223,7 @@
     ts_cell **cell;
     ts_double dcell;
     ts_double shift;
-    ts_double max_occupancy;
+    ts_uint max_occupancy;
 	ts_double dmin_interspecies;
 } ts_cell_list;
 
@@ -251,6 +256,11 @@
 typedef struct ts_poly_list ts_poly_list;
 
 
+typedef struct{
+	ts_float z_max;
+	ts_float z_min;
+	ts_int force_switch;
+} ts_confinement_plane;
 
 typedef struct {
 	long int nshell;
@@ -292,6 +302,9 @@
 	ts_double c0;
 	ts_double w;
 	ts_double F;
+	long int plane_confinement_switch;
+	ts_double plane_d;
+	ts_double plane_F;
 } ts_tape;
 
 
@@ -317,13 +330,14 @@
 	ts_double spring_constant;
 	ts_double pressure;
 	ts_int pswitch;
-    ts_tape *tape;
+ 	ts_tape *tape;
 	ts_double R_nucleus;
 	ts_double R_nucleusX;
 	ts_double R_nucleusY;
 	ts_double R_nucleusZ;
 	ts_double nucleus_center[3];
-    ts_double area;
+	ts_double area;
+	ts_confinement_plane confinement_plane;
 } ts_vesicle;
 
 

--
Gitblit v1.9.3