Numerical Parameters
The &num namelist group defines numerical method parameters; the
input file can contain one or more, but only the last (tag-matching) one is
used. Allowable fields are:
diff_scheme(default'COLLOC_GL2')Difference equation scheme; one of:
'COLLOC_GL2': Second-order Gauss-Legendre collocation'COLLOC_GL4': Fourth-order Gauss-Legendre collocation'COLLOC_GL6': Sixth-order Gauss-Legendre collocation'MAGNUS_GL2': Second-order Gauss-Legendre Magnus'MAGNUS_GL4': Fourth-order Gauss-Legendre Magnus'MAGNUS_GL6': Sixth-order Gauss-Legendre Magnus'MIRK': Fourth-order mono-implicit Runge-Kutta (experimental)'TRAPZ': Trapezoidal, with the prescription by Sugimoto (1970) for non-adiabatic cases
r_root_solver(default'BRENT')Root solver for real arithmetic; one of:
'BRENT': Brent’s method
c_root_solver(default'RIDDERS')Root solver for complex arithmetic; one of
'RIDDERS': Complex Ridders’ method'SECANT': Secant method'SIMPLEX': Simplex method
n_iter_max(default50)Maximum number of iterations in root-finding algorithm
matrix_type(default'BLOCK’)Storage type of system matrix; one of
'BAND': Band-structured'BLOCK': Block-structured
deflate_roots(default.TRUE.)Flag to use root deflation, which can avoid the same eigenfrequency being found multiple times
restrict_roots(default.TRUE.)Flag to check each roots found lies within the bounds of the frequency scan
ad_search(default'BRACKET')Initial search method for adiabatic calculations; one of
'BRACKET': Bracket sign changes in the discriminant function
nad_search(default'AD')Initial search method for non-adiabatic calculations; one of
'AD': Use adiabatic eigenfrequencies'MINMOD': Find minima in the modulus of the discriminant function, along the real-\(\omega\) axis'CONTOUR': Find intersections between real and imaginary zero-contours of the discriminant function
See the Non-Adiabatic Oscillations chapter for more details about these search methods.
tag_list(default'', which matches all)Comma-separated list of
&modetags to match