From 41637ae4c538b5f11b41dacf6f941479e7a5ff00 Mon Sep 17 00:00:00 2001 From: Samo Penic <samo.penic@fe.uni-lj.si> Date: Fri, 06 Jun 2014 09:48:41 +0000 Subject: [PATCH] CV bug fix. Changes in constant volume: verifying the neightbour distances and removed preliminary contraints check. Seems to work now. --- src/tape | 16 +++++++++++----- 1 files changed, 11 insertions(+), 5 deletions(-) diff --git a/src/tape b/src/tape index c1f23e1..cedfc46 100644 --- a/src/tape +++ b/src/tape @@ -12,9 +12,14 @@ # 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=0.0 + +#Constant volume constraint (0 disable constant volume, 1 enable) +constvolswitch=0 +constvolprecision=1e-14 + ####### Polymer (brush) definitions ########### # npoly is a number of polymers attached to npoly distinct vertices on vesicle @@ -26,7 +31,7 @@ ####### 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) @@ -34,7 +39,7 @@ ####### Nucleus (inside the vesicle) ########### # Radius of an impenetrable hard sphere inside the vesicle -R_nucleus=5 +R_nucleus=0 ####### Cell definitions ############ nxmax=60 @@ -44,14 +49,15 @@ ####### Program Control ############ #how many MC sweeps between subsequent records of states to disk -mcsweeps=100 +mcsweeps=20 #how many initial mcsweeps*inititer MC sweeps before recording to disk? inititer=0 #how many records do you want on the disk iteration are there in a run? -iterations=500 +iterations=10 ###### Spherical harmonics ########### +# If 0 then spherical harmonics are not calculated at all. spherical_harmonics_coefficients=21 #shut up if we are using cluster!!! -- Gitblit v1.9.3