Commits January 2011

commits@lists.einsteintoolkit.org

7 participants
57 discussions

[svn:einsteintoolkit] www/ (Rev. 460)
by schnetter＠cct.lsu.edu 31 Jan '11

31 Jan '11

User: eschnett Date: 2011/01/31 04:56 PM Modified: / index.php Log: Add release date to main page File Changes: Directory: / ============ File [modified]: index.php Delta lines: +8 -1 =================================================================== --- index.php 2011-01-31 03:42:04 UTC (rev 459) +++ index.php 2011-01-31 22:56:10 UTC (rev 460) @@ -9,7 +9,14 @@ interdisciplinary research and take advantage of emerging petascale computers and advanced cyberinfrastructure. Please read our pages <a href="about">about</a> the Einstein Toolkit, its <a href="about/governance">governance</a>, and how to <a href="info">get started</a> with the toolkit for more information. -New: We are pleased to announce the second release (code name "<a href="http://en.wikipedia.org/wiki/Subrahmanyan_Chandrasekhar">Chandrasekhar</a>") of the Einstein Toolkit, an open, community developed software infrastructure for relativistic astrophysics. This release is mainly a maintenance release incorporating fixes accumulated since the previous release in June 2010, as well as additional test suites. Read the full announcement <a href="http://einsteintoolkit.org/toolkit/releases/ET_2010_11_announcement.php">here</a>. +November 2010: We are pleased to +<a href="toolkit/releases/ET_2010_11_announcement.php">announce the +second release</a> (code name +"<a href="http://en.wikipedia.org/wiki/Subrahmanyan_Chandrasekhar">Chandrasekhar</a>") +of the Einstein Toolkit, an open, community developed software +infrastructure for relativistic astrophysics. This release is mainly a +maintenance release incorporating fixes accumulated since the previous +release in June 2010, as well as additional test suites. </td> <td valign="top">

1 0

[svn:einsteintoolkit] manifest/trunk/ (Rev. 72)
by knarf＠cct.lsu.edu 31 Jan '11

31 Jan '11

User: knarf Date: 2011/01/31 02:56 PM Modified: /trunk/ einsteintoolkit.th Log: add zlib to thornlist File Changes: Directory: /trunk/ ================== File [modified]: einsteintoolkit.th Delta lines: +1 -0 =================================================================== --- trunk/einsteintoolkit.th 2011-01-24 19:15:35 UTC (rev 71) +++ trunk/einsteintoolkit.th 2011-01-31 20:56:20 UTC (rev 72) @@ -220,6 +220,7 @@ ExternalLibraries/libjpeg ExternalLibraries/LORENE ExternalLibraries/OpenSSL +ExternalLibraries/zlib TAT/TATelliptic

1 0

[svn:einsteintoolkit] Hydro_InitExcision/trunk/test/ (Rev. 113)
by knarf＠cct.lsu.edu 31 Jan '11

31 Jan '11

User: knarf Date: 2011/01/31 09:35 AM Modified: /trunk/test/ diag_flip_pugh_eno.par, diag_flip_pugh_ppm.par, diag_flip_pugh_tvd.par, diag_pugh_eno.par, diag_pugh_ppm.par, diag_pugh_tvd.par, halfsphere_pugh_eno.par.disabled, halfsphere_pugh_ppm.par, halfsphere_pugh_tvd.par, sphere_pugh_eno.par.disabled, sphere_pugh_ppm.par, sphere_pugh_tvd.par, x_flip_pugh_eno.par, x_flip_pugh_ppm.par, x_flip_pugh_tvd.par, x_pugh_eno.par, x_pugh_ppm.par, x_pugh_tvd.par Log: let testsuites use EOS_Omni - which makes them pass again File Changes: Directory: /trunk/test/ ======================= File [modified]: diag_flip_pugh_eno.par Delta lines: +4 -2 =================================================================== --- trunk/test/diag_flip_pugh_eno.par 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/diag_flip_pugh_eno.par 2011-01-31 15:35:54 UTC (rev 113) @@ -54,9 +54,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase" File [modified]: diag_flip_pugh_ppm.par Delta lines: +4 -2 =================================================================== --- trunk/test/diag_flip_pugh_ppm.par 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/diag_flip_pugh_ppm.par 2011-01-31 15:35:54 UTC (rev 113) @@ -54,9 +54,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase" File [modified]: diag_flip_pugh_tvd.par Delta lines: +4 -2 =================================================================== --- trunk/test/diag_flip_pugh_tvd.par 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/diag_flip_pugh_tvd.par 2011-01-31 15:35:54 UTC (rev 113) @@ -54,9 +54,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase" File [modified]: diag_pugh_eno.par Delta lines: +4 -2 =================================================================== --- trunk/test/diag_pugh_eno.par 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/diag_pugh_eno.par 2011-01-31 15:35:54 UTC (rev 113) @@ -54,9 +54,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase" File [modified]: diag_pugh_ppm.par Delta lines: +4 -2 =================================================================== --- trunk/test/diag_pugh_ppm.par 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/diag_pugh_ppm.par 2011-01-31 15:35:54 UTC (rev 113) @@ -54,9 +54,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase" File [modified]: diag_pugh_tvd.par Delta lines: +4 -2 =================================================================== --- trunk/test/diag_pugh_tvd.par 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/diag_pugh_tvd.par 2011-01-31 15:35:54 UTC (rev 113) @@ -54,9 +54,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase" File [modified]: halfsphere_pugh_eno.par.disabled Delta lines: +4 -2 =================================================================== --- trunk/test/halfsphere_pugh_eno.par.disabled 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/halfsphere_pugh_eno.par.disabled 2011-01-31 15:35:54 UTC (rev 113) @@ -58,9 +58,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase" File [modified]: halfsphere_pugh_ppm.par Delta lines: +4 -2 =================================================================== --- trunk/test/halfsphere_pugh_ppm.par 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/halfsphere_pugh_ppm.par 2011-01-31 15:35:54 UTC (rev 113) @@ -58,9 +58,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase" File [modified]: halfsphere_pugh_tvd.par Delta lines: +4 -2 =================================================================== --- trunk/test/halfsphere_pugh_tvd.par 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/halfsphere_pugh_tvd.par 2011-01-31 15:35:54 UTC (rev 113) @@ -58,9 +58,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase" File [modified]: sphere_pugh_eno.par.disabled Delta lines: +4 -2 =================================================================== --- trunk/test/sphere_pugh_eno.par.disabled 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/sphere_pugh_eno.par.disabled 2011-01-31 15:35:54 UTC (rev 113) @@ -58,9 +58,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase" File [modified]: sphere_pugh_ppm.par Delta lines: +4 -2 =================================================================== --- trunk/test/sphere_pugh_ppm.par 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/sphere_pugh_ppm.par 2011-01-31 15:35:54 UTC (rev 113) @@ -58,9 +58,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase" File [modified]: sphere_pugh_tvd.par Delta lines: +4 -2 =================================================================== --- trunk/test/sphere_pugh_tvd.par 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/sphere_pugh_tvd.par 2011-01-31 15:35:54 UTC (rev 113) @@ -58,9 +58,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase" File [modified]: x_flip_pugh_eno.par Delta lines: +4 -2 =================================================================== --- trunk/test/x_flip_pugh_eno.par 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/x_flip_pugh_eno.par 2011-01-31 15:35:54 UTC (rev 113) @@ -53,9 +53,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase" File [modified]: x_flip_pugh_ppm.par Delta lines: +4 -2 =================================================================== --- trunk/test/x_flip_pugh_ppm.par 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/x_flip_pugh_ppm.par 2011-01-31 15:35:54 UTC (rev 113) @@ -53,9 +53,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase" File [modified]: x_flip_pugh_tvd.par Delta lines: +4 -2 =================================================================== --- trunk/test/x_flip_pugh_tvd.par 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/x_flip_pugh_tvd.par 2011-01-31 15:35:54 UTC (rev 113) @@ -53,9 +53,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase" File [modified]: x_pugh_eno.par Delta lines: +4 -2 =================================================================== --- trunk/test/x_pugh_eno.par 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/x_pugh_eno.par 2011-01-31 15:35:54 UTC (rev 113) @@ -53,9 +53,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase" File [modified]: x_pugh_ppm.par Delta lines: +4 -2 =================================================================== --- trunk/test/x_pugh_ppm.par 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/x_pugh_ppm.par 2011-01-31 15:35:54 UTC (rev 113) @@ -53,9 +53,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase" File [modified]: x_pugh_tvd.par Delta lines: +4 -2 =================================================================== --- trunk/test/x_pugh_tvd.par 2010-11-21 03:29:36 UTC (rev 112) +++ trunk/test/x_pugh_tvd.par 2011-01-31 15:35:54 UTC (rev 113) @@ -53,9 +53,11 @@ ADMBase::metric_type = "physical" SpaceMask::use_mask = "yes" -ActiveThorns = "eos_base eos_idealfluid eos_polytrope" +ActiveThorns = "eos_omni" -EOS_IdealFluid::eos_ideal_fluid_gamma = 1.66666666666666 +eos_omni::gl_gamma = 1.66666666666666 +eos_omni::gl_k = 80.0 +eos_omni::poly_k = 80.0 ActiveThorns = "tmunubase hydrobase grhydro initbase"

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/examples/bbh/ (Rev. 23)
by knarf＠cct.lsu.edu 31 Jan '11

31 Jan '11

User: knarf Date: 2011/01/31 08:56 AM Added: /examples/bbh/ bbh.par Log: first version of bbh parameter file File Changes: Directory: /examples/bbh/ ========================= File [added]: bbh.par Delta lines: +502 -0 =================================================================== --- examples/bbh/bbh.par (rev 0) +++ examples/bbh/bbh.par 2011-01-31 14:56:58 UTC (rev 23) @@ -0,0 +1,502 @@ +Cactus::cctk_run_title = "BBH" + +Cactus::cctk_full_warnings = yes +Cactus::highlight_warning_messages = no + +Cactus::terminate = "time" +Cactus::cctk_final_time = 300.0 + + + +ActiveThorns = "IOUtil" + +IO::out_dir = $parfile + + + +ActiveThorns = "AEILocalInterp" + +ActiveThorns = "BLAS LAPACK" + +ActiveThorns = "Fortran" + +ActiveThorns = "GSL" + +ActiveThorns = "GenericFD" + +ActiveThorns = "HDF5" + +ActiveThorns = "LocalInterp" + +ActiveThorns = "LoopControl" + +ActiveThorns = "Slab" + + + +ActiveThorns = "SummationByParts" + +SummationByParts::order = 4 + + + +ActiveThorns = "InitBase" + + + +ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce CarpetSlab" + +Carpet::verbose = no +Carpet::veryverbose = no +Carpet::schedule_barriers = no +Carpet::storage_verbose = no +#Carpet::timers_verbose = no +CarpetLib::output_bboxes = no + +Carpet::domain_from_coordbase = yes +Carpet::max_refinement_levels = 7 + +driver::ghost_size = 3 +Carpet::use_buffer_zones = yes + +Carpet::prolongation_order_space = 5 +Carpet::prolongation_order_time = 2 + +Carpet::convergence_level = 0 + +Carpet::init_fill_timelevels = yes +#Carpet::init_3_timelevels = yes + +Carpet::poison_new_timelevels = yes +CarpetLib::poison_new_memory = yes + +Carpet::output_timers_every = 640 +CarpetLib::print_timestats_every = 640 +CarpetLib::print_memstats_every = 640 + + + +ActiveThorns = "NaNChecker" + +NaNChecker::check_every = 1 # 64 +#NaNChecker::verbose = "all" +#NaNChecker::action_if_found = "just warn" +NaNChecker::action_if_found = "terminate" +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::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift +" + + + +ActiveThorns = "Boundary CartGrid3D CoordBase ReflectionSymmetry RotatingSymmetry180 SymBase" + +CoordBase::domainsize = "minmax" + +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 +CoordBase::dy = 2.00 +CoordBase::dz = 2.00 + +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_z_lower = 1 + +CartGrid3D::type = "coordbase" + +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_z = no + + + +ActiveThorns = "SphericalSurface" + +SphericalSurface::nsurfaces = 5 +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 + +SphericalSurface::ntheta [3] = 39 +SphericalSurface::nphi [3] = 76 +SphericalSurface::nghoststheta[3] = 2 +SphericalSurface::nghostsphi [3] = 2 + +SphericalSurface::ntheta [4] = 39 +SphericalSurface::nphi [4] = 76 +SphericalSurface::nghoststheta[4] = 2 +SphericalSurface::nghostsphi [4] = 2 + + + +ActiveThorns = "CarpetMask" + +CarpetMask::verbose = no + +CarpetMask::excluded_surface [0] = 3 +CarpetMask::excluded_surface_factor[0] = 1.0 + +CarpetMask::excluded_surface [1] = 4 +CarpetMask::excluded_surface_factor[1] = 1.0 + + + +ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal TmunuBase" +ActiveThorns = "CarpetRegrid2 PunctureTracker CarpetTracker" + +CarpetTracker::surface[0] = 0 +CarpetTracker::surface[1] = 1 +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 + + +CarpetRegrid2::regrid_every = 16 +CarpetRegrid2::freeze_unaligned_levels = yes +CarpetRegrid2::symmetry_rotating180 = yes +CarpetRegrid2::verbose = yes + +CarpetRegrid2::num_centres = 2 + +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 + + +ActiveThorns = "MoL Time" + +MoL::ODE_Method = "RK4" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 + +Carpet::time_refinement_factors = "[1, 1, 2, 4, 8, 16, 32, 64, 128, 256]" + +Time::dtfac = 0.25 + + +ADMMacros::spatial_order = 4 + + + +ActiveThorns = "TwoPunctures" + +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" + +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 + +#TODO# TwoPunctures::grid_setup_method = "evaluation" + +TwoPunctures::TP_epsilon = 1.0e-6 +TwoPunctures::TP_Tiny = 1.0e-2 + +TwoPunctures::verbose = yes + + + +ActiveThorns = "ML_BSSN ML_BSSN_Helper NewRad" + +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::harmonicN = 1 # 1+log +ML_BSSN::harmonicF = 2.0 # 1+log +ML_BSSN::ShiftGammaCoeff = 0.75 +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" + + + +ActiveThorns = "Dissipation" + +Dissipation::order = 5 +Dissipation::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 +" + + + +ActiveThorns = "ML_ADMConstraints" + + + +ActiveThorns = "AHFinderDirect" + +AHFinderDirect::find_every = 16 + +AHFinderDirect::run_at_CCTK_ANALYSIS = yes +AHFinderDirect::run_at_CCTK_POSTSTEP = false + +AHFinderDirect::move_origins = yes + +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 = 0 + +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] = 2 +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] = 3 +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] = 4 +AHFinderDirect::reset_horizon_after_not_finding [3] = no +AHFinderDirect::max_allowable_horizon_radius [3] = 6 + + +ActiveThorns = "WeylScal4 Multipole" + +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 = 8 + + +ActiveThorns = "CarpetIOBasic" + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_reductions = "norm2" +IOBasic::outInfo_vars = " + Carpet::physical_time_per_hour + ML_BSSN::ML_trace_curv +" + + + +ActiveThorns = "CarpetIOScalar" + +IOScalar::one_file_per_group = yes + +IOScalar::outScalar_every = 16 +IOScalar::outScalar_vars = " + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift + ML_ADMConstraints::ML_Ham + ML_ADMConstraints::ML_mom + SphericalSurface::sf_radius +" + + + +ActiveThorns = "CarpetIOASCII" + +IOASCII::one_file_per_group = yes + +IOASCII::output_symmetry_points = no +IOASCII::out3D_ghosts = no + +IOASCII::out0D_every = 16 +IOASCII::out0D_vars = " + Carpet::timing + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift + ML_ADMConstraints::ML_Ham + ML_ADMConstraints::ML_mom + SphericalSurface::sf_active + SphericalSurface::sf_valid + SphericalSurface::sf_info + SphericalSurface::sf_radius + SphericalSurface::sf_origin + SphericalSurface::sf_coordinate_descriptors + PunctureTracker::pt_loc +" + +IOASCII::out1D_every = 16 +IOASCII::out1D_vars = " + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift + ML_ADMConstraints::ML_Ham + ML_ADMConstraints::ML_mom + SphericalSurface::sf_radius +" + +IOASCII::out2D_every = 16 +IOASCII::out2D_vars = " + SphericalSurface::sf_radius +" + + + +Activethorns = "CarpetIOHDF5" + +IOHDF5::out_every = 64 +IOHDF5::one_file_per_group = yes +IOHDF5::output_symmetry_points = no +IOHDF5::out3D_ghosts = no +IOHDF5::compression_level = 1 +IOHDF5::use_checksums = yes +IOHDF5::out_vars = " + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ML_ADMConstraints::ML_Ham + ML_ADMConstraints::ML_mom + WeylScal4::Psi4r + WeylScal4::Psi4i +" + +IOHDF5::checkpoint = yes +IO::checkpoint_dir = $parfile +IO::checkpoint_ID = yes +IO::checkpoint_every_walltime_hours = 6.0 +IO::checkpoint_on_terminate = yes + +IO::recover = "autoprobe" +IO::recover_dir = $parfile + + + +ActiveThorns = "Formaline" + + + +ActiveThorns = "TimerReport" + +TimerReport::out_every = 640 +TimerReport::out_filename = "TimerReport" +TimerReport::output_all_timers_together = yes +TimerReport::output_all_timers_readable = yes +TimerReport::n_top_timers = 20

1 0

[svn:einsteintoolkit] www/about/members/ (Rev. 459)
by knarf＠cct.lsu.edu 30 Jan '11

30 Jan '11

User: knarf Date: 2011/01/30 09:42 PM Modified: /about/members/ index.php Log: correct filename File Changes: Directory: /about/members/ ========================== File [modified]: index.php Delta lines: +1 -1 =================================================================== --- about/members/index.php 2011-01-31 03:37:58 UTC (rev 458) +++ about/members/index.php 2011-01-31 03:42:04 UTC (rev 459) @@ -31,7 +31,7 @@ } } - $lines = file('list.txt'); + $lines = file('members.txt'); $groups = 0; $members = 0; $first = true;

1 0

[svn:einsteintoolkit] www/about/members/ (Rev. 458)
by knarf＠cct.lsu.edu 30 Jan '11

30 Jan '11

User: knarf Date: 2011/01/30 09:37 PM Modified: /about/members/ members.txt Log: add Roberto De Pietri File Changes: Directory: /about/members/ ========================== File [modified]: members.txt Delta lines: +3 -0 =================================================================== --- about/members/members.txt 2011-01-28 19:43:16 UTC (rev 457) +++ about/members/members.txt 2011-01-31 03:37:58 UTC (rev 458) @@ -58,6 +58,9 @@ Osaka University Luca Baiotti +Parma University + Roberto De Pietri + http://www.perimeterinstitute.ca/ Perimeter Institute David Rideout

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 22)
by tanja.bode＠physics.gatech.edu 30 Jan '11

30 Jan '11

User: tbode Date: 2011/01/30 04:54 PM Modified: / ET.tex Log: Analysis section draft. File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +32 -4 =================================================================== --- ET.tex 2011-01-24 19:24:58 UTC (rev 21) +++ ET.tex 2011-01-30 22:54:11 UTC (rev 22) @@ -508,14 +508,42 @@ \subsection{Analysis} -\todo{1/2 page Tanja} +%\todo{1/2 page Tanja} -\paragraph{Masses} +Any simulated system may have many applicable quantities for analysis which +are most efficiently calculated during evolution. The Einstein Toolkit +currently includes modules to facilitate the calculation of the most common +analysis quantities: masses, horizons, and gravitational waves. In addition +to these, detailed below, the \codename{ADMConstraints} module calculates the +constraints required, but not imposed by, the formulation for assessing the +simulation's accuracy and stability. Finally, \codename{ADMAnalysis} calculates +several quantities often needed in analyzing the spacetime evolution: the trace +of the extrinsic curvature $K$, the determinant of the 3-metric, the spatial +metric and extrinsic curvature in spherical coordinates, and the spatial Ricci +tensor and trace. -\paragraph{Event and Apparent Horizons} +\subsubsection*{Horizons} For spacetimes which contain a \bh{,} the Toolkit +contains modules to find both the event and apparent horizons and their +properties. The \codename{EHFinder} module traces null surfaces backwards in +time through a pre-evolved spacetime. The Toolkit also includes two separate +apparent horizon finders: Thornburg's \codename{AHFinderDirect} and +\codename{AHFinder}. -\paragraph{Gravitational Waves} +\subsubsection*{Masses} In analyzing spacetimes, it is important to calculate +the masses and momenta of the system. Several Toolkit modules include the +capability of calculating several types of masses. The apparent horizon modules +include in their output the areal masses of the found horizons. The +\codename{Extract} module calculates ADM mass and momenta when extracting +gravitational waves. Additionally, \codename{HydroAnalysis} calculates finds +the center of mass and the point of maximum matter density of a hydronamic +field. +\subsubsection*{Gravitational Waves} Within the Toolkit are modules for extracting +gravitational waves using two formulations: the Moncrief formalism assuming +the spacetime approximates spherical symmetry \codename{Extract}) and the Weyl +scalar $\Psi_4$ constructed on a fiducial tetrad (\codename{WeylScal4}). The +Weyl scalars, as well as any other quantity, can be projected onto a basis +of spin-weighted spherical harmonics using the \codename{Multipole} module. \subsection{Relativity Tools} \todo{1/2 page Peter}

1 0

[svn:einsteintoolkit] www/toolkit/releases/ (Rev. 457)
by knarf＠cct.lsu.edu 28 Jan '11

28 Jan '11

User: knarf Date: 2011/01/28 01:43 PM Modified: /toolkit/releases/ index.php Log: fix typo File Changes: Directory: /toolkit/releases/ ============================= File [modified]: index.php Delta lines: +1 -1 =================================================================== --- toolkit/releases/index.php 2011-01-24 16:20:33 UTC (rev 456) +++ toolkit/releases/index.php 2011-01-28 19:43:16 UTC (rev 457) @@ -2,7 +2,7 @@ include_once($_SERVER['DOCUMENT_ROOT'].'/global/header.php');?> The most recent release of the Einstein Toolkit - is <tt>ET_2010_1</tt>, code named + is <tt>ET_2010_11</tt>, code named "<a href="http://en.wikipedia.org/wiki/Subrahmanyan_Chandrasekhar">Chandrasekhar</a>". We expect this version to remain unchanged, save for correcting serious problems that may be detected. This version should be your

1 0

[svn:einsteintoolkit] EOS_Omni/trunk/src/nuc_eos/ (Rev. 35)
by schnetter＠cct.lsu.edu 26 Jan '11

26 Jan '11

User: eschnett Date: 2011/01/26 10:15 AM Added: /trunk/src/nuc_eos/ linterp.F Removed: /trunk/src/nuc_eos/ temp.f Log: Rename file, stage 2. Change suffix from .f to .F so that Cactus expands macros. File Changes: Directory: /trunk/src/nuc_eos/ ============================== File [added]: linterp.F Delta lines: +134 -0 =================================================================== --- trunk/src/nuc_eos/linterp.F (rev 0) +++ trunk/src/nuc_eos/linterp.F 2011-01-26 16:15:44 UTC (rev 35) @@ -0,0 +1,134 @@ +#include "cctk.h" + + SUBROUTINE intp3d ( x, y, z, f, kt, ft, nx, ny, nz, xt, yt, zt, + . d1, d2, d3 ) +c + implicit none +c +c--------------------------------------------------------------------- +c +c purpose: interpolation of a function of three variables in an +c equidistant(!!!) table. +c +c method: 8-point Lagrange linear interpolation formula +c +c x input vector of first variable +c y input vector of second variable +c z input vector of third variable +c +c f output vector of interpolated function values +c +c kt vector length of input and output vectors +c +c ft 3d array of tabulated function values +c nx x-dimension of table +c ny y-dimension of table +c nz z-dimension of table +c xt vector of x-coordinates of table +c yt vector of y-coordinates of table +c zt vector of z-coordinates of table +c +c d1 centered derivative of ft with respect to x +c d2 centered derivative of ft with respect to y +c d3 centered derivative of ft with respect to z +c Note that d? only make sense when intp3d is called with kt=1 +c--------------------------------------------------------------------- +c +c + +c + integer kt,nx,ny,nz,ktx + double precision x(kt),y(kt),z(kt),f(kt) + double precision xt(nx),yt(ny),zt(nz) + double precision ft(nx,ny,nz) + double precision d1,d2,d3 +c +c + PARAMETER (ktx = 400) + double precision fh(ktx,8), delx(ktx), dely(ktx), delz(ktx), + & a1(ktx), a2(ktx), a3(ktx), a4(ktx), + & a5(ktx), a6(ktx), a7(ktx), a8(ktx) + + double precision dx,dy,dz,dxi,dyi,dzi,dxyi,dxzi,dyzi,dxyzi + integer n,ix,iy,iz + + IF (kt .GT. ktx) call CCTK_WARN (0, '***KTX**') +c +c +c------ determine spacing parameters of (equidistant!!!) table +c + dx = (xt(nx) - xt(1)) / FLOAT(nx-1) + dy = (yt(ny) - yt(1)) / FLOAT(ny-1) + dz = (zt(nz) - zt(1)) / FLOAT(nz-1) +c + dxi = 1. / dx + dyi = 1. / dy + dzi = 1. / dz +c + dxyi = dxi * dyi + dxzi = dxi * dzi + dyzi = dyi * dzi +c + dxyzi = dxi * dyi * dzi +c +c +c------- loop over all points to be interpolated +c + DO n = 1, kt +c +c------- determine location in (equidistant!!!) table +c + ix = 2 + INT( (x(n) - xt(1) - 1.e-10) * dxi ) + iy = 2 + INT( (y(n) - yt(1) - 1.e-10) * dyi ) + iz = 2 + INT( (z(n) - zt(1) - 1.e-10) * dzi ) +c + ix = MAX( 2, MIN( ix, nx ) ) + iy = MAX( 2, MIN( iy, ny ) ) + iz = MAX( 2, MIN( iz, nz ) ) +c +c write(*,*) iy-1,iy,iy+1 +c +c------- set-up auxiliary arrays for Lagrange interpolation +c + delx(n) = xt(ix) - x(n) + dely(n) = yt(iy) - y(n) + delz(n) = zt(iz) - z(n) +c + fh(n,1) = ft(ix , iy , iz ) + fh(n,2) = ft(ix-1, iy , iz ) + fh(n,3) = ft(ix , iy-1, iz ) + fh(n,4) = ft(ix , iy , iz-1) + fh(n,5) = ft(ix-1, iy-1, iz ) + fh(n,6) = ft(ix-1, iy , iz-1) + fh(n,7) = ft(ix , iy-1, iz-1) + fh(n,8) = ft(ix-1, iy-1, iz-1) +c +c------ set up coefficients of the interpolation polynomial and +c evaluate function values +c + a1(n) = fh(n,1) + a2(n) = dxi * ( fh(n,2) - fh(n,1) ) + a3(n) = dyi * ( fh(n,3) - fh(n,1) ) + a4(n) = dzi * ( fh(n,4) - fh(n,1) ) + a5(n) = dxyi * ( fh(n,5) - fh(n,2) - fh(n,3) + fh(n,1) ) + a6(n) = dxzi * ( fh(n,6) - fh(n,2) - fh(n,4) + fh(n,1) ) + a7(n) = dyzi * ( fh(n,7) - fh(n,3) - fh(n,4) + fh(n,1) ) + a8(n) = dxyzi * ( fh(n,8) - fh(n,1) + fh(n,2) + fh(n,3) + + & fh(n,4) - fh(n,5) - fh(n,6) - fh(n,7) ) +c + d1 = -a2(n) + d2 = -a3(n) + d3 = -a4(n) + f(n) = a1(n) + a2(n) * delx(n) + & + a3(n) * dely(n) + & + a4(n) * delz(n) + & + a5(n) * delx(n) * dely(n) + & + a6(n) * delx(n) * delz(n) + & + a7(n) * dely(n) * delz(n) + & + a8(n) * delx(n) * dely(n) * delz(n) +c + ENDDO +c + RETURN + END + File [removed]: temp.f Delta lines: +0 -134 =================================================================== --- trunk/src/nuc_eos/temp.f 2011-01-26 16:15:13 UTC (rev 34) +++ trunk/src/nuc_eos/temp.f 2011-01-26 16:15:44 UTC (rev 35) @@ -1,134 +0,0 @@ -#include "cctk.h" - - SUBROUTINE intp3d ( x, y, z, f, kt, ft, nx, ny, nz, xt, yt, zt, - . d1, d2, d3 ) -c - implicit none -c -c--------------------------------------------------------------------- -c -c purpose: interpolation of a function of three variables in an -c equidistant(!!!) table. -c -c method: 8-point Lagrange linear interpolation formula -c -c x input vector of first variable -c y input vector of second variable -c z input vector of third variable -c -c f output vector of interpolated function values -c -c kt vector length of input and output vectors -c -c ft 3d array of tabulated function values -c nx x-dimension of table -c ny y-dimension of table -c nz z-dimension of table -c xt vector of x-coordinates of table -c yt vector of y-coordinates of table -c zt vector of z-coordinates of table -c -c d1 centered derivative of ft with respect to x -c d2 centered derivative of ft with respect to y -c d3 centered derivative of ft with respect to z -c Note that d? only make sense when intp3d is called with kt=1 -c--------------------------------------------------------------------- -c -c - -c - integer kt,nx,ny,nz,ktx - double precision x(kt),y(kt),z(kt),f(kt) - double precision xt(nx),yt(ny),zt(nz) - double precision ft(nx,ny,nz) - double precision d1,d2,d3 -c -c - PARAMETER (ktx = 400) - double precision fh(ktx,8), delx(ktx), dely(ktx), delz(ktx), - & a1(ktx), a2(ktx), a3(ktx), a4(ktx), - & a5(ktx), a6(ktx), a7(ktx), a8(ktx) - - double precision dx,dy,dz,dxi,dyi,dzi,dxyi,dxzi,dyzi,dxyzi - integer n,ix,iy,iz - - IF (kt .GT. ktx) call CCTK_WARN (0, '***KTX**') -c -c -c------ determine spacing parameters of (equidistant!!!) table -c - dx = (xt(nx) - xt(1)) / FLOAT(nx-1) - dy = (yt(ny) - yt(1)) / FLOAT(ny-1) - dz = (zt(nz) - zt(1)) / FLOAT(nz-1) -c - dxi = 1. / dx - dyi = 1. / dy - dzi = 1. / dz -c - dxyi = dxi * dyi - dxzi = dxi * dzi - dyzi = dyi * dzi -c - dxyzi = dxi * dyi * dzi -c -c -c------- loop over all points to be interpolated -c - DO n = 1, kt -c -c------- determine location in (equidistant!!!) table -c - ix = 2 + INT( (x(n) - xt(1) - 1.e-10) * dxi ) - iy = 2 + INT( (y(n) - yt(1) - 1.e-10) * dyi ) - iz = 2 + INT( (z(n) - zt(1) - 1.e-10) * dzi ) -c - ix = MAX( 2, MIN( ix, nx ) ) - iy = MAX( 2, MIN( iy, ny ) ) - iz = MAX( 2, MIN( iz, nz ) ) -c -c write(*,*) iy-1,iy,iy+1 -c -c------- set-up auxiliary arrays for Lagrange interpolation -c - delx(n) = xt(ix) - x(n) - dely(n) = yt(iy) - y(n) - delz(n) = zt(iz) - z(n) -c - fh(n,1) = ft(ix , iy , iz ) - fh(n,2) = ft(ix-1, iy , iz ) - fh(n,3) = ft(ix , iy-1, iz ) - fh(n,4) = ft(ix , iy , iz-1) - fh(n,5) = ft(ix-1, iy-1, iz ) - fh(n,6) = ft(ix-1, iy , iz-1) - fh(n,7) = ft(ix , iy-1, iz-1) - fh(n,8) = ft(ix-1, iy-1, iz-1) -c -c------ set up coefficients of the interpolation polynomial and -c evaluate function values -c - a1(n) = fh(n,1) - a2(n) = dxi * ( fh(n,2) - fh(n,1) ) - a3(n) = dyi * ( fh(n,3) - fh(n,1) ) - a4(n) = dzi * ( fh(n,4) - fh(n,1) ) - a5(n) = dxyi * ( fh(n,5) - fh(n,2) - fh(n,3) + fh(n,1) ) - a6(n) = dxzi * ( fh(n,6) - fh(n,2) - fh(n,4) + fh(n,1) ) - a7(n) = dyzi * ( fh(n,7) - fh(n,3) - fh(n,4) + fh(n,1) ) - a8(n) = dxyzi * ( fh(n,8) - fh(n,1) + fh(n,2) + fh(n,3) + - & fh(n,4) - fh(n,5) - fh(n,6) - fh(n,7) ) -c - d1 = -a2(n) - d2 = -a3(n) - d3 = -a4(n) - f(n) = a1(n) + a2(n) * delx(n) - & + a3(n) * dely(n) - & + a4(n) * delz(n) - & + a5(n) * delx(n) * dely(n) - & + a6(n) * delx(n) * delz(n) - & + a7(n) * dely(n) * delz(n) - & + a8(n) * delx(n) * dely(n) * delz(n) -c - ENDDO -c - RETURN - END -

1 0

[svn:einsteintoolkit] EOS_Omni/trunk/src/nuc_eos/ (Rev. 34)
by schnetter＠cct.lsu.edu 26 Jan '11

26 Jan '11

User: eschnett Date: 2011/01/26 10:15 AM Added: /trunk/src/nuc_eos/ temp.f Removed: /trunk/src/nuc_eos/ linterp.f Log: Rename file, stage 1 File Changes: Directory: /trunk/src/nuc_eos/ ============================== File [removed]: linterp.f Delta lines: +0 -134 =================================================================== --- trunk/src/nuc_eos/linterp.f 2011-01-26 16:14:37 UTC (rev 33) +++ trunk/src/nuc_eos/linterp.f 2011-01-26 16:15:13 UTC (rev 34) @@ -1,134 +0,0 @@ -#include "cctk.h" - - SUBROUTINE intp3d ( x, y, z, f, kt, ft, nx, ny, nz, xt, yt, zt, - . d1, d2, d3 ) -c - implicit none -c -c--------------------------------------------------------------------- -c -c purpose: interpolation of a function of three variables in an -c equidistant(!!!) table. -c -c method: 8-point Lagrange linear interpolation formula -c -c x input vector of first variable -c y input vector of second variable -c z input vector of third variable -c -c f output vector of interpolated function values -c -c kt vector length of input and output vectors -c -c ft 3d array of tabulated function values -c nx x-dimension of table -c ny y-dimension of table -c nz z-dimension of table -c xt vector of x-coordinates of table -c yt vector of y-coordinates of table -c zt vector of z-coordinates of table -c -c d1 centered derivative of ft with respect to x -c d2 centered derivative of ft with respect to y -c d3 centered derivative of ft with respect to z -c Note that d? only make sense when intp3d is called with kt=1 -c--------------------------------------------------------------------- -c -c - -c - integer kt,nx,ny,nz,ktx - double precision x(kt),y(kt),z(kt),f(kt) - double precision xt(nx),yt(ny),zt(nz) - double precision ft(nx,ny,nz) - double precision d1,d2,d3 -c -c - PARAMETER (ktx = 400) - double precision fh(ktx,8), delx(ktx), dely(ktx), delz(ktx), - & a1(ktx), a2(ktx), a3(ktx), a4(ktx), - & a5(ktx), a6(ktx), a7(ktx), a8(ktx) - - double precision dx,dy,dz,dxi,dyi,dzi,dxyi,dxzi,dyzi,dxyzi - integer n,ix,iy,iz - - IF (kt .GT. ktx) call CCTK_WARN (0, '***KTX**') -c -c -c------ determine spacing parameters of (equidistant!!!) table -c - dx = (xt(nx) - xt(1)) / FLOAT(nx-1) - dy = (yt(ny) - yt(1)) / FLOAT(ny-1) - dz = (zt(nz) - zt(1)) / FLOAT(nz-1) -c - dxi = 1. / dx - dyi = 1. / dy - dzi = 1. / dz -c - dxyi = dxi * dyi - dxzi = dxi * dzi - dyzi = dyi * dzi -c - dxyzi = dxi * dyi * dzi -c -c -c------- loop over all points to be interpolated -c - DO n = 1, kt -c -c------- determine location in (equidistant!!!) table -c - ix = 2 + INT( (x(n) - xt(1) - 1.e-10) * dxi ) - iy = 2 + INT( (y(n) - yt(1) - 1.e-10) * dyi ) - iz = 2 + INT( (z(n) - zt(1) - 1.e-10) * dzi ) -c - ix = MAX( 2, MIN( ix, nx ) ) - iy = MAX( 2, MIN( iy, ny ) ) - iz = MAX( 2, MIN( iz, nz ) ) -c -c write(*,*) iy-1,iy,iy+1 -c -c------- set-up auxiliary arrays for Lagrange interpolation -c - delx(n) = xt(ix) - x(n) - dely(n) = yt(iy) - y(n) - delz(n) = zt(iz) - z(n) -c - fh(n,1) = ft(ix , iy , iz ) - fh(n,2) = ft(ix-1, iy , iz ) - fh(n,3) = ft(ix , iy-1, iz ) - fh(n,4) = ft(ix , iy , iz-1) - fh(n,5) = ft(ix-1, iy-1, iz ) - fh(n,6) = ft(ix-1, iy , iz-1) - fh(n,7) = ft(ix , iy-1, iz-1) - fh(n,8) = ft(ix-1, iy-1, iz-1) -c -c------ set up coefficients of the interpolation polynomial and -c evaluate function values -c - a1(n) = fh(n,1) - a2(n) = dxi * ( fh(n,2) - fh(n,1) ) - a3(n) = dyi * ( fh(n,3) - fh(n,1) ) - a4(n) = dzi * ( fh(n,4) - fh(n,1) ) - a5(n) = dxyi * ( fh(n,5) - fh(n,2) - fh(n,3) + fh(n,1) ) - a6(n) = dxzi * ( fh(n,6) - fh(n,2) - fh(n,4) + fh(n,1) ) - a7(n) = dyzi * ( fh(n,7) - fh(n,3) - fh(n,4) + fh(n,1) ) - a8(n) = dxyzi * ( fh(n,8) - fh(n,1) + fh(n,2) + fh(n,3) + - & fh(n,4) - fh(n,5) - fh(n,6) - fh(n,7) ) -c - d1 = -a2(n) - d2 = -a3(n) - d3 = -a4(n) - f(n) = a1(n) + a2(n) * delx(n) - & + a3(n) * dely(n) - & + a4(n) * delz(n) - & + a5(n) * delx(n) * dely(n) - & + a6(n) * delx(n) * delz(n) - & + a7(n) * dely(n) * delz(n) - & + a8(n) * delx(n) * dely(n) * delz(n) -c - ENDDO -c - RETURN - END - File [added]: temp.f Delta lines: +134 -0 =================================================================== --- trunk/src/nuc_eos/temp.f (rev 0) +++ trunk/src/nuc_eos/temp.f 2011-01-26 16:15:13 UTC (rev 34) @@ -0,0 +1,134 @@ +#include "cctk.h" + + SUBROUTINE intp3d ( x, y, z, f, kt, ft, nx, ny, nz, xt, yt, zt, + . d1, d2, d3 ) +c + implicit none +c +c--------------------------------------------------------------------- +c +c purpose: interpolation of a function of three variables in an +c equidistant(!!!) table. +c +c method: 8-point Lagrange linear interpolation formula +c +c x input vector of first variable +c y input vector of second variable +c z input vector of third variable +c +c f output vector of interpolated function values +c +c kt vector length of input and output vectors +c +c ft 3d array of tabulated function values +c nx x-dimension of table +c ny y-dimension of table +c nz z-dimension of table +c xt vector of x-coordinates of table +c yt vector of y-coordinates of table +c zt vector of z-coordinates of table +c +c d1 centered derivative of ft with respect to x +c d2 centered derivative of ft with respect to y +c d3 centered derivative of ft with respect to z +c Note that d? only make sense when intp3d is called with kt=1 +c--------------------------------------------------------------------- +c +c + +c + integer kt,nx,ny,nz,ktx + double precision x(kt),y(kt),z(kt),f(kt) + double precision xt(nx),yt(ny),zt(nz) + double precision ft(nx,ny,nz) + double precision d1,d2,d3 +c +c + PARAMETER (ktx = 400) + double precision fh(ktx,8), delx(ktx), dely(ktx), delz(ktx), + & a1(ktx), a2(ktx), a3(ktx), a4(ktx), + & a5(ktx), a6(ktx), a7(ktx), a8(ktx) + + double precision dx,dy,dz,dxi,dyi,dzi,dxyi,dxzi,dyzi,dxyzi + integer n,ix,iy,iz + + IF (kt .GT. ktx) call CCTK_WARN (0, '***KTX**') +c +c +c------ determine spacing parameters of (equidistant!!!) table +c + dx = (xt(nx) - xt(1)) / FLOAT(nx-1) + dy = (yt(ny) - yt(1)) / FLOAT(ny-1) + dz = (zt(nz) - zt(1)) / FLOAT(nz-1) +c + dxi = 1. / dx + dyi = 1. / dy + dzi = 1. / dz +c + dxyi = dxi * dyi + dxzi = dxi * dzi + dyzi = dyi * dzi +c + dxyzi = dxi * dyi * dzi +c +c +c------- loop over all points to be interpolated +c + DO n = 1, kt +c +c------- determine location in (equidistant!!!) table +c + ix = 2 + INT( (x(n) - xt(1) - 1.e-10) * dxi ) + iy = 2 + INT( (y(n) - yt(1) - 1.e-10) * dyi ) + iz = 2 + INT( (z(n) - zt(1) - 1.e-10) * dzi ) +c + ix = MAX( 2, MIN( ix, nx ) ) + iy = MAX( 2, MIN( iy, ny ) ) + iz = MAX( 2, MIN( iz, nz ) ) +c +c write(*,*) iy-1,iy,iy+1 +c +c------- set-up auxiliary arrays for Lagrange interpolation +c + delx(n) = xt(ix) - x(n) + dely(n) = yt(iy) - y(n) + delz(n) = zt(iz) - z(n) +c + fh(n,1) = ft(ix , iy , iz ) + fh(n,2) = ft(ix-1, iy , iz ) + fh(n,3) = ft(ix , iy-1, iz ) + fh(n,4) = ft(ix , iy , iz-1) + fh(n,5) = ft(ix-1, iy-1, iz ) + fh(n,6) = ft(ix-1, iy , iz-1) + fh(n,7) = ft(ix , iy-1, iz-1) + fh(n,8) = ft(ix-1, iy-1, iz-1) +c +c------ set up coefficients of the interpolation polynomial and +c evaluate function values +c + a1(n) = fh(n,1) + a2(n) = dxi * ( fh(n,2) - fh(n,1) ) + a3(n) = dyi * ( fh(n,3) - fh(n,1) ) + a4(n) = dzi * ( fh(n,4) - fh(n,1) ) + a5(n) = dxyi * ( fh(n,5) - fh(n,2) - fh(n,3) + fh(n,1) ) + a6(n) = dxzi * ( fh(n,6) - fh(n,2) - fh(n,4) + fh(n,1) ) + a7(n) = dyzi * ( fh(n,7) - fh(n,3) - fh(n,4) + fh(n,1) ) + a8(n) = dxyzi * ( fh(n,8) - fh(n,1) + fh(n,2) + fh(n,3) + + & fh(n,4) - fh(n,5) - fh(n,6) - fh(n,7) ) +c + d1 = -a2(n) + d2 = -a3(n) + d3 = -a4(n) + f(n) = a1(n) + a2(n) * delx(n) + & + a3(n) * dely(n) + & + a4(n) * delz(n) + & + a5(n) * delx(n) * dely(n) + & + a6(n) * delx(n) * delz(n) + & + a7(n) * dely(n) * delz(n) + & + a8(n) * delx(n) * dely(n) * delz(n) +c + ENDDO +c + RETURN + END +

1 0

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits January 2011