<div dir="ltr">Hi Bilal,<div><br></div><div>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:</div><div><a href="https://einsteintoolkit.org/gallery/bbh/index.html">https://einsteintoolkit.org/gallery/bbh/index.html</a><br></div><div><br></div><div>and adding small adjustments until you get it working according to your needs.</div><div><br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div style="font-size:12.8px">-Zach</div><div style="font-size:12.8px"><br></div><span style="font-size:12.8px">*     *     *</span><br style="font-size:12.8px"><span style="font-size:12.8px">Prof. Zachariah Etienne</span></div><div dir="ltr">Physics &amp; Astronomy Dept.<br style="font-size:12.8px"><div style="font-size:12.8px">West Virginia University</div><div><a href="http://astro.phys.wvu.edu/zetienne/" target="_blank">http://astro.phys.wvu.edu/zetienne/</a></div><div><a href="http://blackholesathome.net/" target="_blank">http://blackholesathome.net</a><a href="https://blackholesathome.net" style="font-size:12.8px" target="_blank"><br></a><br></div></div></div></div></div></div></div></div></div></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Aug 27, 2020 at 2:38 PM Bilal Hameed &lt;<a href="mailto:bilalhameedawan962@gmail.com">bilalhameedawan962@gmail.com</a>&gt; wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
  

    
  
  <div>
    <p>Dear Users.</p>
    <p>I am runing --parfile
      repos/ctthorns/CT_MultiLevel/par/twopunctures.par.</p>
    <p>######################################<br>
      # Flesh parameters<br>
      ######################################<br>
      <br>
      Cactus::cctk_run_title = &quot;Solving the Hamiltonian constraint for
      two punctures with CT_MultiLevel (multigrid setup, FMG) and with
      TwoPunctures&quot;<br>
      <br>
      Cactus::cctk_full_warnings         = yes<br>
      Cactus::highlight_warning_messages = no<br>
      Cactus::cctk_timer_output          = &quot;off&quot;<br>
      <br>
      Cactus::terminate       = &quot;any&quot;<br>
      Cactus::cctk_initial_time = 0<br>
      Cactus::cctk_final_time = 30<br>
      Cactus::cctk_itlast     = 200<br>
      <br>
      <br>
#############################################################################################<br>
      ActiveThorns = &quot;Boundary CartGrid3D CoordBase Slab SymBase&quot;<br>
#############################################################################################<br>
      <br>
      CoordBase::domainsize = &quot;minmax&quot;<br>
      <br>
      CoordBase::xmin =  -40.00<br>
      CoordBase::ymin =  -40.00<br>
      CoordBase::zmin =  -40.00<br>
      CoordBase::xmax =  +40.00<br>
      CoordBase::ymax =  +40.00<br>
      CoordBase::zmax =  +40.00<br>
      CoordBase::dx   =  10<br>
      CoordBase::dy   =  10<br>
      CoordBase::dz   =  10<br>
      <br>
      CoordBase::boundary_size_x_lower     = 3<br>
      CoordBase::boundary_size_y_lower     = 3<br>
      CoordBase::boundary_size_z_lower     = 3<br>
      CoordBase::boundary_size_x_upper     = 3<br>
      CoordBase::boundary_size_y_upper     = 3<br>
      CoordBase::boundary_size_z_upper     = 3<br>
      <br>
      CoordBase::boundary_shiftout_x_lower = 1<br>
      CoordBase::boundary_shiftout_y_lower = 1<br>
      CoordBase::boundary_shiftout_z_lower = 1<br>
      CoordBase::boundary_shiftout_x_upper = 1<br>
      CoordBase::boundary_shiftout_y_upper = 1<br>
      CoordBase::boundary_shiftout_z_upper = 1<br>
      <br>
      CartGrid3D::type = &quot;coordbase&quot;<br>
      <br>
##################################################################################<br>
      ActiveThorns = &quot;IOUtil Carpet CarpetLib CarpetInterp LoopControl
      CarpetReduce CarpetSlab&quot;<br>
##################################################################################<br>
      <br>
      Carpet::adaptive_stepsize = &quot;no&quot;<br>
      <br>
      Carpet::schedule_barriers = no<br>
      CarpetLib::output_bboxes  = no<br>
      Carpet::init_fill_timelevels = &quot;yes&quot;<br>
      <br>
      Carpet::domain_from_coordbase = yes<br>
      Carpet::max_refinement_levels = 5<br>
      <br>
      driver::ghost_size       = 3<br>
      Carpet::use_buffer_zones = yes<br>
      <br>
      Carpet::prolongation_order_space = 5<br>
      Carpet::prolongation_order_time  = 2<br>
      <br>
      Carpet::convergence_level = 0<br>
      <br>
      Carpet::poison_new_timelevels = yes<br>
      CarpetLib::poison_new_memory  = yes<br>
      <br>
      Carpet::output_timers_every      = 0<br>
      CarpetLib::print_timestats_every = 0<br>
      CarpetLib::print_memstats_every  = 0<br>
      <br>
      CarpetLib::max_allowed_memory_MB = 12288<br>
      <br>
      ##################################################################<br>
      ActiveThorns = &quot;CarpetRegrid2&quot;<br>
      ##################################################################<br>
      <br>
      Carpet::grid_structure_filename = &quot;grid_structure.asc&quot;<br>
      Carpet::grid_coordinates_filename = &quot;grid_coordinates.asc&quot;<br>
      <br>
      CarpetRegrid2::num_centres  = 1<br>
      CarpetRegrid2::num_levels_1 = 4<br>
      CarpetRegrid2::radius_1[1]  = 50<br>
      CarpetRegrid2::radius_1[2]  = 25<br>
      CarpetRegrid2::radius_1[3]  = 12<br>
      #CarpetRegrid2::radius_1[4]  = 50<br>
      #CarpetRegrid2::radius_1[5]  = 25<br>
      #CarpetRegrid2::radius_1[6]  = 12<br>
      #CarpetRegrid2::radius_1[7]  =  6<br>
      <br>
      #################################<br>
      ActiveThorns = &quot;AEILocalInterp&quot;<br>
      ActiveThorns = &quot;BLAS LAPACK&quot;<br>
      ActiveThorns = &quot;Fortran&quot;<br>
      ActiveThorns = &quot;GenericFD&quot;<br>
      ActiveThorns = &quot;HDF5&quot;<br>
      ActiveThorns = &quot;InitBase&quot;<br>
      ActiveThorns = &quot;LocalInterp&quot;<br>
      #################################<br>
      <br>
      #################################<br>
      ActiveThorns = &quot;MoL Time&quot;<br>
      #################################<br>
      <br>
      MoL::ODE_Method             = &quot;RK4&quot;<br>
      MoL::MoL_Intermediate_Steps = 4<br>
      MoL::MoL_Num_Scratch_Levels = 1<br>
      <br>
      Time::dtfac   = 0.2<br>
      <br>
###################################################################################################<br>
      ActiveThorns = &quot;CT_MultiLevel CT_Analytic ADMBase ML_BSSN
      ML_BSSN_Helper CoordGauge StaticConformal TmunuBase ADMCoupling
      ADMMacros TwoPunctures&quot;<br>
###################################################################################################<br>
      <br>
      ADMBase::metric_timelevels = 3<br>
      ADMBase::lapse_timelevels = 3<br>
      ADMBase::shift_timelevels = 3<br>
      <br>
      ADMBase::metric_type = &quot;physical&quot;<br>
      <br>
      ADMBase::initial_data    = &quot;twopunctures&quot;<br>
      ADMBase::initial_lapse   = &quot;twopunctures-averaged&quot;<br>
      ADMBase::initial_shift   = &quot;zero&quot;<br>
      ADMBase::initial_dtlapse = &quot;zero&quot;<br>
      ADMBase::initial_dtshift = &quot;zero&quot;<br>
      <br>
      ADMBase::evolution_method         = &quot;ML_BSSN&quot;<br>
      ADMBase::lapse_evolution_method   = &quot;ML_BSSN&quot;<br>
      ADMBase::shift_evolution_method   = &quot;ML_BSSN&quot;<br>
      ADMBase::dtlapse_evolution_method = &quot;ML_BSSN&quot;<br>
      ADMBase::dtshift_evolution_method = &quot;ML_BSSN&quot;<br>
      <br>
      ADMBase::admbase_boundary_condition = &quot;none&quot;<br>
      <br>
      TwoPunctures::par_b          =  1.168642873<br>
      TwoPunctures::par_m_plus     =  0.453<br>
      TwoPunctures::par_m_minus    =  0.453<br>
      TwoPunctures::par_P_plus [1] = +0.3331917498<br>
      TwoPunctures::par_P_minus[1] = -0.3331917498<br>
      <br>
      #TwoPunctures::grid_setup_method = &quot;evaluation&quot;<br>
      <br>
      TwoPunctures::TP_epsilon = 1.0e-2<br>
      TwoPunctures::TP_Tiny    = 1.0e-2<br>
      <br>
      TwoPunctures::keep_u_around = &quot;yes&quot;<br>
      <br>
      CT_MultiLevel::cycle_type          = &quot;FMG cycle&quot;<br>
      CT_MultiLevel::tol                 = 1e-12<br>
      CT_MultiLevel::omega               = 1<br>
      CT_MultiLevel::fd_order            = 4<br>
      CT_MultiLevel::output_norms        = &quot;yes&quot;<br>
      CT_MultiLevel::output_walk         = &quot;yes&quot;<br>
      CT_MultiLevel::nrelsteps_down      = 700<br>
      CT_MultiLevel::nrelsteps_up        = 650<br>
      CT_MultiLevel::nrelsteps_bottom    = 700<br>
      CT_MultiLevel::nrelsteps_top       = 650<br>
      CT_MultiLevel::topMGlevel          = 4<br>
      CT_MultiLevel::number_of_equations = 1<br>
      CT_MultiLevel::number_of_auxiliaries = 1<br>
      CT_MultiLevel::fill_ADM              = &quot;yes&quot;<br>
      CT_MultiLevel::fill_Aij              = &quot;Analytic Aij&quot;<br>
      CT_MultiLevel::boundary_conditions   = &quot;TwoPunctures&quot;<br>
      <br>
      CT_MultiLevel::compare_to_exact      = &quot;yes&quot;<br>
      CT_MultiLevel::exact_solution_gfname[0] =
      &quot;TwoPunctures::puncture_u&quot;<br>
      CT_MultiLevel::exact_offset             = -1.0<br>
      <br>
      CT_MultiLevel::inipsi_gfname[0]    = &quot;CT_Analytic::testinipsi&quot;<br>
      CT_MultiLevel::cxx_gfname[0]       = &quot;CT_Analytic::testcxx&quot;<br>
      CT_MultiLevel::cyy_gfname[0]       = &quot;CT_Analytic::testcyy&quot;<br>
      CT_MultiLevel::czz_gfname[0]       = &quot;CT_Analytic::testczz&quot;<br>
      CT_MultiLevel::n0[0]               = -7<br>
      CT_MultiLevel::c0_gfname[0]        =
      &quot;CT_MultiLevel::ct_auxiliary[0]&quot;<br>
      CT_MultiLevel::a0_gfname[0]        = &quot;CT_Analytic::testa0&quot;<br>
      <br>
      CT_MultiLevel::model               = &quot;Bowen-York&quot;<br>
      CT_Analytic::free_data             = &quot;Bowen-York&quot;<br>
      CT_Analytic::ampI                  = 1<br>
      CT_Analytic::ampC                  = 1<br>
      CT_Analytic::massa                 = 0.453<br>
      CT_Analytic::massb                 = 0.453<br>
      CT_Analytic::xa                    = 1.168642873<br>
      CT_Analytic::xb                    = -1.168642873<br>
      CT_Analytic::Pay                   = 0.3331917498<br>
      CT_Analytic::Pby                   = -0.3331917498<br>
      <br>
      CT_Analytic::other_timelevels      = 3<br>
      <br>
      #################################<br>
      ActiveThorns = &quot;CarpetIOBasic&quot;<br>
      #################################<br>
      <br>
      IOBasic::outInfo_every      = 1<br>
      IOBasic::outInfo_reductions = &quot;norm2&quot;<br>
      IOBasic::outInfo_vars       = &quot;<br>
          CT_MultiLevel::psi<br>
              Carpet::physical_time_per_hour<br>
      &quot;<br>
      <br>
      #################################<br>
      ActiveThorns = &quot;CarpetIOScalar&quot;<br>
      #################################<br>
      <br>
      IO::out_dir = $parfile<br>
      <br>
      IOScalar::one_file_per_group = yes<br>
      <br>
      IOScalar::outScalar_every = 0<br>
      IOScalar::outScalar_vars  = &quot;<br>
      &quot;<br>
      <br>
      #################################<br>
      ActiveThorns = &quot;CarpetIOASCII&quot;<br>
      #################################<br>
      <br>
      IOASCII::one_file_per_group = yes<br>
      <br>
      IOASCII::output_symmetry_points = yes<br>
      IOASCII::out3D_ghosts           = yes<br>
      <br>
      IOASCII::out0D_criterion = &quot;divisor&quot;<br>
      IOASCII::out0D_every = 0<br>
      IOASCII::out0D_vars  = &quot;<br>
      &quot;<br>
      <br>
      IOASCII::output_all_timelevels = &quot;no&quot;<br>
      IOASCII::out1D_criterion = &quot;divisor&quot;<br>
      IOASCII::out1D_every = 1<br>
      IOASCII::out1D_vars  = &quot;<br>
          CT_MultiLevel::psi<br>
          CT_MultiLevel::residual<br>
          CT_MultiLevel::err<br>
              ADMBase::metric<br>
              ADMBase::curv<br>
          ML_BSSN::ML_log_confac<br>
          ML_BSSN::ML_ham<br>
          TwoPunctures::puncture_u<br>
      &quot;<br>
      <br>
      IOASCII::out2D_every = 0<br>
      IOASCII::out2D_criterion = &quot;divisor&quot;<br>
      IOASCII::out2D_vars  = &quot;<br>
      &quot;<br>
      <br>
      IOASCII::out3D_every = 0<br>
      IOASCII::out3D_criterion = &quot;divisor&quot;<br>
      IOASCII::out3D_vars  = &quot;<br>
          CT_MultiLevel::psi<br>
          CT_MultiLevel::coeffs<br>
      &quot;<br>
      <br>
      #################################<br>
      Activethorns = &quot;CarpetIOHDF5&quot;<br>
      #################################<br>
      <br>
      IOHDF5::out_criterion = &quot;divisor&quot;<br>
      IOHDF5::out_every          = 1<br>
      IOHDF5::one_file_per_group = yes<br>
      IOHDF5::compression_level  = 0<br>
      IOHDF5::out_vars  = &quot;<br>
          CT_MultiLevel::psi<br>
          CT_MultiLevel::err<br>
          CT_MultiLevel::residual<br>
          ML_BSSN::ML_log_confac<br>
          ML_BSSN::ML_ham<br>
          TwoPunctures::puncture_u<br>
      &quot;<br>
      <br>
      IOHDF5::checkpoint          = &quot;no&quot;<br>
      <br>
      #################################<br>
      ActiveThorns = &quot;Formaline&quot;<br>
      #################################<br>
      <br>
      #################################<br>
      ActiveThorns = &quot;TimerReport&quot;<br>
      #################################<br>
      <br>
      TimerReport::out_every    = 512<br>
      TimerReport::out_filename = &quot;TimerReport&quot;</p>
    <p><b><font color="#000345">My ouput shows -nan after first
          itration.</font> </b></p>
    <pre>   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, 
</pre>
  </div>

_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org" target="_blank">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
</blockquote></div>