# Documentation - Wiki:
# https://gitlab.fizyka.pw.edu.pl/wtools/wslda/-/wikis/home
# -------------- USER DEFINED PARAMETERS --------------
# See: Wiki -> User defined parameters
# Data flow: [Read params from input file] -> [execute process_params( )] -> [pass params to functions]
params[1] = 0.80 # Thomas-Fermi radius for x direction, in units of LX/2
params[2] = 0.90 # Smoothing parameter x1 for harmonic_oscillator_smooth_edges, in units of L/2
params[3] = 0.97 # Smoothing parameter x2 for harmonic_oscillator_smooth_edges, in units of L/2
params[4] = 0.24 # height of the barrier, in eF units
params[5] = 5.0 # width of barrier, in kF^1 units
params[6] = 0.0 # no tilt
# ------------------- INITIALIZATION ------------------
# See: Wiki -> Initialization of the solver
inittype 1 # 0 - create uniform solution and evolve it,
# 10 - read uniform solution from file `inprefix`/uniform.solution
#
# 1 - start from st-wslda-1d solution, inprefix points to folder with s1dpca binary files
# 2 - start from st-wslda-2d solution, inprefix points to folder with s2dpca binary files
# 3 - start from st-wslda-3d solution, inprefix points to folder with s3dpca binary files
#
# 5 - start from td-wslda checkpoint, inprefix points to folder with checkpoint binary files
# ------------------- INPUT/OUTPUT ------------------
outprefix jj-r2 # all output files will start with this prefix
inprefix ../st-jj-example/jj-r1 # prefix of calculation you need to use as input (when inittype!=0)
overwrite 1 # 1 - allow for overwriting existing files, 0 - do not allow for overwriting existing files
iogroups 8 # number of IO groups used for parallel reading of data
# it must be the same as the value used for static calculations.
checkpoint 1 # do checkpoint at the end of calculation: 1-yes, 0-no
checkperiod 0 # periodicity of the checkpoint, meaningful only if checkpoint=1
# if checkperiod>0 then after each checkperiod measurements, the checkpoint will be conducted
# if checkperiod=0 the option is deactivated (default)
walltime 11.9 # after this time force the checkpoint [hr]
# list of variables for writing
# all = rho delta j nu tau V V_ext delta_ext velocity_ext alpha A
# default = rho delta j
writevar default V_ext
writeprec d # precision for data writing: d-double(default), f-float.
# ----------------- SCATTERING LENGTH ---------------
sclgth -10000.0 # scattering length in units of lattice spacing
# meaningful only for FUNCTIONAL=BDG,SLDAE
# in case of FUNCTIONAL=(A)SLDA it is set automatically to infinity
#akF -1.0 # scattering length times Fermi wave vector
# this tag modifies the value of tag sclgth via process_params(...) function
# by default is inactive (akF=0)
# ------------------- TIME EVOLUTION ----------------
dt 0.0035 # integration time step, in units of eF^-1
measurements 200 # number of requested measurements
timesteps 1430 # number of time steps between each measurement
# measurements are written with time resolution: timesteps*dt
# the total trajectory length is: measurements*timesteps*dt
selfstart 1 # use Taylor expansion for first 5 steps?
# See: Wiki -> Starting ABM algorithm
# ------------------ GPUS THREADING SET --------------
# nthreads 512 # number of parallel threads executing by each GPU, default=512
# batch 1000000 # maximal number of wave-functions processed in parallel by each card, default=1000000
# gpuspernode 1 # number of GPUs per node (resource set), default=GPUS_PER_NODE (in machine.h)
# -------------- HIGH K-WAVES FILTER ----------------
# See: Wiki -> Stabilization of the time-dependent code
hkf_mode 1 # 0 - no noise filtering (default)
# 1 - noise filtering for mean-fields only
# 2 - noise filtering of mean-fields and pairing field
hkf_mu 0.9 # mu parameter of the Fermi-Dirac (filtering) function, in Ec units, default=0.9
hkf_T 0.02 # T parameter of the Fermi-Dirac (filtering) function, in Ec units, default=0.02
# ---------------- CONSERVATION LAWS -----------------
# See: Wiki -> Monitoring of conservation laws
# # -> Energy conservation monitoring
Econservation_start 0 # Start to monitor energy conservation from time*eF, default=1e12 (infinity)
Econservation_stop 10000 # Stop to monitor energy conservation at time*eF, default=1e12 (infinity)
Econservation_tol 0.05 # if |[E(t)-E(start_t)]/E(start_t)|>tol then the code will terminate, default=0.05
# # -> Particle number conservation monitoring
Nconservation_start 0 # Start to monitor total particle number conservation from time*eF, default=0
Nconservation_stop 10000 # Stop to monitor total particle number conservation at time*eF, default=1e12 (infinity)
Nconservation_tol 0.05 # if |[N(t)-N(start_t)]/N(start_t)|>tol then the code will terminate, default=0.05
# --------------- QUANTUM FRICTION ------------------
# See: Wiki -> Quantum friction
# qfalpha 0.0 # alpha parameter for quantum friction term (cooling in the normal channel), default=0
# qfbeta 0.0 # beta parameter for quantum friction term (cooling in the pairing channel), default=0
# qfgamma 0.0 # particle control parameter (cooling in the pairing channel), default=0
# qfNreq 0.0 # requested particle number, meaningful only if qfgamma>0
# qfstart 0.0 # start time for evolving with quantum friction [eF], default=0
# qfstop 0.0 # stop time for evolving with quantum friction [eF], default=infinity
# qfswitch 0.0 # the friction will be activated and deactivated gradually over this period of time [eF], default=1
# ---------------- SUBSET TRACKING ------------------
# See: Wiki -> Tracking of selected states
# observables (densities, currents) arising from state in the energy interval En in [subsetMinEn,subsetMaxEn]
# will be computed, if subsetMinEn=subsetMaxEn=0 the functionality is disables
# subsetMinEn 0.0 # in eF units, deafault=0.0
# subsetMaxEn 0.0 # in eF units, deafault=0.0
# subsetShiftDmu 1 # if 1 then apply extra shift of quasiparticle energies by (mu_a-mu_b)/2, default=0
# ----------- IN CASE OF: inittype=0 ----------------
# See: Wiki -> Initialization of the solver
# init0Na 28.0 # requested number of particles of type a, default=Na
# init0Nb 28.0 # requested number of particles of type b, default=Nb
init0na 0.01689 # ...alternatively, you can request the desired density for particles of type a
init0nb 0.01689 # ...alternatively, you can request the desired density for particles of type b
init0muchange 0.1 # Change rate of chemical potential
init0Tstart 0.0 # Start temperature, in units of eF, default 0.2
init0Tstop 0.0 # Stop temperature, in units of eF, default 0.05
init0DeltaT 0.01 # change of temperature in units of eF, default 0.01
init0eps 1.0e-9 # epsilon for convergence, default 1.0e-6
init0scmix 0.5 # mixing parameter in self-consistent process, default 0.25
init0maxiter 2000 # maximum number of iterations, default 100000
# init0kc 3.1415 # momentum cut-off used for uniform solver, default=pi/max(DX,DY,DZ)
init0debug 1 # debug level, default=0 (no debug info),
# init0save 1 # save solution
Last edited by 
Gabriel Wlazłowski
Gabriel Wlazłowski