From 2ae815f58a40d4b86f5183d7e2ec406d5e0f5d65 Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@gmail.com>
Date: Wed, 09 May 2018 14:47:15 +0000
Subject: [PATCH] Prepared for stretching

---
 src/io.c      |    2 ++
 src/tape      |    5 +++++
 src/general.h |    3 +++
 3 files changed, 10 insertions(+), 0 deletions(-)

diff --git a/src/general.h b/src/general.h
index a0e31a5..9a8feab 100644
--- a/src/general.h
+++ b/src/general.h
@@ -195,6 +195,7 @@
 	ts_double znorm;
     ts_double area; // firstly needed for sh.c
     ts_double volume; // firstly needed for sh.c
+	ts_double energy;
 };
 typedef struct ts_triangle ts_triangle;
 
@@ -268,6 +269,8 @@
 	long int pswitch;
     long int constvolswitch;
     long int constareaswitch;
+	long int stretchswitch;
+	ts_double xkA0;
     ts_double constvolprecision;
     	char *multiprocessing;
    	long int brezveze0;
diff --git a/src/io.c b/src/io.c
index 5a5e9a2..f513501 100644
--- a/src/io.c
+++ b/src/io.c
@@ -1165,6 +1165,8 @@
 	CFG_SIMPLE_INT("constvolswitch",&tape->constvolswitch),
 	CFG_SIMPLE_INT("constareaswitch",&tape->constareaswitch),
 	CFG_SIMPLE_FLOAT("constvolprecision",&tape->constvolprecision),
+	CFG_SIMPLE_INT("stretchswitch",&tape->stretchswitch),
+	CFG_SIMPLE_FLOAT("xkA0",&tape->xkA0),	
 	CFG_SIMPLE_FLOAT("pressure",&tape->pressure),
 	CFG_SIMPLE_FLOAT("k_spring",&tape->kspring),
 	CFG_SIMPLE_FLOAT("xi",&tape->xi),
diff --git a/src/tape b/src/tape
index 8f38645..70f526a 100644
--- a/src/tape
+++ b/src/tape
@@ -23,6 +23,11 @@
 #Constant area constraint (0 disable constant area, 2 enable constant area with epsarea)
 constareaswitch=0
 
+
+#Stretching
+stretchswitch=1
+xkA0=10.0
+
 ####### Polymer (brush) definitions ###########
 # npoly is a number of polymers attached to npoly distinct vertices on vesicle
 npoly=0

--
Gitblit v1.9.3