Dear Users.
I am runing --parfile repos/ctthorns/CT_MultiLevel/par/twopunctures.par.
###################################### # Flesh parameters ######################################
Cactus::cctk_run_title = "Solving the Hamiltonian constraint for two punctures with CT_MultiLevel (multigrid setup, FMG) and with TwoPunctures"
Cactus::cctk_full_warnings = yes Cactus::highlight_warning_messages = no Cactus::cctk_timer_output = "off"
Cactus::terminate = "any" Cactus::cctk_initial_time = 0 Cactus::cctk_final_time = 30 Cactus::cctk_itlast = 200
############################################################################################# ActiveThorns = "Boundary CartGrid3D CoordBase Slab SymBase" #############################################################################################
CoordBase::domainsize = "minmax"
CoordBase::xmin = -40.00 CoordBase::ymin = -40.00 CoordBase::zmin = -40.00 CoordBase::xmax = +40.00 CoordBase::ymax = +40.00 CoordBase::zmax = +40.00 CoordBase::dx = 10 CoordBase::dy = 10 CoordBase::dz = 10
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 = 1 CoordBase::boundary_shiftout_y_upper = 1 CoordBase::boundary_shiftout_z_upper = 1
CartGrid3D::type = "coordbase"
################################################################################## ActiveThorns = "IOUtil Carpet CarpetLib CarpetInterp LoopControl CarpetReduce CarpetSlab" ##################################################################################
Carpet::adaptive_stepsize = "no"
Carpet::schedule_barriers = no CarpetLib::output_bboxes = no Carpet::init_fill_timelevels = "yes"
Carpet::domain_from_coordbase = yes Carpet::max_refinement_levels = 5
driver::ghost_size = 3 Carpet::use_buffer_zones = yes
Carpet::prolongation_order_space = 5 Carpet::prolongation_order_time = 2
Carpet::convergence_level = 0
Carpet::poison_new_timelevels = yes CarpetLib::poison_new_memory = yes
Carpet::output_timers_every = 0 CarpetLib::print_timestats_every = 0 CarpetLib::print_memstats_every = 0
CarpetLib::max_allowed_memory_MB = 12288
################################################################## ActiveThorns = "CarpetRegrid2" ##################################################################
Carpet::grid_structure_filename = "grid_structure.asc" Carpet::grid_coordinates_filename = "grid_coordinates.asc"
CarpetRegrid2::num_centres = 1 CarpetRegrid2::num_levels_1 = 4 CarpetRegrid2::radius_1[1] = 50 CarpetRegrid2::radius_1[2] = 25 CarpetRegrid2::radius_1[3] = 12 #CarpetRegrid2::radius_1[4] = 50 #CarpetRegrid2::radius_1[5] = 25 #CarpetRegrid2::radius_1[6] = 12 #CarpetRegrid2::radius_1[7] = 6
################################# ActiveThorns = "AEILocalInterp" ActiveThorns = "BLAS LAPACK" ActiveThorns = "Fortran" ActiveThorns = "GenericFD" ActiveThorns = "HDF5" ActiveThorns = "InitBase" ActiveThorns = "LocalInterp" #################################
################################# ActiveThorns = "MoL Time" #################################
MoL::ODE_Method = "RK4" MoL::MoL_Intermediate_Steps = 4 MoL::MoL_Num_Scratch_Levels = 1
Time::dtfac = 0.2
################################################################################################### ActiveThorns = "CT_MultiLevel CT_Analytic ADMBase ML_BSSN ML_BSSN_Helper CoordGauge StaticConformal TmunuBase ADMCoupling ADMMacros TwoPunctures" ###################################################################################################
ADMBase::metric_timelevels = 3 ADMBase::lapse_timelevels = 3 ADMBase::shift_timelevels = 3
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"
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"
ADMBase::admbase_boundary_condition = "none"
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
#TwoPunctures::grid_setup_method = "evaluation"
TwoPunctures::TP_epsilon = 1.0e-2 TwoPunctures::TP_Tiny = 1.0e-2
TwoPunctures::keep_u_around = "yes"
CT_MultiLevel::cycle_type = "FMG cycle" CT_MultiLevel::tol = 1e-12 CT_MultiLevel::omega = 1 CT_MultiLevel::fd_order = 4 CT_MultiLevel::output_norms = "yes" CT_MultiLevel::output_walk = "yes" CT_MultiLevel::nrelsteps_down = 700 CT_MultiLevel::nrelsteps_up = 650 CT_MultiLevel::nrelsteps_bottom = 700 CT_MultiLevel::nrelsteps_top = 650 CT_MultiLevel::topMGlevel = 4 CT_MultiLevel::number_of_equations = 1 CT_MultiLevel::number_of_auxiliaries = 1 CT_MultiLevel::fill_ADM = "yes" CT_MultiLevel::fill_Aij = "Analytic Aij" CT_MultiLevel::boundary_conditions = "TwoPunctures"
CT_MultiLevel::compare_to_exact = "yes" CT_MultiLevel::exact_solution_gfname[0] = "TwoPunctures::puncture_u" CT_MultiLevel::exact_offset = -1.0
CT_MultiLevel::inipsi_gfname[0] = "CT_Analytic::testinipsi" CT_MultiLevel::cxx_gfname[0] = "CT_Analytic::testcxx" CT_MultiLevel::cyy_gfname[0] = "CT_Analytic::testcyy" CT_MultiLevel::czz_gfname[0] = "CT_Analytic::testczz" CT_MultiLevel::n0[0] = -7 CT_MultiLevel::c0_gfname[0] = "CT_MultiLevel::ct_auxiliary[0]" CT_MultiLevel::a0_gfname[0] = "CT_Analytic::testa0"
CT_MultiLevel::model = "Bowen-York" CT_Analytic::free_data = "Bowen-York" CT_Analytic::ampI = 1 CT_Analytic::ampC = 1 CT_Analytic::massa = 0.453 CT_Analytic::massb = 0.453 CT_Analytic::xa = 1.168642873 CT_Analytic::xb = -1.168642873 CT_Analytic::Pay = 0.3331917498 CT_Analytic::Pby = -0.3331917498
CT_Analytic::other_timelevels = 3
################################# ActiveThorns = "CarpetIOBasic" #################################
IOBasic::outInfo_every = 1 IOBasic::outInfo_reductions = "norm2" IOBasic::outInfo_vars = " CT_MultiLevel::psi Carpet::physical_time_per_hour "
################################# ActiveThorns = "CarpetIOScalar" #################################
IO::out_dir = $parfile
IOScalar::one_file_per_group = yes
IOScalar::outScalar_every = 0 IOScalar::outScalar_vars = " "
################################# ActiveThorns = "CarpetIOASCII" #################################
IOASCII::one_file_per_group = yes
IOASCII::output_symmetry_points = yes IOASCII::out3D_ghosts = yes
IOASCII::out0D_criterion = "divisor" IOASCII::out0D_every = 0 IOASCII::out0D_vars = " "
IOASCII::output_all_timelevels = "no" IOASCII::out1D_criterion = "divisor" IOASCII::out1D_every = 1 IOASCII::out1D_vars = " CT_MultiLevel::psi CT_MultiLevel::residual CT_MultiLevel::err ADMBase::metric ADMBase::curv ML_BSSN::ML_log_confac ML_BSSN::ML_ham TwoPunctures::puncture_u "
IOASCII::out2D_every = 0 IOASCII::out2D_criterion = "divisor" IOASCII::out2D_vars = " "
IOASCII::out3D_every = 0 IOASCII::out3D_criterion = "divisor" IOASCII::out3D_vars = " CT_MultiLevel::psi CT_MultiLevel::coeffs "
################################# Activethorns = "CarpetIOHDF5" #################################
IOHDF5::out_criterion = "divisor" IOHDF5::out_every = 1 IOHDF5::one_file_per_group = yes IOHDF5::compression_level = 0 IOHDF5::out_vars = " CT_MultiLevel::psi CT_MultiLevel::err CT_MultiLevel::residual ML_BSSN::ML_log_confac ML_BSSN::ML_ham TwoPunctures::puncture_u "
IOHDF5::checkpoint = "no"
################################# ActiveThorns = "Formaline" #################################
################################# ActiveThorns = "TimerReport" #################################
TimerReport::out_every = 512 TimerReport::out_filename = "TimerReport"
*My ouput shows -nan after first itration. *
0 0.000 | 1.0013688 | 0.0000000 2 0.250 | -nan | 1.128073e+03 4 0.500 | -nan | 1.382569e+03 6 0.750 | -nan | 1.318078e+03 8 1.000 | -nan | 1.399377e+03 10 1.250 | -nan | 1.336436e+03 12 1.500 | -nan | 1.395822e+03 14 1.750 | -nan | 1.364895e+03 16 2.000 | -nan | 1.398991e+03 18 2.250 | -nan | 1.359115e+03 20 2.500 | -nan | 1.393235e+03 I set refinement levels, radius, topMG. no error found in running.
Kindly,
Hi Bilal,
You are under-resolving the black holes. For stability you will need gridspacing of at least ~1/16 around each black hole (M/16, where M=1 for this binary), and here your gridspacing is 1.25. I would suggest starting by trying to reproduce the binary black hole gallery example: https://einsteintoolkit.org/gallery/bbh/index.html
and adding small adjustments until you get it working according to your needs.
-Zach
* * * Prof. Zachariah Etienne Physics & Astronomy Dept. West Virginia University http://astro.phys.wvu.edu/zetienne/ http://blackholesathome.net https://blackholesathome.net
On Thu, Aug 27, 2020 at 2:38 PM Bilal Hameed bilalhameedawan962@gmail.com wrote:
Dear Users.
I am runing --parfile repos/ctthorns/CT_MultiLevel/par/twopunctures.par.
###################################### # Flesh parameters ######################################
Cactus::cctk_run_title = "Solving the Hamiltonian constraint for two punctures with CT_MultiLevel (multigrid setup, FMG) and with TwoPunctures"
Cactus::cctk_full_warnings = yes Cactus::highlight_warning_messages = no Cactus::cctk_timer_output = "off"
Cactus::terminate = "any" Cactus::cctk_initial_time = 0 Cactus::cctk_final_time = 30 Cactus::cctk_itlast = 200
############################################################################################# ActiveThorns = "Boundary CartGrid3D CoordBase Slab SymBase"
#############################################################################################
CoordBase::domainsize = "minmax"
CoordBase::xmin = -40.00 CoordBase::ymin = -40.00 CoordBase::zmin = -40.00 CoordBase::xmax = +40.00 CoordBase::ymax = +40.00 CoordBase::zmax = +40.00 CoordBase::dx = 10 CoordBase::dy = 10 CoordBase::dz = 10
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 = 1 CoordBase::boundary_shiftout_y_upper = 1 CoordBase::boundary_shiftout_z_upper = 1
CartGrid3D::type = "coordbase"
################################################################################## ActiveThorns = "IOUtil Carpet CarpetLib CarpetInterp LoopControl CarpetReduce CarpetSlab"
##################################################################################
Carpet::adaptive_stepsize = "no"
Carpet::schedule_barriers = no CarpetLib::output_bboxes = no Carpet::init_fill_timelevels = "yes"
Carpet::domain_from_coordbase = yes Carpet::max_refinement_levels = 5
driver::ghost_size = 3 Carpet::use_buffer_zones = yes
Carpet::prolongation_order_space = 5 Carpet::prolongation_order_time = 2
Carpet::convergence_level = 0
Carpet::poison_new_timelevels = yes CarpetLib::poison_new_memory = yes
Carpet::output_timers_every = 0 CarpetLib::print_timestats_every = 0 CarpetLib::print_memstats_every = 0
CarpetLib::max_allowed_memory_MB = 12288
################################################################## ActiveThorns = "CarpetRegrid2" ##################################################################
Carpet::grid_structure_filename = "grid_structure.asc" Carpet::grid_coordinates_filename = "grid_coordinates.asc"
CarpetRegrid2::num_centres = 1 CarpetRegrid2::num_levels_1 = 4 CarpetRegrid2::radius_1[1] = 50 CarpetRegrid2::radius_1[2] = 25 CarpetRegrid2::radius_1[3] = 12 #CarpetRegrid2::radius_1[4] = 50 #CarpetRegrid2::radius_1[5] = 25 #CarpetRegrid2::radius_1[6] = 12 #CarpetRegrid2::radius_1[7] = 6
################################# ActiveThorns = "AEILocalInterp" ActiveThorns = "BLAS LAPACK" ActiveThorns = "Fortran" ActiveThorns = "GenericFD" ActiveThorns = "HDF5" ActiveThorns = "InitBase" ActiveThorns = "LocalInterp" #################################
################################# ActiveThorns = "MoL Time" #################################
MoL::ODE_Method = "RK4" MoL::MoL_Intermediate_Steps = 4 MoL::MoL_Num_Scratch_Levels = 1
Time::dtfac = 0.2
################################################################################################### ActiveThorns = "CT_MultiLevel CT_Analytic ADMBase ML_BSSN ML_BSSN_Helper CoordGauge StaticConformal TmunuBase ADMCoupling ADMMacros TwoPunctures"
###################################################################################################
ADMBase::metric_timelevels = 3 ADMBase::lapse_timelevels = 3 ADMBase::shift_timelevels = 3
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"
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"
ADMBase::admbase_boundary_condition = "none"
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
#TwoPunctures::grid_setup_method = "evaluation"
TwoPunctures::TP_epsilon = 1.0e-2 TwoPunctures::TP_Tiny = 1.0e-2
TwoPunctures::keep_u_around = "yes"
CT_MultiLevel::cycle_type = "FMG cycle" CT_MultiLevel::tol = 1e-12 CT_MultiLevel::omega = 1 CT_MultiLevel::fd_order = 4 CT_MultiLevel::output_norms = "yes" CT_MultiLevel::output_walk = "yes" CT_MultiLevel::nrelsteps_down = 700 CT_MultiLevel::nrelsteps_up = 650 CT_MultiLevel::nrelsteps_bottom = 700 CT_MultiLevel::nrelsteps_top = 650 CT_MultiLevel::topMGlevel = 4 CT_MultiLevel::number_of_equations = 1 CT_MultiLevel::number_of_auxiliaries = 1 CT_MultiLevel::fill_ADM = "yes" CT_MultiLevel::fill_Aij = "Analytic Aij" CT_MultiLevel::boundary_conditions = "TwoPunctures"
CT_MultiLevel::compare_to_exact = "yes" CT_MultiLevel::exact_solution_gfname[0] = "TwoPunctures::puncture_u" CT_MultiLevel::exact_offset = -1.0
CT_MultiLevel::inipsi_gfname[0] = "CT_Analytic::testinipsi" CT_MultiLevel::cxx_gfname[0] = "CT_Analytic::testcxx" CT_MultiLevel::cyy_gfname[0] = "CT_Analytic::testcyy" CT_MultiLevel::czz_gfname[0] = "CT_Analytic::testczz" CT_MultiLevel::n0[0] = -7 CT_MultiLevel::c0_gfname[0] = "CT_MultiLevel::ct_auxiliary[0]" CT_MultiLevel::a0_gfname[0] = "CT_Analytic::testa0"
CT_MultiLevel::model = "Bowen-York" CT_Analytic::free_data = "Bowen-York" CT_Analytic::ampI = 1 CT_Analytic::ampC = 1 CT_Analytic::massa = 0.453 CT_Analytic::massb = 0.453 CT_Analytic::xa = 1.168642873 CT_Analytic::xb = -1.168642873 CT_Analytic::Pay = 0.3331917498 CT_Analytic::Pby = -0.3331917498
CT_Analytic::other_timelevels = 3
################################# ActiveThorns = "CarpetIOBasic" #################################
IOBasic::outInfo_every = 1 IOBasic::outInfo_reductions = "norm2" IOBasic::outInfo_vars = " CT_MultiLevel::psi Carpet::physical_time_per_hour "
################################# ActiveThorns = "CarpetIOScalar" #################################
IO::out_dir = $parfile
IOScalar::one_file_per_group = yes
IOScalar::outScalar_every = 0 IOScalar::outScalar_vars = " "
################################# ActiveThorns = "CarpetIOASCII" #################################
IOASCII::one_file_per_group = yes
IOASCII::output_symmetry_points = yes IOASCII::out3D_ghosts = yes
IOASCII::out0D_criterion = "divisor" IOASCII::out0D_every = 0 IOASCII::out0D_vars = " "
IOASCII::output_all_timelevels = "no" IOASCII::out1D_criterion = "divisor" IOASCII::out1D_every = 1 IOASCII::out1D_vars = " CT_MultiLevel::psi CT_MultiLevel::residual CT_MultiLevel::err ADMBase::metric ADMBase::curv ML_BSSN::ML_log_confac ML_BSSN::ML_ham TwoPunctures::puncture_u "
IOASCII::out2D_every = 0 IOASCII::out2D_criterion = "divisor" IOASCII::out2D_vars = " "
IOASCII::out3D_every = 0 IOASCII::out3D_criterion = "divisor" IOASCII::out3D_vars = " CT_MultiLevel::psi CT_MultiLevel::coeffs "
################################# Activethorns = "CarpetIOHDF5" #################################
IOHDF5::out_criterion = "divisor" IOHDF5::out_every = 1 IOHDF5::one_file_per_group = yes IOHDF5::compression_level = 0 IOHDF5::out_vars = " CT_MultiLevel::psi CT_MultiLevel::err CT_MultiLevel::residual ML_BSSN::ML_log_confac ML_BSSN::ML_ham TwoPunctures::puncture_u "
IOHDF5::checkpoint = "no"
################################# ActiveThorns = "Formaline" #################################
################################# ActiveThorns = "TimerReport" #################################
TimerReport::out_every = 512 TimerReport::out_filename = "TimerReport"
*My ouput shows -nan after first itration. *
0 0.000 | 1.0013688 | 0.0000000 2 0.250 | -nan | 1.128073e+03 4 0.500 | -nan | 1.382569e+03 6 0.750 | -nan | 1.318078e+03 8 1.000 | -nan | 1.399377e+03 10 1.250 | -nan | 1.336436e+03 12 1.500 | -nan | 1.395822e+03 14 1.750 | -nan | 1.364895e+03 16 2.000 | -nan | 1.398991e+03 18 2.250 | -nan | 1.359115e+03 20 2.500 | -nan | 1.393235e+03 I set refinement levels, radius, topMG. no error found in running.
Kindly,
Users mailing list Users@einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/users
Thanks Zach.
On Thu, Aug 27, 2020, 11:50 PM Zach Etienne zachetie@gmail.com wrote:
Hi Bilal,
You are under-resolving the black holes. For stability you will need gridspacing of at least ~1/16 around each black hole (M/16, where M=1 for this binary), and here your gridspacing is 1.25. I would suggest starting by trying to reproduce the binary black hole gallery example: https://einsteintoolkit.org/gallery/bbh/index.html
and adding small adjustments until you get it working according to your needs.
-Zach
Prof. Zachariah Etienne Physics & Astronomy Dept. West Virginia University http://astro.phys.wvu.edu/zetienne/ http://blackholesathome.net https://blackholesathome.net
On Thu, Aug 27, 2020 at 2:38 PM Bilal Hameed bilalhameedawan962@gmail.com wrote:
Dear Users.
I am runing --parfile repos/ctthorns/CT_MultiLevel/par/twopunctures.par.
###################################### # Flesh parameters ######################################
Cactus::cctk_run_title = "Solving the Hamiltonian constraint for two punctures with CT_MultiLevel (multigrid setup, FMG) and with TwoPunctures"
Cactus::cctk_full_warnings = yes Cactus::highlight_warning_messages = no Cactus::cctk_timer_output = "off"
Cactus::terminate = "any" Cactus::cctk_initial_time = 0 Cactus::cctk_final_time = 30 Cactus::cctk_itlast = 200
############################################################################################# ActiveThorns = "Boundary CartGrid3D CoordBase Slab SymBase"
#############################################################################################
CoordBase::domainsize = "minmax"
CoordBase::xmin = -40.00 CoordBase::ymin = -40.00 CoordBase::zmin = -40.00 CoordBase::xmax = +40.00 CoordBase::ymax = +40.00 CoordBase::zmax = +40.00 CoordBase::dx = 10 CoordBase::dy = 10 CoordBase::dz = 10
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 = 1 CoordBase::boundary_shiftout_y_upper = 1 CoordBase::boundary_shiftout_z_upper = 1
CartGrid3D::type = "coordbase"
################################################################################## ActiveThorns = "IOUtil Carpet CarpetLib CarpetInterp LoopControl CarpetReduce CarpetSlab"
##################################################################################
Carpet::adaptive_stepsize = "no"
Carpet::schedule_barriers = no CarpetLib::output_bboxes = no Carpet::init_fill_timelevels = "yes"
Carpet::domain_from_coordbase = yes Carpet::max_refinement_levels = 5
driver::ghost_size = 3 Carpet::use_buffer_zones = yes
Carpet::prolongation_order_space = 5 Carpet::prolongation_order_time = 2
Carpet::convergence_level = 0
Carpet::poison_new_timelevels = yes CarpetLib::poison_new_memory = yes
Carpet::output_timers_every = 0 CarpetLib::print_timestats_every = 0 CarpetLib::print_memstats_every = 0
CarpetLib::max_allowed_memory_MB = 12288
################################################################## ActiveThorns = "CarpetRegrid2" ##################################################################
Carpet::grid_structure_filename = "grid_structure.asc" Carpet::grid_coordinates_filename = "grid_coordinates.asc"
CarpetRegrid2::num_centres = 1 CarpetRegrid2::num_levels_1 = 4 CarpetRegrid2::radius_1[1] = 50 CarpetRegrid2::radius_1[2] = 25 CarpetRegrid2::radius_1[3] = 12 #CarpetRegrid2::radius_1[4] = 50 #CarpetRegrid2::radius_1[5] = 25 #CarpetRegrid2::radius_1[6] = 12 #CarpetRegrid2::radius_1[7] = 6
################################# ActiveThorns = "AEILocalInterp" ActiveThorns = "BLAS LAPACK" ActiveThorns = "Fortran" ActiveThorns = "GenericFD" ActiveThorns = "HDF5" ActiveThorns = "InitBase" ActiveThorns = "LocalInterp" #################################
################################# ActiveThorns = "MoL Time" #################################
MoL::ODE_Method = "RK4" MoL::MoL_Intermediate_Steps = 4 MoL::MoL_Num_Scratch_Levels = 1
Time::dtfac = 0.2
################################################################################################### ActiveThorns = "CT_MultiLevel CT_Analytic ADMBase ML_BSSN ML_BSSN_Helper CoordGauge StaticConformal TmunuBase ADMCoupling ADMMacros TwoPunctures"
###################################################################################################
ADMBase::metric_timelevels = 3 ADMBase::lapse_timelevels = 3 ADMBase::shift_timelevels = 3
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"
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"
ADMBase::admbase_boundary_condition = "none"
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
#TwoPunctures::grid_setup_method = "evaluation"
TwoPunctures::TP_epsilon = 1.0e-2 TwoPunctures::TP_Tiny = 1.0e-2
TwoPunctures::keep_u_around = "yes"
CT_MultiLevel::cycle_type = "FMG cycle" CT_MultiLevel::tol = 1e-12 CT_MultiLevel::omega = 1 CT_MultiLevel::fd_order = 4 CT_MultiLevel::output_norms = "yes" CT_MultiLevel::output_walk = "yes" CT_MultiLevel::nrelsteps_down = 700 CT_MultiLevel::nrelsteps_up = 650 CT_MultiLevel::nrelsteps_bottom = 700 CT_MultiLevel::nrelsteps_top = 650 CT_MultiLevel::topMGlevel = 4 CT_MultiLevel::number_of_equations = 1 CT_MultiLevel::number_of_auxiliaries = 1 CT_MultiLevel::fill_ADM = "yes" CT_MultiLevel::fill_Aij = "Analytic Aij" CT_MultiLevel::boundary_conditions = "TwoPunctures"
CT_MultiLevel::compare_to_exact = "yes" CT_MultiLevel::exact_solution_gfname[0] = "TwoPunctures::puncture_u" CT_MultiLevel::exact_offset = -1.0
CT_MultiLevel::inipsi_gfname[0] = "CT_Analytic::testinipsi" CT_MultiLevel::cxx_gfname[0] = "CT_Analytic::testcxx" CT_MultiLevel::cyy_gfname[0] = "CT_Analytic::testcyy" CT_MultiLevel::czz_gfname[0] = "CT_Analytic::testczz" CT_MultiLevel::n0[0] = -7 CT_MultiLevel::c0_gfname[0] = "CT_MultiLevel::ct_auxiliary[0]" CT_MultiLevel::a0_gfname[0] = "CT_Analytic::testa0"
CT_MultiLevel::model = "Bowen-York" CT_Analytic::free_data = "Bowen-York" CT_Analytic::ampI = 1 CT_Analytic::ampC = 1 CT_Analytic::massa = 0.453 CT_Analytic::massb = 0.453 CT_Analytic::xa = 1.168642873 CT_Analytic::xb = -1.168642873 CT_Analytic::Pay = 0.3331917498 CT_Analytic::Pby = -0.3331917498
CT_Analytic::other_timelevels = 3
################################# ActiveThorns = "CarpetIOBasic" #################################
IOBasic::outInfo_every = 1 IOBasic::outInfo_reductions = "norm2" IOBasic::outInfo_vars = " CT_MultiLevel::psi Carpet::physical_time_per_hour "
################################# ActiveThorns = "CarpetIOScalar" #################################
IO::out_dir = $parfile
IOScalar::one_file_per_group = yes
IOScalar::outScalar_every = 0 IOScalar::outScalar_vars = " "
################################# ActiveThorns = "CarpetIOASCII" #################################
IOASCII::one_file_per_group = yes
IOASCII::output_symmetry_points = yes IOASCII::out3D_ghosts = yes
IOASCII::out0D_criterion = "divisor" IOASCII::out0D_every = 0 IOASCII::out0D_vars = " "
IOASCII::output_all_timelevels = "no" IOASCII::out1D_criterion = "divisor" IOASCII::out1D_every = 1 IOASCII::out1D_vars = " CT_MultiLevel::psi CT_MultiLevel::residual CT_MultiLevel::err ADMBase::metric ADMBase::curv ML_BSSN::ML_log_confac ML_BSSN::ML_ham TwoPunctures::puncture_u "
IOASCII::out2D_every = 0 IOASCII::out2D_criterion = "divisor" IOASCII::out2D_vars = " "
IOASCII::out3D_every = 0 IOASCII::out3D_criterion = "divisor" IOASCII::out3D_vars = " CT_MultiLevel::psi CT_MultiLevel::coeffs "
################################# Activethorns = "CarpetIOHDF5" #################################
IOHDF5::out_criterion = "divisor" IOHDF5::out_every = 1 IOHDF5::one_file_per_group = yes IOHDF5::compression_level = 0 IOHDF5::out_vars = " CT_MultiLevel::psi CT_MultiLevel::err CT_MultiLevel::residual ML_BSSN::ML_log_confac ML_BSSN::ML_ham TwoPunctures::puncture_u "
IOHDF5::checkpoint = "no"
################################# ActiveThorns = "Formaline" #################################
################################# ActiveThorns = "TimerReport" #################################
TimerReport::out_every = 512 TimerReport::out_filename = "TimerReport"
*My ouput shows -nan after first itration. *
0 0.000 | 1.0013688 | 0.0000000 2 0.250 | -nan | 1.128073e+03 4 0.500 | -nan | 1.382569e+03 6 0.750 | -nan | 1.318078e+03 8 1.000 | -nan | 1.399377e+03 10 1.250 | -nan | 1.336436e+03 12 1.500 | -nan | 1.395822e+03 14 1.750 | -nan | 1.364895e+03 16 2.000 | -nan | 1.398991e+03 18 2.250 | -nan | 1.359115e+03 20 2.500 | -nan | 1.393235e+03 I set refinement levels, radius, topMG. no error found in running.
Kindly,
Users mailing list Users@einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/users
Hello Bilal,
following up on the private conversation we had:
hmm.
So you may not want to change itlast.
Note that CT_MultiLevel is an elliptic solver, so you are constructing initial data.
Looking at:
repos/ctthorns/CT_MultiLevel/schedule.cc
there is only:
schedule CT_MultiLevel AT CCTK_INITIAL after CT_ScalarField_setConfRho { LANG: C OPTIONS: GLOBAL-LATE } "Main MultiLevel function"
ie only in INITIAL which is for itlast=0.
To the actually run an evolution (the black holes actually moving) you would have to read in the data produced by CT_MultiLevel into an evolution simulation. This would involve using a parameter file similar to the one that Zach pointed you to and a thorn like the ReadInterpolate thorn to read in the results data from CT_MultiLevel into the evolution simulation.
How to do that is not fully trivial, unfortunately. You would have to ask someone on the mailing list if they have an example parameter file. Possible the CT_MultiLevel author, who had used my thorn before.
Yours, Roland
Thanks Zach.
On Thu, Aug 27, 2020, 11:50 PM Zach Etienne zachetie@gmail.com wrote:
Hi Bilal,
You are under-resolving the black holes. For stability you will need gridspacing of at least ~1/16 around each black hole (M/16, where M=1 for this binary), and here your gridspacing is 1.25. I would suggest starting by trying to reproduce the binary black hole gallery example: https://einsteintoolkit.org/gallery/bbh/index.html
and adding small adjustments until you get it working according to your needs.
-Zach
Prof. Zachariah Etienne Physics & Astronomy Dept. West Virginia University http://astro.phys.wvu.edu/zetienne/ http://blackholesathome.net https://blackholesathome.net
On Thu, Aug 27, 2020 at 2:38 PM Bilal Hameed bilalhameedawan962@gmail.com wrote:
Dear Users.
I am runing --parfile repos/ctthorns/CT_MultiLevel/par/twopunctures.par.
###################################### # Flesh parameters ######################################
Cactus::cctk_run_title = "Solving the Hamiltonian constraint for two punctures with CT_MultiLevel (multigrid setup, FMG) and with TwoPunctures"
Cactus::cctk_full_warnings = yes Cactus::highlight_warning_messages = no Cactus::cctk_timer_output = "off"
Cactus::terminate = "any" Cactus::cctk_initial_time = 0 Cactus::cctk_final_time = 30 Cactus::cctk_itlast = 200
############################################################################################# ActiveThorns = "Boundary CartGrid3D CoordBase Slab SymBase"
#############################################################################################
CoordBase::domainsize = "minmax"
CoordBase::xmin = -40.00 CoordBase::ymin = -40.00 CoordBase::zmin = -40.00 CoordBase::xmax = +40.00 CoordBase::ymax = +40.00 CoordBase::zmax = +40.00 CoordBase::dx = 10 CoordBase::dy = 10 CoordBase::dz = 10
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 = 1 CoordBase::boundary_shiftout_y_upper = 1 CoordBase::boundary_shiftout_z_upper = 1
CartGrid3D::type = "coordbase"
################################################################################## ActiveThorns = "IOUtil Carpet CarpetLib CarpetInterp LoopControl CarpetReduce CarpetSlab"
##################################################################################
Carpet::adaptive_stepsize = "no"
Carpet::schedule_barriers = no CarpetLib::output_bboxes = no Carpet::init_fill_timelevels = "yes"
Carpet::domain_from_coordbase = yes Carpet::max_refinement_levels = 5
driver::ghost_size = 3 Carpet::use_buffer_zones = yes
Carpet::prolongation_order_space = 5 Carpet::prolongation_order_time = 2
Carpet::convergence_level = 0
Carpet::poison_new_timelevels = yes CarpetLib::poison_new_memory = yes
Carpet::output_timers_every = 0 CarpetLib::print_timestats_every = 0 CarpetLib::print_memstats_every = 0
CarpetLib::max_allowed_memory_MB = 12288
################################################################## ActiveThorns = "CarpetRegrid2" ##################################################################
Carpet::grid_structure_filename = "grid_structure.asc" Carpet::grid_coordinates_filename = "grid_coordinates.asc"
CarpetRegrid2::num_centres = 1 CarpetRegrid2::num_levels_1 = 4 CarpetRegrid2::radius_1[1] = 50 CarpetRegrid2::radius_1[2] = 25 CarpetRegrid2::radius_1[3] = 12 #CarpetRegrid2::radius_1[4] = 50 #CarpetRegrid2::radius_1[5] = 25 #CarpetRegrid2::radius_1[6] = 12 #CarpetRegrid2::radius_1[7] = 6
################################# ActiveThorns = "AEILocalInterp" ActiveThorns = "BLAS LAPACK" ActiveThorns = "Fortran" ActiveThorns = "GenericFD" ActiveThorns = "HDF5" ActiveThorns = "InitBase" ActiveThorns = "LocalInterp" #################################
################################# ActiveThorns = "MoL Time" #################################
MoL::ODE_Method = "RK4" MoL::MoL_Intermediate_Steps = 4 MoL::MoL_Num_Scratch_Levels = 1
Time::dtfac = 0.2
################################################################################################### ActiveThorns = "CT_MultiLevel CT_Analytic ADMBase ML_BSSN ML_BSSN_Helper CoordGauge StaticConformal TmunuBase ADMCoupling ADMMacros TwoPunctures"
###################################################################################################
ADMBase::metric_timelevels = 3 ADMBase::lapse_timelevels = 3 ADMBase::shift_timelevels = 3
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"
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"
ADMBase::admbase_boundary_condition = "none"
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
#TwoPunctures::grid_setup_method = "evaluation"
TwoPunctures::TP_epsilon = 1.0e-2 TwoPunctures::TP_Tiny = 1.0e-2
TwoPunctures::keep_u_around = "yes"
CT_MultiLevel::cycle_type = "FMG cycle" CT_MultiLevel::tol = 1e-12 CT_MultiLevel::omega = 1 CT_MultiLevel::fd_order = 4 CT_MultiLevel::output_norms = "yes" CT_MultiLevel::output_walk = "yes" CT_MultiLevel::nrelsteps_down = 700 CT_MultiLevel::nrelsteps_up = 650 CT_MultiLevel::nrelsteps_bottom = 700 CT_MultiLevel::nrelsteps_top = 650 CT_MultiLevel::topMGlevel = 4 CT_MultiLevel::number_of_equations = 1 CT_MultiLevel::number_of_auxiliaries = 1 CT_MultiLevel::fill_ADM = "yes" CT_MultiLevel::fill_Aij = "Analytic Aij" CT_MultiLevel::boundary_conditions = "TwoPunctures"
CT_MultiLevel::compare_to_exact = "yes" CT_MultiLevel::exact_solution_gfname[0] = "TwoPunctures::puncture_u" CT_MultiLevel::exact_offset = -1.0
CT_MultiLevel::inipsi_gfname[0] = "CT_Analytic::testinipsi" CT_MultiLevel::cxx_gfname[0] = "CT_Analytic::testcxx" CT_MultiLevel::cyy_gfname[0] = "CT_Analytic::testcyy" CT_MultiLevel::czz_gfname[0] = "CT_Analytic::testczz" CT_MultiLevel::n0[0] = -7 CT_MultiLevel::c0_gfname[0] = "CT_MultiLevel::ct_auxiliary[0]" CT_MultiLevel::a0_gfname[0] = "CT_Analytic::testa0"
CT_MultiLevel::model = "Bowen-York" CT_Analytic::free_data = "Bowen-York" CT_Analytic::ampI = 1 CT_Analytic::ampC = 1 CT_Analytic::massa = 0.453 CT_Analytic::massb = 0.453 CT_Analytic::xa = 1.168642873 CT_Analytic::xb = -1.168642873 CT_Analytic::Pay = 0.3331917498 CT_Analytic::Pby = -0.3331917498
CT_Analytic::other_timelevels = 3
################################# ActiveThorns = "CarpetIOBasic" #################################
IOBasic::outInfo_every = 1 IOBasic::outInfo_reductions = "norm2" IOBasic::outInfo_vars = " CT_MultiLevel::psi Carpet::physical_time_per_hour "
################################# ActiveThorns = "CarpetIOScalar" #################################
IO::out_dir = $parfile
IOScalar::one_file_per_group = yes
IOScalar::outScalar_every = 0 IOScalar::outScalar_vars = " "
################################# ActiveThorns = "CarpetIOASCII" #################################
IOASCII::one_file_per_group = yes
IOASCII::output_symmetry_points = yes IOASCII::out3D_ghosts = yes
IOASCII::out0D_criterion = "divisor" IOASCII::out0D_every = 0 IOASCII::out0D_vars = " "
IOASCII::output_all_timelevels = "no" IOASCII::out1D_criterion = "divisor" IOASCII::out1D_every = 1 IOASCII::out1D_vars = " CT_MultiLevel::psi CT_MultiLevel::residual CT_MultiLevel::err ADMBase::metric ADMBase::curv ML_BSSN::ML_log_confac ML_BSSN::ML_ham TwoPunctures::puncture_u "
IOASCII::out2D_every = 0 IOASCII::out2D_criterion = "divisor" IOASCII::out2D_vars = " "
IOASCII::out3D_every = 0 IOASCII::out3D_criterion = "divisor" IOASCII::out3D_vars = " CT_MultiLevel::psi CT_MultiLevel::coeffs "
################################# Activethorns = "CarpetIOHDF5" #################################
IOHDF5::out_criterion = "divisor" IOHDF5::out_every = 1 IOHDF5::one_file_per_group = yes IOHDF5::compression_level = 0 IOHDF5::out_vars = " CT_MultiLevel::psi CT_MultiLevel::err CT_MultiLevel::residual ML_BSSN::ML_log_confac ML_BSSN::ML_ham TwoPunctures::puncture_u "
IOHDF5::checkpoint = "no"
################################# ActiveThorns = "Formaline" #################################
################################# ActiveThorns = "TimerReport" #################################
TimerReport::out_every = 512 TimerReport::out_filename = "TimerReport"
*My ouput shows -nan after first itration. *
0 0.000 | 1.0013688 | 0.0000000 2 0.250 | -nan | 1.128073e+03 4 0.500 | -nan | 1.382569e+03 6 0.750 | -nan | 1.318078e+03 8 1.000 | -nan | 1.399377e+03 10 1.250 | -nan | 1.336436e+03 12 1.500 | -nan | 1.395822e+03 14 1.750 | -nan | 1.364895e+03 16 2.000 | -nan | 1.398991e+03 18 2.250 | -nan | 1.359115e+03 20 2.500 | -nan | 1.393235e+03 I set refinement levels, radius, topMG. no error found in running.
Kindly,
Users mailing list Users@einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/users
Thanks Sir. But I am still confused about, why changing "itlast" generated that issue as it only controls "last iteration"?
On Fri, Aug 28, 2020, 12:28 AM Roland Haas rhaas@illinois.edu wrote:
Hello Bilal,
following up on the private conversation we had:
hmm.
So you may not want to change itlast.
Note that CT_MultiLevel is an elliptic solver, so you are constructing initial data.
Looking at:
repos/ctthorns/CT_MultiLevel/schedule.cc
there is only:
schedule CT_MultiLevel AT CCTK_INITIAL after CT_ScalarField_setConfRho { LANG: C OPTIONS: GLOBAL-LATE } "Main MultiLevel function"
ie only in INITIAL which is for itlast=0.
To the actually run an evolution (the black holes actually moving) you would have to read in the data produced by CT_MultiLevel into an evolution simulation. This would involve using a parameter file similar to the one that Zach pointed you to and a thorn like the ReadInterpolate thorn to read in the results data from CT_MultiLevel into the evolution simulation.
How to do that is not fully trivial, unfortunately. You would have to ask someone on the mailing list if they have an example parameter file. Possible the CT_MultiLevel author, who had used my thorn before.
Yours, Roland
Thanks Zach.
On Thu, Aug 27, 2020, 11:50 PM Zach Etienne zachetie@gmail.com wrote:
Hi Bilal,
You are under-resolving the black holes. For stability you will need gridspacing of at least ~1/16 around each black hole (M/16, where
M=1
for this binary), and here your gridspacing is 1.25. I would suggest starting by trying to reproduce the binary black hole gallery example: https://einsteintoolkit.org/gallery/bbh/index.html
and adding small adjustments until you get it working according to your needs.
-Zach
Prof. Zachariah Etienne Physics & Astronomy Dept. West Virginia University http://astro.phys.wvu.edu/zetienne/ http://blackholesathome.net https://blackholesathome.net
On Thu, Aug 27, 2020 at 2:38 PM Bilal Hameed <
bilalhameedawan962@gmail.com>
wrote:
Dear Users.
I am runing --parfile
repos/ctthorns/CT_MultiLevel/par/twopunctures.par.
###################################### # Flesh parameters ######################################
Cactus::cctk_run_title = "Solving the Hamiltonian constraint for two punctures with CT_MultiLevel (multigrid setup, FMG) and with
TwoPunctures"
Cactus::cctk_full_warnings = yes Cactus::highlight_warning_messages = no Cactus::cctk_timer_output = "off"
Cactus::terminate = "any" Cactus::cctk_initial_time = 0 Cactus::cctk_final_time = 30 Cactus::cctk_itlast = 200
#############################################################################################
ActiveThorns = "Boundary CartGrid3D CoordBase Slab SymBase"
#############################################################################################
CoordBase::domainsize = "minmax"
CoordBase::xmin = -40.00 CoordBase::ymin = -40.00 CoordBase::zmin = -40.00 CoordBase::xmax = +40.00 CoordBase::ymax = +40.00 CoordBase::zmax = +40.00 CoordBase::dx = 10 CoordBase::dy = 10 CoordBase::dz = 10
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 = 1 CoordBase::boundary_shiftout_y_upper = 1 CoordBase::boundary_shiftout_z_upper = 1
CartGrid3D::type = "coordbase"
##################################################################################
ActiveThorns = "IOUtil Carpet CarpetLib CarpetInterp LoopControl CarpetReduce CarpetSlab"
##################################################################################
Carpet::adaptive_stepsize = "no"
Carpet::schedule_barriers = no CarpetLib::output_bboxes = no Carpet::init_fill_timelevels = "yes"
Carpet::domain_from_coordbase = yes Carpet::max_refinement_levels = 5
driver::ghost_size = 3 Carpet::use_buffer_zones = yes
Carpet::prolongation_order_space = 5 Carpet::prolongation_order_time = 2
Carpet::convergence_level = 0
Carpet::poison_new_timelevels = yes CarpetLib::poison_new_memory = yes
Carpet::output_timers_every = 0 CarpetLib::print_timestats_every = 0 CarpetLib::print_memstats_every = 0
CarpetLib::max_allowed_memory_MB = 12288
################################################################## ActiveThorns = "CarpetRegrid2" ##################################################################
Carpet::grid_structure_filename = "grid_structure.asc" Carpet::grid_coordinates_filename = "grid_coordinates.asc"
CarpetRegrid2::num_centres = 1 CarpetRegrid2::num_levels_1 = 4 CarpetRegrid2::radius_1[1] = 50 CarpetRegrid2::radius_1[2] = 25 CarpetRegrid2::radius_1[3] = 12 #CarpetRegrid2::radius_1[4] = 50 #CarpetRegrid2::radius_1[5] = 25 #CarpetRegrid2::radius_1[6] = 12 #CarpetRegrid2::radius_1[7] = 6
################################# ActiveThorns = "AEILocalInterp" ActiveThorns = "BLAS LAPACK" ActiveThorns = "Fortran" ActiveThorns = "GenericFD" ActiveThorns = "HDF5" ActiveThorns = "InitBase" ActiveThorns = "LocalInterp" #################################
################################# ActiveThorns = "MoL Time" #################################
MoL::ODE_Method = "RK4" MoL::MoL_Intermediate_Steps = 4 MoL::MoL_Num_Scratch_Levels = 1
Time::dtfac = 0.2
###################################################################################################
ActiveThorns = "CT_MultiLevel CT_Analytic ADMBase ML_BSSN
ML_BSSN_Helper
CoordGauge StaticConformal TmunuBase ADMCoupling ADMMacros
TwoPunctures"
###################################################################################################
ADMBase::metric_timelevels = 3 ADMBase::lapse_timelevels = 3 ADMBase::shift_timelevels = 3
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"
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"
ADMBase::admbase_boundary_condition = "none"
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
#TwoPunctures::grid_setup_method = "evaluation"
TwoPunctures::TP_epsilon = 1.0e-2 TwoPunctures::TP_Tiny = 1.0e-2
TwoPunctures::keep_u_around = "yes"
CT_MultiLevel::cycle_type = "FMG cycle" CT_MultiLevel::tol = 1e-12 CT_MultiLevel::omega = 1 CT_MultiLevel::fd_order = 4 CT_MultiLevel::output_norms = "yes" CT_MultiLevel::output_walk = "yes" CT_MultiLevel::nrelsteps_down = 700 CT_MultiLevel::nrelsteps_up = 650 CT_MultiLevel::nrelsteps_bottom = 700 CT_MultiLevel::nrelsteps_top = 650 CT_MultiLevel::topMGlevel = 4 CT_MultiLevel::number_of_equations = 1 CT_MultiLevel::number_of_auxiliaries = 1 CT_MultiLevel::fill_ADM = "yes" CT_MultiLevel::fill_Aij = "Analytic Aij" CT_MultiLevel::boundary_conditions = "TwoPunctures"
CT_MultiLevel::compare_to_exact = "yes" CT_MultiLevel::exact_solution_gfname[0] = "TwoPunctures::puncture_u" CT_MultiLevel::exact_offset = -1.0
CT_MultiLevel::inipsi_gfname[0] = "CT_Analytic::testinipsi" CT_MultiLevel::cxx_gfname[0] = "CT_Analytic::testcxx" CT_MultiLevel::cyy_gfname[0] = "CT_Analytic::testcyy" CT_MultiLevel::czz_gfname[0] = "CT_Analytic::testczz" CT_MultiLevel::n0[0] = -7 CT_MultiLevel::c0_gfname[0] = "CT_MultiLevel::ct_auxiliary[0]" CT_MultiLevel::a0_gfname[0] = "CT_Analytic::testa0"
CT_MultiLevel::model = "Bowen-York" CT_Analytic::free_data = "Bowen-York" CT_Analytic::ampI = 1 CT_Analytic::ampC = 1 CT_Analytic::massa = 0.453 CT_Analytic::massb = 0.453 CT_Analytic::xa = 1.168642873 CT_Analytic::xb = -1.168642873 CT_Analytic::Pay = 0.3331917498 CT_Analytic::Pby = -0.3331917498
CT_Analytic::other_timelevels = 3
################################# ActiveThorns = "CarpetIOBasic" #################################
IOBasic::outInfo_every = 1 IOBasic::outInfo_reductions = "norm2" IOBasic::outInfo_vars = " CT_MultiLevel::psi Carpet::physical_time_per_hour "
################################# ActiveThorns = "CarpetIOScalar" #################################
IO::out_dir = $parfile
IOScalar::one_file_per_group = yes
IOScalar::outScalar_every = 0 IOScalar::outScalar_vars = " "
################################# ActiveThorns = "CarpetIOASCII" #################################
IOASCII::one_file_per_group = yes
IOASCII::output_symmetry_points = yes IOASCII::out3D_ghosts = yes
IOASCII::out0D_criterion = "divisor" IOASCII::out0D_every = 0 IOASCII::out0D_vars = " "
IOASCII::output_all_timelevels = "no" IOASCII::out1D_criterion = "divisor" IOASCII::out1D_every = 1 IOASCII::out1D_vars = " CT_MultiLevel::psi CT_MultiLevel::residual CT_MultiLevel::err ADMBase::metric ADMBase::curv ML_BSSN::ML_log_confac ML_BSSN::ML_ham TwoPunctures::puncture_u "
IOASCII::out2D_every = 0 IOASCII::out2D_criterion = "divisor" IOASCII::out2D_vars = " "
IOASCII::out3D_every = 0 IOASCII::out3D_criterion = "divisor" IOASCII::out3D_vars = " CT_MultiLevel::psi CT_MultiLevel::coeffs "
################################# Activethorns = "CarpetIOHDF5" #################################
IOHDF5::out_criterion = "divisor" IOHDF5::out_every = 1 IOHDF5::one_file_per_group = yes IOHDF5::compression_level = 0 IOHDF5::out_vars = " CT_MultiLevel::psi CT_MultiLevel::err CT_MultiLevel::residual ML_BSSN::ML_log_confac ML_BSSN::ML_ham TwoPunctures::puncture_u "
IOHDF5::checkpoint = "no"
################################# ActiveThorns = "Formaline" #################################
################################# ActiveThorns = "TimerReport" #################################
TimerReport::out_every = 512 TimerReport::out_filename = "TimerReport"
*My ouput shows -nan after first itration. *
0 0.000 | 1.0013688 | 0.0000000 2 0.250 | -nan | 1.128073e+03 4 0.500 | -nan | 1.382569e+03 6 0.750 | -nan | 1.318078e+03 8 1.000 | -nan | 1.399377e+03 10 1.250 | -nan | 1.336436e+03 12 1.500 | -nan | 1.395822e+03 14 1.750 | -nan | 1.364895e+03 16 2.000 | -nan | 1.398991e+03 18 2.250 | -nan | 1.359115e+03 20 2.500 | -nan | 1.393235e+03 I set refinement levels, radius, topMG. no error found in running.
Kindly,
Users mailing list Users@einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/users
-- My email is as private as my paper mail. I therefore support encrypting and signing email messages. Get my PGP key from http://pgp.mit.edu .
Hello Bilal,
Thanks Sir. But I am still confused about, why changing "itlast" generated that issue as it only controls "last iteration"?
Changing itlast to > 0 makes Cactus enter an evolution loop where it will cycle time levels bringing fresh, unitialized time levels to the top of the set of time levels. Since there is nothing that would evolve the data or re-initialize the timelevels your are left with unitialized time levels, i.e. NaNs.
Just to be explicit: iteration is not some sort of iteration in a non-linear solver. It is (more or less) the time step of a hyperbolic evolution.
If you want more multigrid (W, V, whatnot) cycles then you have to change CT_MultiLevel's parameters. I would suggest having a look at param.ccl of CT_MultiLevel to see which ones are available.
Yours, Roland
Thanks!
On Fri, Aug 28, 2020, 6:34 AM Roland Haas rhaas@illinois.edu wrote:
Hello Bilal,
Thanks Sir. But I am still confused about, why changing "itlast"
generated
that issue as it only controls "last iteration"?
Changing itlast to > 0 makes Cactus enter an evolution loop where it will cycle time levels bringing fresh, unitialized time levels to the top of the set of time levels. Since there is nothing that would evolve the data or re-initialize the timelevels your are left with unitialized time levels, i.e. NaNs.
Just to be explicit: iteration is not some sort of iteration in a non-linear solver. It is (more or less) the time step of a hyperbolic evolution.
If you want more multigrid (W, V, whatnot) cycles then you have to change CT_MultiLevel's parameters. I would suggest having a look at param.ccl of CT_MultiLevel to see which ones are available.
Yours, Roland
-- My email is as private as my paper mail. I therefore support encrypting and signing email messages. Get my PGP key from http://pgp.mit.edu .
users@lists.einsteintoolkit.org
-
Bilal Hameed -
Roland Haas -
Zach Etienne