#
# AUTO constants file
#
# dictionary (mapping) of U(*) to user-defined names
unames = {1:'u', 2:'v'}
#
# dictionary (mapping) of PAR(*) to user-defined names
parnames = {1:'delta', 2:'lambda'}
#
NDIM = 1 # problem dimension
IPS = 1 # problem type: 0 = algebraic bifurcation problem
# 1 = stationary solutions of ODEs with detection of Hopf bifurcations
# -1 = fixed points of discrete dynamical system
# -2 = time integration using implicit Euler
# 2 = computation of periodic solutions
# 4 = boundary value problem
# 5 = algebraic optimization problems
# 7 = boundary value problem with computation of Floquet multipliers
# 9 = used in connection with HomCont algorithms
# 11 = spatially uniform solutions of a system of parabolic PDEs,
# with detection of traveling wave bifurcations
# 12 = Continuation of traveling wave solutions to a system of
# parabolic PDEs (reaction diffusion form).
# 14 = Time evolution for a system of parabolic PDEs subject to
# periodic boundary conditions.
# 15 = Optimization of periodic solutions.
# 16 = This option is similar to IPS=14, except that the user supplies
# the boundary conditions.
#
#
IRS = 0 # start solution label
#TY # start solution type
ILP = 1 # fold detection: 0 = off
# 1 = on
#
ICP = [1,2] # continuation parameters: list of indices j for par(j)
#
NTST = 50 # number of mesh intervals
NCOL = 4 # number of collocation points
IAD = 3 # mesh adaption every IAD steps: 0 = off
# > 0 = on
ISP = 2 # bifurcation detection: 0 = off
# 1 = BP for alg. probl. and comp. of Floq. mult.
# 2 = all (note one Floq. mult. must be near 1. Sec. 10.7)
# 3 = BP, no Floq. mult.
#
ISW = 1 # branch switching: 1 = normal
# -1 = switch branch (BP, HB, PD)
# 2 = switch to two-parameter continuation (LP, symm. BP, HB, TR)
# 3 = switch to three-parameter continuation (BP)
#
IPLT = 0 # select principal solution measure. Sec. 10.9.10
NBC = 0 # number of boundary conditions
NINT = 0 # number of integral conditions
#
NMX = 100 # maximum number of steps
RL0 = -10 # parameter interval: RL0 <= lambda <= RL1
RL1 = 100
A0 = 0 # interval of principal solution measure: A0 <= |.| <= A1
A1 = 1000
#
NPR = 20 # print and save restart data every NPR steps
MXBF = 2 # if IPS = 0,1: automatic branch switching for the first MXBF bifurcation points
IID = 3 # control diagnostic output; 0 = none
# 1 = little
# 2 = normal
# 3 = additional for algebraic problem
# 4 = additional for ODE
# 5 = extensive
#
ITMX = 8 # maximum number of iterations for locating special solutions/points
ITNW = 5 # maximum number of correction steps
NWTN = 3 # corrector uses full Newton for NWTN steps
JAC = 0 # user defines derivatives: 0 = no
# 1 = yes
#
EPSL = 1.e-10 # convergence criterion: parameters
EPSU = 1.e-10 # solution components
EPSS = 1.e-8 # special points
#
DS = -1.e-1 # start step size
DSMIN = 1.e-9 # step size interval DSMIN <= h <= DSMAX
DSMAX = 5.e-1
IADS = 1 # Step size adaption every IADS steps: 0=off
# > 0 = on
#
THL={11:0.} # list of parameter weights
THU={} # list of solution weights
#
UZR = {} # list of parameter index:values for runtime output
UZSTOP = {} # list of parameter index:values for stopping condition
SP = {} # list of bifurcations to check and bifurcation stop conditions. Sec. 10.8.2
STOP = {} # list of bifurcation/UZ stop conditions. Sec. 10.6.1
#
# Define file names: e, s, dat, sv
# other switches: IBR (sec. 10.9.6), LAB (sec. 10.9.7), IIS (sec. 10.9.8)
# HOMCONT: NUNSTAB, NSTAB, IEQUIB, ITWIST, ISTART, IREV, IFIXED, IPSI