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