From 1121faa13a2038facad22073f0fc610903d98691 Mon Sep 17 00:00:00 2001
From: Samo Penic <samo@CAE-linux.(none)>
Date: Fri, 05 Sep 2014 20:18:05 +0000
Subject: [PATCH] First variant of constant volume constrant (a new one proposed by Miha after reading his article). It seems to work, however there are still some things to be done, such as Miha's derivation of the epsvol (0.003% is used at the moment) and solving the problem of additional global variables.

---
 src/tape |   28 ++++++++++++++++++++--------
 1 files changed, 20 insertions(+), 8 deletions(-)

diff --git a/src/tape b/src/tape
index dbaabdb..b0ac65c 100644
--- a/src/tape
+++ b/src/tape
@@ -12,13 +12,18 @@
 
 # Pressure calculations
 # (pswitch=1: calc. p*dV energy contribution)
-pswitch = 1
+pswitch = 0
 # pressure difference: p_inside - p_outside (in units kT/l_min^3):
-pressure=-1.0
+pressure=0.0
+
+#Constant volume constraint (0 disable constant volume, 1 enable)
+constvolswitch=1
+constvolprecision=1e-14
+
 
 ####### Polymer (brush) definitions ###########
 # npoly is a number of polymers attached to npoly distinct vertices on vesicle
-npoly=0
+npoly=500
 # nmono is a number of monomers in each polymer
 nmono=10
 # Spring constant between monomers of the polymer
@@ -26,15 +31,15 @@
 
 ####### Filament (inside the vesicle) definitions ###########
 # nfil is a number of filaments inside the vesicle
-nfil=2
+nfil=0
 # nfono is a number of monomers in each filament
 nfono=300
 # Persistence lenght of the filaments (in units l_min)
-xi=5000
+xi=0
 
 ####### Nucleus (inside the vesicle) ###########
 # Radius of an impenetrable hard sphere inside the vesicle
-R_nucleus=5
+R_nucleus=0
 
 #######  Cell definitions ############
 nxmax=60
@@ -44,13 +49,20 @@
 
 ####### Program Control ############
 #how many MC sweeps between subsequent records of states to disk
-mcsweeps=10
+#200000
+mcsweeps=20
 #how many initial mcsweeps*inititer MC sweeps before recording to disk?
+#2
 inititer=0
 #how many records do you want on the disk iteration are there in a run?
-iterations=20
+#10000
+iterations=100
 
 
+###### Spherical harmonics ###########
+# If 0 then spherical harmonics are not calculated at all.
+spherical_harmonics_coefficients=21
+
 #shut up if we are using cluster!!!
 quiet=false
 

--
Gitblit v1.9.3