[Users] -nan output.
Roland Haas
rhaas at illinois.edu
Thu Aug 27 14:28:30 CDT 2020
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 at 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 at 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 at 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 .
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 833 bytes
Desc: OpenPGP digital signature
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20200827/ffd57c8d/attachment.bin
More information about the Users
mailing list