Commits

commits@lists.einsteintoolkit.org

2 participants
6997 discussions

[svn:einsteintoolkit] incoming/PITTNullCode/ (Rev. 45)
by yosef＠astro.rit.edu 30 Jun '11

30 Jun '11

User: zlochower Date: 2011/06/30 10:43 AM Added: /PITTNullCode/SphericalHarmonicDecomp/par/ qc0-mclachlan-CCE_Cauchy.par /PITTNullCode/SphericalHarmonicRecon/par/ qc0-mclachlan-CCE_Null.par Log: paramter files (based on the qc0 parfile in simfactory) for CCE. qc0-mclachlan-CCE_Cauchy.par uses McLachlan to produce the worldtube data. qc0-mclachlan-CCE_Null.par then evolves these data to produce the news. File Changes: Directory: /PITTNullCode/SphericalHarmonicDecomp/par/ ===================================================== File [added]: qc0-mclachlan-CCE_Cauchy.par Delta lines: +547 -0 =================================================================== --- PITTNullCode/SphericalHarmonicDecomp/par/qc0-mclachlan-CCE_Cauchy.par (rev 0) +++ PITTNullCode/SphericalHarmonicDecomp/par/qc0-mclachlan-CCE_Cauchy.par 2011-06-30 15:43:02 UTC (rev 45) @@ -0,0 +1,547 @@ +Cactus::cctk_run_title = "QC-0-CCE" + +Cactus::cctk_full_warnings = yes +Cactus::highlight_warning_messages = no + +Cactus::terminate = "time" +Cactus::cctk_final_time = 100.0 + + + +ActiveThorns = "IOUtil" + +IO::out_dir = $parfile + + + +ActiveThorns = "TerminationTrigger" + +TerminationTrigger::max_walltime = 0.0 # hours +TerminationTrigger::on_remaining_walltime = 15.0 # minutes + +TerminationTrigger::termination_from_file = yes +TerminationTrigger::create_termination_file = yes +TerminationTrigger::termination_file = "../TERMINATE" + + + +ActiveThorns = "AEILocalInterp" + +ActiveThorns = "BLAS LAPACK" + +ActiveThorns = "Fortran" + +ActiveThorns = "GSL" + +ActiveThorns = "GenericFD" + +ActiveThorns = "HDF5" + +ActiveThorns = "LocalInterp" + +ActiveThorns = "LoopControl" + +ActiveThorns = "Slab" + +ActiveThorns = "TGRtensor" + + + +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 = 10 + +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_3_timelevels = yes + +Carpet::poison_new_timelevels = yes +CarpetLib::poison_new_memory = yes + +Carpet::output_timers_every = 512 +CarpetLib::print_timestats_every = 512 +CarpetLib::print_memstats_every = 512 + + + +ActiveThorns = "NaNChecker" + +NaNChecker::check_every = 1 # 512 +NaNChecker::action_if_found = "terminate" +NaNChecker::check_vars = " + ML_BSSN::ML_Gamma + ML_BSSN::ML_lapse + ML_BSSN::ML_shift + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_trace_curv + ML_BSSN::ML_curv + 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 = 4 +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 + + + +ActiveThorns = "CarpetMask" + +CarpetMask::verbose = no + +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 + + + +ActiveThorns = "CarpetRegrid2 CarpetTracker" + +CarpetTracker::surface[0] = 0 +CarpetTracker::surface[1] = 1 +CarpetTracker::surface[2] = 2 + +CarpetRegrid2::regrid_every = 128 +CarpetRegrid2::freeze_unaligned_levels = yes +CarpetRegrid2::symmetry_rotating180 = yes +CarpetRegrid2::verbose = yes + +CarpetRegrid2::num_centres = 3 + +CarpetRegrid2::num_levels_1 = 7 +CarpetRegrid2::position_x_1 = +1.168642873 +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 = -1.168642873 +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 + +CarpetRegrid2::num_levels_3 = 7 +CarpetRegrid2::active_3 = no +CarpetRegrid2::radius_3[ 1] = 64.0 +CarpetRegrid2::radius_3[ 2] = 16.0 +CarpetRegrid2::radius_3[ 3] = 8.0 +CarpetRegrid2::radius_3[ 4] = 4.0 +CarpetRegrid2::radius_3[ 5] = 2.0 +CarpetRegrid2::radius_3[ 6] = 1.0 +CarpetRegrid2::movement_threshold_3 = 0.16 + + + +ActiveThorns = "MoL Time" + +MoL::ODE_Method = "RK4" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 + +Time::dtfac = 0.40 + + + +ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal TmunuBase" + +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 = 1.168642873 +TwoPunctures::par_m_plus = 0.453 +TwoPunctures::par_m_minus = 0.453 +TwoPunctures::par_P_plus [1] = +0.3331917498 +TwoPunctures::par_P_minus[1] = -0.3331917498 + +#TODO# TwoPunctures::grid_setup_method = "evaluation" + +TwoPunctures::TP_epsilon = 1.0e-2 +TwoPunctures::TP_Tiny = 1.0e-2 + +TwoPunctures::verbose = yes + + + +ActiveThorns = "ML_BSSN ML_BSSN_Helper" + +ADMBase::evolution_method = "ML_BSSN" +ADMBase::lapse_evolution_method = "ML_BSSN" +ADMBase::shift_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::BetaDriver = 1.0 +ML_BSSN::LapseAdvectionCoeff = 1 +ML_BSSN::ShiftAdvectionCoeff = 1 + +ML_BSSN::my_boundary_condition = "Minkowski" + +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_Gamma + ML_BSSN::ML_lapse + ML_BSSN::ML_shift + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_trace_curv + ML_BSSN::ML_curv +" + + + +ActiveThorns = "ML_ADMConstraints" + + + +ActiveThorns = "AHFinderDirect" + +AHFinderDirect::find_every = 128 + +AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = no + +AHFinderDirect::move_origins = yes +AHFinderDirect::reshape_while_moving = yes +AHFinderDirect::predict_origin_movement = yes + +AHFinderDirect::geometry_interpolator_name = "Lagrange polynomial interpolation" +AHFinderDirect::geometry_interpolator_pars = " + order=4 + boundary_off_centering_tolerance={1.0e-10 1.0e-10 1.0e-10 1.0e-10 1.0e-10 1.0e-10} + boundary_extrapolation_tolerance={0.0 0.0 0.0 0.0 0.0 0.0} +" +AHFinderDirect::surface_interpolator_name = "Lagrange polynomial interpolation" +AHFinderDirect::surface_interpolator_pars = " + order=4 + boundary_off_centering_tolerance={1.0e-10 1.0e-10 1.0e-10 1.0e-10} + boundary_extrapolation_tolerance={0.0 0.0 0.0 0.0} +" + +AHFinderDirect::output_h_every = 0 + +AHFinderDirect::N_horizons = 4 + +AHFinderDirect::origin_x [1] = +1.168642873 +AHFinderDirect::initial_guess__coord_sphere__x_center[1] = +1.168642873 +AHFinderDirect::initial_guess__coord_sphere__radius [1] = 0.25 +AHFinderDirect::which_surface_to_store_info [1] = 0 +AHFinderDirect::reset_horizon_after_not_finding [1] = no +AHFinderDirect::dont_find_after_individual_time [1] = 30.0 + +AHFinderDirect::origin_x [2] = -1.168642873 +AHFinderDirect::initial_guess__coord_sphere__x_center[2] = -1.168642873 +AHFinderDirect::initial_guess__coord_sphere__radius [2] = 0.25 +AHFinderDirect::which_surface_to_store_info [2] = 1 +AHFinderDirect::reset_horizon_after_not_finding [2] = no +AHFinderDirect::dont_find_after_individual_time [2] = 30.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::find_after_individual_time [3] = 15.0 + +AHFinderDirect::surface_definition [4] = "expansion product" +AHFinderDirect::surface_selection [4] = "areal radius" +AHFinderDirect::desired_value [4] = 50.0 +AHFinderDirect::initial_guess__coord_sphere__radius [4] = 50.0 +AHFinderDirect::which_surface_to_store_info [4] = 3 +AHFinderDirect::reset_horizon_after_not_finding [4] = no + + + +ActiveThorns = "QuasiLocalMeasures" + +QuasiLocalMeasures::verbose = yes +QuasiLocalMeasures::interpolator = "Lagrange polynomial interpolation" +QuasiLocalMeasures::interpolator_options = " + order=4 + boundary_off_centering_tolerance={1.0e-10 1.0e-10 1.0e-10 1.0e-10 1.0e-10 1.0e-10} + boundary_extrapolation_tolerance={0.0 0.0 0.0 0.0 0.0 0.0} +" +QuasiLocalMeasures::spatial_order = 4 + +QuasiLocalMeasures::num_surfaces = 4 +QuasiLocalMeasures::surface_index[0] = 0 +QuasiLocalMeasures::surface_index[1] = 1 +QuasiLocalMeasures::surface_index[2] = 2 +QuasiLocalMeasures::surface_index[3] = 3 + + + +ActiveThorns = "Multipole WeylScal4" + +#WeylScal4::Psi4r_group_bound = "scalar" +#WeylScal4::Psi4i_group_bound = "scalar" + +Multipole::l_max = 8 +Multipole::nradii = 4 +Multipole::radius[0] = 100.0 +Multipole::radius[1] = 80.0 +Multipole::radius[2] = 60.0 +Multipole::radius[3] = 40.0 +Multipole::out_1d_every = 128 +Multipole::variables = " + WeylScal4::psi4r{sw=-2 cmplx='WeylScal4::psi4i' name='WeylScal4::psi4'} +" + +ActiveThorns="SphericalHarmonicDecomp" + +SphericalHarmonicDecomp::extract_spacetime_metric_every=32 +SphericalHarmonicDecomp::num_radii=3 +SphericalHarmonicDecomp::EM_Rin[0]=18 +SphericalHarmonicDecomp::EM_Rout[0]=22 +SphericalHarmonicDecomp::EM_Rin[1]=47 +SphericalHarmonicDecomp::EM_Rout[1]=53 +SphericalHarmonicDecomp::EM_Rin[2]=94 +SphericalHarmonicDecomp::EM_Rout[2]=106 +SphericalHarmonicDecomp::num_l_modes=7 +SphericalHarmonicDecomp::num_n_modes=7 +SphericalHarmonicDecomp::num_mu_points=41 +SphericalHarmonicDecomp::num_phi_points=82 +SphericalHarmonicDecomp::num_x_points=28 + + + +ActiveThorns = "CarpetIOBasic" + +IOBasic::outInfo_every = 128 +IOBasic::outInfo_reductions = "norm2" +IOBasic::outInfo_vars = " + Carpet::physical_time_per_hour + Carpet::local_grid_points_per_second + ML_ADMConstraints::H + WeylScal4::psi4r +" + + + +ActiveThorns = "CarpetIOScalar" + +IOScalar::one_file_per_group = yes + +IOScalar::outScalar_every = 128 +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 + WeylScal4::psi4r + WeylScal4::psi4i + QuasiLocalMeasures::qlm_newman_penrose + QuasiLocalMeasures::qlm_weyl_scalars + QuasiLocalMeasures::qlm_ricci_scalars +" + + + +ActiveThorns = "CarpetIOASCII" + +IOASCII::one_file_per_group = yes + +#IOASCII::output_symmetry_points = no +#IOASCII::output_ghost_points = no + +IOASCII::out0D_every = 128 +IOASCII::out0D_vars = " + Carpet::timing + QuasiLocalMeasures::qlm_state + QuasiLocalMeasures::qlm_grid_int + QuasiLocalMeasures::qlm_grid_real + QuasiLocalMeasures::qlm_scalars + QuasiLocalMeasures::qlm_multipole_moments +" + +IOASCII::out1D_every = 128 +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 + WeylScal4::psi4r + WeylScal4::psi4i + Multipole::harmonics + QuasiLocalMeasures::qlm_newman_penrose + QuasiLocalMeasures::qlm_weyl_scalars + QuasiLocalMeasures::qlm_ricci_scalars + QuasiLocalMeasures::qlm_killing_vector + QuasiLocalMeasures::qlm_invariant_coordinates + QuasiLocalMeasures::qlm_3determinant +" + +IOASCII::out2D_every = 128 +IOASCII::out2D_vars = " + SphericalSurface::sf_radius +" + + + +Activethorns = "CarpetIOHDF5" + +IOHDF5::out_every = 512 +IOHDF5::one_file_per_group = yes +#IOHDF5::output_symmetry_points = no +#IOHDF5::output_ghost_points = no +IOHDF5::compression_level = 1 +IOHDF5::use_checksums = yes +IOHDF5::out_vars = " + CarpetReduce::weight + ADMBase::alp + ML_ADMConstraints::ML_Ham + WeylScal4::psi4r + WeylScal4::psi4i + QuasiLocalMeasures::qlm_newman_penrose + QuasiLocalMeasures::qlm_weyl_scalars + QuasiLocalMeasures::qlm_ricci_scalars +" + +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 = 512 +TimerReport::out_filename = "TimerReport" +TimerReport::output_all_timers_together = yes +TimerReport::output_all_timers_readable = yes +TimerReport::n_top_timers = 20 +ADMBase::metric_timelevels=3 Directory: /PITTNullCode/SphericalHarmonicRecon/par/ ==================================================== File [added]: qc0-mclachlan-CCE_Null.par Delta lines: +109 -0 =================================================================== --- PITTNullCode/SphericalHarmonicRecon/par/qc0-mclachlan-CCE_Null.par (rev 0) +++ PITTNullCode/SphericalHarmonicRecon/par/qc0-mclachlan-CCE_Null.par 2011-06-30 15:43:02 UTC (rev 45) @@ -0,0 +1,109 @@ +#--------------------------------------- +# PARAM FILE FOR THE DYNAMIC MINKOWSKY +#--------------------------------------- + +ActiveThorns = "Fortran MoL SymBase CoordBase CartGrid3D Time" +ActiveThorns = "IOASCII IOBasic IOUtil" +ActiveThorns = "Boundary LocalReduce SpaceMask" +ActiveThorns = "Pugh PughReduce PughSlab PUGHInterp" +ActiveThorns = "AEILocalInterp NullGrid NullInterp NullVars" +ActiveThorns = "NullEvolve NullDecomp NullSHRExtract SphericalHarmonicRecon" +ActiveThorns = "NullNews NullConstr" +ActiveThorns = "LocalInterp HDF5 IOHDF5 IOHDF5Util" + +#Time------------------------------ +cactus::cctk_initial_time = 0.0 +cactus::cctk_final_time = 82.6 +cactus::terminate = "time" +#time::dtfac = 1 + +#Null Time----------------------------- +Time::timestep_method = "given" +Time::timestep = 0.05 + +#Initial data---------------------- +NullEvolve::boundary_data = "SHRE" +NullEvolve::initial_J_data = "smooth_J" +NullEvolve::first_order_scheme = "yes" + +#Extraction radius--------------- +SphericalHarmonicRecon::r_extract=20.0 #50.0 #100.0 +NullSHRExtract::cr = 20.0 #50.0 #100.0 +NullGrid::null_rwt = 20.0 #50.0 #100.0 + +#Override extraction parameters--------------- +#SphericalHarmonicRecon::time_derivative_in_file = "yes" +SphericalHarmonicRecon::metric_data_filename = "metric_obs_0_Decomp_ft.h5" +#SphericalHarmonicRecon::metric_data_filename = "metric_obs_1_Decomp.h5" +#SphericalHarmonicRecon::metric_data_filename = "metric_obs_2_Decomp.h5" + +#Extraction data------------------------ +NullSHRExtract::mass = 1 +NullSHRExtract::WT_metric = "Full" +NullSHRExtract::l_max = 6 + +#Null Grid --------------------------- +NullGrid::null_xin = 0.45 +NullGrid::N_ang_ghost_pts = 3 +NullGrid::N_ang_stencil_size = 3 +NullGrid::N_ang_ev_outside_eq = 20 +NullGrid::N_ang_pts_inside_eq = 101 +NullGrid::N_radial_pts = 113 + +#DEBUG: +NullEvolve::Diagnostics_Coord_x = 0.8 # 0.50877754829751409 + +#Interpolation-------------------------- +NullInterp::interpolation_order = 4 +NullInterp::stereo_patch_type = "circle" +NullInterp::deriv_accuracy = 4 + +#Dissipation------------------------------------------------------------------ +#NullEvolve::dissip_W = 0.000001 +#NullEvolve::dissip_Q = 0.000001 +NullEvolve::dissip_J = 0.0001 +NullEvolve::dissip_Jx = 0.001 +NullEvolve::dissip_mask_type = "zero at rD0, one at rD1" #0r01r1 +NullEvolve::N_dissip_zero_outside_eq = 11 +NullEvolve::N_dissip_one_outside_eq = 10 + +#Nan check-------------------------------- +ActiveThorns = "NaNChecker" +NaNChecker::action_if_found = "terminate" +NaNChecker::check_every = 1 +NaNChecker::check_vars = "nullvars::realcharfuncs" +NaNChecker::report_max = 1 + +#Output data and specifications----------- +IO::out_dir = NT100_BR20_ft +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = "NullNews::NewsB[0]" +IOBasic::outInfo_reductions = "norm2" +IOBasic::outScalar_every = 100 +IOBasic::outScalar_vars = "NullNews::NewsB[0] NullNews::Psi4[0] NullNews::Jn[0] NullNews::Jn_l[0] NullNews::Jn[1] NullNews::Jn_l[1]" +IOBasic::outScalar_reductions = "norm2" +IOBasic::out_format = ".16E" +IOASCII::out2D_every = 200 +IOASCII::out2D_vars = "NullNews::NewsB[0] NullNews::Psi4[0] NullNews::Jn[0] NullNews::Jn[1] NullNews::Jn_l[0] NullNews::Jn_l[1]" +IOASCII::out2D_style = "gnuplot f(x,y)" +IOASCII::out_format = ".16E" + +#Worltube Spherical Harmonics profiling--------------------------- +NullDecomp::use_rsYlm = no +NullDecomp::l_max = 4 +NullSHRExtract::WT_spherical_harmonics = yes + +#News profiling--------------------------- +NullNews::write_spherical_harmonics = yes +NullNews::debug_output = yes +NullNews::interp_to_constant_uBondi = yes +NullNews::news_interp_order = 4 +NullNews::mask_Psi4 = yes +NullNews::mask_NewsB = yes + +#Cactus::terminate = time +# +#Cactus::cctk_final_time = 0.0 +IO::checkpoint_on_terminate = yes +IO::checkpoint_ID = no +

1 0

[svn:einsteintoolkit] EinsteinExamples/trunk/par/ (Rev. 25)
by knarf＠cct.lsu.edu 28 Jun '11

28 Jun '11

User: knarf Date: 2011/06/28 06:01 PM Modified: /trunk/par/ static_tov.par Log: let the parfile use EOS_Omni, as this is required by GRHydro now File Changes: Directory: /trunk/par/ ====================== File [modified]: static_tov.par Delta lines: +1 -4 =================================================================== --- trunk/par/static_tov.par 2011-06-28 23:00:40 UTC (rev 24) +++ trunk/par/static_tov.par 2011-06-28 23:01:04 UTC (rev 25) @@ -114,7 +114,7 @@ ## Shift Condition: \partial_t beta^i = 0 # Hydro paramters -ActiveThorns = "EOS_Base EOS_IdealFluid EOS_Polytrope" +ActiveThorns = "EOS_Omni" ActiveThorns = "GRHydro" HydroBase::evolution_method = "GRHydro" @@ -129,10 +129,7 @@ #GRHydro::GRHydro = 18 # Tmunu(10), rho,press,eps,w_lorentz,vel, tau #GRHydro::GRHydro = 10 # gij(6), alpha, beta(3) -EOS_Polytrope::eos_gamma = 2.0 -EOS_Polytrope::eos_k = 100.0 - ActiveThorns = "GenericFD NewRad" ActiveThorns = "ML_BSSN ML_BSSN_Helper" ADMBase::evolution_method = "ML_BSSN"

1 0

[svn:einsteintoolkit] EinsteinExamples/branches/ET_2011_05/par/ (Rev. 24)
by knarf＠cct.lsu.edu 28 Jun '11

28 Jun '11

User: knarf Date: 2011/06/28 06:00 PM Modified: /branches/ET_2011_05/par/ static_tov.par Log: let the parfile use EOS_Omni, as this is required by GRHydro now File Changes: Directory: /branches/ET_2011_05/par/ ==================================== File [modified]: static_tov.par Delta lines: +1 -4 =================================================================== --- branches/ET_2011_05/par/static_tov.par 2011-04-21 04:44:59 UTC (rev 23) +++ branches/ET_2011_05/par/static_tov.par 2011-06-28 23:00:40 UTC (rev 24) @@ -114,7 +114,7 @@ ## Shift Condition: \partial_t beta^i = 0 # Hydro paramters -ActiveThorns = "EOS_Base EOS_IdealFluid EOS_Polytrope" +ActiveThorns = "EOS_Omni" ActiveThorns = "GRHydro" HydroBase::evolution_method = "GRHydro" @@ -129,10 +129,7 @@ #GRHydro::GRHydro = 18 # Tmunu(10), rho,press,eps,w_lorentz,vel, tau #GRHydro::GRHydro = 10 # gij(6), alpha, beta(3) -EOS_Polytrope::eos_gamma = 2.0 -EOS_Polytrope::eos_k = 100.0 - ActiveThorns = "GenericFD NewRad" ActiveThorns = "ML_BSSN ML_BSSN_Helper" ADMBase::evolution_method = "ML_BSSN"

1 0

[svn:einsteintoolkit] pyGWAnalysis/trunk/DataAnalysis/ (Rev. 5)
by reisswig＠tapir.caltech.edu 28 Jun '11

28 Jun '11

User: reisswig Date: 2011/06/28 09:44 AM Modified: /trunk/DataAnalysis/ DataAnalysis.py Log: Fix same bug in a related routine. File Changes: Directory: /trunk/DataAnalysis/ =============================== File [modified]: DataAnalysis.py Delta lines: +2 -2 =================================================================== --- trunk/DataAnalysis/DataAnalysis.py 2011-06-28 14:31:10 UTC (rev 4) +++ trunk/DataAnalysis/DataAnalysis.py 2011-06-28 14:44:54 UTC (rev 5) @@ -239,9 +239,9 @@ for ii in range(0, fr.Length()): if (fr.x[ii] != 0): - if fr.x[ii] < 0: fr.f[ii] = 0 + #if fr.x[ii] < 0: fr.f[ii] = 0 div = 2*pi*fr.x[ii] - if fr.x[ii] < f0: div = 2*pi*f0 + if abs(fr.x[ii]) < abs(f0): div = 2*pi*abs(f0)*sign(fr.x[ii]) fr.f[ii] = complex(0.,1.)**k * fr.f[ii]/(div**k) else:

1 0

[svn:einsteintoolkit] pyGWAnalysis/trunk/DataAnalysis/ (Rev. 4)
by reisswig＠tapir.caltech.edu 28 Jun '11

28 Jun '11

User: reisswig Date: 2011/06/28 09:31 AM Modified: /trunk/DataAnalysis/ DataAnalysis.py Log: Fixed a bug in integration from news->h File Changes: Directory: /trunk/DataAnalysis/ =============================== File [modified]: DataAnalysis.py Delta lines: +3 -2 =================================================================== --- trunk/DataAnalysis/DataAnalysis.py 2011-06-14 23:19:15 UTC (rev 3) +++ trunk/DataAnalysis/DataAnalysis.py 2011-06-28 14:31:10 UTC (rev 4) @@ -298,9 +298,9 @@ for ii in range(0, fr.Length()): if (fr.x[ii] != 0): - if fr.x[ii] < 0: fr.f[ii] = 0 + #if fr.x[ii] < 0: fr.f[ii] = 0 div = 2*pi*fr.x[ii] - if abs(fr.x[ii]) < abs(frange[0]): div = 2*pi*frange[0]*sign(fr.x[ii]) + if abs(fr.x[ii]) < abs(frange[0]): div = 2*pi*abs(frange[0])*sign(fr.x[ii]) #if fr.x[ii] > frange[1]: div = 4*pi**2*frange[1]**2 #if DF1fr.x[ii] > frange[0] and DF1fr.x[ii] < frange[1]: #div = (1.0-GaussBlend(DF1fr.x[ii], frange[0], (frange[1]-frange[0])/100)) * div + GaussBlend(DF1fr.x[ii], frange[0], (frange[1]-frange[0])/100) * 4*pi**2*frange[0]**2 # div = (1 - (DF1fr.x[ii]-frange[0])/(frange[1]-frange[0])) * 4*pi**2*frange[0]**2 + (DF1fr.x[ii]-frange[0])/(frange[1]-frange[0]) * div @@ -317,6 +317,7 @@ return res + def PolynomialFilterDriver(WFfr, x0, x1, blend=False): """Given a waveform in Fourier domain, we fit a polynomial between frequncy range frange as described in Reisswig&Pollney 2010"""

1 0

[svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 93)
by knarf＠cct.lsu.edu 27 Jun '11

27 Jun '11

User: knarf Date: 2011/06/27 10:54 AM Modified: / ET.tex Log: define a command to mark sections which are done or almost done, for an easier overview File Changes: Directory: / ============ File [modified]: ET.tex Delta lines: +8 -4 =================================================================== --- ET.tex 2011-06-10 22:52:35 UTC (rev 92) +++ ET.tex 2011-06-27 15:54:01 UTC (rev 93) @@ -35,6 +35,11 @@ \newcommand{\comment}[1]{{\color{blue}$\blacksquare$~\textsf{[Comment: #1]}}} \newcommand{\pages}[1]{{\color{blue}$\blacksquare$~\textsf{[#1]}}} +% Use this macro for sections which are done, so that those comments +% can easily be disabled. +\newcommand{\pagesdone}[1]{{\color{green}$\blacksquare$~\textsf{[#1]}}} +%\newcommand{\pagesdone}[1]{} + \newcommand{\BCM}[1]{{\bf \color{blue} [BCM: #1] }} %BCM comments % Don't use tt font for urls @@ -98,7 +103,6 @@ \tableofcontents \section{Introduction\pages{1 Frank}} -\comment{About 1 page} Scientific progress in the field of numerical relativity has always been closely tied with the availability and ease-of-use of enabling @@ -332,7 +336,7 @@ \section{Core Technologies\pages{0.5 Frank}} -\subsection{Cactus Framework\pages{1 Gab/Frank}} +\subsection{Cactus Framework\pages{1 Frank}} The \textbf{Cactus Framework}~\cite{CS_cactus_web,CS_Goodale02a,CS_cactususersguide} is @@ -375,7 +379,7 @@ The Einstein Toolkit now collects the widely used parts of CactusEinstein, combined with contributions from the community. -\subsection{Adaptive Mesh Refinement\pages{1 Erik}} +\subsection{Adaptive Mesh Refinement\pagesdone{1 Erik}} In Cactus, infrastructure capabilities such as memory management, parallelisation, time evolution, mesh refinement, and I/O are @@ -454,7 +458,7 @@ world-wide. To date, more than 90 peer-reviewed publication and more than 15 student theses are based on Carpet \cite{CarpetCode:web}. -\subsection{Simulation Factory\pages{1 Erik}} +\subsection{Simulation Factory\pagesdone{1 Erik}} Today's supercomputers differ significantly in their hardware configuration, available software, directory structure,

1 0

[svn:einsteintoolkit] incoming/ (Rev. 44)
by knarf＠cct.lsu.edu 21 Jun '11

21 Jun '11

User: knarf Date: 2011/06/21 12:25 PM Removed: / test Log: remove test file File Changes: Directory: / ============ File [removed]: test Delta lines: +0 -0 Deleted: test

1 0

[svn:einsteintoolkit] incoming/ (Rev. 43)
by knarf＠cct.lsu.edu 21 Jun '11

21 Jun '11

User: knarf Date: 2011/06/21 12:24 PM Added: / test Log: testing File Changes: Directory: / ============ File [added]: test Delta lines: +0 -0 Added: test

1 0

[svn:einsteintoolkit] incoming/PITTNullCode/ (Rev. 42)
by reisswig＠tapir.caltech.edu 15 Jun '11

15 Jun '11

User: reisswig Date: 2011/06/15 03:28 PM Added: /PITTNullCode/SphericalHarmonicRecon/ README, configuration.ccl, interface.ccl, param.ccl, schedule.ccl /PITTNullCode/SphericalHarmonicRecon/codes/ Makefile, README, ascii_output.c, dampfilter.c, fftwfilter.c, findlast.c, printtime.c, readmeta.c, setmeta.c, util.c, util.h /PITTNullCode/SphericalHarmonicRecon/doc/ documentation.tex /PITTNullCode/SphericalHarmonicRecon/par/ /PITTNullCode/SphericalHarmonicRecon/src/ Chebyshev.cc, Chebyshev.hh, Recon.cc, make.code.defn /PITTNullCode/SphericalHarmonicRecon/test/ Modified: /PITTNullCode/SphericalHarmonicDecomp/src/ ExtractMetric.cc, matrix.cc Log: Latest version. File Changes: Directory: /PITTNullCode/SphericalHarmonicRecon/ ================================================ File [added]: README Delta lines: <Unavailable> File [added]: configuration.ccl Delta lines: <Unavailable> File [added]: interface.ccl Delta lines: <Unavailable> File [added]: param.ccl Delta lines: <Unavailable> File [added]: schedule.ccl Delta lines: <Unavailable> Directory: /PITTNullCode/SphericalHarmonicRecon/codes/ ====================================================== File [added]: Makefile Delta lines: <Unavailable> File [added]: README Delta lines: <Unavailable> File [added]: ascii_output.c Delta lines: <Unavailable> File [added]: dampfilter.c Delta lines: <Unavailable> File [added]: fftwfilter.c Delta lines: <Unavailable> File [added]: findlast.c Delta lines: <Unavailable> File [added]: printtime.c Delta lines: <Unavailable> File [added]: readmeta.c Delta lines: <Unavailable> File [added]: setmeta.c Delta lines: <Unavailable> File [added]: util.c Delta lines: <Unavailable> File [added]: util.h Delta lines: <Unavailable> Directory: /PITTNullCode/SphericalHarmonicRecon/doc/ ==================================================== File [added]: documentation.tex Delta lines: <Unavailable> Directory: /PITTNullCode/SphericalHarmonicDecomp/src/ ===================================================== File [modified]: ExtractMetric.cc Delta lines: <Unavailable> File [modified]: matrix.cc Delta lines: <Unavailable> Directory: /PITTNullCode/SphericalHarmonicRecon/src/ ==================================================== File [added]: Chebyshev.cc Delta lines: <Unavailable> File [added]: Chebyshev.hh Delta lines: <Unavailable> File [added]: Recon.cc Delta lines: <Unavailable> File [added]: make.code.defn Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] incoming/PITTNullCode/SphericalHarmonicRecon/ (Rev. 41)
by reisswig＠tapir.caltech.edu 15 Jun '11

15 Jun '11

User: reisswig Date: 2011/06/15 03:26 PM Removed: /PITTNullCode/SphericalHarmonicRecon/ Log: Clean version.

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits