[Users] -nan output.

Bilal Hameed bilalhameedawan962 at gmail.com
Thu Aug 27 20:05:28 CDT 2020


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 at 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 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 --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/users/attachments/20200828/d4df24b4/attachment-0001.html 


More information about the Users mailing list