Commits May 2011

commits@lists.einsteintoolkit.org

6 participants
21 discussions

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 248)
by bcmsma＠astro.rit.edu 23 May '11

23 May '11

User: bmundim Date: 2011/05/23 05:31 PM Modified: /trunk/src/ GRHydro_Con2PrimM_pt.c Log: Use CCTK_REAL instead of double. File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2PrimM_pt.c Delta lines: +2 -2 =================================================================== --- trunk/src/GRHydro_Con2PrimM_pt.c 2011-05-09 01:03:12 UTC (rev 247) +++ trunk/src/GRHydro_Con2PrimM_pt.c 2011-05-23 22:31:18 UTC (rev 248) @@ -680,12 +680,12 @@ **********************************************************************/ static CCTK_REAL eos_info(CCTK_REAL W, CCTK_REAL vsq, CCTK_REAL *dpdw, CCTK_REAL *dpdvsq, CCTK_REAL gammaeos, struct LocGlob *lgp) { - register double ftmp,gtmp; + register CCTK_REAL ftmp,gtmp; ftmp = 1. - vsq; gtmp = sqrt(ftmp); - double gam_m1_o_gam = ((gammaeos-1.)/gammaeos); + CCTK_REAL gam_m1_o_gam = ((gammaeos-1.)/gammaeos); *dpdw = gam_m1_o_gam * ftmp ; *dpdvsq = gam_m1_o_gam * ( 0.5 * lgp->D/gtmp - W ) ;

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/examples/bbh/ (Rev. 91)
by bcmsma＠astro.rit.edu 22 May '11

22 May '11

User: bmundim Date: 2011/05/22 05:05 PM Added: /examples/bbh/ BBHHighRes.par, BBHMedRes.par Modified: /examples/bbh/ BBHLowRes.par Log: Increase number of collocation points in TwoPunctures. Add higher resolutions parameter files. File Changes: Directory: /examples/bbh/ ========================= File [added]: BBHHighRes.par Delta lines: <Unavailable> File [modified]: BBHLowRes.par Delta lines: <Unavailable> File [added]: BBHMedRes.par Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] www/about/ (Rev. 590)
by bcmsma＠astro.rit.edu 21 May '11

21 May '11

User: bmundim Date: 2011/05/21 07:31 PM Modified: /about/ index.php Log: Oops, I missed one GetComponent link... File Changes: Directory: /about/ ================== File [modified]: index.php Delta lines: +1 -1 =================================================================== --- about/index.php 2011-05-22 00:29:02 UTC (rev 589) +++ about/index.php 2011-05-22 00:31:07 UTC (rev 590) @@ -62,7 +62,7 @@ the source code and associated tools for simulations, including Cactus thorns. The component list is written using the Component Retrieval Language and can be checked out using - the <a href="https://svn.cactuscode.org/Utilities/branches/ET_2011_05/Scripts/GetCompone…" + the <a href="https://github.com/gridaphobe/CRL/raw/ET_2011_05/GetComponents" class="nostyle">GetComponents</a> tools. Distribution: <a href="https://github.com/gridaphobe/CRL/raw/ET_2011_05/GetComponents">https://github.com/gridaphobe/CRL/raw/ET_2011_05/GetComponents</a>

1 0

[svn:einsteintoolkit] www/about/ (Rev. 589)
by bcmsma＠astro.rit.edu 21 May '11

21 May '11

User: bmundim Date: 2011/05/21 07:29 PM Modified: /about/ index.php Log: Fix link to GetComponents. File Changes: Directory: /about/ ================== File [modified]: index.php Delta lines: +1 -2 =================================================================== --- about/index.php 2011-05-12 19:06:04 UTC (rev 588) +++ about/index.php 2011-05-22 00:29:02 UTC (rev 589) @@ -64,8 +64,7 @@ Retrieval Language and can be checked out using the <a href="https://svn.cactuscode.org/Utilities/branches/ET_2011_05/Scripts/GetCompone…" class="nostyle">GetComponents</a> tools. -Distribution: <a href="https://svn.einsteintoolkit.org/manifest/branches/ET_2011_05/einsteintoolki…">https://svn.einsteintoolkit.org/manifest/branches/ET_2011_05/einsteintoolki…</a> -<a href="https://svn.cactuscode.org/Utilities/branches/ET_2011_05/Scripts/GetCompone…">https://svn.cactuscode.org/Utilities/branches/ET_2011_05/Scripts/GetCompone…</a> +Distribution: <a href="https://github.com/gridaphobe/CRL/raw/ET_2011_05/GetComponents">https://github.com/gridaphobe/CRL/raw/ET_2011_05/GetComponents</a> </td> </tr>

1 0

[svn:einsteintoolkit] TwoPunctures/trunk/ (Rev. 123)
by knarf＠cct.lsu.edu 20 May '11

20 May '11

User: knarf Date: 2011/05/20 06:02 AM Modified: /trunk/ interface.ccl Log: don't prolongate puncture_u File Changes: Directory: /trunk/ ================== File [modified]: interface.ccl Delta lines: +1 -1 =================================================================== --- trunk/interface.ccl 2011-04-21 04:47:02 UTC (rev 122) +++ trunk/interface.ccl 2011-05-20 11:02:32 UTC (rev 123) @@ -4,7 +4,7 @@ INHERITS: ADMBase StaticConformal grid -REAL puncture_u TYPE=gf +REAL puncture_u TYPE=gf tags='prolongation="none"' CCTK_INT FUNCTION Set_Rho_ADM( \ CCTK_POINTER_TO_CONST IN cctkGH, \

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/examples/bbh/ (Rev. 90)
by bcmsma＠astro.rit.edu 19 May '11

19 May '11

User: bmundim Date: 2011/05/19 02:02 PM Modified: /examples/bbh/ BBHLowRes.par Log: WeylScal4::fd_order = "4th" and more points in Multipole. File Changes: Directory: /examples/bbh/ ========================= File [modified]: BBHLowRes.par Delta lines: +6 -5 =================================================================== --- examples/bbh/BBHLowRes.par 2011-05-18 14:22:33 UTC (rev 89) +++ examples/bbh/BBHLowRes.par 2011-05-19 19:02:03 UTC (rev 90) @@ -98,8 +98,8 @@ # Utility parameters: #------------------------------------------------------------------------------ -#NaNChecker::check_every = 16 -NaNChecker::check_every = 1 +NaNChecker::check_every = 16 +#NaNChecker::check_every = 1 NaNChecker::check_vars = " ML_BSSN::ML_log_confac ML_BSSN::ML_metric @@ -505,6 +505,7 @@ IOHDF5::out2D_yz = "no" IOHDF5::out2D_xyplane_z = 0.0 IOHDF5::out2D_vars = " + CarpetReduce::weight Grid::coordinates ADMBase::lapse ADMBase::metric @@ -514,7 +515,6 @@ WeylScal4::Psi4r WeylScal4::Psi4i " -# CarpetReduce::weight #IOHDF5::out3D_ghosts = "no" @@ -613,12 +613,13 @@ Multipole::radius[1] = 40 Multipole::radius[2] = 50 Multipole::radius[3] = 60 -Multipole::ntheta = 60 -Multipole::nphi = 120 +Multipole::ntheta = 120 +Multipole::nphi = 240 Multipole::variables = "WeylScal4::Psi4r{sw=-2 cmplx='WeylScal4::Psi4i' name='psi4'}" Multipole::out_every = 4 Multipole::l_max = 4 +WeylScal4::fd_order = "4th" #------------------------------------------------------------------------------ # Checkpoint/Recovery:

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/examples/tov/ (Rev. 89)
by knarf＠cct.lsu.edu 18 May '11

18 May '11

User: knarf Date: 2011/05/18 09:22 AM Modified: /examples/tov/ tov.rpar Log: rename parameter File Changes: Directory: /examples/tov/ ========================= File [modified]: tov.rpar Delta lines: +1 -1 =================================================================== --- examples/tov/tov.rpar 2011-05-18 14:07:46 UTC (rev 88) +++ examples/tov/tov.rpar 2011-05-18 14:22:33 UTC (rev 89) @@ -192,7 +192,7 @@ ADMBase::initial_dtshift = "zero" #TOVSolver::TOV_Rho_Central[0] = 1.28e-3 -TOVSolver::Mass[0] = 1.4 +TOVSolver::TOV_Mass[0] = 1.4 TOVSolver::TOV_Gamma[0] = 2.0 TOVSolver::TOV_K[0] = $K_TOV $K_change

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/examples/tov/ (Rev. 88)
by knarf＠cct.lsu.edu 18 May '11

18 May '11

User: knarf Date: 2011/05/18 09:07 AM Removed: /examples/tov/ tov_pert.rpar Modified: /examples/tov/ tov.rpar Log: combine rpars File Changes: Directory: /examples/tov/ ========================= File [modified]: tov.rpar Delta lines: +3 -2 =================================================================== --- examples/tov/tov.rpar 2011-05-18 13:46:36 UTC (rev 87) +++ examples/tov/tov.rpar 2011-05-18 14:07:46 UTC (rev 88) @@ -8,7 +8,7 @@ my $K_TOV = 100; my $K_change = ""; if ($type eq "_pert") { - $K_TOV = 99; + $K_TOV = 90; $K_change = "TOVSolver::TOV_Transform_to_K = 100"; } @@ -191,7 +191,8 @@ ADMBase::initial_dtlapse = "zero" ADMBase::initial_dtshift = "zero" -TOVSolver::TOV_Rho_Central[0] = 1.28e-3 +#TOVSolver::TOV_Rho_Central[0] = 1.28e-3 +TOVSolver::Mass[0] = 1.4 TOVSolver::TOV_Gamma[0] = 2.0 TOVSolver::TOV_K[0] = $K_TOV $K_change File [removed]: tov_pert.rpar Delta lines: +0 -242 =================================================================== --- examples/tov/tov_pert.rpar 2011-05-18 13:46:36 UTC (rev 87) +++ examples/tov/tov_pert.rpar 2011-05-18 14:07:46 UTC (rev 88) @@ -1,242 +0,0 @@ -#!/usr/bin/perl -w -use strict; - -foreach my $res ((8, 4, 2)) { - - my $out_every = 256/$res; - my $par = <<PAR_EOF; -## mclachlan tov_static -ActiveThorns = "Time MoL" -ActiveThorns = "Coordbase CartGrid3d Boundary StaticConformal" -ActiveThorns = "SymBase ADMBase TmunuBase HydroBase InitBase ADMCoupling ADMMacros" -ActiveThorns = "IOUtil" -#ActiveThorns = "Formaline" -ActiveThorns = "SpaceMask CoordGauge Constants LocalReduce aeilocalinterp LoopControl" -ActiveThorns = "Carpet CarpetLib CarpetReduce CarpetRegrid2 CarpetInterp" -ActiveThorns = "CarpetIOASCII CarpetIOScalar CarpetIOHDF5 CarpetIOBasic" -ActiveThorns = "SystemStatistics" - -ActiveThorns = "ADMConstraints" - -# grid parameters -CartGrid3D::type = "coordbase" -CartGrid3D::domain = "full" -CartGrid3D::avoid_origin = "no" -CoordBase::xmin = 0.0 -CoordBase::ymin = 0.0 -CoordBase::zmin = 0.0 -CoordBase::xmax = 240.0 -CoordBase::ymax = 240.0 -CoordBase::zmax = 240.0 -CoordBase::dx = $res -CoordBase::dy = $res -CoordBase::dz = $res - -CoordBase::boundary_size_x_lower = 3 -CoordBase::boundary_size_y_lower = 3 -CoordBase::boundary_size_z_lower = 3 -CoordBase::boundary_size_x_upper = 3 -CoordBase::boundary_size_y_upper = 3 -CoordBase::boundary_size_z_upper = 3 -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_y_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 -CoordBase::boundary_shiftout_x_upper = 0 -CoordBase::boundary_shiftout_y_upper = 0 -CoordBase::boundary_shiftout_z_upper = 0 - - - - -ActiveThorns = "ReflectionSymmetry" - -ReflectionSymmetry::reflection_x = "yes" -ReflectionSymmetry::reflection_y = "yes" -ReflectionSymmetry::reflection_z = "yes" -ReflectionSymmetry::avoid_origin_x = "no" -ReflectionSymmetry::avoid_origin_y = "no" -ReflectionSymmetry::avoid_origin_z = "no" - - - -TmunuBase::stress_energy_storage = yes -TmunuBase::stress_energy_at_RHS = yes -TmunuBase::timelevels = 1 -TmunuBase::prolongation_type = none - -HydroBase::timelevels = 3 - -ADMMacros::spatial_order = 4 - -ADMBase::metric_type = "physical" - -ADMConstraints::bound = "static" -ADMConstraints::constraints_timelevels = 3 -ADMConstraints::constraints_persist = yes - -SpaceMask::use_mask = "yes" - -Cactus::terminate = "time" -Cactus::cctk_final_time = 1000 - -Carpet::domain_from_coordbase = "yes" -Carpet::enable_all_storage = no -Carpet::use_buffer_zones = "yes" - -Carpet::poison_new_timelevels = "yes" -Carpet::check_for_poison = "no" -Carpet::poison_value = 113 - -Carpet::init_3_timelevels = no -Carpet::init_fill_timelevels = "yes" - -CarpetLib::poison_new_memory = "yes" -CarpetLib::poison_value = 114 - -# system specific Carpet paramters -Carpet::max_refinement_levels = 10 -driver::ghost_size = 3 -Carpet::prolongation_order_space = 3 -Carpet::prolongation_order_time = 2 - -CarpetRegrid2::regrid_every = 0 -CarpetRegrid2::num_centres = 1 - -CarpetRegrid2::num_levels_1 = 5 -CarpetRegrid2::radius_1[1] =120.0 -CarpetRegrid2::radius_1[2] = 60.0 -CarpetRegrid2::radius_1[3] = 30.0 -CarpetRegrid2::radius_1[4] = 15.0 - - -time::dtfac = 0.25 -MoL::ODE_Method = "rk4" -MoL::MoL_Intermediate_Steps = 4 -MoL::MoL_Num_Scratch_Levels = 1 - -# check all physical variables for NaNs -ActiveThorns = "NaNChecker" -#NaNChecker::check_every = 1 -#NaNChecker::action_if_found = "just warn" #"terminate", "just warn", "abort" -#NaNChecker::check_vars = "ADMBase::metric ADMBase::lapse ADMBase::shift HydroBase::rho HydroBase::eps HydroBase::press HydroBase::vel" - -## Lapse Condition: \\partial_t alpha = - alpha K -## Shift Condition: \\partial_t beta^i = 0 - -# Hydro paramters -ActiveThorns = "EOS_Omni" -ActiveThorns = "GRHydro" - -HydroBase::evolution_method = "GRHydro" - -GRHydro::riemann_solver = "Marquina" -GRHydro::GRHydro_eos_type = "Polytype" -GRHydro::GRHydro_eos_table = "2D_Polytrope" -GRHydro::recon_method = "ppm" -GRHydro::GRHydro_stencil = 3 -GRHydro::bound = "none" -GRHydro::rho_abs_min = 1.e-10 -#GRHydro::GRHydro = 18 # Tmunu(10), rho,press,eps,w_lorentz,vel, tau -#GRHydro::GRHydro = 10 # gij(6), alpha, beta(3) - - -ActiveThorns = "GenericFD NewRad" -ActiveThorns = "ML_BSSN ML_BSSN_Helper" -ADMBase::evolution_method = "ML_BSSN" -ADMBase::lapse_evolution_method = "ML_BSSN" -ADMBase::shift_evolution_method = "ML_BSSN" -ADMBase::dtlapse_evolution_method= "ML_BSSN" -ADMBase::dtshift_evolution_method= "ML_BSSN" - -TmunuBase::support_old_CalcTmunu_mechanism = "no" - -ML_BSSN::timelevels = 3 - -ML_BSSN::harmonicN = 1 # 1+log -ML_BSSN::harmonicF = 1.0 # 1+log -ML_BSSN::LapseACoeff = 1.0 -ML_BSSN::ShiftBCoeff = 1.0 -ML_BSSN::ShiftGammaCoeff = 0.0 -ML_BSSN::AlphaDriver = 0.0 -ML_BSSN::BetaDriver = 0.0 -ML_BSSN::LapseAdvectionCoeff = 0.0 -ML_BSSN::ShiftAdvectionCoeff = 0.0 - -ML_BSSN::MinimumLapse = 1.0e-8 - -ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas" -ML_BSSN::my_rhs_boundary_condition = "NewRad" - -ML_BSSN::ML_log_confac_bound = "none" -ML_BSSN::ML_metric_bound = "none" -ML_BSSN::ML_Gamma_bound = "none" -ML_BSSN::ML_trace_curv_bound = "none" -ML_BSSN::ML_curv_bound = "none" -ML_BSSN::ML_lapse_bound = "none" -ML_BSSN::ML_dtlapse_bound = "none" -ML_BSSN::ML_shift_bound = "none" -ML_BSSN::ML_dtshift_bound = "none" - - - -# init parameters -InitBase::initial_data_setup_method = "init_some_levels" - -ActiveThorns = "TOVSolver" - -ADMBase::initial_data = "tov" -ADMBase::initial_lapse = "tov" -ADMBase::initial_shift = "tov" -ADMBase::initial_dtlapse = "zero" -ADMBase::initial_dtshift = "zero" - -#TOVSolver::TOV_Rho_Central[0] = 1.28e-3 -TOVSolver::TOV_Mass[0] = 1.4 -TOVSolver::TOV_Gamma[0] = 2.0 -TOVSolver::TOV_K[0] = 100.0 - -# I/O - -IO::out_dir = \$parfile - -IOBasic::outInfo_every = $out_every -IOBasic::outInfo_vars = "HydroBase::rho ADMBase::lapse" - -IOScalar::outScalar_every = $out_every -IOScalar::one_file_per_group = yes -IOScalar::outScalar_vars = " - HydroBase::rho - HydroBase::press - HydroBase::eps - HydroBase::vel - ADMBase::lapse - ADMBase::metric - ADMBase::curv - ADMConstraints::ham - ADMConstraints::momentum - SystemStatistics::process_memory -" - -IOASCII::out1D_every = $out_every -IOASCII::one_file_per_group = yes -IOASCII::output_symmetry_points = no -IOASCII::out3D_ghosts = no -IOASCII::out3D_outer_ghosts = no -IOASCII::out1D_vars = " - HydroBase::rho - HydroBase::press - HydroBase::eps - HydroBase::vel - ADMBase::lapse - ADMBase::metric - ADMBase::curv - ADMConstraints::ham - ADMConstraints::momentum -" -PAR_EOF - - open(OUT, ">tov_${res}.par"); - print OUT "$par"; - close(OUT); -} -

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 87)
by knarf＠cct.lsu.edu 18 May '11

18 May '11

User: knarf Date: 2011/05/18 08:46 AM Added: /examples/tov/ tov_pert.rpar Log: parameter file for perturbed TOV star - note, this so far uses uncommited code within the TOVSolver Directory Changes: Directory: /svn:executable/ =========================== + * File Changes: Directory: /examples/tov/ ========================= File [added]: tov_pert.rpar Delta lines: +242 -0 =================================================================== --- examples/tov/tov_pert.rpar (rev 0) +++ examples/tov/tov_pert.rpar 2011-05-18 13:46:36 UTC (rev 87) @@ -0,0 +1,242 @@ +#!/usr/bin/perl -w +use strict; + +foreach my $res ((8, 4, 2)) { + + my $out_every = 256/$res; + my $par = <<PAR_EOF; +## mclachlan tov_static +ActiveThorns = "Time MoL" +ActiveThorns = "Coordbase CartGrid3d Boundary StaticConformal" +ActiveThorns = "SymBase ADMBase TmunuBase HydroBase InitBase ADMCoupling ADMMacros" +ActiveThorns = "IOUtil" +#ActiveThorns = "Formaline" +ActiveThorns = "SpaceMask CoordGauge Constants LocalReduce aeilocalinterp LoopControl" +ActiveThorns = "Carpet CarpetLib CarpetReduce CarpetRegrid2 CarpetInterp" +ActiveThorns = "CarpetIOASCII CarpetIOScalar CarpetIOHDF5 CarpetIOBasic" +ActiveThorns = "SystemStatistics" + +ActiveThorns = "ADMConstraints" + +# grid parameters +CartGrid3D::type = "coordbase" +CartGrid3D::domain = "full" +CartGrid3D::avoid_origin = "no" +CoordBase::xmin = 0.0 +CoordBase::ymin = 0.0 +CoordBase::zmin = 0.0 +CoordBase::xmax = 240.0 +CoordBase::ymax = 240.0 +CoordBase::zmax = 240.0 +CoordBase::dx = $res +CoordBase::dy = $res +CoordBase::dz = $res + +CoordBase::boundary_size_x_lower = 3 +CoordBase::boundary_size_y_lower = 3 +CoordBase::boundary_size_z_lower = 3 +CoordBase::boundary_size_x_upper = 3 +CoordBase::boundary_size_y_upper = 3 +CoordBase::boundary_size_z_upper = 3 +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 +CoordBase::boundary_shiftout_x_upper = 0 +CoordBase::boundary_shiftout_y_upper = 0 +CoordBase::boundary_shiftout_z_upper = 0 + + + + +ActiveThorns = "ReflectionSymmetry" + +ReflectionSymmetry::reflection_x = "yes" +ReflectionSymmetry::reflection_y = "yes" +ReflectionSymmetry::reflection_z = "yes" +ReflectionSymmetry::avoid_origin_x = "no" +ReflectionSymmetry::avoid_origin_y = "no" +ReflectionSymmetry::avoid_origin_z = "no" + + + +TmunuBase::stress_energy_storage = yes +TmunuBase::stress_energy_at_RHS = yes +TmunuBase::timelevels = 1 +TmunuBase::prolongation_type = none + +HydroBase::timelevels = 3 + +ADMMacros::spatial_order = 4 + +ADMBase::metric_type = "physical" + +ADMConstraints::bound = "static" +ADMConstraints::constraints_timelevels = 3 +ADMConstraints::constraints_persist = yes + +SpaceMask::use_mask = "yes" + +Cactus::terminate = "time" +Cactus::cctk_final_time = 1000 + +Carpet::domain_from_coordbase = "yes" +Carpet::enable_all_storage = no +Carpet::use_buffer_zones = "yes" + +Carpet::poison_new_timelevels = "yes" +Carpet::check_for_poison = "no" +Carpet::poison_value = 113 + +Carpet::init_3_timelevels = no +Carpet::init_fill_timelevels = "yes" + +CarpetLib::poison_new_memory = "yes" +CarpetLib::poison_value = 114 + +# system specific Carpet paramters +Carpet::max_refinement_levels = 10 +driver::ghost_size = 3 +Carpet::prolongation_order_space = 3 +Carpet::prolongation_order_time = 2 + +CarpetRegrid2::regrid_every = 0 +CarpetRegrid2::num_centres = 1 + +CarpetRegrid2::num_levels_1 = 5 +CarpetRegrid2::radius_1[1] =120.0 +CarpetRegrid2::radius_1[2] = 60.0 +CarpetRegrid2::radius_1[3] = 30.0 +CarpetRegrid2::radius_1[4] = 15.0 + + +time::dtfac = 0.25 +MoL::ODE_Method = "rk4" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 + +# check all physical variables for NaNs +ActiveThorns = "NaNChecker" +#NaNChecker::check_every = 1 +#NaNChecker::action_if_found = "just warn" #"terminate", "just warn", "abort" +#NaNChecker::check_vars = "ADMBase::metric ADMBase::lapse ADMBase::shift HydroBase::rho HydroBase::eps HydroBase::press HydroBase::vel" + +## Lapse Condition: \\partial_t alpha = - alpha K +## Shift Condition: \\partial_t beta^i = 0 + +# Hydro paramters +ActiveThorns = "EOS_Omni" +ActiveThorns = "GRHydro" + +HydroBase::evolution_method = "GRHydro" + +GRHydro::riemann_solver = "Marquina" +GRHydro::GRHydro_eos_type = "Polytype" +GRHydro::GRHydro_eos_table = "2D_Polytrope" +GRHydro::recon_method = "ppm" +GRHydro::GRHydro_stencil = 3 +GRHydro::bound = "none" +GRHydro::rho_abs_min = 1.e-10 +#GRHydro::GRHydro = 18 # Tmunu(10), rho,press,eps,w_lorentz,vel, tau +#GRHydro::GRHydro = 10 # gij(6), alpha, beta(3) + + +ActiveThorns = "GenericFD NewRad" +ActiveThorns = "ML_BSSN ML_BSSN_Helper" +ADMBase::evolution_method = "ML_BSSN" +ADMBase::lapse_evolution_method = "ML_BSSN" +ADMBase::shift_evolution_method = "ML_BSSN" +ADMBase::dtlapse_evolution_method= "ML_BSSN" +ADMBase::dtshift_evolution_method= "ML_BSSN" + +TmunuBase::support_old_CalcTmunu_mechanism = "no" + +ML_BSSN::timelevels = 3 + +ML_BSSN::harmonicN = 1 # 1+log +ML_BSSN::harmonicF = 1.0 # 1+log +ML_BSSN::LapseACoeff = 1.0 +ML_BSSN::ShiftBCoeff = 1.0 +ML_BSSN::ShiftGammaCoeff = 0.0 +ML_BSSN::AlphaDriver = 0.0 +ML_BSSN::BetaDriver = 0.0 +ML_BSSN::LapseAdvectionCoeff = 0.0 +ML_BSSN::ShiftAdvectionCoeff = 0.0 + +ML_BSSN::MinimumLapse = 1.0e-8 + +ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas" +ML_BSSN::my_rhs_boundary_condition = "NewRad" + +ML_BSSN::ML_log_confac_bound = "none" +ML_BSSN::ML_metric_bound = "none" +ML_BSSN::ML_Gamma_bound = "none" +ML_BSSN::ML_trace_curv_bound = "none" +ML_BSSN::ML_curv_bound = "none" +ML_BSSN::ML_lapse_bound = "none" +ML_BSSN::ML_dtlapse_bound = "none" +ML_BSSN::ML_shift_bound = "none" +ML_BSSN::ML_dtshift_bound = "none" + + + +# init parameters +InitBase::initial_data_setup_method = "init_some_levels" + +ActiveThorns = "TOVSolver" + +ADMBase::initial_data = "tov" +ADMBase::initial_lapse = "tov" +ADMBase::initial_shift = "tov" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +#TOVSolver::TOV_Rho_Central[0] = 1.28e-3 +TOVSolver::TOV_Mass[0] = 1.4 +TOVSolver::TOV_Gamma[0] = 2.0 +TOVSolver::TOV_K[0] = 100.0 + +# I/O + +IO::out_dir = \$parfile + +IOBasic::outInfo_every = $out_every +IOBasic::outInfo_vars = "HydroBase::rho ADMBase::lapse" + +IOScalar::outScalar_every = $out_every +IOScalar::one_file_per_group = yes +IOScalar::outScalar_vars = " + HydroBase::rho + HydroBase::press + HydroBase::eps + HydroBase::vel + ADMBase::lapse + ADMBase::metric + ADMBase::curv + ADMConstraints::ham + ADMConstraints::momentum + SystemStatistics::process_memory +" + +IOASCII::out1D_every = $out_every +IOASCII::one_file_per_group = yes +IOASCII::output_symmetry_points = no +IOASCII::out3D_ghosts = no +IOASCII::out3D_outer_ghosts = no +IOASCII::out1D_vars = " + HydroBase::rho + HydroBase::press + HydroBase::eps + HydroBase::vel + ADMBase::lapse + ADMBase::metric + ADMBase::curv + ADMConstraints::ham + ADMConstraints::momentum +" +PAR_EOF + + open(OUT, ">tov_${res}.par"); + print OUT "$par"; + close(OUT); +} + Property changes on: examples/tov/tov_pert.rpar ___________________________________________________________________

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/examples/bbh/ (Rev. 86)
by bcmsma＠astro.rit.edu 17 May '11

17 May '11

User: bmundim Date: 2011/05/17 12:27 PM Modified: /examples/bbh/ BBHLowRes.par Log: Parameter file used for this first bbh run. Comments: 1) According to Peter and Ian comments, I will have to change the order in WeylScal4::fd_order and maybe increase the angular resolution in Multipole. 2) This parameter file runs smoothly with the ET curie release version. When the ET development version is used, it develops Nans in the AMR boundary closest to Pi symmetry boundary around 10M of evolution for the second finest grid level. File Changes: Directory: /examples/bbh/ ========================= File [modified]: BBHLowRes.par Delta lines: +74 -69 =================================================================== --- examples/bbh/BBHLowRes.par 2011-05-11 15:34:34 UTC (rev 85) +++ examples/bbh/BBHLowRes.par 2011-05-17 17:27:28 UTC (rev 86) @@ -11,12 +11,12 @@ #Cactus::terminate = "never" #Cactus::terminate = "time" #Cactus::cctk_final_time = 1500.0 -Cactus::terminate = "iteration" +#Cactus::terminate = "iteration" #Cactus::cctk_itlast = 0 -Cactus::cctk_itlast = 256 -#Cactus::terminate="runtime" +#Cactus::cctk_itlast = 256 +Cactus::terminate="runtime" #Cactus::max_runtime=28800 # 46 days -#Cactus::max_runtime=1410 # 23.5 hours +Cactus::max_runtime=1410 # 23.5 hours #Cactus::max_runtime=1380 # 23 hours #Cactus::max_runtime=690 # 11.5 hours #Cactus::max_runtime=450 # 7.5 hours @@ -54,6 +54,8 @@ #------------------------------------------------------------------------------ # Diagnostic parameters: #------------------------------------------------------------------------------ +AHFinderDirect::verbose_level = "physics details" + Carpet::output_timers_every = 1024 Carpet::storage_verbose = "no" #Carpet::timers_verbose = "no" @@ -68,11 +70,16 @@ #CarpetLib::print_timestats_every = 1024 CarpetMask::verbose = "yes" - +CarpetReduce::verbose = "no" CarpetRegrid2::verbose = "yes" +#CarpetRegrid2::veryverbose = "yes" -#NaNChecker::verbose = "all" +CarpetTracker::verbose = "yes" +Dissipation::verbose = "no" + +NaNChecker::verbose = "all" + TimerReport::out_every = 1024 TimerReport::out_filename = "TimerReport" #TimerReport::output_all_timers = "yes" @@ -83,18 +90,16 @@ TwoPunctures::verbose = "yes" -dissipation::verbose = "no" -#CarpetTracker::verbose = "yes" - -AHFinderDirect::verbose_level = "physics details" +PunctureTracker::verbose = "yes" QuasiLocalMeasures::verbose = "yes" +SphericalSurface::verbose = "yes" #------------------------------------------------------------------------------ # Utility parameters: #------------------------------------------------------------------------------ -NaNChecker::check_every = 16 -#NaNChecker::check_every = 1 +#NaNChecker::check_every = 16 +NaNChecker::check_every = 1 NaNChecker::check_vars = " ML_BSSN::ML_log_confac ML_BSSN::ML_metric @@ -215,7 +220,7 @@ CarpetInterp::check_tree_search = "no" CarpetRegrid2::freeze_unaligned_levels = "yes" -CarpetRegrid2::snap_to_coarse = "yes" +#CarpetRegrid2::snap_to_coarse = "yes" CarpetRegrid2::symmetry_rotating180 = "yes" # System specific Carpet parameters: @@ -298,8 +303,8 @@ TwoPunctures::npoints_B = 16 TwoPunctures::npoints_phi = 14 TwoPunctures::Newton_tol = 1.e-10 -#TwoPunctures::grid_setup_method = "evaluation" -TwoPunctures::grid_setup_method = "Taylor expansion" +TwoPunctures::grid_setup_method = "evaluation" +#TwoPunctures::grid_setup_method = "Taylor expansion" TwoPunctures::par_b = 3.0 TwoPunctures::par_m_plus = 0.47656 @@ -316,6 +321,7 @@ InitBase::initial_data_setup_method = "init_all_levels" ADMMacros::spatial_order = 4 +SummationByParts::order = 4 #SpaceMask::use_mask = "yes" @@ -374,7 +380,7 @@ ML_BSSN::harmonicN = 1 # 1+log ML_BSSN::harmonicF = 2.0 # 1+log ML_BSSN::ShiftGammaCoeff = 0.75 -ML_BSSN::AlphaDriver = 1.0 +#ML_BSSN::AlphaDriver = 1.0 ML_BSSN::BetaDriver = 1.0 ML_BSSN::LapseAdvectionCoeff = 1.0 ML_BSSN::ShiftAdvectionCoeff = 1.0 @@ -432,8 +438,8 @@ IOASCII::one_file_per_group = "yes" IOASCII::output_symmetry_points = "no" -IOASCII::output_ghost_points = "no" -#IOASCII::out3D_ghosts = "no" +#IOASCII::output_ghost_points = "no" #hg version only! +IOASCII::out3D_ghosts = "no" #IOASCII::out3D_outer_ghosts = "no" @@ -451,32 +457,34 @@ " #Set these IOASCII options for initial data only: -#IOASCII::out1D_every = 512000000 +IOASCII::out1D_every = 512000000 +IOASCII::out1D_d = "no" +IOASCII::out1D_vars = " + CarpetReduce::weight + ADMBase::lapse + ADMBase::metric + ADMBase::curv + ADMBase::shift + ML_BSSN::ML_Ham + ML_BSSN::ML_mom + ML_ADMConstraints::ML_Ham + ML_ADMConstraints::ML_mom +" + +# To draw the AH circunference in the z=0 plane +# set these IOASCII options (for evolution only): +# NOTE: it is not possible yet to output array lines +# as set below. +#IOASCII::out1D_every = 32 +#IOASCII::out1D_x = "no" +#IOASCII::out1D_z = "no" #IOASCII::out1D_d = "no" +#IOASCII::out1D_yline_xi = 19 +#IOASCII::out1D_yline_zi = 0 #IOASCII::out1D_vars = " -# CarpetReduce::weight -# ADMBase::lapse -# ADMBase::metric -# ADMBase::curv -# ADMBase::shift -# ML_BSSN::ML_Ham -# ML_BSSN::ML_mom -# ML_ADMConstraints::ML_Ham -# ML_ADMConstraints::ML_mom +# SphericalSurface::sf_radius #" -# To draw the AH circunference in the z=0 plane -# set these IOASCII options (for evolution only): -IOASCII::out1D_every = 32 -IOASCII::out1D_x = "no" -IOASCII::out1D_z = "no" -IOASCII::out1D_d = "no" -IOASCII::out1D_yline_xi = 19 -IOASCII::out1D_yline_zi = 0 -IOASCII::out1D_vars = " - SphericalSurface::sf_radius -" - IOASCII::out2D_every = 512000000 IOASCII::out2D_vars = " ML_ADMConstraints::ML_Ham @@ -500,6 +508,7 @@ Grid::coordinates ADMBase::lapse ADMBase::metric + ADMBase::curv ML_ADMConstraints::ML_Ham ML_ADMConstraints::ML_Mom WeylScal4::Psi4r @@ -512,13 +521,14 @@ #IOHDF5::out3D_outer_ghosts = "no" IOHDF5::out_every = 512000000 +#IOHDF5::out_every = 256 IOHDF5::out_vars = " + CarpetReduce::weight ADMBase::lapse ML_ADMConstraints::ML_Ham ML_ADMConstraints::ML_Mom " # Grid::coordinates -# CarpetReduce::weight # ADMBase::metric # ML_BSSN::ML_Ham # ML_BSSN::ML_mom @@ -528,22 +538,25 @@ #------------------------------------------------------------------------------ AHFinderDirect::find_every = 32 -AHFinderDirect::run_at_CCTK_ANALYSIS = "yes" -AHFinderDirect::run_at_CCTK_POSTSTEP = "no" +#AHFinderDirect::run_at_CCTK_ANALYSIS = "yes" +#AHFinderDirect::run_at_CCTK_POSTSTEP = "no" +AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = "no" -AHFinderDirect::move_origins = yes +AHFinderDirect::move_origins = "yes" +AHFinderDirect::reshape_while_moving = "yes" +AHFinderDirect::predict_origin_movement = "yes" # Hermite to order 3 to avoid discontinuities in the metric spatial derivatives: -AHFinderDirect::geometry_interpolator_name = "Hermite polynomial interpolation" -AHFinderDirect::geometry_interpolator_pars = "order=3" -AHFinderDirect::surface_interpolator_name = "Hermite polynomial interpolation" -AHFinderDirect::surface_interpolator_pars = "order=3" -#AHFinderDirect::geometry_interpolator_name = "Lagrange polynomial interpolation" -#AHFinderDirect::geometry_interpolator_pars = "order=4" -#AHFinderDirect::surface_interpolator_name = "Lagrange polynomial interpolation" -#AHFinderDirect::surface_interpolator_pars = "order=4" +#AHFinderDirect::geometry_interpolator_name = "Hermite polynomial interpolation" +#AHFinderDirect::geometry_interpolator_pars = "order=3" +#AHFinderDirect::surface_interpolator_name = "Hermite polynomial interpolation" +#AHFinderDirect::surface_interpolator_pars = "order=3" +AHFinderDirect::geometry_interpolator_name = "Lagrange polynomial interpolation" +AHFinderDirect::geometry_interpolator_pars = "order=4" +AHFinderDirect::surface_interpolator_name = "Lagrange polynomial interpolation" +AHFinderDirect::surface_interpolator_pars = "order=4" -AHFinderDirect::output_h_every = 512 +AHFinderDirect::output_h_every = 0 AHFinderDirect::N_horizons = 3 @@ -610,23 +623,15 @@ #------------------------------------------------------------------------------ # Checkpoint/Recovery: #------------------------------------------------------------------------------ -CarpetIOHDF5::checkpoint = "yes" -#CarpetIOHDF5::checkpoint = "no" - -IO::checkpoint_ID = "yes" #may save you from calculating the initial data again -IO::checkpoint_every = 6144 -#IO::checkpoint_every = 6912 -##IO::checkpoint_every = 9750 +IOHDF5::checkpoint = "yes" +IO::checkpoint_dir = $parfile +IO::checkpoint_ID = "yes" +IO::checkpoint_every = 6144 +#IO::checkpoint_every_walltime_hours = 6.0 IO::checkpoint_keep=3 -IO::checkpoint_on_terminate = "yes" -#IO::checkpoint_dir = $parfile -#IO::checkpoint_file = "checkpoint" -#IO::verbose = "full" +IO::checkpoint_on_terminate = "yes" -#IO::recover_dir = $parfile -#IO::recover_file = "checkpoint" -IO::recover = "autoprobe" -#IO::recover = "auto" +IO::recover = "autoprobe" +IO::recover_dir = $parfile -

3 2

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits May 2011