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()