[Users] -nan output.
Bilal Hameed
bilalhameedawan962 at gmail.com
Thu Aug 27 13:52:18 CDT 2020
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
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/users/attachments/20200827/6835c28e/attachment-0001.html
More information about the Users
mailing list