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 := 16000 output := 0-*/20 # # Matter Parameters # amp1 := 0.45 r01 := 0.0 delta1 := 3.0 # +1 --- in-going # 0 --- time-symmetric # -1 --- out-going signum1 := 20 rho_center1 := 0 z_center1 := 20.0 cyl1 := 1 ic_v1 := -0.05 #FP -- second piece of initial profile amp1_b := 0.45 r01_b := 0.0 delta1_b := 3.0 # +1 --- in-going # 0 --- time-symmetric # -1 --- out-going signum1_b := 20 rho_center1_b := 0 z_center1_b := -20.0 cyl1_b := 1 ic_v1_b := 0.05 # # 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 f_a1 := 1.0 f_a2 := 1.0 f_m1 := 1 f_m2 := 1 f_v1 := 0.0 f_v2 := 0.0 f_z01 := 21 f_z02 := -21 # # 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 := 250 Nrho0 := 256 Nz0 := 512 rhomin := 0 rhomax := 80 zmin := -80 zmax := 80 ser := 0 fout := 1 epsiter := 1.0e-4 ic_epsiter := 1.0e-8 tag := "bf_" lambda := .30 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 := 50000 nout := 1000 # # Multigrid parameters # nvcycle := 1 preswp := 3 pstswp := 3 # # Excision parameters # do_excise := 1 exp_out_thresh := -0.0 ex_buffer := 30 ex_method := 32 ex_method_am := 32 ex_mg_buffer := 3 ex_test := 0 ex_in_buffer := 2.0 ex_start_t := 15.0 repop_extrap := 2 ex_fs_psi := 1 ex_fs_phi := 1 ex_fs_bar := 0 can_reduce_ex_zone := 1 ex_alpha_dbc := 0.2 ex_ri_psi := 1 ex_bz_dt := 40 ex_brho_dt := 0 ex_bz_max_s := 1.5 ex_brho_max_s := 1 ex_bz_s_am := 0.5 ex_brho_s_am := 0.5 ex_mf_all := 1 ex_ri_iters_pcs := 2 ex_ri_csteps := 10 ex_ri_byrelax := 0 sph_omegabar := 0 extrapolate := 0 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 apph_mxstep := 10 apph_search_dt := 0.25 apph_find_z := -1 apph_find_a := 1 apph_range := 10 apph_delta_f_z := 0.5 apph_num_c := 0 # 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