parameters for hlc_graxi ######################################################################## # # # id Standard Initial Data File # # # # To construct otherwise identical ID files, run # # "mi id" which will create files idd0, idd1, etc # # which then differ only by the resolution. # # # ######################################################################## iter := 4000 output := 0-*/10 # # Matter Parameters # amp1 := 0.25 r01 := 0.0 delta1 := 1 # +1 --- in-going # 0 --- time-symmetric # -1 --- out-going signum1 := 0 rho_center1 := 0 z_center1 := 0.0 cyl1 := .005 #FP -- second piece of initial profile amp1_b := 0.0 r01_b := 0.0 delta1_b := 3.0 # +1 --- in-going # 0 --- time-symmetric # -1 --- out-going signum1_b := 0 rho_center1_b := 8.0 z_center1_b := 0.0 cyl1_b := 1 # # Geometry Parameters # # sigma ampsi := 0.00000 r0si := 8.0 deltasi := 2.0 signumsi := 0 rho_centersi := 0.0 z_centersi := 0.0 cylsi := 1 # omega ampo := 0.00000 r0o := 7.0 deltao := 3.0 signumo := 0 rho_centero := 0.0 z_centero := 0.0 cylo := 1 # # Tracing level: # # > 0 ---- Canonical CN loop output # > 1 ---- Canonical MG tracing # > 2 ---- Output relativized errors in CN loop # > 3 ---- parameter dumping # > 4 ---- level output in MG # trace_level := 1 # # General Parameters # MAXITER := 50 maxstep := 500 Nrho0 := 128 Nz0 := 256 rhomin := 0 rhomax := 40 zmin := -40 zmax := 40 ser := 0 fout := 1 epsiter := 1.0e-3 ic_epsiter := 1.0e-6 tag := "cyl_0_005_a025_ex_more_" lambda := .3 epsdis := 0.5 # # Specify form of Initial Data: # IDATA: # 0 ---- general Gaussians with parameters set below # 1 ---- specified ansatz, and computed source terms # 2 ---- load in scalar field profiles # IDATA := 0 # # Restrict model to Flat Space? # (if one, then metric restricted to flat metric) # flat_space := 0 # # force_cont_off == 1 do *not* use continuation methods in ID # != 1 do use them # force_cont_off := 0 # # Output initial time sdf's? (if non-zero) # output_init := 0 # # Solve by relaxation ? (if equal to one) # nsweeps -- max number of relaxation sweeps # nout -- sweep interval for output (according to trace_level) # solve_by_relax := 0 nsweeps := 5000 nout := 1000 # # Multigrid parameters # nvcycle := 1 preswp := 6 pstswp := 6 # # Excision parameters # do_excise := 1 exp_out_thresh := -0.0 ex_buffer := 2.00 ex_method := 2 ex_method_am := 2 ex_mg_buffer := 3 ex_test := 0 ex_in_buffer := 2.00 ex_start_t := 0 repop_extrap := 2 ex_fs_psi := 3 ex_fs_bar := 0 ex_fs_phi := 1 can_reduce_ex_zone := 2 ex_alpha_dbc := -1 ex_ri_psi := 0 ex_bz_dt := 0 ex_brho_dt := 0 ex_bz_max_s := 1 ex_brho_max_s := 1 ex_bz_s_am := 1 ex_brho_s_am := 1 ex_ri_iters_pcs := 2 ex_ri_csteps := 10 ex_ri_byrelax := 0 sph_omegabar := 0 extrapolate := 2 dbc_alpha := 1 dbc_betarho := 1 dbc_betaz := 1 dbc_psi := 0 ic_dbc_alpha := 1 ic_dbc_betarho := 1 ic_dbc_betaz := 1 ic_dbc_psi := 0 dbc_hh := 0 dbc_on_all4 := 0 dbc_on_n_finest := 0 b_sweep_x := 1 sym_diff := 1 # The following set of alternative constraint solver # parameters are (currently) only used in constraints.inc, # i.e. the initial conditions are still MG with relaxation # on coarsest. #bicg instead of MG on finest level use_bicg := 0 bicg_max_iter := 5000 bicg_eps_iter := 1.0e-5 bicg_newton_iters := 1 use_evo_for_psi := 1 compute_smallest_eig := 0 apph_search_dt := 0.5 apph_range := 1 apph_mxstep := 10 apph_find_z := 0 apph_find_a := 1 apph_num_c := 0 apph_min_l := 0 apph_apph_delta_A := 0.005 apph_apph_delta_f_a := 0.1 MG_w := 1.0 MG_w_scale := 0.5 MG_w_rescale := 1.025 MG_w_min := 0.05