jg.runtime.rexi_ci_max_imag = 2
jg.runtime.rexi_ci_mu = 0
jg.runtime.rexi_ci_primitive = 'circle'
elif 'etdrk' in jg.runtime.timestepping_method:
# CI REXI via file input
jg.runtime.rexi_method = 'file'
cirexi = CIREXI()
coeffs_phi0 = cirexi.setup("phi0", N=32, R=2).toFloat()
coeffs_phi1 = cirexi.setup("phi1", N=32, R=2).toFloat()
coeffs_phi2 = cirexi.setup("phi2", N=32, R=2).toFloat()
coeffs_phi3 = cirexi.setup("phi3", N=32, R=2).toFloat()
jg.runtime.rexi_files_coefficients = [
coeffs_phi0, coeffs_phi1, coeffs_phi2, coeffs_phi3
]
else:
# B REXI via file
jg.runtime.rexi_method = 'file'
brexi = BREXI()
coeffs = brexi.setup(N=8, quadrature_method='gauss').toFloat()
jg.runtime.rexi_files_coefficients = [coeffs]
else:
jg.runtime.rexi_method = None
jg.gen_jobscript_directory()
# Validate with C-implementation of CI-REXI method
jg.runtime.rexi_method = "ci"
jg.runtime.rexi_ci_n = N
jg.runtime.rexi_ci_max_real = lambda_max_real
jg.runtime.rexi_ci_max_imag = lambda_include_imag
jg.gen_jobscript_directory()
# Back to original version
jg.runtime.rexi_method = "file"
elif rexi_file_method == "brexi":
if function_name != "phi0":
continue
brexi = BREXI(efloat_mode=efloat_mode)
# CI-REXI: Number of quadrature poles
N_list = [8, 10, 12, 14, 16]
#N_list = [8, 10]
quadrature_method_list = ["gauss", "radau", "chebyshev"]
for (N, quadrature_method) in product(N_list,
quadrature_method_list):
coeffs = brexi.setup(
N=N, quadrature_method=quadrature_method).toFloat()
jg.runtime.rexi_files_coefficients = [coeffs]
jg.gen_jobscript_directory()
else: