jg.runtime.viscosity = 0.0
jg.runtime.rexi_method = ''
# Update TIME parallelization
ptime = JobParallelizationDimOptions('time')
ptime.num_cores_per_rank = 1
ptime.num_threads_per_rank = 1 #pspace.num_cores_per_rank
ptime.num_ranks = 1
pspace = JobParallelizationDimOptions('space')
pspace.num_cores_per_rank = 1
pspace.num_threads_per_rank = params_pspace_num_cores_per_rank[-1] / 4
pspace.num_ranks = 1
jg.reference_job = False
# Setup parallelization
jg.setup_parallelization([pspace, ptime])
if verbose:
pspace.print()
ptime.print()
jg.parallelization.print()
unique_id_filter = []
#unique_id_filter.append('runtime.disc_space')
unique_id_filter.append('runtime.simparams')
unique_id_filter.append('runtime.reuse_plans')
#unique_id_filter.append('runtime.benchmark')
#unique_id_filter.append('runtime.timestepping')
print("Reference")
tsm = ts_methods[0]
jg.parallelization.max_wallclock_seconds = ref_max_wallclock_seconds
SetupSpectralMethods(jg)
jg.runtime.timestep_size = 2 # second #jg.runtime.output_timestep_size/100.0
jg.runtime.timestepping_method = tsm[0]
jg.runtime.timestepping_order = tsm[1]
jg.runtime.timestepping_order2 = tsm[2]
jg.runtime.space_res_physical = -1
#jg.runtime.space_res_spectral = 1024
jg.runtime.space_res_spectral = phys_res_list[0]
# Tag this as a reference job
jg.reference_job = True
jg.gen_jobscript_directory()
jg.reference_job = False
# Use this one as the reference solution!
jg.reference_job_unique_id = jg.job_unique_id
#
# Use only 2 iterations for Semi-Lagrangian methods
#
unique_id_filter.append('runtime.semi_lagrangian')
jg.runtime.semi_lagrangian_iterations = 2
jg.runtime.semi_lagrangian_convergence_threshold = -1
jg.parallelization.max_wallclock_seconds = max_wallclock_seconds
jg.compile.plane_spectral_space = 'enable' jg.compile.plane_spectral_dealiasing = 'enable' #Setup method jg.runtime.timestepping_method = 'ln_erk' jg.runtime.timestepping_order = 4 jg.runtime.timestepping_order2 = 4 jg.runtime.space_res_physical = -1 #Setup space info jg.runtime.space_res_spectral = 64 #Viscosity jg.runtime.viscosity = 0.0 #Banchmark to be used naming: waves_N_k0_k1_d0_dwest_deast_k0_k1_d0_dwest_deast_k0_k1_d0_dwest_deast #jg.runtime.benchmark_normal_modes_case ="waves_1_3_2_0_100_0" #jg.runtime.benchmark_normal_modes_case ="waves_2_1_3_0_100_0_7_2_100_0_0" #jg.runtime.benchmark_normal_modes_case ="waves_1_1_3_100_0_0" jg.runtime.benchmark_normal_modes_case = "waves_1_1_2_0_1_0" # Tag this as a reference job jg.reference_job = True jg.gen_jobscript_directory() # Use this one as the reference solution! jg.reference_job_unique_id = jg.job_unique_id # Write compile script jg.write_compilecommands()
