Commits June 2016

commits@lists.einsteintoolkit.org

7 participants
129 discussions

[svn:einsteintoolkit] www/about/gallery/ (Rev. 1527)
by barry.wardell＠gmail.com 05 Jun '16

05 Jun '16

User: barry.wardell Date: 2016/06/05 05:04 PM Added: /about/gallery/NsNsToHMNS/ index.php Modified: /about/gallery/ index.php Log: gallery: Move BNS to a separate page File Changes: Directory: /about/gallery/NsNsToHMNS/ ===================================== File [added]: index.php Delta lines: +30 -0 =================================================================== --- about/gallery/NsNsToHMNS/index.php (rev 0) +++ about/gallery/NsNsToHMNS/index.php 2016-06-05 22:04:53 UTC (rev 1527) @@ -0,0 +1,30 @@ +<?php $title="Gallery: Binary, inspiraling neutron stars forming a hypermassive neutron star"; +include_once($_SERVER['DOCUMENT_ROOT'].'/global/header.php');?> +<script type="text/javascript" src="/global/jwplayer/jwplayer.js"></script> + + <p>This simulation shows how to evolve a pair of neutron stars, reading in + initial data provided by the <a href="http://www.lorene.obspm.fr/">LORENE</a> code. The video shows + the rest mass density in an equatorial slice of the simulation. The + emitted gravitational waves of the inspiral and due to the oscillations in + the formed hypermassive neutron star are shown as well. </p> + <table> + <tr><th>Parameter file</th><td><a href="NsNsToHMNS.par">NsNsToHMNS.par</a></td></tr> + <tr><th>Thornlist</th><td><a href="NsNsToHMNS.th">NsNsToHMNS.th</a> (Thorn NSTracker added to release)</td></tr> + <tr><th>Initial data file</th><td><a href="G2_I12vs12_D4R33T21_45km.resu.xz">G2_I12vs12_D4R33T21_45km.resu.xz</a> (uncompress using <a href="http://en.wikipedia.org/wiki/XZ_Utils">unxz</a>)</td></tr> + <tr><th>Support scripts</th><td><a href="scripts.tar.gz">scripts.tar.gz</a></td></tr> + <tr><th>approx. memory</th><td>8.8 GB</td></tr> + <tr><th>approx. runtime</th><td>12*48 SU</td></tr> + </table> + <a href="mp_Psi4_l2_m2_r300.00.asc"><img src="mp_Psi4_l2_m2_r300.00.png" alt="Psi_4^{2,2} at r=300M over time"></a> + <div id="NsNsToHMNS_rho">Loading the player...</div> + Download movie: <a href="rho.mp4">rho.mp4</a> + <script type="text/javascript"> + jwplayer("NsNsToHMNS_rho").setup({ + file: "rho.mp4", + image: "rho.png", + width: 640, height: 480 + }); + </script> + + +<?php include_once($_SERVER['DOCUMENT_ROOT'].'/global/footer.php');?> Directory: /about/gallery/ ========================== File [modified]: index.php Delta lines: +3 -25 =================================================================== --- about/gallery/index.php 2016-06-05 22:04:49 UTC (rev 1526) +++ about/gallery/index.php 2016-06-05 22:04:53 UTC (rev 1527) @@ -22,38 +22,16 @@ <p> <ul> <li><a href="tov">Single, stable neutron star</a></li> +<li><a href="NsNsToHMNS">Binary, inspiraling neutron stars forming a hypermassive neutron star</a></li> + </ul> </p>  - <h4> Binary, inspiraling neutron stars forming a hypermassive neutron star</h4> - <p>This simulation shows how to evolve a pair of neutron stars, reading in - initial data provided by the <a href="http://www.lorene.obspm.fr/">LORENE</a> code. The video shows - the rest mass density in an equatorial slice of the simulation. The - emitted gravitational waves of the inspiral and due to the oscillations in - the formed hypermassive neutron star are shown as well. </p> - <table> - <tr><th>Parameter file</th><td><a href="NsNsToHMNS/NsNsToHMNS.par">NsNsToHMNS.par</a></td></tr> - <tr><th>Thornlist</th><td><a href="NsNsToHMNS/NsNsToHMNS.th">NsNsToHMNS.th</a> (Thorn NSTracker added to release)</td></tr> - <tr><th>Initial data file</th><td><a href="NsNsToHMNS/G2_I12vs12_D4R33T21_45km.resu.xz">G2_I12vs12_D4R33T21_45km.resu.xz</a> (uncompress using <a href="http://en.wikipedia.org/wiki/XZ_Utils">unxz</a>)</td></tr> - <tr><th>Support scripts</th><td><a href="NsNsToHMNS/scripts.tar.gz">scripts.tar.gz</a></td></tr> - <tr><th>approx. memory</th><td>8.8 GB</td></tr> - <tr><th>approx. runtime</th><td>12*48 SU</td></tr> - </table> - <a href="NsNsToHMNS/mp_Psi4_l2_m2_r300.00.asc"><img src="NsNsToHMNS/mp_Psi4_l2_m2_r300.00.png" alt="Psi_4^{2,2} at r=300M over time"></a> - <div id="NsNsToHMNS_rho">Loading the player...</div> - Download movie: <a href="NsNsToHMNS/rho.mp4">rho.mp4</a> - <script type="text/javascript"> - jwplayer("NsNsToHMNS_rho").setup({ - file: "NsNsToHMNS/rho.mp4", - image: "NsNsToHMNS/rho.png", - width: 640, height: 480 - }); - </script> - <h4>Poisson equation</h4> +<h4>Poisson equation</h4> <p>This example shows how to solve the Poisson equation with a generic source term using

1 0

[svn:einsteintoolkit] www/about/gallery/tov/ (Rev. 1526)
by barry.wardell＠gmail.com 05 Jun '16

05 Jun '16

User: barry.wardell Date: 2016/06/05 05:04 PM Modified: /about/gallery/tov/ index.php Log: tov: Correct title File Changes: Directory: /about/gallery/tov/ ============================== File [modified]: index.php Delta lines: +1 -2 =================================================================== --- about/gallery/tov/index.php 2016-06-05 22:04:43 UTC (rev 1525) +++ about/gallery/tov/index.php 2016-06-05 22:04:49 UTC (rev 1526) @@ -1,8 +1,7 @@ -<?php $title="Gallery: Binary black hole GW150914"; +<?php $title="Gallery: Single stable neutron star"; include_once($_SERVER['DOCUMENT_ROOT'].'/global/header.php');?> <script type="text/javascript" src="/global/jwplayer/jwplayer.js"></script> - <h4> Single, stable neutron star</h4> <p>This simulation shows how to evolve a stable, single neutron star. It uses a TOV solution as initial data and follows the initial (numerical) perturbations, writing the maximum density to disk. What can be seen is that the oscillates for a few cycles (only very

1 0

[svn:einsteintoolkit] www/about/gallery/ (Rev. 1525)
by barry.wardell＠gmail.com 05 Jun '16

05 Jun '16

User: barry.wardell Date: 2016/06/05 05:04 PM Added: /about/gallery/tov/ index.php, tov.par Removed: /about/gallery/ tov.par Modified: /about/gallery/ index.php Log: gallery: Move TOV into separate page File Changes: Directory: /about/gallery/ ========================== File [modified]: index.php Delta lines: +7 -17 =================================================================== --- about/gallery/index.php 2016-06-05 22:04:37 UTC (rev 1524) +++ about/gallery/index.php 2016-06-05 22:04:43 UTC (rev 1525) @@ -18,23 +18,13 @@  - <h4> Single, stable neutron star</h4> - <p>This simulation shows how to evolve a stable, single neutron star. It uses a TOV solution - as initial data and follows the initial (numerical) perturbations, writing the maximum - density to disk. What can be seen is that the oscillates for a few cycles (only very - slitely because the solution is stationary and the induced numerical perturbation is small), - before numerical errors due to the lack of resolution spoil the results. The size of this - simulation is intentionally kept as small as possible, thus the resolution is about as low - as youcan go without loosing any confidence in at least quantitatively correct results. - It is intended to be used as quick examples in lectures, workshops or summer schools where - both time and computational power are typically scarce.</p> - <table> - <tr><th>Parameter file</th><td><a href="tov.par">tov.par</a></td></tr> - <tr><th>Thornlist</th><td><a href="https://bitbucket.org/einsteintoolkit/manifest/raw/ET_2015_11/einsteintoolk…">ET Somerville release</a></td></tr> - <tr><th>approx. memory</th><td>0.25 GB</td></tr> - <tr><th>approx. runtime</th><td>2.5 SU</td></tr> - </table> - <a href="tov/rho_max.pdf"><img src="tov/rho_max.png" alt="Maximum density over time"></a> + +<p> +<ul> +<li><a href="tov">Single, stable neutron star</a></li> +</ul> +</p> +  File [removed]: tov.par Delta lines: +0 -287 =================================================================== --- about/gallery/tov.par 2016-06-05 22:04:37 UTC (rev 1524) +++ about/gallery/tov.par 2016-06-05 22:04:43 UTC (rev 1525) @@ -1,287 +0,0 @@ -# Example parameter file for a static TOV star. Everything is evolved, but -# because this is a solution to the GR and hydro equations, nothing changes -# much. What can be seen is the initial perturbation (due to numerical errors) -# ringing down (look at the density maximum), and later numerical errors -# governing the solution. Try higher resolutions to decrease this error. - -# Some basic stuff -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" - -# Finalize -Cactus::terminate = "time" -Cactus::cctk_final_time = 800 # devide by ~203 to get ms - -# Termination Trigger -ActiveThorns = "TerminationTrigger" -TerminationTrigger::max_walltime = 24 # hours -TerminationTrigger::on_remaining_walltime = 15 # minutes -TerminationTrigger::check_file_every = 512 -TerminationTrigger::termination_file = "TerminationTrigger.txt" -TerminationTrigger::termination_from_file = "yes" -TerminationTrigger::create_termination_file = "yes" - -# grid parameters -Carpet::domain_from_coordbase = "yes" -CartGrid3D::type = "coordbase" -CartGrid3D::domain = "full" -CartGrid3D::avoid_origin = "no" -CoordBase::xmin = 0.0 -CoordBase::ymin = 0.0 -CoordBase::zmin = 0.0 -CoordBase::xmax = 48.0 -CoordBase::ymax = 48.0 -CoordBase::zmax = 48.0 -# Change these parameters to change resolution. The ?max settings above -# have to be multiples of these. 'dx' is the size of one cell in x-direction. -# Making this smaller means using higher resolution, because more points will -# be used to cover the same space. -CoordBase::dx = 4.0 -CoordBase::dy = 4.0 -CoordBase::dz = 4.0 - -CarpetRegrid2::regrid_every = 0 -CarpetRegrid2::num_centres = 1 -CarpetRegrid2::num_levels_1 = 3 -CarpetRegrid2::radius_1[1] = 24.0 -CarpetRegrid2::radius_1[2] = 12.0 - -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" - -# storage and coupling -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 - -SpaceMask::use_mask = "yes" - -Carpet::enable_all_storage = no -Carpet::use_buffer_zones = "yes" - -Carpet::poison_new_timelevels = "yes" -Carpet::check_for_poison = "no" - -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 - -# Time integration -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 -# This can save you computing time, so it's not a bad idea to do this -# once in a whioe. -ActiveThorns = "NaNChecker" -NaNChecker::check_every = 16384 -NaNChecker::action_if_found = "terminate" #"terminate", "just warn", "abort" -NaNChecker::check_vars = "ADMBase::metric ADMBase::lapse ADMBase::shift HydroBase::rho HydroBase::eps HydroBase::press HydroBase::vel" - -# Hydro paramters - -ActiveThorns = "EOS_Omni GRHydro" - -HydroBase::evolution_method = "GRHydro" - -GRHydro::riemann_solver = "HLLE" -GRHydro::GRHydro_eos_type = "General" -GRHydro::GRHydro_eos_table = "Ideal_Fluid" -GRHydro::recon_method = "ppm" -GRHydro::GRHydro_stencil = 3 -GRHydro::bound = "none" -GRHydro::rho_abs_min = 1.e-10 -GRHydro::GRHydro_atmo_tolerance = 1.e-3 - -# Curvature evolution parameters - -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" - -ML_BSSN::timelevels = 3 - -ML_BSSN::harmonicN = 1 # 1+log -ML_BSSN::harmonicF = 2.0 # 1+log -ML_BSSN::ShiftBCoeff = 1 -ML_BSSN::ShiftGammaCoeff = 0.75 -ML_BSSN::BetaDriver = 2.66 -ML_BSSN::LapseAdvectionCoeff = 0.0 -ML_BSSN::ShiftAdvectionCoeff = 0.0 - -ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas" -ML_BSSN::my_rhs_boundary_condition = "NewRad" - -# Some dissipation to get rid of high-frequency noise -ActiveThorns = "SphericalSurface Dissipation" -Dissipation::verbose = "no" -Dissipation::epsdis = 0.01 -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 -" - - -# init parameters -InitBase::initial_data_setup_method = "init_some_levels" - -# Use TOV as initial data -ActiveThorns = "TOVSolver" - -HydroBase::initial_hydro = "tov" -ADMBase::initial_data = "tov" -ADMBase::initial_lapse = "tov" -ADMBase::initial_shift = "tov" -ADMBase::initial_dtlapse = "zero" -ADMBase::initial_dtshift = "zero" - -# Parameters for initial star -TOVSolver::TOV_Rho_Central[0] = 1.28e-3 -TOVSolver::TOV_Gamma = 2 -TOVSolver::TOV_K = 100 - -# Set equation of state for evolution -EOS_Omni::poly_gamma = 2 -EOS_Omni::poly_k = 100 -EOS_Omni::gl_gamma = 2 -EOS_Omni::gl_k = 100 - -# I/O - -cactus::cctk_timer_output = "full" - -# Use (create if necessary) an output directory named like the -# parameter file (minus the .par) -IO::out_dir = ${parfile} - -# Write one file overall per output (variable/group) -# In production runs, comment this or set to "proc" to get one file -# per MPI process -IO::out_mode = "onefile" - -# Some screen output -IOBasic::outInfo_every = 16 -IOBasic::outInfo_vars = "Carpet::physical_time_per_hour HydroBase::rho{reductions='maximum'}" - -# Scalar output -IOScalar::outScalar_every = 512 -IOScalar::one_file_per_group = "yes" -IOScalar::outScalar_reductions = "norm1 norm2 norm_inf sum maximum minimum" -IOScalar::outScalar_vars = " - HydroBase::rho{reductions='maximum'} - HydroBase::press{reductions='maximum'} - HydroBase::eps{reductions='minimum maximum'} - HydroBase::vel{reductions='minimum maximum'} - HydroBase::w_lorentz{reductions='minimum maximum'} - ADMBase::lapse{reductions='minimum maximum'} - ADMBase::shift{reductions='minimum maximum'} - ML_BSSN::ML_Ham{reductions='norm1 norm2 maximum minimum norm_inf'} - ML_BSSN::ML_mom{reductions='norm1 norm2 maximum minimum norm_inf'} - GRHydro::dens{reductions='minimum maximum sum'} - Carpet::timing{reductions='average'} -" - -# 1D ASCII output. Disable for production runs! -IOASCII::out1D_every = 512 -IOASCII::one_file_per_group = yes -IOASCII::output_symmetry_points = no -IOASCII::out1D_vars = " - HydroBase::rho - HydroBase::press - HydroBase::eps - HydroBase::vel - ADMBase::lapse - ADMBase::metric - ADMBase::curv - ML_BSSN::ML_Ham - ML_BSSN::ML_mom -" - -# 2D HDF5 output -CarpetIOHDF5::output_buffer_points = "no" - -CarpetIOHDF5::out2D_every = 32768 -CarpetIOHDF5::out2D_vars = " - HydroBase::rho - HydroBase::eps - HydroBase::vel - HydroBase::w_lorentz - ADMBase::lapse - ADMBase::shift - ADMBase::metric - ML_BSSN::ML_Ham - ML_BSSN::ML_mom - " - -# Checkpointing options -IOHDF5::checkpoint = "yes" -IO::checkpoint_dir = $parfile -IO::recover_dir = $parfile -IO::recover = "autoprobe" -IO::checkpoint_ID = "yes" -IO::checkpoint_every = 1048576 -IO::checkpoint_keep = 3 -IO::checkpoint_on_terminate = "yes" - -# Enable to get detailed timing information -#ActiveThorns = "TimerReport" -#TimerReport::out_every = 1024 -#TimerReport::out_filename = "TimerReport" -#TimerReport::output_all_timers_readable = "yes" -#TimerReport::output_all_timers = "yes" - -# Enable for profiling -#Carpet::output_timers_every = 1024 - Directory: /about/gallery/tov/ ============================== File [added]: index.php Delta lines: +23 -0 =================================================================== --- about/gallery/tov/index.php (rev 0) +++ about/gallery/tov/index.php 2016-06-05 22:04:43 UTC (rev 1525) @@ -0,0 +1,23 @@ +<?php $title="Gallery: Binary black hole GW150914"; +include_once($_SERVER['DOCUMENT_ROOT'].'/global/header.php');?> +<script type="text/javascript" src="/global/jwplayer/jwplayer.js"></script> + + <h4> Single, stable neutron star</h4> + <p>This simulation shows how to evolve a stable, single neutron star. It uses a TOV solution + as initial data and follows the initial (numerical) perturbations, writing the maximum + density to disk. What can be seen is that the oscillates for a few cycles (only very + slitely because the solution is stationary and the induced numerical perturbation is small), + before numerical errors due to the lack of resolution spoil the results. The size of this + simulation is intentionally kept as small as possible, thus the resolution is about as low + as youcan go without loosing any confidence in at least quantitatively correct results. + It is intended to be used as quick examples in lectures, workshops or summer schools where + both time and computational power are typically scarce.</p> + <table> + <tr><th>Parameter file</th><td><a href="tov.par">tov.par</a></td></tr> + <tr><th>Thornlist</th><td><a href="https://bitbucket.org/einsteintoolkit/manifest/raw/ET_2015_11/einsteintoolk…">ET Somerville release</a></td></tr> + <tr><th>approx. memory</th><td>0.25 GB</td></tr> + <tr><th>approx. runtime</th><td>2.5 SU</td></tr> + </table> + <a href="rho_max.pdf"><img src="rho_max.png" alt="Maximum density over time"></a> + +<?php include_once($_SERVER['DOCUMENT_ROOT'].'/global/footer.php');?> File [added]: tov.par Delta lines: +287 -0 =================================================================== --- about/gallery/tov/tov.par (rev 0) +++ about/gallery/tov/tov.par 2016-06-05 22:04:43 UTC (rev 1525) @@ -0,0 +1,287 @@ +# Example parameter file for a static TOV star. Everything is evolved, but +# because this is a solution to the GR and hydro equations, nothing changes +# much. What can be seen is the initial perturbation (due to numerical errors) +# ringing down (look at the density maximum), and later numerical errors +# governing the solution. Try higher resolutions to decrease this error. + +# Some basic stuff +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" + +# Finalize +Cactus::terminate = "time" +Cactus::cctk_final_time = 800 # devide by ~203 to get ms + +# Termination Trigger +ActiveThorns = "TerminationTrigger" +TerminationTrigger::max_walltime = 24 # hours +TerminationTrigger::on_remaining_walltime = 15 # minutes +TerminationTrigger::check_file_every = 512 +TerminationTrigger::termination_file = "TerminationTrigger.txt" +TerminationTrigger::termination_from_file = "yes" +TerminationTrigger::create_termination_file = "yes" + +# grid parameters +Carpet::domain_from_coordbase = "yes" +CartGrid3D::type = "coordbase" +CartGrid3D::domain = "full" +CartGrid3D::avoid_origin = "no" +CoordBase::xmin = 0.0 +CoordBase::ymin = 0.0 +CoordBase::zmin = 0.0 +CoordBase::xmax = 48.0 +CoordBase::ymax = 48.0 +CoordBase::zmax = 48.0 +# Change these parameters to change resolution. The ?max settings above +# have to be multiples of these. 'dx' is the size of one cell in x-direction. +# Making this smaller means using higher resolution, because more points will +# be used to cover the same space. +CoordBase::dx = 4.0 +CoordBase::dy = 4.0 +CoordBase::dz = 4.0 + +CarpetRegrid2::regrid_every = 0 +CarpetRegrid2::num_centres = 1 +CarpetRegrid2::num_levels_1 = 3 +CarpetRegrid2::radius_1[1] = 24.0 +CarpetRegrid2::radius_1[2] = 12.0 + +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" + +# storage and coupling +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 + +SpaceMask::use_mask = "yes" + +Carpet::enable_all_storage = no +Carpet::use_buffer_zones = "yes" + +Carpet::poison_new_timelevels = "yes" +Carpet::check_for_poison = "no" + +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 + +# Time integration +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 +# This can save you computing time, so it's not a bad idea to do this +# once in a whioe. +ActiveThorns = "NaNChecker" +NaNChecker::check_every = 16384 +NaNChecker::action_if_found = "terminate" #"terminate", "just warn", "abort" +NaNChecker::check_vars = "ADMBase::metric ADMBase::lapse ADMBase::shift HydroBase::rho HydroBase::eps HydroBase::press HydroBase::vel" + +# Hydro paramters + +ActiveThorns = "EOS_Omni GRHydro" + +HydroBase::evolution_method = "GRHydro" + +GRHydro::riemann_solver = "HLLE" +GRHydro::GRHydro_eos_type = "General" +GRHydro::GRHydro_eos_table = "Ideal_Fluid" +GRHydro::recon_method = "ppm" +GRHydro::GRHydro_stencil = 3 +GRHydro::bound = "none" +GRHydro::rho_abs_min = 1.e-10 +GRHydro::GRHydro_atmo_tolerance = 1.e-3 + +# Curvature evolution parameters + +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" + +ML_BSSN::timelevels = 3 + +ML_BSSN::harmonicN = 1 # 1+log +ML_BSSN::harmonicF = 2.0 # 1+log +ML_BSSN::ShiftBCoeff = 1 +ML_BSSN::ShiftGammaCoeff = 0.75 +ML_BSSN::BetaDriver = 2.66 +ML_BSSN::LapseAdvectionCoeff = 0.0 +ML_BSSN::ShiftAdvectionCoeff = 0.0 + +ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas" +ML_BSSN::my_rhs_boundary_condition = "NewRad" + +# Some dissipation to get rid of high-frequency noise +ActiveThorns = "SphericalSurface Dissipation" +Dissipation::verbose = "no" +Dissipation::epsdis = 0.01 +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 +" + + +# init parameters +InitBase::initial_data_setup_method = "init_some_levels" + +# Use TOV as initial data +ActiveThorns = "TOVSolver" + +HydroBase::initial_hydro = "tov" +ADMBase::initial_data = "tov" +ADMBase::initial_lapse = "tov" +ADMBase::initial_shift = "tov" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +# Parameters for initial star +TOVSolver::TOV_Rho_Central[0] = 1.28e-3 +TOVSolver::TOV_Gamma = 2 +TOVSolver::TOV_K = 100 + +# Set equation of state for evolution +EOS_Omni::poly_gamma = 2 +EOS_Omni::poly_k = 100 +EOS_Omni::gl_gamma = 2 +EOS_Omni::gl_k = 100 + +# I/O + +cactus::cctk_timer_output = "full" + +# Use (create if necessary) an output directory named like the +# parameter file (minus the .par) +IO::out_dir = ${parfile} + +# Write one file overall per output (variable/group) +# In production runs, comment this or set to "proc" to get one file +# per MPI process +IO::out_mode = "onefile" + +# Some screen output +IOBasic::outInfo_every = 16 +IOBasic::outInfo_vars = "Carpet::physical_time_per_hour HydroBase::rho{reductions='maximum'}" + +# Scalar output +IOScalar::outScalar_every = 512 +IOScalar::one_file_per_group = "yes" +IOScalar::outScalar_reductions = "norm1 norm2 norm_inf sum maximum minimum" +IOScalar::outScalar_vars = " + HydroBase::rho{reductions='maximum'} + HydroBase::press{reductions='maximum'} + HydroBase::eps{reductions='minimum maximum'} + HydroBase::vel{reductions='minimum maximum'} + HydroBase::w_lorentz{reductions='minimum maximum'} + ADMBase::lapse{reductions='minimum maximum'} + ADMBase::shift{reductions='minimum maximum'} + ML_BSSN::ML_Ham{reductions='norm1 norm2 maximum minimum norm_inf'} + ML_BSSN::ML_mom{reductions='norm1 norm2 maximum minimum norm_inf'} + GRHydro::dens{reductions='minimum maximum sum'} + Carpet::timing{reductions='average'} +" + +# 1D ASCII output. Disable for production runs! +IOASCII::out1D_every = 512 +IOASCII::one_file_per_group = yes +IOASCII::output_symmetry_points = no +IOASCII::out1D_vars = " + HydroBase::rho + HydroBase::press + HydroBase::eps + HydroBase::vel + ADMBase::lapse + ADMBase::metric + ADMBase::curv + ML_BSSN::ML_Ham + ML_BSSN::ML_mom +" + +# 2D HDF5 output +CarpetIOHDF5::output_buffer_points = "no" + +CarpetIOHDF5::out2D_every = 32768 +CarpetIOHDF5::out2D_vars = " + HydroBase::rho + HydroBase::eps + HydroBase::vel + HydroBase::w_lorentz + ADMBase::lapse + ADMBase::shift + ADMBase::metric + ML_BSSN::ML_Ham + ML_BSSN::ML_mom + " + +# Checkpointing options +IOHDF5::checkpoint = "yes" +IO::checkpoint_dir = $parfile +IO::recover_dir = $parfile +IO::recover = "autoprobe" +IO::checkpoint_ID = "yes" +IO::checkpoint_every = 1048576 +IO::checkpoint_keep = 3 +IO::checkpoint_on_terminate = "yes" + +# Enable to get detailed timing information +#ActiveThorns = "TimerReport" +#TimerReport::out_every = 1024 +#TimerReport::out_filename = "TimerReport" +#TimerReport::output_all_timers_readable = "yes" +#TimerReport::output_all_timers = "yes" + +# Enable for profiling +#Carpet::output_timers_every = 1024 +

1 0

[svn:einsteintoolkit] www/about/gallery/ (Rev. 1524)
by barry.wardell＠gmail.com 05 Jun '16

05 Jun '16

User: barry.wardell Date: 2016/06/05 05:04 PM Added: /about/gallery/gw150914/ index.php Modified: /about/gallery/ index.php Log: gallery: Move gw150914 to a separate page File Changes: Directory: /about/gallery/ ========================== File [modified]: index.php Delta lines: +4 -47 =================================================================== --- about/gallery/index.php 2016-06-05 22:04:33 UTC (rev 1523) +++ about/gallery/index.php 2016-06-05 22:04:37 UTC (rev 1524) @@ -90,55 +90,12 @@ <p><img src="poisson/err.png" alt="Error" width="40%" align="top"></p> <h3>Using external but public code</h3> - - <h4> Binary black hole GW150914, as observed by LIGO; the first detection of gravitational waves</h4> - <p>This simulation shows how to evolve a binary black hole system for the last 6 orbits through to merger. </p> - - <h5>Physical parameters</h5> - <table> - <tr><th>Initial separation D</th><td>10 M [TODO: give in Schwarzschild radii?]</td></tr> - <tr><th>Mass ratio q = m1/m2</th><td>36/29 ~ 1.24</td></tr> - <tr><th>Spin chi1 = a1/m1</th><td>0.31</td></tr> - <tr><th>Spin chi2 = a2/m2</th><td>-0.46</td></tr> - </table> - - <h5>Physical properties</h5> - <table> - <tr><th>Number of orbits</th><td>6</td></tr> - <tr><th>Time to merger</th><td>900 M [TODO: give in physical units for GW150914?]</td></tr> - <tr><th>Mass of final BH</th><td></td></tr> - <tr><th>Spin of final BH</th><td></td></tr> - </table> - - <h5>Computational details</h5> - <?php - $sim_stats = json_decode(file_get_contents("gw150914/stats.json"), true); - ?></p> - <table> - <tr><th>Parameter file</th><td><a href="GW150914/GW150914.rpar">GW150914.rpar</a></td></tr> - <tr><th>Thornlist</th><td><a href="GW150914/GW150914.th">GW150914.th</a> (Llama multi-block code added to ET release version)</td></tr> - <tr><th>Submission command</th><td>simfactory/bin/sim create-submit GW150914_28 --define N 28 --parfile repos/GW150914/ParameterFiles/GW150914.rpar --procs 128 --walltime 24:00:00</td></tr> - <tr><th>Total memory</th><td><?php echo $sim_stats['TotalMemory']?> GB</td></tr> - <tr><th>Run time</th><td><?php echo $sim_stats['Runtime']?> days on <?php echo $sim_stats['Cores']?> cores (Intel(R) Xeon(R) CPU E5-2630 v3 @ 2.40GHz)</td></tr> - <tr><th>Cost</th><td><?php echo $sim_stats['CostInCoreHours']?> core hours</td></tr> - </table> - <p> - - <h5>Horizon centroid coordinate trajectories</h5> - <a href="gw150914/h22.asc"><img src="gw150914/trajectories.png" alt="BH trajectories" width="242" height="237"></a> +<ul> +<li><a href="gw150914">Binary black hole GW150914</a></li> +</ul> +</p> - <h5>Gravitational waveform</h5> - <a href="gw150914/h22.asc"><img src="gw150914/waveform.png" alt="h_+^{2,2} at r=500M as a function of retarded time" width="542" height="191"></a> - - <h5>Tutorials</h5> - <ul> - <li><a href="https://bitbucket.org/barrywardell/gw150914sim/wiki/Simulation%20tutorial">Compile and run</a>: - Compile the code and run the simulation</li> - <li><a href="https://bitbucket.org/barrywardell/gw150914sim/wiki/Analysis%20tutorial">Analysis</a>: Analyse the results using SimulationTools for Mathematica</li> - <li><a href="https://bitbucket.org/barrywardell/gw150914sim/wiki/Visualising%20the%20spa…">Visualisation</a>: Visualise the data using VisIt</li> - </ul> -  Directory: /about/gallery/gw150914/ =================================== File [added]: index.php Delta lines: +59 -0 =================================================================== --- about/gallery/gw150914/index.php (rev 0) +++ about/gallery/gw150914/index.php 2016-06-05 22:04:37 UTC (rev 1524) @@ -0,0 +1,59 @@ +<?php $title="Gallery: Binary black hole GW150914"; +include_once($_SERVER['DOCUMENT_ROOT'].'/global/header.php');?> +<script type="text/javascript" src="/global/jwplayer/jwplayer.js"></script> + +<h4>As observed by LIGO; the first detection of gravitational waves</h4> + + <p>This simulation shows how to evolve a binary black hole system for the last 6 orbits through to merger. </p> + + <div> + <h5>Physical parameters</h5> + <table> + <tr><th>Initial separation D</th><td>10 M</td></tr> + <tr><th>Mass ratio q = m1/m2</th><td>36/29 ~ 1.24</td></tr> + <tr><th>Spin chi1 = a1/m1</th><td>0.31</td></tr> + <tr><th>Spin chi2 = a2/m2</th><td>-0.46</td></tr> + </table> + + <h5>Physical properties</h5> + <table> + <tr><th>Number of orbits</th><td>6</td></tr> + <tr><th>Time to merger</th><td>900 M</td></tr> + <tr><th>Mass of final BH</th><td></td></tr> + <tr><th>Spin of final BH</th><td></td></tr> + </table> + + <h5>Computational details</h5> + <?php + $sim_stats = json_decode(file_get_contents("stats.json"), true); + ?></p> + <table> + <tr><th>Parameter file</th><td><a href="GW150914.rpar">GW150914.rpar</a></td></tr> + <tr><th>Thornlist</th><td><a href="GW150914.th">GW150914.th</a> (Llama multi-block code added to ET release version)</td></tr> + <tr><th>Submission command</th><td>simfactory/bin/sim create-submit GW150914_28 --define N 28 --parfile repos/GW150914/ParameterFiles/GW150914.rpar --procs 128 --walltime 24:00:00</td></tr> + <tr><th>Total memory</th><td><?php echo $sim_stats['TotalMemory']?> GB</td></tr> + <tr><th>Run time</th><td><?php echo $sim_stats['Runtime']?> days on <?php echo $sim_stats['Cores']?> cores (Intel(R) Xeon(R) CPU E5-2630 v3 @ 2.40GHz)</td></tr> + <tr><th>Cost</th><td><?php echo $sim_stats['CostInCoreHours']?> core hours</td></tr> + </table> + +<p> + + <h5>Horizon centroid coordinate trajectories</h5> + <a href="h22.asc"><img src="trajectories.png" alt="BH trajectories" width="242" height="237"></a> + + <h5>Gravitational waveform</h5> + <a href="h22.asc"><img src="waveform.png" alt="h_+^{2,2} at r=500M as a function of retarded time" width="542" height="191"></a> + + <h5>Tutorials</h5> + <ul> + <li><a href="https://bitbucket.org/barrywardell/gw150914sim/wiki/Simulation%20tutorial">Compile and run</a>: + Compile the code and run the simulation</li> + <li><a href="https://bitbucket.org/barrywardell/gw150914sim/wiki/Analysis%20tutorial">Analysis</a>: Analyse the results using SimulationTools for Mathematica</li> + <li><a href="https://bitbucket.org/barrywardell/gw150914sim/wiki/Visualising%20the%20spa…">Visualisation</a>: Visualise the data using VisIt</li> + </ul> + + + +<?php include_once($_SERVER['DOCUMENT_ROOT'].'/global/footer.php');?>

1 0

[svn:einsteintoolkit] www/about/gallery/gw150914/ (Rev. 1523)
by barry.wardell＠gmail.com 05 Jun '16

05 Jun '16

User: barry.wardell Date: 2016/06/05 05:04 PM Added: /about/gallery/gw150914/ stats.json Log: gw150914: Generate statistics File Changes: Directory: /about/gallery/gw150914/ =================================== File [added]: stats.json Delta lines: +6 -0 =================================================================== --- about/gallery/gw150914/stats.json (rev 0) +++ about/gallery/gw150914/stats.json 2016-06-05 22:04:33 UTC (rev 1523) @@ -0,0 +1,6 @@ +{ + "CostInCoreHours": 6995, + "Runtime": "2.3", + "Cores": 128, + "TotalMemory": 82 +} \ No newline at end of file

1 0

[svn:einsteintoolkit] www/about/gallery/ (Rev. 1522)
by barry.wardell＠gmail.com 05 Jun '16

05 Jun '16

User: barry.wardell Date: 2016/06/05 05:04 PM Modified: /about/gallery/ index.php Log: gw150914: Use statistics in page File Changes: Directory: /about/gallery/ ========================== File [modified]: index.php Delta lines: +8 -2 =================================================================== --- about/gallery/index.php 2016-06-05 22:04:19 UTC (rev 1521) +++ about/gallery/index.php 2016-06-05 22:04:27 UTC (rev 1522) @@ -111,14 +111,20 @@ </table> <h5>Computational details</h5> + <?php + $sim_stats = json_decode(file_get_contents("gw150914/stats.json"), true); + ?></p> <table> <tr><th>Parameter file</th><td><a href="GW150914/GW150914.rpar">GW150914.rpar</a></td></tr> <tr><th>Thornlist</th><td><a href="GW150914/GW150914.th">GW150914.th</a> (Llama multi-block code added to ET release version)</td></tr> <tr><th>Submission command</th><td>simfactory/bin/sim create-submit GW150914_28 --define N 28 --parfile repos/GW150914/ParameterFiles/GW150914.rpar --procs 128 --walltime 24:00:00</td></tr> - <tr><th>Approx. memory (N=28 resolution)</th><td>64 GB</td></tr> - <tr><th>Approx. runtime (N=28 resolution)</th><td>XX days on 128 cores (Intel(R) Xeon(R) CPU E5-2630 v3 @ 2.40GHz)</td></tr> + <tr><th>Total memory</th><td><?php echo $sim_stats['TotalMemory']?> GB</td></tr> + <tr><th>Run time</th><td><?php echo $sim_stats['Runtime']?> days on <?php echo $sim_stats['Cores']?> cores (Intel(R) Xeon(R) CPU E5-2630 v3 @ 2.40GHz)</td></tr> + <tr><th>Cost</th><td><?php echo $sim_stats['CostInCoreHours']?> core hours</td></tr> </table> +<p> + <h5>Horizon centroid coordinate trajectories</h5> <a href="gw150914/h22.asc"><img src="gw150914/trajectories.png" alt="BH trajectories" width="242" height="237"></a>

1 0

[svn:einsteintoolkit] www/about/gallery/gw150914/ (Rev. 1521)
by barry.wardell＠gmail.com 05 Jun '16

05 Jun '16

User: barry.wardell Date: 2016/06/05 05:04 PM Modified: /about/gallery/gw150914/ Makefile Log: gw150914: Make statistics from Makefile File Changes: Directory: /about/gallery/gw150914/ =================================== File [modified]: Makefile Delta lines: +7 -4 =================================================================== --- about/gallery/gw150914/Makefile 2016-06-05 22:04:14 UTC (rev 1520) +++ about/gallery/gw150914/Makefile 2016-06-05 22:04:19 UTC (rev 1521) @@ -1,5 +1,5 @@ -PLOTS = $(patsubst %.m,%.png,$(wildcard *.m)) +PLOTS = trajectories.png waveform.png -include make.conf @@ -7,10 +7,13 @@ PLOTFLAGS = --simulations-directory $(SIMULATIONS_DIRECTORY) endif -all: $(PLOTS) +all: $(PLOTS) stats.json -$(PLOTS): %.png : %.m +$(PLOTS): %.png: %.m MathematicaScript -script $< $(PLOTFLAGS) +stats.json: stats.m + MathematicaScript -script $< + clean: - rm -f $(PLOTS) + rm -f $(PLOTS) stats.json

1 0

[svn:einsteintoolkit] www/about/gallery/gw150914/ (Rev. 1520)
by barry.wardell＠gmail.com 05 Jun '16

05 Jun '16

User: barry.wardell Date: 2016/06/05 05:04 PM Added: /about/gallery/gw150914/ stats.m Log: gw150914: Add script for computing statistics of simulation File Changes: Directory: /about/gallery/gw150914/ =================================== File [added]: stats.m Delta lines: +26 -0 =================================================================== --- about/gallery/gw150914/stats.m (rev 0) +++ about/gallery/gw150914/stats.m 2016-06-05 22:04:14 UTC (rev 1520) @@ -0,0 +1,26 @@ + +<< SimulationTools`; + +simDirs = ReplaceList[$ScriptCommandLine, + {___, "--simulations-directory", val_, ___} :> val]; + +$SimulationPath = Join[simDirs, $SimulationPath]; + +outDir = Directory[]; + +sim = "GW150914_28"; + +walltimeHours = ReadWalltimeHours[sim]; +cores = 128; +walltimeDays = walltimeHours/24; +cost = cores walltimeHours; +processes = cores/8; +totalMemory = processes*Max[ReadSimulationAverageMemoryUsage[sim]]; + +stats = + {"CostInCoreHours" -> Round[cost], + "Runtime" -> ToString@NumberForm[walltimeDays,{Infinity,1}], + "Cores" -> cores, + "TotalMemory" -> Round[totalMemory/1024]}; + +Export["stats.json", stats];

1 0

[svn:einsteintoolkit] www/about/gallery/gw150914/ (Rev. 1519)
by barry.wardell＠gmail.com 05 Jun '16

05 Jun '16

User: barry.wardell Date: 2016/06/05 05:04 PM Modified: /about/gallery/gw150914/ trajectories.png, waveform.png Log: gw150914: Regenerate plots File Changes: Directory: /about/gallery/gw150914/ =================================== File [modified]: trajectories.png Delta lines: <Unavailable> File [modified]: waveform.png Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] www/about/gallery/gw150914/ (Rev. 1518)
by barry.wardell＠gmail.com 05 Jun '16

05 Jun '16

User: barry.wardell Date: 2016/06/05 05:04 PM Modified: /about/gallery/gw150914/ Makefile Log: gw150914: Fix logic in Makefile File Changes: Directory: /about/gallery/gw150914/ =================================== File [modified]: Makefile Delta lines: +5 -1 =================================================================== --- about/gallery/gw150914/Makefile 2016-06-05 22:03:59 UTC (rev 1517) +++ about/gallery/gw150914/Makefile 2016-06-05 22:04:03 UTC (rev 1518) @@ -3,10 +3,14 @@ -include make.conf +ifdef SIMULATIONS_DIRECTORY +PLOTFLAGS = --simulations-directory $(SIMULATIONS_DIRECTORY) +endif + all: $(PLOTS) $(PLOTS): %.png : %.m - MathematicaScript -script $< --simulations-directory $(SIMULATIONS_DIRECTORY) + MathematicaScript -script $< $(PLOTFLAGS) clean: rm -f $(PLOTS)

1 0

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits June 2016