Commits May 2011

commits@lists.einsteintoolkit.org

6 participants
21 discussions

[svn:einsteintoolkit] TOVSolver/trunk/src/ (Rev. 128)
by cott＠tapir.caltech.edu 17 May '11

17 May '11

User: cott Date: 2011/05/17 09:53 PM Modified: /trunk/src/ tov.c Log: * set up consistent values of press and eps in the atmosphere. File Changes: Directory: /trunk/src/ ====================== File [modified]: tov.c Delta lines: +10 -4 =================================================================== --- trunk/src/tov.c 2011-04-21 04:47:01 UTC (rev 127) +++ trunk/src/tov.c 2011-05-18 02:53:28 UTC (rev 128) @@ -773,11 +773,14 @@ my_psi4); rho[i3D] = max_rho; + eps[i3D] = eps_point[star]; + press[i3D] = press_point[star]; /* Set atmosphere according to chosen star */ - if(rho[i3D] <= TOV_Atmosphere[star]) + if(rho[i3D] <= TOV_Atmosphere[star]) { rho[i3D] = TOV_Atmosphere[star]; - eps[i3D] = eps_point[star]; - press[i3D] = press_point[star]; + press[i3D] = TOV_K[star] * pow(rho[i3D],TOV_Gamma[star]); + eps[i3D] = press[i3D]/(TOV_Gamma[star]-1.0)/rho[i3D]; + } } else if (CCTK_EQUALS(TOV_Combine_Method, "average")) @@ -823,8 +826,11 @@ } /* Reset atmosphere according to chosen star */ - if(rho[i3D] <= TOV_Atmosphere[star_i]) + if(rho[i3D] <= TOV_Atmosphere[star_i]) { rho[i3D] = TOV_Atmosphere[star_i]; + press[i3D] = TOV_K[star_i] * pow(rho[i3D],TOV_Gamma[star_i]); + eps[i3D] = press[i3D]/(TOV_Gamma[star_i]-1.0)/rho[i3D]; + } } if (TOV_Conformal_Flat_Three_Metric)

1 0

[svn:einsteintoolkit] www/global/ (Rev. 588)
by knarf＠cct.lsu.edu 12 May '11

12 May '11

User: knarf Date: 2011/05/12 02:06 PM Modified: /global/ navigation.php Log: correct typo File Changes: Directory: /global/ =================== File [modified]: navigation.php Delta lines: +1 -1 =================================================================== --- global/navigation.php 2011-05-07 04:38:08 UTC (rev 587) +++ global/navigation.php 2011-05-12 19:06:04 UTC (rev 588) @@ -23,7 +23,7 @@ <li><a href="/documentation/">Documentation</a> <ul class="subnav"> - <li><a href="/documentation/licenses">Citating</a> + <li><a href="/documentation/licenses">Citing</a> <li><a href="/documentation/UsersGuide/UsersGuide.html">User's Guide</a></li> <li><a href="/documentation/ThornGuide/ThornGuide.html">Thorn Guide</a></li> <li><a href="/documentation/ReferenceManual/ReferenceManual.html">Reference Manual</a></li>

1 0

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

11 May '11

User: bmundim Date: 2011/05/11 10:34 AM Added: /examples/bbh/ BBHLowRes.par Log: Updated BBH parameter file. Main differences, besides reorganization, consist in: 1) reducing the total number of spherical surfaces from 5 to 3. 2) Change the AHFinderDirect interpolator to Hermite polynomial. 3) Add a few parameters to CarpetLib, CarpetInterp and CarpetRegrid2. Issues this parameter file is facing so far: a) Its evolution is stopped with Nans around 8.6M. Nans appear to come from the BSSN fields, but still needs further investigation. b) Either CarpetMask or CarpetReduce don't seem to do their job. CarpetMask uses the AH surfaces to mask their interior. At least it says it does so. However the reduction operations still don't mask out the AH interior, at least when performed by CarpetIOScalar or CarpetIOBasic. c) When CarpetReduce::weight is output, its values seem strange to me. Most of the domain has 1/2 set. A few points close to the boundary 1/4 and the boundaries 0 (physical and ghost). No mask out values inside AH is shown, ie they still are set to 1/2. These values seem to make more sense to AMR restriction operations. d) CarpetIOASCII 1D output is not working for grid arrays as expected for grid functions. The line output for SphericalSurface::sf_radius shows that. File Changes: Directory: /examples/bbh/ ========================= File [added]: BBHLowRes.par Delta lines: +632 -0 =================================================================== --- examples/bbh/BBHLowRes.par (rev 0) +++ examples/bbh/BBHLowRes.par 2011-05-11 15:34:34 UTC (rev 85) @@ -0,0 +1,632 @@ +#============================================================================== +# BBH: TwoPunctures - Mclachlan +#============================================================================== + +#------------------------------------------------------------------------------ +# Cactus parameters: +#------------------------------------------------------------------------------ +Cactus::cctk_run_title = "BBH" +Cactus::cctk_full_warnings = "yes" + +#Cactus::terminate = "never" +#Cactus::terminate = "time" +#Cactus::cctk_final_time = 1500.0 +Cactus::terminate = "iteration" +#Cactus::cctk_itlast = 0 +Cactus::cctk_itlast = 256 +#Cactus::terminate="runtime" +#Cactus::max_runtime=28800 # 46 days +#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 +#Cactus::max_runtime=210 # 3.5 hours +#Cactus::max_runtime=110 # 1.83 hours +#Cactus::max_runtime=55 # 55min + +#------------------------------------------------------------------------------ +# Activate all necessary thorns: +#------------------------------------------------------------------------------ + +ActiveThorns = "Boundary CartGrid3D CoordBase Fortran InitBase IOUtil LocalReduce SymBase Time" +ActiveThorns = "AEILocalInterp LocalInterp" +#ActiveThorns = "MoL ReflectionSymmetry RotatingSymmetry90 RotatingSymmetry180 Slab SpaceMask SphericalSurface" +ActiveThorns = "MoL ReflectionSymmetry RotatingSymmetry180 Slab SpaceMask SphericalSurface" +#ActiveThorns = "MoL ReflectionSymmetry Slab SpaceMask SphericalSurface" +ActiveThorns = "Carpet CarpetInterp CarpetIOASCII CarpetIOHDF5 CarpetIOScalar CarpetLib CarpetIOBasic CarpetReduce CarpetRegrid2 CarpetSlab CarpetTracker CarpetMask LoopControl" +#ActiveThorns = "Formaline NaNChecker TerminationTrigger TimerReport" +ActiveThorns = "NaNChecker TerminationTrigger TimerReport" +ActiveThorns = "ADMbase ADMcoupling ADMmacros CoordGauge StaticConformal" +ActiveThorns = "PunctureTracker" +ActiveThorns = "TmunuBase" +ActiveThorns = "QuasiLocalMeasures" +#ActiveThorns = "ADMConstraints" +ActiveThorns = "BLAS LAPACK GSL HDF5" +ActiveThorns = "TwoPunctures" +ActiveThorns = "SummationByParts" +ActiveThorns = "GenericFD NewRad" +ActiveThorns = "ML_BSSN ML_BSSN_Helper ML_ADMConstraints" +ActiveThorns = "Dissipation" +ActiveThorns = "AHFinderDirect" +ActiveThorns = "WeylScal4 Multipole" +ActiveThorns = "SystemStatistics" + +#------------------------------------------------------------------------------ +# Diagnostic parameters: +#------------------------------------------------------------------------------ +Carpet::output_timers_every = 1024 +Carpet::storage_verbose = "no" +#Carpet::timers_verbose = "no" +Carpet::verbose = "no" +#Carpet::verbose = "yes" # only for Initial Data +Carpet::veryverbose = "no" +Carpet::grid_structure_filename = "carpet-grid-structure" +Carpet::grid_coordinates_filename = "carpet-grid-coordinates" + +CarpetLib::output_bboxes = "no" +CarpetLib::print_memstats_every = 1024 +#CarpetLib::print_timestats_every = 1024 + +CarpetMask::verbose = "yes" + +CarpetRegrid2::verbose = "yes" + +#NaNChecker::verbose = "all" + +TimerReport::out_every = 1024 +TimerReport::out_filename = "TimerReport" +#TimerReport::output_all_timers = "yes" +#TimerReport::output_all_timers_together = "yes" +#TimerReport::output_all_timers_readable = "yes" +#TimerReport::before_checkpoint = "yes" +TimerReport::n_top_timers = 40 + +TwoPunctures::verbose = "yes" + +dissipation::verbose = "no" +#CarpetTracker::verbose = "yes" + +AHFinderDirect::verbose_level = "physics details" +QuasiLocalMeasures::verbose = "yes" + +#------------------------------------------------------------------------------ +# Utility parameters: +#------------------------------------------------------------------------------ + +NaNChecker::check_every = 16 +#NaNChecker::check_every = 1 +NaNChecker::check_vars = " + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_trace_curv + ML_BSSN::ML_curv + ML_BSSN::ML_Gamma + ML_BSSN::ML_lapse + ML_BSSN::ML_shift + ML_BSSN::ML_dtlapse + ML_BSSN::ML_dtshift + ADMBase::curv + ADMBase::metric + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift +" +NaNChecker::action_if_found = "terminate" +#NaNChecker::action_if_found = "abort" +#NaNChecker::action_if_found = "just warn" #"terminate", "just warn", "abort" +#NaNChecker::check_after=0 + +#TerminationTrigger::max_walltime = 0.0 # hours +TerminationTrigger::max_walltime = 24.0 # hours +TerminationTrigger::on_remaining_walltime = 10.0 # minutes +TerminationTrigger::termination_from_file = "yes" +TerminationTrigger::create_termination_file = "yes" +TerminationTrigger::termination_file = "../TERMINATE" + + +#------------------------------------------------------------------------------ +# Run parameters: +#------------------------------------------------------------------------------ + +#------ +# Grid: +#------ + +Time::dtfac = 0.25 +MoL::ODE_Method = "rk4" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 + +CartGrid3D::type = "coordbase" +CartGrid3D::domain = "full" +CartGrid3D::avoid_origin = "no" + +CoordBase::domainsize = "minmax" +CoordBase::spacing = "gridspacing" # "gridspacing" or "numcells" + +# On a workstation (Total required memory: ~1.9GB): +#CoordBase::xmin = 0.00 +#CoordBase::ymin =-12.00 +#CoordBase::zmin = 0.00 +#CoordBase::xmax = 12.00 +#CoordBase::ymax = 12.00 +#CoordBase::zmax = 12.00 +#CoordBase::dx = 0.40 # dx or ncells_x +#CoordBase::dy = 0.40 # dy or ncells_y +#CoordBase::dz = 0.40 # dz or ncells_z + +# On a HPC (Total required memory: ~6.5GB): +CoordBase::xmin = 0.00 +CoordBase::ymin =-120.00 +CoordBase::zmin = 0.00 +CoordBase::xmax = 120.00 +CoordBase::ymax = 120.00 +CoordBase::zmax = 120.00 +CoordBase::dx = 2.00 # dx or ncells_x +CoordBase::dy = 2.00 # dy or ncells_y +CoordBase::dz = 2.00 # dz or ncells_z + +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 = 0 +CoordBase::boundary_shiftout_z_lower = 1 +CoordBase::boundary_shiftout_x_upper = 0 +CoordBase::boundary_shiftout_y_upper = 0 +CoordBase::boundary_shiftout_z_upper = 0 + +Driver::ghost_size = 3 + +ReflectionSymmetry::reflection_x = "no" +ReflectionSymmetry::reflection_y = "no" +ReflectionSymmetry::reflection_z = "yes" +ReflectionSymmetry::avoid_origin_x = "no" +ReflectionSymmetry::avoid_origin_y = "no" +ReflectionSymmetry::avoid_origin_z = "no" + + +# General Carpet parameters: +Carpet::domain_from_coordbase = "yes" +Carpet::enable_all_storage = "no" +Carpet::use_buffer_zones = "yes" +Carpet::schedule_barriers = "no" + +Carpet::poison_new_timelevels = "yes" +Carpet::check_for_poison = "no" +#Carpet::poison_value = 113 # git version only! + +Carpet::init_3_timelevels = "no" +Carpet::init_fill_timelevels = "yes" + +CarpetLib::poison_new_memory = "yes" +CarpetLib::poison_value = 114 +CarpetLib::check_bboxes = "no" +CarpetLib::interleave_communications = "yes" +#CarpetLib::vary_tags = "yes" # git version only! +CarpetLib::combine_sends = "yes" +#CarpetLib::reduce_mpi_waitall = "yes" # git version only! + +CarpetInterp::tree_search = "yes" +CarpetInterp::check_tree_search = "no" + +CarpetRegrid2::freeze_unaligned_levels = "yes" +CarpetRegrid2::snap_to_coarse = "yes" +CarpetRegrid2::symmetry_rotating180 = "yes" + +# System specific Carpet parameters: +Carpet::time_refinement_factors = "[1, 1, 2, 4, 8, 16, 32, 64, 128, 256]" +Carpet::max_refinement_levels = 9 +Carpet::prolongation_order_space = 5 +Carpet::prolongation_order_time = 2 + +CarpetRegrid2::regrid_every = 32 +CarpetRegrid2::num_centres = 2 + +# On a workstation (Total required memory: ~1.9GB): +#CarpetRegrid2::num_levels_1 = 6 +#CarpetRegrid2::position_x_1 = +3.0 +#CarpetRegrid2::radius_1[ 1] = 2.8 +#CarpetRegrid2::radius_1[ 2] = 2.0 +#CarpetRegrid2::radius_1[ 3] = 1.2 +#CarpetRegrid2::radius_1[ 4] = 0.6 +#CarpetRegrid2::radius_1[ 5] = 0.3 +#CarpetRegrid2::movement_threshold_1 = 0.16 + +#CarpetRegrid2::num_levels_2 = 6 +#CarpetRegrid2::position_x_2 = -3.0 +#CarpetRegrid2::radius_2[ 1] = 2.8 +#CarpetRegrid2::radius_2[ 2] = 2.0 +#CarpetRegrid2::radius_2[ 3] = 1.2 +#CarpetRegrid2::radius_2[ 4] = 0.6 +#CarpetRegrid2::radius_2[ 5] = 0.3 +#CarpetRegrid2::movement_threshold_2 = 0.16 + +# On a HPC (Total required memory: ~6.5GB): +CarpetRegrid2::num_levels_1 = 7 +CarpetRegrid2::position_x_1 = +3.0 +CarpetRegrid2::radius_1[ 1] = 64.0 +CarpetRegrid2::radius_1[ 2] = 16.0 +CarpetRegrid2::radius_1[ 3] = 8.0 +CarpetRegrid2::radius_1[ 4] = 4.0 +CarpetRegrid2::radius_1[ 5] = 2.0 +CarpetRegrid2::radius_1[ 6] = 1.0 +CarpetRegrid2::movement_threshold_1 = 0.16 + +CarpetRegrid2::num_levels_2 = 7 +CarpetRegrid2::position_x_2 = -3.0 +CarpetRegrid2::radius_2[ 1] = 64.0 +CarpetRegrid2::radius_2[ 2] = 16.0 +CarpetRegrid2::radius_2[ 3] = 8.0 +CarpetRegrid2::radius_2[ 4] = 4.0 +CarpetRegrid2::radius_2[ 5] = 2.0 +CarpetRegrid2::radius_2[ 6] = 1.0 +CarpetRegrid2::movement_threshold_2 = 0.16 + +CarpetMask::excluded_surface [0] = 0 +CarpetMask::excluded_surface_factor[0] = 1.0 +CarpetMask::excluded_surface [1] = 1 +CarpetMask::excluded_surface_factor[1] = 1.0 +CarpetMask::excluded_surface [2] = 2 +CarpetMask::excluded_surface_factor[2] = 1.0 + +CarpetTracker::surface[0] = 0 +CarpetTracker::surface[1] = 1 + +#------ +# MODEL: +#------ + +ADMBase::metric_type = "physical" +ADMBase::initial_data = "twopunctures" +ADMBase::initial_lapse = "twopunctures-averaged" +ADMBase::initial_shift = "zero" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" +# needed for AHFinderDirect: +ADMBase::metric_timelevels = 3 + + + +#TwoPunctures::keep_u_around = "yes" + +TwoPunctures::npoints_A = 16 +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::par_b = 3.0 +TwoPunctures::par_m_plus = 0.47656 +TwoPunctures::par_m_minus = 0.47656 +TwoPunctures::par_P_plus [1] = +0.13808 +TwoPunctures::par_P_minus[1] = -0.13808 + +TwoPunctures::TP_Tiny = 1.e-10 + +#---------- +# Numerics: +#---------- + +InitBase::initial_data_setup_method = "init_all_levels" + +ADMMacros::spatial_order = 4 + +#SpaceMask::use_mask = "yes" + +#Dissipation::epsdis = 0.01 +#Dissipation::extra_dissipation_in_horizons = "yes" +#Dissipation::update_ah_every = 16 +#Dissipation::ah_slope = 0.2 +#Dissipation::ah_max_epsdis=0.1 +#Dissipation::surface_number[0] = 1 +#Dissipation::horizon_number[0] = 1 +Dissipation::order = 5 +Dissipation::vars = " + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_curv + ML_BSSN::ML_trace_curv + ML_BSSN::ML_Gamma + ML_BSSN::ML_lapse + ML_BSSN::ML_shift + ML_BSSN::ML_dtlapse + ML_BSSN::ML_dtshift +" + +SphericalSurface::nsurfaces = 3 +SphericalSurface::maxntheta = 39 +SphericalSurface::maxnphi = 76 + +SphericalSurface::ntheta [0] = 39 +SphericalSurface::nphi [0] = 76 +SphericalSurface::nghoststheta[0] = 2 +SphericalSurface::nghostsphi [0] = 2 + +SphericalSurface::ntheta [1] = 39 +SphericalSurface::nphi [1] = 76 +SphericalSurface::nghoststheta[1] = 2 +SphericalSurface::nghostsphi [1] = 2 + +SphericalSurface::ntheta [2] = 39 +SphericalSurface::nphi [2] = 76 +SphericalSurface::nghoststheta[2] = 2 +SphericalSurface::nghostsphi [2] = 2 + + +#----------- +# Evolution: +#----------- + +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" + +ML_BSSN::timelevels = 3 + +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::BetaDriver = 1.0 +ML_BSSN::LapseAdvectionCoeff = 1.0 +ML_BSSN::ShiftAdvectionCoeff = 1.0 + +ML_BSSN::MinimumLapse = 1.0e-8 + +ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas" +ML_BSSN::my_rhs_boundary_condition = "NewRad" +Boundary::radpower = 2 + +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" + + +#------------------------------------------------------------------------------ +# Output: +#------------------------------------------------------------------------------ + +IO::out_dir = $parfile +IO::out_fileinfo = "none" + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_reductions = "maximum" +IOBasic::outInfo_vars = " + Carpet::physical_time_per_hour + SystemStatistics::maxrss_mb + ML_ADMConstraints::ML_Ham + ML_BSSN::ML_trace_curv +" + +IOScalar::outScalar_every = 32 +IOScalar::one_file_per_group = "yes" +IOScalar::outScalar_reductions = "minimum maximum average norm1 norm2" +IOScalar::outScalar_vars = " + ADMBase::lapse + ADMBase::shift + ADMBase::metric + ADMBase::curv + ADMBase::dtlapse + ADMBase::dtshift + ML_ADMConstraints::ML_Ham + ML_ADMConstraints::ML_mom + ML_BSSN::ML_Ham + ML_BSSN::ML_mom + SystemStatistics::process_memory_mb{out_every=4} + SphericalSurface::sf_radius +" + +IOASCII::one_file_per_group = "yes" +IOASCII::output_symmetry_points = "no" +IOASCII::output_ghost_points = "no" +#IOASCII::out3D_ghosts = "no" +#IOASCII::out3D_outer_ghosts = "no" + + +IOASCII::out0D_every = 32 +IOASCII::out0D_vars = " + Carpet::timing + PunctureTracker::pt_loc + QuasiLocalMeasures::qlm_scalars + SphericalSurface::sf_active + SphericalSurface::sf_valid + SphericalSurface::sf_info + SphericalSurface::sf_radius + SphericalSurface::sf_origin + SphericalSurface::sf_coordinate_descriptors +" + +#Set these IOASCII options for initial data only: +#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): +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 + ML_ADMConstraints::ML_mom + SphericalSurface::sf_radius{out_every=32} +" +# ADMBase::lapse +# ADMBase::metric +# ADMBase::curv +# ML_BSSN::ML_Ham +# ML_BSSN::ML_mom + +IOHDF5::one_file_per_group = "yes" +IOHDF5::open_one_input_file_at_a_time = "yes" +IOHDF5::out2D_every = 512 +IOHDF5::out2D_xy = "yes" +IOHDF5::out2D_xz = "no" +IOHDF5::out2D_yz = "no" +IOHDF5::out2D_xyplane_z = 0.0 +IOHDF5::out2D_vars = " + Grid::coordinates + ADMBase::lapse + ADMBase::metric + ML_ADMConstraints::ML_Ham + ML_ADMConstraints::ML_Mom + WeylScal4::Psi4r + WeylScal4::Psi4i + " +# CarpetReduce::weight + + +#IOHDF5::out3D_ghosts = "no" +#IOHDF5::out3D_outer_ghosts = "no" + +IOHDF5::out_every = 512000000 +IOHDF5::out_vars = " + ADMBase::lapse + ML_ADMConstraints::ML_Ham + ML_ADMConstraints::ML_Mom +" +# Grid::coordinates +# CarpetReduce::weight +# ADMBase::metric +# ML_BSSN::ML_Ham +# ML_BSSN::ML_mom + +#------------------------------------------------------------------------------ +# Analysis: +#------------------------------------------------------------------------------ +AHFinderDirect::find_every = 32 + +AHFinderDirect::run_at_CCTK_ANALYSIS = "yes" +AHFinderDirect::run_at_CCTK_POSTSTEP = "no" + +AHFinderDirect::move_origins = 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::output_h_every = 512 + +AHFinderDirect::N_horizons = 3 + +AHFinderDirect::origin_x [1] = +3.0 +AHFinderDirect::initial_guess__coord_sphere__x_center[1] = +3.0 +AHFinderDirect::initial_guess__coord_sphere__radius [1] = 0.25 +AHFinderDirect::which_surface_to_store_info [1] = 0 +AHFinderDirect::set_mask_for_individual_horizon [1] = no +AHFinderDirect::reset_horizon_after_not_finding [1] = no +AHFinderDirect::track_origin_from_grid_scalar [1] = yes +AHFinderDirect::track_origin_source_x [1] = "PunctureTracker::pt_loc_x[0]" +AHFinderDirect::track_origin_source_y [1] = "PunctureTracker::pt_loc_y[0]" +AHFinderDirect::track_origin_source_z [1] = "PunctureTracker::pt_loc_z[0]" +AHFinderDirect::max_allowable_horizon_radius [1] = 3 + +AHFinderDirect::origin_x [2] = -3.0 +AHFinderDirect::initial_guess__coord_sphere__x_center[2] = -3.0 +AHFinderDirect::initial_guess__coord_sphere__radius [2] = 0.25 +AHFinderDirect::which_surface_to_store_info [2] = 1 +AHFinderDirect::set_mask_for_individual_horizon [2] = no +AHFinderDirect::reset_horizon_after_not_finding [2] = no +AHFinderDirect::track_origin_from_grid_scalar [2] = yes +AHFinderDirect::track_origin_source_x [2] = "PunctureTracker::pt_loc_x[1]" +AHFinderDirect::track_origin_source_y [2] = "PunctureTracker::pt_loc_y[1]" +AHFinderDirect::track_origin_source_z [2] = "PunctureTracker::pt_loc_z[1]" +AHFinderDirect::max_allowable_horizon_radius [2] = 3 + +AHFinderDirect::origin_x [3] = 0 +AHFinderDirect::find_after_individual_time [3] = 100.0 +AHFinderDirect::initial_guess__coord_sphere__x_center[3] = 0 +AHFinderDirect::initial_guess__coord_sphere__radius [3] = 1.0 +AHFinderDirect::which_surface_to_store_info [3] = 2 +AHFinderDirect::reset_horizon_after_not_finding [3] = no +AHFinderDirect::max_allowable_horizon_radius [3] = 6 + + +PunctureTracker::track [0] = "yes" +PunctureTracker::initial_x [0] = 3.0 +PunctureTracker::which_surface_to_store_info[0] = 0 +PunctureTracker::track [1] = "yes" +PunctureTracker::initial_x [1] = -3.0 +PunctureTracker::which_surface_to_store_info[1] = 1 + +QuasiLocalMeasures::num_surfaces = 3 +QuasiLocalMeasures::spatial_order = 4 +QuasiLocalMeasures::interpolator = "Lagrange polynomial interpolation" +QuasiLocalMeasures::interpolator_options = "order=4" +QuasiLocalMeasures::surface_index [0] = 0 +QuasiLocalMeasures::surface_index [1] = 1 +QuasiLocalMeasures::surface_index [2] = 2 + +Multipole::nradii = 4 +Multipole::radius[0] = 30 +Multipole::radius[1] = 40 +Multipole::radius[2] = 50 +Multipole::radius[3] = 60 +Multipole::ntheta = 60 +Multipole::nphi = 120 +Multipole::variables = "WeylScal4::Psi4r{sw=-2 cmplx='WeylScal4::Psi4i' name='psi4'}" +Multipole::out_every = 4 +Multipole::l_max = 4 + + +#------------------------------------------------------------------------------ +# 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 +IO::checkpoint_keep=3 +IO::checkpoint_on_terminate = "yes" +#IO::checkpoint_dir = $parfile +#IO::checkpoint_file = "checkpoint" +#IO::verbose = "full" + +#IO::recover_dir = $parfile +#IO::recover_file = "checkpoint" +IO::recover = "autoprobe" +#IO::recover = "auto" + + +

1 0

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

09 May '11

User: knarf Date: 2011/05/09 10:48 AM Modified: / ET.tex Log: enable kasner example pics in tex File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +11 -11 =================================================================== --- ET.tex 2011-05-09 15:11:56 UTC (rev 83) +++ ET.tex 2011-05-09 15:48:54 UTC (rev 84) @@ -1816,8 +1816,8 @@ \label{sec:bbh-example} \subsection{TOV\pages{2 Frank}} -Show stable TOV evolution + \subsection{Collapse\pages{2 Christian}} Show TOV collapse and BH formation @@ -1858,16 +1858,16 @@ of spacetime~\ref{eq:kasner}, along with the associated error for a sequence of time resolutions. -%\begin{figure} -%\includegraphics[width=0.45\textwidth]{kasner.png} -%\includegraphics[width=0.45\textwidth]{err.png} -%\caption{Left: the evolution of a vacuum spacetime of the type~\ref{eq:gowdyT3}, -%with $P=Q=\lambda=0$; the initial data are chosen as -%$\gamma_{ij}=\delta_{ij}$ and $K_{ij}={\rm diag}(-2/3,-2/3,1/3)$. -%Right: the numerical error for a sequence of four time resolutions $dt=[0.0125,0.025,0.05,0.1]$; -%the errors are scaled according to the expectation for fourth-order convergence. -%\label{fig:kasner}} -%\end{figure} +\begin{figure} + \includegraphics[width=0.45\textwidth]{kasner.png} + \includegraphics[width=0.45\textwidth]{err.png} + \caption{Left: the evolution of a vacuum spacetime of the type~\ref{eq:gowdyT3}, + with $P=Q=\lambda=0$; the initial data are chosen as + $\gamma_{ij}=\delta_{ij}$ and $K_{ij}={\rm diag}(-2/3,-2/3,1/3)$. + Right: the numerical error for a sequence of four time resolutions $dt=[0.0125,0.025,0.05,0.1]$; + the errors are scaled according to the expectation for fourth-order convergence. + \label{fig:kasner}} +\end{figure} \section{Future Work\pages{1 Frank}}

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 83)
by eloisa.bentivegna＠aei.mpg.de 09 May '11

09 May '11

User: bentivegna Date: 2011/05/09 10:11 AM Added: /figures/ err.png, kasner.png Log: Figures for the cosmology example section. Directory Changes: Directory: /svn:mime-type/ ========================== + application/octet-stream File Changes: Directory: /figures/ ==================== File [added]: err.png Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: figures/err.png ___________________________________________________________________ File [added]: kasner.png Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: figures/kasner.png ___________________________________________________________________

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 247)
by schnetter＠cct.lsu.edu 08 May '11

08 May '11

User: eschnett Date: 2011/05/08 08:03 PM Modified: /trunk/ interface.ccl Log: Make Y_e fluxes publicly visible File Changes: Directory: /trunk/ ================== File [modified]: interface.ccl Delta lines: +2 -0 =================================================================== --- trunk/interface.ccl 2011-04-28 21:33:20 UTC (rev 246) +++ trunk/interface.ccl 2011-05-09 01:03:12 UTC (rev 247) @@ -529,7 +529,9 @@ } "Maximum characteristic speeds for the flux splitting" real Y_e_con_rhs TYPE=GF tags='Prolongation="None" checkpoint="no"' "RHS for the electron fraction" +public: real Y_e_con_flux TYPE=GF tags='Prolongation="None" checkpoint="no"' "Flux for the electron fraction" +private: real Y_e_plus TYPE=GF tags='Prolongation="None" checkpoint="no"' "Plus state for the electron fraction" real Y_e_minus TYPE=GF tags='Prolongation="None" checkpoint="no"' "Minus state for the electron fraction"

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 82)
by roland.haas＠physics.gatech.edu 07 May '11

07 May '11

User: rhaas Date: 2011/05/08 12:23 AM Modified: / ET.tex Log: work a bit on language in simulation domain section File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +57 -36 =================================================================== --- ET.tex 2011-04-28 19:13:37 UTC (rev 81) +++ ET.tex 2011-05-08 05:23:51 UTC (rev 82) @@ -1597,11 +1597,12 @@ \subsection{Simulation Domain, Symmetries, Boundaries\pages{3 Roland}} \subsubsection{Domains and Coordinates.} -The Einstein Toolkit provides a module \codename{CoordBase} to facilate +The Einstein Toolkit provided thorn \codename{CoordBase} facilates the specification of the simulation domain independend from the actual -evolution module used. -The simulation domain is specified via the parameter file \todo{ES: - ensure this has been defined above}. Cactus distinguishes between +evolution thorn used. +The simulation domain is specified at run time via a parameter file at the +same time as parameters describing the physical system are specified. +Cactus distinguishes between the \emph{physical} domain, which lives in the continuum, and \emph{discrete} domain, which consists of a discrete set of grid points. The physical domain is defined via its coordinate extent and @@ -1618,33 +1619,39 @@ the numerical resolution do not affect the extent of the physical domain, i.e.\ that the discrete domains converge to the physical domain in the limint of infinite resolution. -\codename{CoordBase} exposes a public interface that allows other -modules to query the domain and description in a uniform way. -\codename{Carpet} queries \codename{CoordBase} for the discrete -grid when creating the hierarchy of grids. Evolution modules use this +\codename{CoordBase} exposes a public runtime interface that allows +other +thorns to query the domain description in a uniform way. This is used +by +\codename{Carpet} to query \codename{CoordBase} for the discrete +grid when creating the hierarchy of grids, automatically ensuring a +consistent grid description between the two thorns. +Evolution thorns such as \codename{McLachlan} use the domain +information information to decide which points are evolved and therefore require the evaluation of the righ-hand-side expression and which ones are set via boundary or symmetry conditions. \subsubsection{Symmetries, Boundary Conditions.} -The Einstein Toolkit includes a set of modules, \codename{Boundary} -and \codename{SymBase}, to provide a generic interface to boundary and +The Einstein Toolkit includes a set of thorns, \codename{Boundary} +and \codename{SymBase}, to provide a generic interface to specify and +implement boundary and symmetry conditions. -The toolkit includes built-in suppport for a set of reflecting or +For symmetry conditions the toolkit includes built-in suppport for a +set of reflecting or rotating symmetry conditions that can be used to reduce the size of the simulation domain. These symmetries include periodicity in any of the coordinate directions -(via the \codename{Periodic} module), reflections about the +(via the \codename{Periodic} module), reflections across the coordinate planes (via the \codename{Reflection} module), -and $90^{\circ}$ and $180^{\circ}$ rotational +$90^{\circ}$ and $180^{\circ}$ rotational (via the \codename{RotatingSymmetry90} and \codename{RotatingSymmetry180} modules respectively) -symmetries about the $z$ axis. -\codename{Cartoon2D} provides a continuous -rotational symmetry rather than a discrete -symmetry~\cite{Alcubierre1999ab}. \codename{Cartoon2D} allows fully +symmetries about the $z$ axis, and a continuous rotational symmetry (via +the \codename{Cartoon2D} thorn)~\cite{Alcubierre1999ab}. +\codename{Cartoon2D} allows fully three dimensional codes to be used in axissymmtric problems by evolving a slice in the $y=0$ plane and using the rotational symmetry to populate ghost points off the plane (see Figure~\ref{fig:cartoon-plane}). @@ -1658,15 +1665,16 @@ \label{fig:cartoon-plane} \end{figure} -In applying symmetries, the +In applying symmetries to populate ghost zones, the transformation behaviour of tensorial quantities (including tensor densities and non-tensors such as Christoffel symbols) is correctly taken into account. -The interpolation routines present in \codename{Cactus} automatically +Similarly the interpolation routines present in \codename{Cactus} +automatically take the existing symmetries into account, transparently mapping points to the numerical domain and respecting the transformation behaviour of -tensorial quantities. Symmetries are handled transparently +tensorial quantities. Thus symmetries are handled transparently from the point of view of user modules (see Figure~\ref{fig:faces} for an illustration). \begin{figure}[htbp] @@ -1682,16 +1690,20 @@ \label{fig:faces} \end{figure} -Thorn \codename{Boundary} provides basic boundary conditions. A boundary +Thorn \codename{Boundary} provides basic physical boundary conditions. +A boundary condition suitable for matter fields which approach a constant -(atmosphere) value is provided via either the ``flat'' or ``scalar'' +(atmosphere) value on the boundary is provided via either the ``flat'' +or ``scalar'' (Dirichlet) boundary conditions. These boundary conditions set the values on the boundary to either those close by in the direction of the -normal to the boundary or to a fixed given value, respectively. +normal to the boundary (``flat'') or to a fixed given value +(``scalar''). The ``copy'' boundary condition uses a secondary grid variable to fill in the boundary values, allowing arbitrary boundary conditions to be imposed. -A fall-off (``Robin'') boundary condition is provided suitable for field +\codename{Boundary} provides a fall-off (``Robin'') boundary condition +suitable for fields whose behaviour near the boundary is described by \begin{equation} f(r) = f_0 + \frac{k}{r^n}\text{,} @@ -1709,25 +1721,31 @@ the paper is almost ready)?} \subsubsection{Adaptive Mesh Refinement} -The \codename{Carpet} driver module provides adaptive mesh refinement and -Multi-Block capabilities to modules in the Einstein Toolkit. As in most +The \codename{Carpet} driver thorn provides adaptive mesh refinement +optionally using subcycling in time and level-dependend spatial +refinement factors as well as +Multi-Block capabilities to thorns in the Einstein Toolkit. As in most adaptive mesh refinement frameworks, care has been taken to separate driver related code from physics codes. A typical evolution or analysis code like \codename{McLachlan} or \codename{WeylScal4} is agnostic to the mesh refinement. -Carpet supports feature-based mesh refinement, which is based on +\codename{Carpet} supports feature-based mesh refinement, which is +based on extracting e.g.\ the locations of black holes or neutron stars, and then constructing a mesh hierarchy (stacks of refined regions) based on the locations, sizes, and speeds of these objects. This allows tracking objects as they move through the domain. One can also add or -remove stacks if e.g.\ the number of objects changes. If initial +remove stacks if e.g.\ the number of objects changes. Fully AMR based on +a local error estimate is supported by \codename{Carpet} however the +Einstein Toolkit does not provide a suitable redgridding thorn to create +such a grid. If initial conditions are constructed outside of Carpet (which is often the case), then the initial mesh hierarchy has to be defined manually. -The Einstein toolkit currently provides two ``regridding'' modules in +The Einstein toolkit currently provides two regridding modules in the \codename{CarpetRegrid} and \codename{CarpetRegrid2} thorns. Both thorns primarily support box-in-box type refined meshes, which is well suited to current binary black hole simulations in which the -high-resolution region is centered on the individual black holes. +high-resolution regions are centered on the individual black holes. Figure~\ref{fig:bbh-boxes} shows a typical set of nested boxes during the inspiral phase of a binary black hole merger simulation. \begin{figure}[htbp] @@ -1742,18 +1760,20 @@ \codename{RotatingSymmetry180} to reduce the computational domain.} \label{fig:bbh-boxes} \end{figure} + \codename{CarpetRegrid} provides a number of different ways to specify the refined regions, either as a set of boxes centered around the origin or as an explicit list of boxes that make up the grid hierarchy. Traditionally groups using \codename{CarpetRegrid} have employed auxiliary thorns, that are not part of the Einstein Toolkit, to create -the list of boxes based on information obtained e.g. from apparent +this list of boxes based on information obtained e.g.\ from apparent horizon tracking. \codename{CarpetRegrid2} provides a user friendly interface to define sets of nested boxes that follow black holes or other tracked objects. \codename{CarpetRegrid2} supports up to 10 sets -of nested regions which can be either moving or stationary. The number -of refinement levels in each region as well as the radii of each nested -box are allowed to change during runtime, making it possible to e.g. +of nested boxes which can be either moving or stationary. The number +of refinement levels in each set of boxes as well as the radii of each +nested +box are allowed to change during runtime, making it possible to e.g.\ adapt the shape of the refined region to the surface of a star. \codename{CarpetRegrid2} contains code to handle the $\pi$-symmetry provided by \codename{RotatingSymmetry180}, enforcing the symmetry on @@ -1777,12 +1797,13 @@ only as a source for the boundary condition.} \label{fig:rot180-grid} \end{figure} + \codename{CarpetTracker} provides a simple interface to slave \codename{CarpetRegrid2}'s regions to the object trackers \codename{PunctureTracker} and \codename{NSTracker} (see -section~\ref{sec:object-tracking}) by copying the position of tracker +section~\ref{sec:object-tracking}) by copying the position of tracked objects -positions out of a \codename{SphericalSurface} into +from a \codename{SphericalSurface} into \codename{CarpetRegrid2}.

1 0

[svn:einsteintoolkit] www/about/members/ (Rev. 587)
by knarf＠cct.lsu.edu 06 May '11

06 May '11

User: knarf Date: 2011/05/06 11:38 PM Modified: /about/members/ members.txt Log: add David Garrison File Changes: Directory: /about/members/ ========================== File [modified]: members.txt Delta lines: +3 -0 =================================================================== --- about/members/members.txt 2011-05-07 04:31:31 UTC (rev 586) +++ about/members/members.txt 2011-05-07 04:38:08 UTC (rev 587) @@ -94,6 +94,9 @@ Universidad Nacional AutÃ³noma de MÃ©xico Jose Manuel Torres +University of Houston Clear Lake + David Garrison + University of Maryland Bruno Giacomazzo

1 0

[svn:einsteintoolkit] www/documentation/ (Rev. 586)
by knarf＠cct.lsu.edu 06 May '11

06 May '11

User: knarf Date: 2011/05/06 11:31 PM Modified: /documentation/ request_email Log: reformat File Changes: Directory: /documentation/ ========================== File [modified]: request_email Delta lines: +14 -13 =================================================================== --- documentation/request_email 2011-05-07 04:20:25 UTC (rev 585) +++ documentation/request_email 2011-05-07 04:31:31 UTC (rev 586) @@ -1,11 +1,13 @@ Thank you for your request for an Einstein Toolkit demonstration account. Please go to <https://allocations.loni.org/login_request.php> -and request a LONI account. Enter "Frank Loeffler" as your LONI -Contact/Collaborator. Enter "Einstein Toolkit" as your research area. -You do not need to enter a Globus DN. +and request a LONI account. Use an email address connected to your +afficilation, independent providers like yahoo or google will not be +accepted. Enter "Frank Loeffler" as your LONI Contact/Collaborator. +Enter "Einstein Toolkit" as your research area. You do not need to +enter a Globus DN. -(If you read XX here, this email was not properly created, it should - be two numbers!) +(If you read XX here, this email was not properly created, it should be two +numbers!) Please use the user name 'ettestXX'. This account is only meant for testing the Einstein Toolkit and will be removed within one to two months, without warning, and possibly including all stored data. @@ -13,9 +15,9 @@ Once you have received your account (hopefully within one or two days), go to <https://allocations.loni.org/balances.php> and log in. Select "Request Allocation" on the right hand side, and then "Join -Allocation". Enter "Frank Loeffler" and then select the -"loni_cactus06" allocation. This will be approved within one or two -days, and you will receive an email confirming the allocation. +Allocation". Enter "Frank Loeffler" and then select the "loni_cactus06" +allocation. This will be approved within one or two days, and you will +receive an email confirming the allocation. Note that this allocation is reserved for running the examples for the Einstein Toolkit, please limit your total usage to 2 kSU. If you need @@ -26,9 +28,8 @@ and files. The user guide for Queen Bee is at <http://www.loni.org/teragrid/users_guide.php>. -If you continue to use the Einstein Toolkit, please join our -consortium by signing up at -<http://einsteintoolkit.org/about/members/join>. Please also write to -<users(a)einsteintoolkit.org> if you have problems or questions about -the examples. +If you continue to use the Einstein Toolkit, please join our consortium +by signing up at <http://einsteintoolkit.org/about/members/join>. +Please also write to <users(a)einsteintoolkit.org> if you have problems +or questions about the examples.

1 0

[svn:einsteintoolkit] www/documentation/ (Rev. 585)
by knarf＠cct.lsu.edu 06 May '11

06 May '11

User: knarf Date: 2011/05/06 11:20 PM Added: /documentation/ request_email Log: template for request answer for loni account File Changes: Directory: /documentation/ ========================== File [added]: request_email Delta lines: +34 -0 =================================================================== --- documentation/request_email (rev 0) +++ documentation/request_email 2011-05-07 04:20:25 UTC (rev 585) @@ -0,0 +1,34 @@ +Thank you for your request for an Einstein Toolkit demonstration +account. Please go to <https://allocations.loni.org/login_request.php> +and request a LONI account. Enter "Frank Loeffler" as your LONI +Contact/Collaborator. Enter "Einstein Toolkit" as your research area. +You do not need to enter a Globus DN. + +(If you read XX here, this email was not properly created, it should + be two numbers!) +Please use the user name 'ettestXX'. This account is only meant for +testing the Einstein Toolkit and will be removed within one to two +months, without warning, and possibly including all stored data. + +Once you have received your account (hopefully within one or two +days), go to <https://allocations.loni.org/balances.php> and log in. +Select "Request Allocation" on the right hand side, and then "Join +Allocation". Enter "Frank Loeffler" and then select the +"loni_cactus06" allocation. This will be approved within one or two +days, and you will receive an email confirming the allocation. + +Note that this allocation is reserved for running the examples for the +Einstein Toolkit, please limit your total usage to 2 kSU. If you need +more resources once you reach 2 kSU, feel free to contact the toolkit +maintainers for suggestions on where you might find a real allocation +for your work. We also may remove you from the allocation if you have +used over 2 kSU, but you will still be able to access all your data +and files. The user guide for Queen Bee is at +<http://www.loni.org/teragrid/users_guide.php>. + +If you continue to use the Einstein Toolkit, please join our +consortium by signing up at +<http://einsteintoolkit.org/about/members/join>. Please also write to +<users(a)einsteintoolkit.org> if you have problems or questions about +the examples. +

1 0

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits May 2011