From 02d65c3c9baa0e76744ce001c55ec699ab3b956f Mon Sep 17 00:00:00 2001 From: Miha Fošnarič <miha.fosnaric@fe.uni-lj.si> Date: Mon, 11 Jul 2016 13:09:15 +0000 Subject: [PATCH] A fix in energy.c. Normals are summed now. --- src/tape | 77 ++++++++++++++++++++++++++++++++++---- 1 files changed, 69 insertions(+), 8 deletions(-) diff --git a/src/tape b/src/tape index 64f529a..16a5480 100644 --- a/src/tape +++ b/src/tape @@ -1,14 +1,50 @@ ####### Vesicle definitions ########### # nshell is a number of divisions of dipyramid -nshell=20 -# dmax is the square of the bond length -dmax=1.67 -# bending rigidity of the membrane -xk0=25.0 -# max step size +nshell=10 +# dmax is the max. bond length (in units l_min) +dmax=1.7 +# dmin_interspecies in the min. dist. between different vertex species (in units l_min) +dmin_interspecies=1.2 +# bending rigidity of the membrane (in units kT) +xk0=10.0 +# max step size (in units l_min) stepsize=0.15 +# Pressure calculations +# (pswitch=1: calc. p*dV energy contribution) +pswitch = 0 +# pressure difference: p_inside - p_outside (in units kT/l_min^3): +pressure=-10.0 +#Constant volume constraint (0 disable constant volume, 1 enable wiht additional vertex move, 2 enable with epsvol) +constvolswitch=0 +constvolprecision=1e-14 + +#Constant area constraint (0 disable constant area, 2 enable constant area with epsarea) +constareaswitch=0 + +####### Polymer (brush) definitions ########### +# npoly is a number of polymers attached to npoly distinct vertices on vesicle +npoly=0 +# nmono is a number of monomers in each polymer +nmono=2 +# Spring constant between monomers of the polymer +k_spring=800 + +####### Filament (inside the vesicle) definitions ########### +# nfil is a number of filaments inside the vesicle +nfil=0 +# nfono is a number of monomers in each filament +nfono=3 +# Persistence lenght of the filaments (in units l_min) +xi=100 + +####### Nucleus (inside the vesicle) ########### +# Radius of an impenetrable hard sphere inside the vesicle +R_nucleus=0 +R_nucleusX=0 +R_nucleusY=0 +R_nucleusZ=0 ####### Cell definitions ############ nxmax=60 nymax=60 @@ -16,11 +52,25 @@ ####### Program Control ############ -#how many iteration are there in a run? -iterations=20000 +#how many MC sweeps between subsequent records of states to disk +#200000 +mcsweeps=2000 +#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? +#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 + + #what type of multiprocessing? (*none, smp, cluster, distributed, cuda, auto) #currently only none makes sense. @@ -32,3 +82,14 @@ #max number of processes in distributed (voluntary) environment distributed_processes=50 #cuda options??? + + +#NirGov branch only! +#number of vertices with spontaneous curvature (integer) +number_of_vertices_with_c0=100 +#spontaneous curvature (float) +c0=0.5 +#energy of attraction of vertices with spontaneous curvature (float, positive value for attraction) +w=10.0 +#direct force on vesicles with spontaneous curvature (float) +F=2.0 -- Gitblit v1.9.3