Simulation name: nsns_09_063_illinois16 Running simulation nsns_09_063_illinois16 Preparing: Checking: /scratch1/scratchdirs/chlin/simulations/nsns_09_063_illinois16/output-0000-active nid02984 Thu Jul 19 11:24:40 PDT 2018 Environment: Starting: INFO (Cactus): Increasing logging level from 0 to 3 INFO (Cactus): Increasing logging level from 0 to 3 -------------------------------------------------------------------------------- 10 1 0101 ************************ 01 1010 10 The Cactus Code V4.2.3 1010 1101 011 www.cactuscode.org 1001 100101 ************************ 00010101 100011 (c) Copyright The Authors 0100 GNU Licensed. No Warranty 0101 -------------------------------------------------------------------------------- Cactus version: 4.2.3 Compile date: Jul 19 2018 (11:14:34) Run date: Jul 19 2018 (11:24:44-0700) Run host: nid02984 (pid=34463) Working directory: /scratch1/scratchdirs/chlin/simulations/nsns_09_063_illinois16/output-0000 Executable: /scratch1/scratchdirs/chlin/simulations/nsns_09_063_illinois16/SIMFACTORY/exe/cactus_sim Parameter file: /scratch1/scratchdirs/chlin/simulations/nsns_09_063_illinois16/output-0000/nsns_test.par -------------------------------------------------------------------------------- Activating thorn Cactus...Success -> active implementation Cactus Activation requested for --->ADMBase ADMCoupling ADMMacros AEILocalInterp AHFinderDirect Boundary Carpet CarpetInterp CarpetIOASCII CarpetIOBasic CarpetIOHDF5 CarpetIOScalar CarpetLib CarpetReduce CarpetRegrid2 CarpetSlab CartGrid3D CoordBase CoordGauge CycleClock Dissipation EOS_Base EOS_Omni Fortran GenericFD GSL HDF5 HydroBase InitBase IOUtil LocalInterp LocalReduce LoopControl ML_ADMConstraints ML_BSSN ML_BSSN_Helper MoL MPI NaNChecker NewRad Slab SpaceMask SphericalSurface StaticConformal SymBase SystemStatistics TerminationTrigger Time Timers TmunuBase Vectors zlib Meudon_Bin_NS IllinoisGRMHD ID_converter_ILGRMHD Lorene SummationByParts VolumeIntegrals_GRMHD particle_tracerET smallbPoynET convert_to_HydroBase Seed_Magnetic_Fields_BNS <--- Thorn LoopControl requests automatic activation of hwloc Thorn Lorene requests automatic activation of OpenBLAS Thorn hwloc requests automatic activation of pciutils Activating thorn ADMBase...Success -> active implementation ADMBase Activating thorn ADMCoupling...Success -> active implementation ADMCoupling Activating thorn ADMMacros...Success -> active implementation ADMMacros Activating thorn AEILocalInterp...Success -> active implementation AEILocalInterp Activating thorn AHFinderDirect...Success -> active implementation AHFinderDirect Activating thorn Boundary...Success -> active implementation boundary Activating thorn Carpet...Success -> active implementation Driver Activating thorn CarpetInterp...Success -> active implementation interp Activating thorn CarpetIOASCII...Success -> active implementation IOASCII Activating thorn CarpetIOBasic...Success -> active implementation IOBasic Activating thorn CarpetIOHDF5...Success -> active implementation IOHDF5 Activating thorn CarpetIOScalar...Success -> active implementation IOScalar Activating thorn CarpetLib...Success -> active implementation CarpetLib Activating thorn CarpetReduce...Success -> active implementation reduce Activating thorn CarpetRegrid2...Success -> active implementation CarpetRegrid2 Activating thorn CarpetSlab...Success -> active implementation Hyperslab Activating thorn CartGrid3D...Success -> active implementation grid Activating thorn convert_to_HydroBase...Success -> active implementation Convert_to_HydroBase Activating thorn CoordBase...Success -> active implementation CoordBase Activating thorn CoordGauge...Success -> active implementation CoordGauge Activating thorn CycleClock...Success -> active implementation CycleClock Activating thorn Dissipation...Success -> active implementation Dissipation Activating thorn EOS_Base...Success -> active implementation EOS_Base Activating thorn EOS_Omni...Success -> active implementation EOS_Omni Activating thorn Fortran...Success -> active implementation Fortran Activating thorn GenericFD...Success -> active implementation GenericFD Activating thorn GSL...Success -> active implementation GSL Activating thorn HDF5...Success -> active implementation HDF5 Activating thorn hwloc...Success -> active implementation hwloc Activating thorn HydroBase...Success -> active implementation HydroBase Activating thorn ID_converter_ILGRMHD...Success -> active implementation ID_converter_ILGRMHD Activating thorn IllinoisGRMHD...Success -> active implementation IllinoisGRMHD Activating thorn InitBase...Success -> active implementation InitBase Activating thorn IOUtil...Success -> active implementation IO Activating thorn LocalInterp...Success -> active implementation LocalInterp Activating thorn LocalReduce...Success -> active implementation LocalReduce Activating thorn LoopControl...Success -> active implementation LoopControl Activating thorn Lorene...Success -> active implementation LORENE Activating thorn Meudon_Bin_NS...Success -> active implementation Meudon_Bin_NS Activating thorn ML_ADMConstraints...Success -> active implementation ML_ADMConstraints Activating thorn ML_BSSN...Success -> active implementation ML_BSSN Activating thorn ML_BSSN_Helper...Success -> active implementation ML_BSSN_Helper Activating thorn MoL...Success -> active implementation MethodOfLines Activating thorn MPI...Success -> active implementation MPI Activating thorn NaNChecker...Success -> active implementation NaNChecker Activating thorn NewRad...Success -> active implementation NewRad Activating thorn OpenBLAS...Success -> active implementation OpenBLAS Activating thorn particle_tracerET...Success -> active implementation particle_tracerET Activating thorn pciutils...Success -> active implementation pciutils Activating thorn Seed_Magnetic_Fields_BNS...Success -> active implementation Seed_Magnetic_Fields_Privt Activating thorn Slab...Success -> active implementation Slab Activating thorn smallbPoynET...Success -> active implementation smallbPoynET Activating thorn SpaceMask...Success -> active implementation SpaceMask Activating thorn SphericalSurface...Success -> active implementation SphericalSurface Activating thorn StaticConformal...Success -> active implementation StaticConformal Activating thorn SummationByParts...Success -> active implementation SummationByParts Activating thorn SymBase...Success -> active implementation SymBase Activating thorn SystemStatistics...Success -> active implementation SystemStatistics Activating thorn TerminationTrigger...Success -> active implementation TerminationTrigger Activating thorn Time...Success -> active implementation time Activating thorn Timers...Success -> active implementation Timers Activating thorn TmunuBase...Success -> active implementation TmunuBase Activating thorn Vectors...Success -> active implementation Vectors Activating thorn VolumeIntegrals_GRMHD...Success -> active implementation VolumeIntegrals_GRMHD Activating thorn zlib...Success -> active implementation zlib WARNING level 1 from host nid02984 process 0 while executing schedule bin (none), routine (no thorn)::(no routine) in thorn IOUtil, file nsns_test.par:393: -> Unknown add operator: - Interpreting as literal string with value 'ABEoutput_directory-nsns_test' WARNING level 1 from host nid02984 process 0 while executing schedule bin (none), routine (no thorn)::(no routine) in thorn IOUtil, file nsns_test.par:397: -> Unknown add operator: - Interpreting as literal string with value 'ABEoutput_directory-nsns_test' WARNING level 1 from host nid02984 process 0 while executing schedule bin (none), routine (no thorn)::(no routine) in thorn IOUtil, file nsns_test.par:401: -> Unknown add operator: - Interpreting as literal string with value 'ABEoutput_directory-nsns_test' INFO (IOUtil): Recovery directory 'ABEoutput_directory-nsns_test' doesn't exist -------------------------------------------------------------------------------- if (recover initial data) Recover parameters endif Startup routines [CCTK_STARTUP] Carpet::MultiModel_Startup: Multi-model Startup routine CycleClock::CycleClock_Setup: Set up CycleClock LoopControl::lc_setup: Set up LoopControl ML_BSSN_Helper::ML_BSSN_SetGroupTags: [meta] Set checkpointing and prolongation group tags Timers::Timer_Startup: Prepare hierarchical timers Carpet::Driver_Startup: Startup routine IOUtil::IOUtil_Startup: Startup routine CarpetInterp::CarpetInterpStartup: Startup routine CarpetReduce::CarpetReduceStartup: Startup routine CartGrid3D::SymmetryStartup: Register GH Extension for GridSymmetry CoordBase::CoordBase_Startup: Register a GH extension to store the coordinate system handles AEILocalInterp::AEILocalInterp_U_Startup: register CCTK_InterpLocalUniform() interpolation operators HydroBase::HydroBase_StartUp: Startup banner CarpetIOASCII::CarpetIOASCIIStartup: [global] Startup routine LocalInterp::LocalInterp_Startup: register LocalInterp's interpolation operators LocalReduce::LocalReduce_Startup: Startup routine CarpetIOBasic::CarpetIOBasicStartup: [global] Startup routine ML_ADMConstraints::ML_ADMConstraints_Startup: [meta] create banner ML_BSSN::ML_BSSN_Startup: [meta] create banner ML_BSSN_Helper::ML_BSSN_RegisterSlicing: [meta] Register slicing CarpetIOHDF5::CarpetIOHDF5_Startup: Startup routine MoL::MoL_Startup: Startup banner Slab::Slab_InitMPIDatatypes: Create MPI datatypes for complex variables in C SymBase::SymBase_Startup: Register GH Extension for SymBase CarpetIOScalar::CarpetIOScalarStartup: [global] Startup routine Vectors::Vectors_Startup: Print startup message Startup routines which need an existing grid hierarchy [CCTK_WRAGH] ADMBase::Einstein_InitSymBound: [global] Set up GF symmetries Boundary::Boundary_RegisterBCs: [global] Register boundary conditions that this thorn provides CarpetRegrid2::CarpetRegrid2_Initialise: [global] Initialise locations of refined regions CartGrid3D::RegisterCartGrid3DCoords: [meta] Register coordinates for the Cartesian grid CoordGauge::Einstein_ActivateSlicing: Initialize slicing, setup priorities for mixed slicings CoordGauge::Einstein_SetNextSlicing: Identify the slicing for the next iteration EOS_Omni::EOS_Omni_Startup: [global] Set up conversion factors and other fun stuff ML_BSSN_Helper::ML_BSSN_ParamCompat: [meta] Handle parameter backward compatibility MoL::MoL_SetupIndexArrays: Set up the MoL bookkeeping index arrays MoL::MoL_SetScheduleStatus: [global] Set the flag so it is ok to register with MoL TmunuBase::TmunuBase_SetStressEnergyState: [global] Set the stress_energy_state variable GROUP MoL_Register: The group where physics thorns register variables with MoL IllinoisGRMHD::IllinoisGRMHD_RegisterVars: [meta] Register evolved, rhs variables in IllinoisGRMHD for MoL ML_ADMConstraints::ML_ADMConstraints_RegisterVars: [meta] Register Variables for MoL ML_BSSN::ML_BSSN_RegisterVars: [meta] Register Variables for MoL ML_BSSN_Helper::ML_BSSN_RegisterConstrained: [meta] Register ADMBase variables as constrained Slab::Slab_InitTimers: Initialise timers SpaceMask::MaskSym: [global] Set grid symmetries for mask SpaceMask::MaskSym_emask: [global] Set grid symmetries for emask (compatibility mode) GROUP SymBase_Wrapper: Wrapper group for SymBase GROUP SymmetryRegister: Register your symmetries here CartGrid3D::RegisterSymmetryBoundaries: [meta] Register symmetry boundaries ML_ADMConstraints::ML_ADMConstraints_RegisterSymmetries: [meta] register symmetries ML_BSSN::ML_BSSN_RegisterSymmetries: [meta] register symmetries SymBase::SymBase_Statistics: Print symmetry boundary face descriptions TerminationTrigger::TerminationTrigger_StartTimer: Start timer TerminationTrigger::TerminationTrigger_CreateFile: Create termination file MoL::MoL_ReportNumberVariables: [meta] Report how many of each type of variable there are Parameter checking routines [CCTK_PARAMCHECK] ADMBase::ADMBase_ParamCheck: [global] Check consistency of parameters Boundary::Boundary_Check: Check dimension of grid variables Carpet::CarpetParamCheck: Parameter checking routine CarpetLib::CarpetLib_test_prolongate_3d_rf2: [global] Test prolongation operators CarpetRegrid2::CarpetRegrid2_ParamCheck: Check parameters CartGrid3D::ParamCheck_CartGrid3D: Check coordinates for CartGrid3D Dissipation::dissipation_paramcheck: Check dissipation parameters for consistency Fortran::CheckFortranParameters: Test whether Fortran parameters work correctly ML_BSSN_Helper::ML_BSSN_ParamCheck: [meta] Check parameters MoL::MoL_ParamCheck: Basic parameter checking Vectors::Vectors_Test: Run correctness tests. Initialisation if (NOT (recover initial data AND recovery_mode is 'strict')) [CCTK_PREREGRIDINITIAL] Set up grid hierarchy [CCTK_POSTREGRIDINITIAL] SpaceMask::MaskZero: Initialise mask to zero CartGrid3D::SpatialCoordinates: Set Coordinates after regridding GROUP MaskBase_SetupMask: Set up the weight function GROUP MaskBase_SetupMaskAll: Set up the weight function CarpetReduce::MaskBase_AllocateMask: [global] Allocate the weight function CarpetReduce::MaskBase_InitMask: [global] [loop-local] Initialise the weight function GROUP SetupIMaskInternal: Set up the integer weight function (schedule other routines in here) CarpetReduce::CoordBase_SetupMask: [global] [loop-local] Set up the outer boundaries of the weight function CarpetReduce::CarpetMaskSetup: [global] [loop-singlemap] Set up the weight function for the restriction regions GROUP SetupIMask: Set up the integer weight function (schedule other routines in here) CarpetReduce::MaskBase_SetMask: [global] [loop-local] Set the weight function GROUP SetupMask: Set up the real weight function (schedule other routines in here) CarpetReduce::MaskBase_TestMask: [global] Test the weight function Dissipation::setup_epsdis: Setup spatially varying dissipation GROUP MoL_PseudoEvolutionBoundaries: Apply boundary conditions to pseudo-evolved quantities GROUP ML_ADMConstraints_evaluate_bc_group: ML_ADMConstraints_evaluate ML_ADMConstraints::ML_ADMConstraints_evaluate_SelectBCs: [level] ML_ADMConstraints_evaluate_SelectBCs GROUP ML_ADMConstraints_evaluate_ApplyBCs: Apply BCs for groups set in ML_ADMConstraints_evaluate GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ML_BSSN_ConstraintsEverywhere_bc_group: ML_BSSN_ConstraintsEverywhere ML_BSSN::ML_BSSN_ConstraintsEverywhere_SelectBCs: [level] ML_BSSN_ConstraintsEverywhere_SelectBCs GROUP ML_BSSN_ConstraintsEverywhere_ApplyBCs: Apply BCs for groups set in ML_BSSN_ConstraintsEverywhere GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ML_BSSN_ConstraintsInterior_bc_group: ML_BSSN_ConstraintsInterior ML_BSSN::ML_BSSN_ConstraintsInterior_SelectBCs: [level] ML_BSSN_ConstraintsInterior_SelectBCs GROUP ML_BSSN_ConstraintsInterior_ApplyBCs: Apply BCs for groups set in ML_BSSN_ConstraintsInterior GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions SpaceMask::MaskOne: Set mask to one AHFinderDirect::AHFinderDirect_maybe_do_masks: set mask(s) based on apparent horizon position(s) [CCTK_BASEGRID] ADMBase::ADMBase_SetShiftStateOn: Set the shift_state variable to 1 ADMBase::ADMBase_SetDtLapseStateOn: Set the dtlapse_state variable to 1 ADMBase::ADMBase_SetDtShiftStateOn: Set the dtshift_state variable to 1 ADMMacros::ADMMacros_SetLocalSpatialOrder: Initialize the local_spatial_order CartGrid3D::SpatialSpacings: Set up ranges for spatial 3D Cartesian coordinates (on all grids) CartGrid3D::SpatialCoordinates: Set up spatial 3D Cartesian coordinates on the GH SphericalSurface::SphericalSurface_SetupRes: [global] [loop-local] Set surface resolution automatically AHFinderDirect::AHFinderDirect_setup: [global] setup data structures Dissipation::dissipation_basegrid: Ensure that there are enough ghost zones IllinoisGRMHD::IllinoisGRMHD_InitSymBound: Schedule symmetries ML_ADMConstraints::ML_ADMConstraints_CheckBoundaries: [meta] check boundaries treatment ML_BSSN::ML_BSSN_CheckBoundaries: [meta] check boundaries treatment NaNChecker::NaNChecker_ResetCounter: [global] Reset the NaNChecker::NaNsFound counter SpaceMask::MaskZero: Initialise mask to zero SpaceMask::MaskOne: Set old style mask to one SphericalSurface::SphericalSurface_Setup: [global] Calculate surface coordinate descriptors GROUP MaskBase_SetupMask: Set up the weight function GROUP MaskBase_SetupMaskAll: Set up the weight function CarpetReduce::MaskBase_AllocateMask: [global] Allocate the weight function CarpetReduce::MaskBase_InitMask: [global] [loop-local] Initialise the weight function GROUP SetupIMaskInternal: Set up the integer weight function (schedule other routines in here) CarpetReduce::CoordBase_SetupMask: [global] [loop-local] Set up the outer boundaries of the weight function CarpetReduce::CarpetMaskSetup: [global] [loop-singlemap] Set up the weight function for the restriction regions GROUP SetupIMask: Set up the integer weight function (schedule other routines in here) CarpetReduce::MaskBase_SetMask: [global] [loop-local] Set the weight function GROUP SetupMask: Set up the real weight function (schedule other routines in here) CarpetReduce::MaskBase_TestMask: [global] Test the weight function SphericalSurface::SphericalSurface_Set: [global] Set surface radii to be used for initial setup in other thorns GROUP SphericalSurface_HasBeenSet: Set the spherical surfaces before this group, and use it afterwards SphericalSurface::SphericalSurface_CheckState: [global] Test the state of the spherical surfaces SummationByParts::SBP_SetNormMask: Setup the mask for the calculation of the norm SymBase::SymBase_Check: Check whether the driver set up the grid consistently Time::Time_Initialise: [global] Initialise Time variables Time::TemporalSpacings: [singlemap] Set timestep based on Courant condition (courant_static) [CCTK_INITIAL] StaticConformal::StaticConformal_InitialiseState: Set the conformal_state variable to 0 GROUP ADMBase_InitialData: Schedule group for calculating ADM initial data GROUP ADMBase_InitialGauge: Schedule group for the ADM initial gauge condition ADMBase::ADMBase_ShiftZero: Set the shift to 0 at all points ADMBase::ADMBase_DtLapseZero: Set the dtlapse to 0 at all points ADMBase::ADMBase_DtShiftZero: Set the dtshift to 0 at all points CarpetIOASCII::CarpetIOASCIIInit: [global] Initialisation routine CarpetIOBasic::CarpetIOBasicInit: [global] Initialisation routine CarpetIOHDF5::CarpetIOHDF5_Init: [global] Initialisation routine CarpetIOScalar::CarpetIOScalarInit: [global] Initialisation routine GROUP HydroBase_Initial: HydroBase initial data group HydroBase::HydroBase_Zero: Set up vacuum hydro initial data HydroBase::HydroBase_Avec_zero: Set vector potential to 0 HydroBase::HydroBase_Aphi_zero: Set vector potential Phi to 0 Meudon_Bin_NS::Meudon_Bin_NS_initialise: Set up binary neutron star initial data Seed_Magnetic_Fields_BNS::Seed_Magnetic_Fields_Privt: Set up binary neutron star seed magnetic fields. GROUP ADMBase_PostInitial: Schedule group for modifying the ADM initial data, such as e.g. adding noise GROUP HydroBase_Prim2ConInitial: Recover the conservative variables from the primitive variables GROUP IllinoisGRMHD_ID_Converter: Translate ET-generated, HydroBase-compatible initial data and convert into variables used by IllinoisGRMHD ID_converter_ILGRMHD::first_initialdata: [local] Convert HydroBase initial data (ID) to ID that IllinoisGRMHD can read. ID_converter_ILGRMHD::third_initialdata: Schedule symmetries -- Actually just a placeholder function to ensure prolongation / processor syncs are done BEFORE the primitives solver. ID_converter_ILGRMHD::fourth_initialdata: Compute B and B_stagger from A ID_converter_ILGRMHD::fifth_initialdata: Compute primitive variables from conservatives. This is non-trivial, requiring a Newton-Raphson root-finder. ML_BSSN::ML_BSSN_InitialADMBase1Everywhere: ML_BSSN_InitialADMBase1Everywhere ML_BSSN::ML_BSSN_InitialADMBase2Interior: ML_BSSN_InitialADMBase2Interior ML_BSSN::ML_BSSN_InitialADMBase2BoundaryScalar: ML_BSSN_InitialADMBase2BoundaryScalar ML_BSSN_Helper::ML_BSSN_ExtrapolateGammas: Extrapolate Gammas and time derivatives of lapse and shift MoL::MoL_StartLoop: [level] Initialise the step size control GROUP SetTmunu: Calculate the stress-energy tensor TmunuBase::TmunuBase_ZeroTmunu: Initialise the stress-energy tensor to zero GROUP AddToTmunu: Add to the stress-energy tensor here IllinoisGRMHD::IllinoisGRMHD_conserv_to_prims: Compute primitive variables from conservatives. This is non-trivial, requiring a Newton-Raphson root-finder. IllinoisGRMHD::IllinoisGRMHD_outer_boundaries_on_P_rho_b_vx_vy_vz: Apply outflow-only, flat BCs on {P,rho_b,vx,vy,vz}. Outflow only == velocities pointed inward from outer boundary are set to zero. Convert_to_HydroBase::Convert_to_HydroBase: Convert IllinoisGRMHD-native variables to HydroBase VolumeIntegrals_GRMHD::VI_GRMHD_file_output_routine_Startup: Create directory for VI GRMHD file output. VolumeIntegrals_GRMHD::VI_GRMHD_InitializeIntegralCounterToZero: [global] Initialize IntegralCounter variable to zero particle_tracerET::particle_tracerET_file_output_routine_Startup: Create directory for file output. [CCTK_POSTINITIAL] CarpetIOHDF5::CarpetIOHDF5_CloseFiles: [global] Close all filereader input files GROUP MoL_PostStepModify: The group for physics thorns to schedule enforcing constraints ML_BSSN::ML_BSSN_EnforceEverywhere: ML_BSSN_EnforceEverywhere GROUP MoL_PostStep: Ensure that everything is correct after the initial data have been set up ML_BSSN::ML_BSSN_SelectBoundConds: [level] select boundary conditions GROUP ML_BSSN_ApplyBCs: Apply boundary conditions controlled by thorn Boundary GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions ML_ADMConstraints::ML_ADMConstraints_SelectBoundConds: [level] select boundary conditions GROUP ML_ADMConstraints_ApplyBCs: Apply boundary conditions controlled by thorn Boundary GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions ML_BSSN::ML_BSSN_ADMBaseInterior: ML_BSSN_ADMBaseInterior ML_BSSN::ML_BSSN_ADMBaseBoundaryScalar: ML_BSSN_ADMBaseBoundaryScalar ML_BSSN::ML_BSSN_ADMBaseEverywhere: ML_BSSN_ADMBaseEverywhere ML_BSSN_Helper::ML_BSSN_ADMBase_SelectBCs: [level] Select boundary conditions for ADMBase variables GROUP ML_BSSN_ADMBase_ApplyBCs: Apply boundary conditions to ADMBase variables GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ADMBase_SetADMVars: Set the ADM variables before this group, and use them afterwards GROUP HydroBase_PostStep: Post step tasks for hydro thorns GROUP HydroBase_Boundaries: HydroBase-internal Boundary conditions group GROUP HydroBase_Select_Boundaries: Group to schedule the boundary condition functions GROUP HydroBase_ApplyBCs: Apply the boundary conditions of HydroBase GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions IllinoisGRMHD::Boundary_SYNCs: Schedule symmetries -- Actually just a placeholder function to ensure prolongations / processor syncs are done BEFORE outer boundaries are updated. IllinoisGRMHD::IllinoisGRMHD_outer_boundaries_on_A_mu: Apply linear extrapolation BCs on A_{mu}, so that BCs are flat on B^i. IllinoisGRMHD::IllinoisGRMHD_compute_B_and_Bstagger_from_A: Compute B and B_stagger from A, SYNC: grmhd_primitives_Bi,grmhd_primitives_Bi_stagger GROUP HydroBase_Con2Prim: Convert from conservative to primitive variables GROUP SetTmunu: Group for calculating the stress-energy tensor TmunuBase::TmunuBase_ZeroTmunu: Initialise the stress-energy tensor to zero GROUP AddToTmunu: Add to the stress-energy tensor here IllinoisGRMHD::IllinoisGRMHD_conserv_to_prims: Compute primitive variables from conservatives. This is non-trivial, requiring a Newton-Raphson root-finder. IllinoisGRMHD::IllinoisGRMHD_outer_boundaries_on_P_rho_b_vx_vy_vz: Apply outflow-only, flat BCs on {P,rho_b,vx,vy,vz}. Outflow only == velocities pointed inward from outer boundary are set to zero. GROUP MoL_PseudoEvolution: Calculate pseudo-evolved quantities GROUP ADMBase_SetADMVars: Set the ADM variables before this group, and use them afterwards GROUP ML_ADMConstraints_evaluate_group: ML_ADMConstraints_evaluate ML_ADMConstraints::ML_ADMConstraints_evaluate: ML_ADMConstraints_evaluate GROUP ML_ADMConstraints_evaluate_bc_group: ML_ADMConstraints_evaluate ML_ADMConstraints::ML_ADMConstraints_evaluate_SelectBCs: [level] ML_ADMConstraints_evaluate_SelectBCs GROUP ML_ADMConstraints_evaluate_ApplyBCs: Apply BCs for groups set in ML_ADMConstraints_evaluate GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ML_BSSN_ConstraintsEverywhere_group: ML_BSSN_ConstraintsEverywhere ML_BSSN::ML_BSSN_ConstraintsEverywhere: ML_BSSN_ConstraintsEverywhere GROUP ML_BSSN_ConstraintsEverywhere_bc_group: ML_BSSN_ConstraintsEverywhere ML_BSSN::ML_BSSN_ConstraintsEverywhere_SelectBCs: [level] ML_BSSN_ConstraintsEverywhere_SelectBCs GROUP ML_BSSN_ConstraintsEverywhere_ApplyBCs: Apply BCs for groups set in ML_BSSN_ConstraintsEverywhere GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ML_BSSN_ConstraintsInterior_group: ML_BSSN_ConstraintsInterior ML_BSSN::ML_BSSN_ConstraintsInterior: ML_BSSN_ConstraintsInterior GROUP ML_BSSN_ConstraintsInterior_bc_group: ML_BSSN_ConstraintsInterior ML_BSSN::ML_BSSN_ConstraintsInterior_SelectBCs: [level] ML_BSSN_ConstraintsInterior_SelectBCs GROUP ML_BSSN_ConstraintsInterior_ApplyBCs: Apply BCs for groups set in ML_BSSN_ConstraintsInterior GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions SummationByParts::SBP_CheckGridSizes: Check grid sizes and ghost zones Initialise finer grids recursively Restrict from finer grids [CCTK_POSTRESTRICTINITIAL] GROUP MoL_PostStep: Ensure that everything is correct after restriction ML_BSSN::ML_BSSN_SelectBoundConds: [level] select boundary conditions GROUP ML_BSSN_ApplyBCs: Apply boundary conditions controlled by thorn Boundary GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions ML_ADMConstraints::ML_ADMConstraints_SelectBoundConds: [level] select boundary conditions GROUP ML_ADMConstraints_ApplyBCs: Apply boundary conditions controlled by thorn Boundary GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions ML_BSSN::ML_BSSN_ADMBaseInterior: ML_BSSN_ADMBaseInterior ML_BSSN::ML_BSSN_ADMBaseBoundaryScalar: ML_BSSN_ADMBaseBoundaryScalar ML_BSSN::ML_BSSN_ADMBaseEverywhere: ML_BSSN_ADMBaseEverywhere ML_BSSN_Helper::ML_BSSN_ADMBase_SelectBCs: [level] Select boundary conditions for ADMBase variables GROUP ML_BSSN_ADMBase_ApplyBCs: Apply boundary conditions to ADMBase variables GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ADMBase_SetADMVars: Set the ADM variables before this group, and use them afterwards GROUP HydroBase_PostStep: Post step tasks for hydro thorns GROUP HydroBase_Boundaries: HydroBase-internal Boundary conditions group GROUP HydroBase_Select_Boundaries: Group to schedule the boundary condition functions GROUP HydroBase_ApplyBCs: Apply the boundary conditions of HydroBase GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions IllinoisGRMHD::Boundary_SYNCs: Schedule symmetries -- Actually just a placeholder function to ensure prolongations / processor syncs are done BEFORE outer boundaries are updated. IllinoisGRMHD::IllinoisGRMHD_outer_boundaries_on_A_mu: Apply linear extrapolation BCs on A_{mu}, so that BCs are flat on B^i. IllinoisGRMHD::IllinoisGRMHD_compute_B_and_Bstagger_from_A: Compute B and B_stagger from A, SYNC: grmhd_primitives_Bi,grmhd_primitives_Bi_stagger GROUP HydroBase_Con2Prim: Convert from conservative to primitive variables GROUP SetTmunu: Group for calculating the stress-energy tensor TmunuBase::TmunuBase_ZeroTmunu: Initialise the stress-energy tensor to zero GROUP AddToTmunu: Add to the stress-energy tensor here IllinoisGRMHD::IllinoisGRMHD_conserv_to_prims: Compute primitive variables from conservatives. This is non-trivial, requiring a Newton-Raphson root-finder. IllinoisGRMHD::IllinoisGRMHD_outer_boundaries_on_P_rho_b_vx_vy_vz: Apply outflow-only, flat BCs on {P,rho_b,vx,vy,vz}. Outflow only == velocities pointed inward from outer boundary are set to zero. GROUP MoL_PseudoEvolutionBoundaries: Apply boundary conditions to pseudo-evolved quantities GROUP ML_ADMConstraints_evaluate_bc_group: ML_ADMConstraints_evaluate ML_ADMConstraints::ML_ADMConstraints_evaluate_SelectBCs: [level] ML_ADMConstraints_evaluate_SelectBCs GROUP ML_ADMConstraints_evaluate_ApplyBCs: Apply BCs for groups set in ML_ADMConstraints_evaluate GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ML_BSSN_ConstraintsEverywhere_bc_group: ML_BSSN_ConstraintsEverywhere ML_BSSN::ML_BSSN_ConstraintsEverywhere_SelectBCs: [level] ML_BSSN_ConstraintsEverywhere_SelectBCs GROUP ML_BSSN_ConstraintsEverywhere_ApplyBCs: Apply BCs for groups set in ML_BSSN_ConstraintsEverywhere GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ML_BSSN_ConstraintsInterior_bc_group: ML_BSSN_ConstraintsInterior ML_BSSN::ML_BSSN_ConstraintsInterior_SelectBCs: [level] ML_BSSN_ConstraintsInterior_SelectBCs GROUP ML_BSSN_ConstraintsInterior_ApplyBCs: Apply BCs for groups set in ML_BSSN_ConstraintsInterior GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions [CCTK_POSTPOSTINITIAL] AHFinderDirect::AHFinderDirect_import_mask: [global] [loop-local] import the excision mask AHFinderDirect::AHFinderDirect_find_horizons: [global] find apparent horizon(s) after this time step AHFinderDirect::AHFinderDirect_store: [global] store apparent horizon(s) into spherical surface(s) AHFinderDirect::AHFinderDirect_save: [global] save apparent horizon(s) into Cactus variables AHFinderDirect::AHFinderDirect_maybe_do_masks: set mask(s) based on apparent horizon position(s) GROUP Con2Prim: Convert from conservative to primitive variables (might be redundant) GROUP IllinoisGRMHD_PostPostInitial: HydroBase_Con2Prim in CCTK_POSTPOSTINITIAL sets conserv to prim then outer boundaries (OBs, which are technically disabled). The post OB SYNCs actually reprolongate the conservative variables, making cons and prims INCONSISTENT. So here we redo the con2prim, avoiding the SYNC afterward, then copy the result to other timelevels IllinoisGRMHD::postid: Schedule symmetries -- Actually just a placeholder function to ensure prolongations / processor syncs are done BEFORE outer boundaries are updated. IllinoisGRMHD::compute_b: Compute B and B_stagger from A SYNC: grmhd_primitives_Bi,grmhd_primitives_Bi_stagger IllinoisGRMHD::IllinoisGRMHD_conserv_to_prims: Compute primitive variables from conservatives. This is non-trivial, requiring a Newton-Raphson root-finder. IllinoisGRMHD::mhdpostid: Compute post-initialdata quantities [CCTK_POSTSTEP] AHFinderDirect::AHFinderDirect_import_mask: [global] [loop-local] import the excision mask AHFinderDirect::AHFinderDirect_find_horizons: [global] find apparent horizon(s) after this time step AHFinderDirect::AHFinderDirect_store: [global] store apparent horizon(s) into spherical surface(s) AHFinderDirect::AHFinderDirect_save: [global] save apparent horizon(s) into Cactus variables AHFinderDirect::AHFinderDirect_maybe_do_masks: set mask(s) based on apparent horizon position(s) SphericalSurface::SphericalSurface_Set: [global] Set surface radii GROUP zzz_NaNChecker_NaNCheck: Check for NaNs and count them in NaNChecker::NaNsFound NaNChecker::NaNChecker_NaNCheck_Prepare: [level] Prepare data structures to check for NaNs NaNChecker::NaNChecker_NaNCheck_Check: [local] Check for NaNs NaNChecker::NaNChecker_NaNCheck_Finish: [level] Count NaNs in NaNChecker::NaNsFound NaNChecker::NaNChecker_TakeAction: [global] [loop-level] Output NaNChecker::NaNmask and take action according to NaNChecker::action_if_found SpaceMask::CheckMask: Ensure that all mask values are legal GROUP SphericalSurface_HasBeenSet: Set the spherical surfaces before this group, and use it afterwards SphericalSurface::SphericalSurface_CheckState: [global] Test the state of the spherical surfaces Dissipation::setup_epsdis: Setup spatially varying dissipation endif if (recover initial data) [CCTK_BASEGRID] ADMBase::ADMBase_SetShiftStateOn: Set the shift_state variable to 1 ADMBase::ADMBase_SetDtLapseStateOn: Set the dtlapse_state variable to 1 ADMBase::ADMBase_SetDtShiftStateOn: Set the dtshift_state variable to 1 ADMMacros::ADMMacros_SetLocalSpatialOrder: Initialize the local_spatial_order CartGrid3D::SpatialSpacings: Set up ranges for spatial 3D Cartesian coordinates (on all grids) CartGrid3D::SpatialCoordinates: Set up spatial 3D Cartesian coordinates on the GH SphericalSurface::SphericalSurface_SetupRes: [global] [loop-local] Set surface resolution automatically AHFinderDirect::AHFinderDirect_setup: [global] setup data structures Dissipation::dissipation_basegrid: Ensure that there are enough ghost zones IllinoisGRMHD::IllinoisGRMHD_InitSymBound: Schedule symmetries ML_ADMConstraints::ML_ADMConstraints_CheckBoundaries: [meta] check boundaries treatment ML_BSSN::ML_BSSN_CheckBoundaries: [meta] check boundaries treatment NaNChecker::NaNChecker_ResetCounter: [global] Reset the NaNChecker::NaNsFound counter SpaceMask::MaskZero: Initialise mask to zero SpaceMask::MaskOne: Set old style mask to one SphericalSurface::SphericalSurface_Setup: [global] Calculate surface coordinate descriptors GROUP MaskBase_SetupMask: Set up the weight function GROUP MaskBase_SetupMaskAll: Set up the weight function CarpetReduce::MaskBase_AllocateMask: [global] Allocate the weight function CarpetReduce::MaskBase_InitMask: [global] [loop-local] Initialise the weight function GROUP SetupIMaskInternal: Set up the integer weight function (schedule other routines in here) CarpetReduce::CoordBase_SetupMask: [global] [loop-local] Set up the outer boundaries of the weight function CarpetReduce::CarpetMaskSetup: [global] [loop-singlemap] Set up the weight function for the restriction regions GROUP SetupIMask: Set up the integer weight function (schedule other routines in here) CarpetReduce::MaskBase_SetMask: [global] [loop-local] Set the weight function GROUP SetupMask: Set up the real weight function (schedule other routines in here) CarpetReduce::MaskBase_TestMask: [global] Test the weight function SphericalSurface::SphericalSurface_Set: [global] Set surface radii to be used for initial setup in other thorns GROUP SphericalSurface_HasBeenSet: Set the spherical surfaces before this group, and use it afterwards SphericalSurface::SphericalSurface_CheckState: [global] Test the state of the spherical surfaces SummationByParts::SBP_SetNormMask: Setup the mask for the calculation of the norm SymBase::SymBase_Check: Check whether the driver set up the grid consistently Time::Time_Initialise: [global] Initialise Time variables Time::TemporalSpacings: [singlemap] Set timestep based on Courant condition (courant_static) [CCTK_RECOVER_VARIABLES] IOUtil::IOUtil_RecoverGH: [level] Checkpoint recovery routine [CCTK_POST_RECOVER_VARIABLES] AHFinderDirect::AHFinderDirect_recover: [global] import horizon data from Cactus variables CarpetIOHDF5::CarpetIOHDF5_InitCheckpointingIntervals: [global] Initialisation of checkpointing intervals after recovery CarpetIOHDF5::CarpetIOHDF5_CloseFiles: [meta] Close all initial data checkpoint files after recovery GROUP MaskBase_SetupMask: Set up the weight function GROUP MaskBase_SetupMaskAll: Set up the weight function CarpetReduce::MaskBase_AllocateMask: [global] Allocate the weight function CarpetReduce::MaskBase_InitMask: [global] [loop-local] Initialise the weight function GROUP SetupIMaskInternal: Set up the integer weight function (schedule other routines in here) CarpetReduce::CoordBase_SetupMask: [global] [loop-local] Set up the outer boundaries of the weight function CarpetReduce::CarpetMaskSetup: [global] [loop-singlemap] Set up the weight function for the restriction regions GROUP SetupIMask: Set up the integer weight function (schedule other routines in here) CarpetReduce::MaskBase_SetMask: [global] [loop-local] Set the weight function GROUP SetupMask: Set up the real weight function (schedule other routines in here) CarpetReduce::MaskBase_TestMask: [global] Test the weight function GROUP zzz_NaNChecker_NaNCheck: Check for NaNs and count them in NaNChecker::NaNsFound NaNChecker::NaNChecker_NaNCheck_Prepare: [level] Prepare data structures to check for NaNs NaNChecker::NaNChecker_NaNCheck_Check: [local] Check for NaNs NaNChecker::NaNChecker_NaNCheck_Finish: [level] Count NaNs in NaNChecker::NaNsFound NaNChecker::NaNChecker_TakeAction: [global] [loop-level] Output NaNChecker::NaNmask and take action according to NaNChecker::action_if_found TerminationTrigger::TerminationTrigger_ResetMinutes: [global] Reset Watchtime VolumeIntegrals_GRMHD::VI_GRMHD_file_output_routine_Startup: Create directory for VI GRMHD file output. VolumeIntegrals_GRMHD::VI_GRMHD_InitializeIntegralCounterToZero: [global] Initialize IntegralCounter variable to zero endif if (checkpoint initial data) [CCTK_CPINITIAL] CarpetIOHDF5::CarpetIOHDF5_InitialDataCheckpoint: [meta] Initial data checkpoint routine endif if (analysis) [CCTK_ANALYSIS] VolumeIntegrals_GRMHD::VI_GRMHD_InitializeIntegralCounter: [global] Initialize IntegralCounter variable while (VolumeIntegrals_GRMHD::IntegralCounter) GROUP VI_GRMHD_VolumeIntegralGroup: Evaluate all volume integrals VolumeIntegrals_GRMHD::VI_GRMHD_ComputeIntegrand: [global] [loop-local] Compute Integrand VolumeIntegrals_GRMHD::VI_GRMHD_DoSum: [global] Do Sum VolumeIntegrals_GRMHD::VI_GRMHD_DecrementIntegralCounter: [global] Decrement IntegralCounter variable end while Convert_to_HydroBase::Convert_to_HydroBase: [global] [loop-local] Convert IllinoisGRMHD-native variables to HydroBase LoopControl::lc_statistics_analysis: [meta] Output LoopControl statistics GROUP ML_BSSN_EvolutionAnalysis: Calculate RHS at analysis ML_BSSN::ML_BSSN_EvolutionAnalysisInit: ML_BSSN_EvolutionAnalysisInit ML_BSSN::ML_BSSN_EvolutionAnalysisInterior: ML_BSSN_EvolutionAnalysisInterior ML_BSSN_Helper::ML_BSSN_NewRad: Apply NewRad boundary conditions to RHS SystemStatistics::SystemStatistics_Collect: [global] Collect system statistics TerminationTrigger::TerminationTrigger_CheckWalltime: Check elapsed job walltime TerminationTrigger::TerminationTrigger_CheckFile: Check termination file CarpetLib::CarpetLib_printmemstats: [global] Print memory statistics if desired CarpetLib::CarpetLib_printtimestats: [global] Print timing statistics if desired VolumeIntegrals_GRMHD::VI_GRMHD_file_output: [global] Output VolumeIntegral results to disk particle_tracerET::convert_to_MHD_3velocity: [global] [loop-local] Convert to MHD 3 velocity GROUP particle_tracerET: Particle tracerET subroutines particle_tracerET::InitializeRK4IterationCounterToOne: [global] Initialize RK4IterationCounter variable to one particle_tracerET::InitializeParticlePositions: [global] Initialize particle positions particle_tracerET::DoOneRK4StepForParticleTracerET: [global] Do one RK4 step for particle tracerET particle_tracerET::particle_tracerET_file_output: [global] Output particle tracerET data to disk smallbPoynET::compute_bi_b2_Poyn_fluxET: [global] [loop-local] Set b^mu, b^2, and Poynting flux gridfunctions. endif Output grid variables do loop over timesteps [CCTK_PREREGRID] Change grid hierarchy [CCTK_POSTREGRID] SpaceMask::MaskZero: Initialise mask to zero CartGrid3D::SpatialCoordinates: Set Coordinates after regridding GROUP MaskBase_SetupMask: Set up the weight function GROUP MaskBase_SetupMaskAll: Set up the weight function CarpetReduce::MaskBase_AllocateMask: [global] Allocate the weight function CarpetReduce::MaskBase_InitMask: [global] [loop-local] Initialise the weight function GROUP SetupIMaskInternal: Set up the integer weight function (schedule other routines in here) CarpetReduce::CoordBase_SetupMask: [global] [loop-local] Set up the outer boundaries of the weight function CarpetReduce::CarpetMaskSetup: [global] [loop-singlemap] Set up the weight function for the restriction regions GROUP SetupIMask: Set up the integer weight function (schedule other routines in here) CarpetReduce::MaskBase_SetMask: [global] [loop-local] Set the weight function GROUP SetupMask: Set up the real weight function (schedule other routines in here) CarpetReduce::MaskBase_TestMask: [global] Test the weight function Dissipation::setup_epsdis: Setup spatially varying dissipation GROUP MoL_PostStep: Ensure that everything is correct after regridding ML_BSSN::ML_BSSN_SelectBoundConds: [level] select boundary conditions GROUP ML_BSSN_ApplyBCs: Apply boundary conditions controlled by thorn Boundary GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions ML_ADMConstraints::ML_ADMConstraints_SelectBoundConds: [level] select boundary conditions GROUP ML_ADMConstraints_ApplyBCs: Apply boundary conditions controlled by thorn Boundary GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions ML_BSSN::ML_BSSN_ADMBaseInterior: ML_BSSN_ADMBaseInterior ML_BSSN::ML_BSSN_ADMBaseBoundaryScalar: ML_BSSN_ADMBaseBoundaryScalar ML_BSSN::ML_BSSN_ADMBaseEverywhere: ML_BSSN_ADMBaseEverywhere ML_BSSN_Helper::ML_BSSN_ADMBase_SelectBCs: [level] Select boundary conditions for ADMBase variables GROUP ML_BSSN_ADMBase_ApplyBCs: Apply boundary conditions to ADMBase variables GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ADMBase_SetADMVars: Set the ADM variables before this group, and use them afterwards GROUP HydroBase_PostStep: Post step tasks for hydro thorns GROUP HydroBase_Boundaries: HydroBase-internal Boundary conditions group GROUP HydroBase_Select_Boundaries: Group to schedule the boundary condition functions GROUP HydroBase_ApplyBCs: Apply the boundary conditions of HydroBase GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions IllinoisGRMHD::Boundary_SYNCs: Schedule symmetries -- Actually just a placeholder function to ensure prolongations / processor syncs are done BEFORE outer boundaries are updated. IllinoisGRMHD::IllinoisGRMHD_outer_boundaries_on_A_mu: Apply linear extrapolation BCs on A_{mu}, so that BCs are flat on B^i. IllinoisGRMHD::IllinoisGRMHD_compute_B_and_Bstagger_from_A: Compute B and B_stagger from A, SYNC: grmhd_primitives_Bi,grmhd_primitives_Bi_stagger GROUP HydroBase_Con2Prim: Convert from conservative to primitive variables GROUP SetTmunu: Group for calculating the stress-energy tensor TmunuBase::TmunuBase_ZeroTmunu: Initialise the stress-energy tensor to zero GROUP AddToTmunu: Add to the stress-energy tensor here IllinoisGRMHD::IllinoisGRMHD_conserv_to_prims: Compute primitive variables from conservatives. This is non-trivial, requiring a Newton-Raphson root-finder. IllinoisGRMHD::IllinoisGRMHD_outer_boundaries_on_P_rho_b_vx_vy_vz: Apply outflow-only, flat BCs on {P,rho_b,vx,vy,vz}. Outflow only == velocities pointed inward from outer boundary are set to zero. GROUP MoL_PseudoEvolutionBoundaries: Apply boundary conditions to pseudo-evolved quantities GROUP ML_ADMConstraints_evaluate_bc_group: ML_ADMConstraints_evaluate ML_ADMConstraints::ML_ADMConstraints_evaluate_SelectBCs: [level] ML_ADMConstraints_evaluate_SelectBCs GROUP ML_ADMConstraints_evaluate_ApplyBCs: Apply BCs for groups set in ML_ADMConstraints_evaluate GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ML_BSSN_ConstraintsEverywhere_bc_group: ML_BSSN_ConstraintsEverywhere ML_BSSN::ML_BSSN_ConstraintsEverywhere_SelectBCs: [level] ML_BSSN_ConstraintsEverywhere_SelectBCs GROUP ML_BSSN_ConstraintsEverywhere_ApplyBCs: Apply BCs for groups set in ML_BSSN_ConstraintsEverywhere GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ML_BSSN_ConstraintsInterior_bc_group: ML_BSSN_ConstraintsInterior ML_BSSN::ML_BSSN_ConstraintsInterior_SelectBCs: [level] ML_BSSN_ConstraintsInterior_SelectBCs GROUP ML_BSSN_ConstraintsInterior_ApplyBCs: Apply BCs for groups set in ML_BSSN_ConstraintsInterior GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions SpaceMask::MaskOne: Set mask to one AHFinderDirect::AHFinderDirect_maybe_do_masks: set mask(s) based on apparent horizon position(s) Rotate timelevels iteration = iteration+1 t = t+dt [CCTK_PRESTEP] CoordGauge::Einstein_SetNextSlicing: Identify the slicing for the next iteration LoopControl::lc_steer: [meta] Update LoopControl algorithm preferences NaNChecker::NaNChecker_ResetCounter: [global] Reset the NaNChecker::NaNsFound counter [CCTK_EVOL] MoL::MoL_StartLoop: [level] Initialise the step size control while (MoL::MoL_Stepsize_Bad) GROUP MoL_Evolution: A single Cactus evolution step using MoL GROUP MoL_StartStep: MoL internal setup for the evolution step MoL::MoL_SetCounter: [level] Set the counter for the ODE method to loop over MoL::MoL_SetTime: [level] Ensure the correct time and timestep are used MoL::MoL_AllocateScratchSpace: [level] Allocate storage for scratch levels GROUP MoL_PreStep: Physics thorns can schedule preloop setup routines in here MoL::MoL_AllocateScratch: Allocate sufficient space for array scratch variables MoL::MoL_InitialCopy: Ensure the data is in the correct timelevel while (MoL::MoL_Intermediate_Step) GROUP MoL_Step: The loop over the intermediate steps for the ODE integrator MoL::MoL_InitRHS: Initialise the RHS functions GROUP MoL_CalcRHS: Physics thorns schedule the calculation of the discrete spatial operator in here GROUP HydroBase_RHS: Groups for scheduling tasks for calculating RHS of hydro variables IllinoisGRMHD::IllinoisGRMHD_RHS_eval: Evaluate RHSs of GR Hydro & GRMHD equations ML_BSSN::ML_BSSN_EvolutionBoundaryScalar: ML_BSSN_EvolutionBoundaryScalar GROUP ML_BSSN_EvolutionInteriorSplitBy: ML_BSSN::ML_BSSN_EvolutionInteriorSplitBy1: ML_BSSN_EvolutionInteriorSplitBy1 ML_BSSN::ML_BSSN_EvolutionInteriorSplitBy2: ML_BSSN_EvolutionInteriorSplitBy2 ML_BSSN::ML_BSSN_EvolutionInteriorSplitBy3: ML_BSSN_EvolutionInteriorSplitBy3 ML_BSSN_Helper::ML_BSSN_NewRad: Apply NewRad boundary conditions to RHS GROUP MoL_PostRHS: Modify RHS functions Dissipation::dissipation_add: Add Kreiss-Oliger dissipation to the right hand sides GROUP MoL_RHSBoundaries: Any 'final' modifications to the RHS functions (boundaries etc.) MoL::MoL_Add: Updates calculated with the efficient Runge-Kutta 4 method MoL::MoL_DecrementCounter: [level] Alter the counter number MoL::MoL_ResetTime: [level] If necessary, change the time GROUP MoL_PostStepModify: The group for physics thorns to schedule enforcing constraints ML_BSSN::ML_BSSN_EnforceEverywhere: ML_BSSN_EnforceEverywhere GROUP MoL_PostStep: The group for physics thorns to schedule boundary calls etc. ML_BSSN::ML_BSSN_SelectBoundConds: [level] select boundary conditions GROUP ML_BSSN_ApplyBCs: Apply boundary conditions controlled by thorn Boundary GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions ML_ADMConstraints::ML_ADMConstraints_SelectBoundConds: [level] select boundary conditions GROUP ML_ADMConstraints_ApplyBCs: Apply boundary conditions controlled by thorn Boundary GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions ML_BSSN::ML_BSSN_ADMBaseInterior: ML_BSSN_ADMBaseInterior ML_BSSN::ML_BSSN_ADMBaseBoundaryScalar: ML_BSSN_ADMBaseBoundaryScalar ML_BSSN::ML_BSSN_ADMBaseEverywhere: ML_BSSN_ADMBaseEverywhere ML_BSSN_Helper::ML_BSSN_ADMBase_SelectBCs: [level] Select boundary conditions for ADMBase variables GROUP ML_BSSN_ADMBase_ApplyBCs: Apply boundary conditions to ADMBase variables GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ADMBase_SetADMVars: Set the ADM variables before this group, and use them afterwards GROUP HydroBase_PostStep: Post step tasks for hydro thorns GROUP HydroBase_Boundaries: HydroBase-internal Boundary conditions group GROUP HydroBase_Select_Boundaries: Group to schedule the boundary condition functions GROUP HydroBase_ApplyBCs: Apply the boundary conditions of HydroBase GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions IllinoisGRMHD::Boundary_SYNCs: Schedule symmetries -- Actually just a placeholder function to ensure prolongations / processor syncs are done BEFORE outer boundaries are updated. IllinoisGRMHD::IllinoisGRMHD_outer_boundaries_on_A_mu: Apply linear extrapolation BCs on A_{mu}, so that BCs are flat on B^i. IllinoisGRMHD::IllinoisGRMHD_compute_B_and_Bstagger_from_A: Compute B and B_stagger from A, SYNC: grmhd_primitives_Bi,grmhd_primitives_Bi_stagger GROUP HydroBase_Con2Prim: Convert from conservative to primitive variables GROUP SetTmunu: Group for calculating the stress-energy tensor TmunuBase::TmunuBase_ZeroTmunu: Initialise the stress-energy tensor to zero GROUP AddToTmunu: Add to the stress-energy tensor here IllinoisGRMHD::IllinoisGRMHD_conserv_to_prims: Compute primitive variables from conservatives. This is non-trivial, requiring a Newton-Raphson root-finder. IllinoisGRMHD::IllinoisGRMHD_outer_boundaries_on_P_rho_b_vx_vy_vz: Apply outflow-only, flat BCs on {P,rho_b,vx,vy,vz}. Outflow only == velocities pointed inward from outer boundary are set to zero. MoL::MoL_ResetDeltaTime: [level] If necessary, change the timestep end while MoL::MoL_FinishLoop: [level] Control the step size MoL::MoL_RestoreSandR: Restoring the Save and Restore variables to the original state MoL::MoL_FreeScratchSpace: [level] Free storage for scratch levels end while GROUP MoL_PseudoEvolution: Calculate pseudo-evolved quantities GROUP ADMBase_SetADMVars: Set the ADM variables before this group, and use them afterwards GROUP ML_ADMConstraints_evaluate_group: ML_ADMConstraints_evaluate ML_ADMConstraints::ML_ADMConstraints_evaluate: ML_ADMConstraints_evaluate GROUP ML_ADMConstraints_evaluate_bc_group: ML_ADMConstraints_evaluate ML_ADMConstraints::ML_ADMConstraints_evaluate_SelectBCs: [level] ML_ADMConstraints_evaluate_SelectBCs GROUP ML_ADMConstraints_evaluate_ApplyBCs: Apply BCs for groups set in ML_ADMConstraints_evaluate GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ML_BSSN_ConstraintsEverywhere_group: ML_BSSN_ConstraintsEverywhere ML_BSSN::ML_BSSN_ConstraintsEverywhere: ML_BSSN_ConstraintsEverywhere GROUP ML_BSSN_ConstraintsEverywhere_bc_group: ML_BSSN_ConstraintsEverywhere ML_BSSN::ML_BSSN_ConstraintsEverywhere_SelectBCs: [level] ML_BSSN_ConstraintsEverywhere_SelectBCs GROUP ML_BSSN_ConstraintsEverywhere_ApplyBCs: Apply BCs for groups set in ML_BSSN_ConstraintsEverywhere GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ML_BSSN_ConstraintsInterior_group: ML_BSSN_ConstraintsInterior ML_BSSN::ML_BSSN_ConstraintsInterior: ML_BSSN_ConstraintsInterior GROUP ML_BSSN_ConstraintsInterior_bc_group: ML_BSSN_ConstraintsInterior ML_BSSN::ML_BSSN_ConstraintsInterior_SelectBCs: [level] ML_BSSN_ConstraintsInterior_SelectBCs GROUP ML_BSSN_ConstraintsInterior_ApplyBCs: Apply BCs for groups set in ML_BSSN_ConstraintsInterior GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions Evolve finer grids recursively Restrict from finer grids [CCTK_POSTRESTRICT] GROUP MoL_PostStep: Ensure that everything is correct after restriction ML_BSSN::ML_BSSN_SelectBoundConds: [level] select boundary conditions GROUP ML_BSSN_ApplyBCs: Apply boundary conditions controlled by thorn Boundary GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions ML_ADMConstraints::ML_ADMConstraints_SelectBoundConds: [level] select boundary conditions GROUP ML_ADMConstraints_ApplyBCs: Apply boundary conditions controlled by thorn Boundary GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions ML_BSSN::ML_BSSN_ADMBaseInterior: ML_BSSN_ADMBaseInterior ML_BSSN::ML_BSSN_ADMBaseBoundaryScalar: ML_BSSN_ADMBaseBoundaryScalar ML_BSSN::ML_BSSN_ADMBaseEverywhere: ML_BSSN_ADMBaseEverywhere ML_BSSN_Helper::ML_BSSN_ADMBase_SelectBCs: [level] Select boundary conditions for ADMBase variables GROUP ML_BSSN_ADMBase_ApplyBCs: Apply boundary conditions to ADMBase variables GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ADMBase_SetADMVars: Set the ADM variables before this group, and use them afterwards GROUP HydroBase_PostStep: Post step tasks for hydro thorns GROUP HydroBase_Boundaries: HydroBase-internal Boundary conditions group GROUP HydroBase_Select_Boundaries: Group to schedule the boundary condition functions GROUP HydroBase_ApplyBCs: Apply the boundary conditions of HydroBase GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions IllinoisGRMHD::Boundary_SYNCs: Schedule symmetries -- Actually just a placeholder function to ensure prolongations / processor syncs are done BEFORE outer boundaries are updated. IllinoisGRMHD::IllinoisGRMHD_outer_boundaries_on_A_mu: Apply linear extrapolation BCs on A_{mu}, so that BCs are flat on B^i. IllinoisGRMHD::IllinoisGRMHD_compute_B_and_Bstagger_from_A: Compute B and B_stagger from A, SYNC: grmhd_primitives_Bi,grmhd_primitives_Bi_stagger GROUP HydroBase_Con2Prim: Convert from conservative to primitive variables GROUP SetTmunu: Group for calculating the stress-energy tensor TmunuBase::TmunuBase_ZeroTmunu: Initialise the stress-energy tensor to zero GROUP AddToTmunu: Add to the stress-energy tensor here IllinoisGRMHD::IllinoisGRMHD_conserv_to_prims: Compute primitive variables from conservatives. This is non-trivial, requiring a Newton-Raphson root-finder. IllinoisGRMHD::IllinoisGRMHD_outer_boundaries_on_P_rho_b_vx_vy_vz: Apply outflow-only, flat BCs on {P,rho_b,vx,vy,vz}. Outflow only == velocities pointed inward from outer boundary are set to zero. GROUP MoL_PseudoEvolutionBoundaries: Apply boundary conditions to pseudo-evolved quantities GROUP ML_ADMConstraints_evaluate_bc_group: ML_ADMConstraints_evaluate ML_ADMConstraints::ML_ADMConstraints_evaluate_SelectBCs: [level] ML_ADMConstraints_evaluate_SelectBCs GROUP ML_ADMConstraints_evaluate_ApplyBCs: Apply BCs for groups set in ML_ADMConstraints_evaluate GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ML_BSSN_ConstraintsEverywhere_bc_group: ML_BSSN_ConstraintsEverywhere ML_BSSN::ML_BSSN_ConstraintsEverywhere_SelectBCs: [level] ML_BSSN_ConstraintsEverywhere_SelectBCs GROUP ML_BSSN_ConstraintsEverywhere_ApplyBCs: Apply BCs for groups set in ML_BSSN_ConstraintsEverywhere GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ML_BSSN_ConstraintsInterior_bc_group: ML_BSSN_ConstraintsInterior ML_BSSN::ML_BSSN_ConstraintsInterior_SelectBCs: [level] ML_BSSN_ConstraintsInterior_SelectBCs GROUP ML_BSSN_ConstraintsInterior_ApplyBCs: Apply BCs for groups set in ML_BSSN_ConstraintsInterior GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions [CCTK_POSTSTEP] AHFinderDirect::AHFinderDirect_import_mask: [global] [loop-local] import the excision mask AHFinderDirect::AHFinderDirect_find_horizons: [global] find apparent horizon(s) after this time step AHFinderDirect::AHFinderDirect_store: [global] store apparent horizon(s) into spherical surface(s) AHFinderDirect::AHFinderDirect_save: [global] save apparent horizon(s) into Cactus variables AHFinderDirect::AHFinderDirect_maybe_do_masks: set mask(s) based on apparent horizon position(s) SphericalSurface::SphericalSurface_Set: [global] Set surface radii GROUP zzz_NaNChecker_NaNCheck: Check for NaNs and count them in NaNChecker::NaNsFound NaNChecker::NaNChecker_NaNCheck_Prepare: [level] Prepare data structures to check for NaNs NaNChecker::NaNChecker_NaNCheck_Check: [local] Check for NaNs NaNChecker::NaNChecker_NaNCheck_Finish: [level] Count NaNs in NaNChecker::NaNsFound NaNChecker::NaNChecker_TakeAction: [global] [loop-level] Output NaNChecker::NaNmask and take action according to NaNChecker::action_if_found SpaceMask::CheckMask: Ensure that all mask values are legal GROUP SphericalSurface_HasBeenSet: Set the spherical surfaces before this group, and use it afterwards SphericalSurface::SphericalSurface_CheckState: [global] Test the state of the spherical surfaces Dissipation::setup_epsdis: Setup spatially varying dissipation if (checkpoint) [CCTK_CHECKPOINT] CarpetIOHDF5::CarpetIOHDF5_EvolutionCheckpoint: [meta] Evolution checkpoint routine endif if (analysis) [CCTK_ANALYSIS] VolumeIntegrals_GRMHD::VI_GRMHD_InitializeIntegralCounter: [global] Initialize IntegralCounter variable while (VolumeIntegrals_GRMHD::IntegralCounter) GROUP VI_GRMHD_VolumeIntegralGroup: Evaluate all volume integrals VolumeIntegrals_GRMHD::VI_GRMHD_ComputeIntegrand: [global] [loop-local] Compute Integrand VolumeIntegrals_GRMHD::VI_GRMHD_DoSum: [global] Do Sum VolumeIntegrals_GRMHD::VI_GRMHD_DecrementIntegralCounter: [global] Decrement IntegralCounter variable end while Convert_to_HydroBase::Convert_to_HydroBase: [global] [loop-local] Convert IllinoisGRMHD-native variables to HydroBase LoopControl::lc_statistics_analysis: [meta] Output LoopControl statistics GROUP ML_BSSN_EvolutionAnalysis: Calculate RHS at analysis ML_BSSN::ML_BSSN_EvolutionAnalysisInit: ML_BSSN_EvolutionAnalysisInit ML_BSSN::ML_BSSN_EvolutionAnalysisInterior: ML_BSSN_EvolutionAnalysisInterior ML_BSSN_Helper::ML_BSSN_NewRad: Apply NewRad boundary conditions to RHS SystemStatistics::SystemStatistics_Collect: [global] Collect system statistics TerminationTrigger::TerminationTrigger_CheckWalltime: Check elapsed job walltime TerminationTrigger::TerminationTrigger_CheckFile: Check termination file CarpetLib::CarpetLib_printmemstats: [global] Print memory statistics if desired CarpetLib::CarpetLib_printtimestats: [global] Print timing statistics if desired VolumeIntegrals_GRMHD::VI_GRMHD_file_output: [global] Output VolumeIntegral results to disk particle_tracerET::convert_to_MHD_3velocity: [global] [loop-local] Convert to MHD 3 velocity GROUP particle_tracerET: Particle tracerET subroutines particle_tracerET::InitializeRK4IterationCounterToOne: [global] Initialize RK4IterationCounter variable to one particle_tracerET::InitializeParticlePositions: [global] Initialize particle positions particle_tracerET::DoOneRK4StepForParticleTracerET: [global] Do one RK4 step for particle tracerET particle_tracerET::particle_tracerET_file_output: [global] Output particle tracerET data to disk smallbPoynET::compute_bi_b2_Poyn_fluxET: [global] [loop-local] Set b^mu, b^2, and Poynting flux gridfunctions. endif Output grid variables enddo Termination routines [CCTK_TERMINATE] CarpetIOHDF5::CarpetIOHDF5_TerminationCheckpoint: [meta] Termination checkpoint routine LoopControl::lc_statistics_terminate: [meta] Output LoopControl statistics MoL::MoL_FreeIndexArrays: Free the MoL bookkeeping index arrays Shutdown routines [CCTK_SHUTDOWN] Timers::Timer_Shutdown: Prepare hierarchical timers Routines run after changing the grid hierarchy: [CCTK_POSTREGRID] SpaceMask::MaskZero: Initialise mask to zero CartGrid3D::SpatialCoordinates: Set Coordinates after regridding GROUP MaskBase_SetupMask: Set up the weight function GROUP MaskBase_SetupMaskAll: Set up the weight function CarpetReduce::MaskBase_AllocateMask: [global] Allocate the weight function CarpetReduce::MaskBase_InitMask: [global] [loop-local] Initialise the weight function GROUP SetupIMaskInternal: Set up the integer weight function (schedule other routines in here) CarpetReduce::CoordBase_SetupMask: [global] [loop-local] Set up the outer boundaries of the weight function CarpetReduce::CarpetMaskSetup: [global] [loop-singlemap] Set up the weight function for the restriction regions GROUP SetupIMask: Set up the integer weight function (schedule other routines in here) CarpetReduce::MaskBase_SetMask: [global] [loop-local] Set the weight function GROUP SetupMask: Set up the real weight function (schedule other routines in here) CarpetReduce::MaskBase_TestMask: [global] Test the weight function Dissipation::setup_epsdis: Setup spatially varying dissipation GROUP MoL_PostStep: Ensure that everything is correct after regridding ML_BSSN::ML_BSSN_SelectBoundConds: [level] select boundary conditions GROUP ML_BSSN_ApplyBCs: Apply boundary conditions controlled by thorn Boundary GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions ML_ADMConstraints::ML_ADMConstraints_SelectBoundConds: [level] select boundary conditions GROUP ML_ADMConstraints_ApplyBCs: Apply boundary conditions controlled by thorn Boundary GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions ML_BSSN::ML_BSSN_ADMBaseInterior: ML_BSSN_ADMBaseInterior ML_BSSN::ML_BSSN_ADMBaseBoundaryScalar: ML_BSSN_ADMBaseBoundaryScalar ML_BSSN::ML_BSSN_ADMBaseEverywhere: ML_BSSN_ADMBaseEverywhere ML_BSSN_Helper::ML_BSSN_ADMBase_SelectBCs: [level] Select boundary conditions for ADMBase variables GROUP ML_BSSN_ADMBase_ApplyBCs: Apply boundary conditions to ADMBase variables GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ADMBase_SetADMVars: Set the ADM variables before this group, and use them afterwards GROUP HydroBase_PostStep: Post step tasks for hydro thorns GROUP HydroBase_Boundaries: HydroBase-internal Boundary conditions group GROUP HydroBase_Select_Boundaries: Group to schedule the boundary condition functions GROUP HydroBase_ApplyBCs: Apply the boundary conditions of HydroBase GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions IllinoisGRMHD::Boundary_SYNCs: Schedule symmetries -- Actually just a placeholder function to ensure prolongations / processor syncs are done BEFORE outer boundaries are updated. IllinoisGRMHD::IllinoisGRMHD_outer_boundaries_on_A_mu: Apply linear extrapolation BCs on A_{mu}, so that BCs are flat on B^i. IllinoisGRMHD::IllinoisGRMHD_compute_B_and_Bstagger_from_A: Compute B and B_stagger from A, SYNC: grmhd_primitives_Bi,grmhd_primitives_Bi_stagger GROUP HydroBase_Con2Prim: Convert from conservative to primitive variables GROUP SetTmunu: Group for calculating the stress-energy tensor TmunuBase::TmunuBase_ZeroTmunu: Initialise the stress-energy tensor to zero GROUP AddToTmunu: Add to the stress-energy tensor here IllinoisGRMHD::IllinoisGRMHD_conserv_to_prims: Compute primitive variables from conservatives. This is non-trivial, requiring a Newton-Raphson root-finder. IllinoisGRMHD::IllinoisGRMHD_outer_boundaries_on_P_rho_b_vx_vy_vz: Apply outflow-only, flat BCs on {P,rho_b,vx,vy,vz}. Outflow only == velocities pointed inward from outer boundary are set to zero. GROUP MoL_PseudoEvolutionBoundaries: Apply boundary conditions to pseudo-evolved quantities GROUP ML_ADMConstraints_evaluate_bc_group: ML_ADMConstraints_evaluate ML_ADMConstraints::ML_ADMConstraints_evaluate_SelectBCs: [level] ML_ADMConstraints_evaluate_SelectBCs GROUP ML_ADMConstraints_evaluate_ApplyBCs: Apply BCs for groups set in ML_ADMConstraints_evaluate GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ML_BSSN_ConstraintsEverywhere_bc_group: ML_BSSN_ConstraintsEverywhere ML_BSSN::ML_BSSN_ConstraintsEverywhere_SelectBCs: [level] ML_BSSN_ConstraintsEverywhere_SelectBCs GROUP ML_BSSN_ConstraintsEverywhere_ApplyBCs: Apply BCs for groups set in ML_BSSN_ConstraintsEverywhere GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions GROUP ML_BSSN_ConstraintsInterior_bc_group: ML_BSSN_ConstraintsInterior ML_BSSN::ML_BSSN_ConstraintsInterior_SelectBCs: [level] ML_BSSN_ConstraintsInterior_SelectBCs GROUP ML_BSSN_ConstraintsInterior_ApplyBCs: Apply BCs for groups set in ML_BSSN_ConstraintsInterior GROUP BoundaryConditions: Execute all boundary conditions Boundary::Boundary_ApplyPhysicalBCs: Apply all requested local physical boundary conditions CartGrid3D::CartGrid3D_ApplyBC: Apply symmetry boundary conditions Boundary::Boundary_ClearSelection: [level] Unselect all grid variables for boundary conditions SpaceMask::MaskOne: Set mask to one AHFinderDirect::AHFinderDirect_maybe_do_masks: set mask(s) based on apparent horizon position(s) -------------------------------------------------------------------------------- INFO (Carpet): Multi-Model listing: model 0: "world" INFO (Carpet): Multi-Model process distribution: processes 0-1: model 0 "world" INFO (Carpet): Multi-Model: This is process 0, model 0 "world" INFO (CarpetLib): Process startup time was 4.74 seconds Current core file size limit: hard=[unlimited], soft=0 MB Current memory size limit: hard=[unlimited], soft=[unlimited] INFO (CycleClock): Measuring CycleClock tick via OpenMP... INFO (CycleClock): Calibrated CycleClock: 0.416645 ns per clock tick (2.40013 GHz) INFO (Vectors): Using vector size 4 for architecture AVX (64-bit precision) -------------------------------------------------------------------------------- AMR driver provided by Carpet -------------------------------------------------------------------------------- HydroBase: Let it flow. -------------------------------------------------------------------------------- AMR 0D ASCII I/O provided by CarpetIOASCII -------------------------------------------------------------------------------- AMR 1D ASCII I/O provided by CarpetIOASCII -------------------------------------------------------------------------------- AMR 2D ASCII I/O provided by CarpetIOASCII -------------------------------------------------------------------------------- AMR 3D ASCII I/O provided by CarpetIOASCII -------------------------------------------------------------------------------- AMR info I/O provided by CarpetIOBasic -------------------------------------------------------------------------------- ML_ADMConstraints -------------------------------------------------------------------------------- ML_BSSN -------------------------------------------------------------------------------- AMR HDF5 I/O provided by CarpetIOHDF5 -------------------------------------------------------------------------------- AMR 0D HDF5 I/O provided by CarpetIOHDF5 -------------------------------------------------------------------------------- AMR 1D HDF5 I/O provided by CarpetIOHDF5 -------------------------------------------------------------------------------- AMR 2D HDF5 I/O provided by CarpetIOHDF5 -------------------------------------------------------------------------------- AMR 3D HDF5 I/O provided by CarpetIOHDF5 -------------------------------------------------------------------------------- MoL: Generalized time integration. -------------------------------------------------------------------------------- AMR scalar I/O provided by CarpetIOScalar -------------------------------------------------------------------------------- INFO (Carpet): MPI is enabled INFO (Carpet): Carpet is running on 2 processes INFO (Carpet): This is process 0 INFO (Carpet): OpenMP is enabled INFO (Carpet): This process contains 12 threads, this is thread 0 INFO (Carpet): There are 24 threads in total INFO (Carpet): There are 12 threads per process INFO (Carpet): This process runs on host nid02984, pid=34463 INFO (Carpet): This process runs on 12 cores: 0-5, 24-29 INFO (Carpet): Thread 0 runs on 12 cores: 0-5, 24-29 INFO (Carpet): Thread 1 runs on 12 cores: 0-5, 24-29 INFO (Carpet): Thread 2 runs on 12 cores: 0-5, 24-29 INFO (Carpet): Thread 3 runs on 12 cores: 0-5, 24-29 INFO (Carpet): Thread 4 runs on 12 cores: 0-5, 24-29 INFO (Carpet): Thread 5 runs on 12 cores: 0-5, 24-29 INFO (Carpet): Thread 6 runs on 12 cores: 0-5, 24-29 INFO (Carpet): Thread 7 runs on 12 cores: 0-5, 24-29 INFO (Carpet): Thread 8 runs on 12 cores: 0-5, 24-29 INFO (Carpet): Thread 9 runs on 12 cores: 0-5, 24-29 INFO (Carpet): Thread 10 runs on 12 cores: 0-5, 24-29 INFO (Carpet): Thread 11 runs on 12 cores: 0-5, 24-29 INFO (Carpet): This simulation is running in 3 dimensions INFO (Carpet): Boundary specification for map 0: nboundaryzones: [[3,3,3],[3,3,3]] is_internal : [[0,0,0],[0,0,0]] is_staggered : [[0,0,0],[0,0,0]] shiftout : [[0,0,0],[0,0,0]] INFO (Carpet): CoordBase domain specification for map 0: physical extent: [-540,-540,-540] : [540,540,540] ([1080,1080,1080]) interior extent: [-522,-522,-522] : [522,522,522] ([1044,1044,1044]) exterior extent: [-576,-576,-576] : [576,576,576] ([1152,1152,1152]) base_spacing : [18,18,18] INFO (Carpet): Adapted domain specification for map 0: convergence factor: 2 convergence level : 0 physical extent : [-540,-540,-540] : [540,540,540] ([1080,1080,1080]) interior extent : [-522,-522,-522] : [522,522,522] ([1044,1044,1044]) exterior extent : [-576,-576,-576] : [576,576,576] ([1152,1152,1152]) spacing : [18,18,18] INFO (Carpet): Base grid specification for map 0: number of grid points : [65,65,65] number of coarse grid ghost points: [[3,3,3],[3,3,3]] INFO (Carpet): Buffer zone counts (excluding ghosts): [0]: [[0,0,0],[0,0,0]] [1]: [[9,9,9],[9,9,9]] [2]: [[9,9,9],[9,9,9]] [3]: [[9,9,9],[9,9,9]] [4]: [[9,9,9],[9,9,9]] [5]: [[9,9,9],[9,9,9]] [6]: [[9,9,9],[9,9,9]] INFO (Carpet): Overlap zone counts: [0]: [[0,0,0],[0,0,0]] [1]: [[0,0,0],[0,0,0]] [2]: [[0,0,0],[0,0,0]] [3]: [[0,0,0],[0,0,0]] [4]: [[0,0,0],[0,0,0]] [5]: [[0,0,0],[0,0,0]] [6]: [[0,0,0],[0,0,0]] INFO (Carpet): Group and variable statistics: INFO (Carpet): There are 994 grid functions in 96 groups INFO (Carpet): There are 322 grid scalars in 58 groups INFO (Carpet): There are 140 1-dimensional grid arrays in 14 groups INFO (Carpet): There are 5 2-dimensional grid arrays in 3 groups INFO (Carpet): There are 1 3-dimensional grid arrays in 1 groups INFO (Carpet): (The number of variables counts all time levels) INFO (CarpetIOASCII): I/O Method 'IOASCII_0D' registered: 0D AMR output of grid variables to ASCII files INFO (CarpetIOASCII): Periodic 0D AMR output requested for: CARPET::physical_time_per_hour SPHERICALSURFACE::sf_active[0] SPHERICALSURFACE::sf_active[1] SPHERICALSURFACE::sf_active[2] INFO (CarpetIOASCII): I/O Method 'IOASCII_1D' registered: 1D AMR output of grid variables to ASCII files INFO (CarpetIOASCII): Periodic 1D AMR output requested for: ADMBASE::gxx ADMBASE::gxy ADMBASE::gxz ADMBASE::gyy ADMBASE::gyz ADMBASE::gzz ADMBASE::kxx ADMBASE::kxy ADMBASE::kxz ADMBASE::kyy ADMBASE::kyz ADMBASE::kzz ADMBASE::alp ADMBASE::betax ADMBASE::betay ADMBASE::betaz CARPETREDUCE::weight ILLINOISGRMHD::rho_b ILLINOISGRMHD::P ILLINOISGRMHD::vx ILLINOISGRMHD::vy ILLINOISGRMHD::vz ILLINOISGRMHD::Bx ILLINOISGRMHD::By ILLINOISGRMHD::Bz ML_BSSN::H ML_BSSN::M1 ML_BSSN::M2 ML_BSSN::M3 INFO (CarpetIOASCII): I/O Method 'IOASCII_2D' registered: 2D AMR output of grid variables to ASCII files INFO (CarpetIOASCII): Periodic 2D AMR output requested for: ADMBASE::gxx ADMBASE::gxy ADMBASE::gxz ADMBASE::gyy ADMBASE::gyz ADMBASE::gzz ADMBASE::kxx ADMBASE::kxy ADMBASE::kxz ADMBASE::kyy ADMBASE::kyz ADMBASE::kzz ADMBASE::alp ADMBASE::betax ADMBASE::betay ADMBASE::betaz HYDROBASE::rho HYDROBASE::eps HYDROBASE::vel[0] HYDROBASE::vel[1] HYDROBASE::vel[2] HYDROBASE::w_lorentz ILLINOISGRMHD::rho_b ILLINOISGRMHD::P ILLINOISGRMHD::vx ILLINOISGRMHD::vy ILLINOISGRMHD::vz ILLINOISGRMHD::Bx ILLINOISGRMHD::By ILLINOISGRMHD::Bz ILLINOISGRMHD::failure_checker ML_BSSN::H ML_BSSN::M1 ML_BSSN::M2 ML_BSSN::M3 SPHERICALSURFACE::sf_radius[0] SPHERICALSURFACE::sf_radius[1] SPHERICALSURFACE::sf_radius[2] SMALLBPOYNET::smallbx SMALLBPOYNET::smallby SMALLBPOYNET::smallbz SMALLBPOYNET::smallb2 SMALLBPOYNET::minus_one_minus_u_0 INFO (CarpetIOASCII): I/O Method 'IOASCII_3D' registered: 3D AMR output of grid variables to ASCII files INFO (CarpetIOHDF5): I/O Method 'IOHDF5' registered: AMR output of grid variables to HDF5 files INFO (CarpetIOHDF5): Periodic AMR output requested for: ILLINOISGRMHD::rho_b ILLINOISGRMHD::P ILLINOISGRMHD::vx ILLINOISGRMHD::vy ILLINOISGRMHD::vz ILLINOISGRMHD::Bx ILLINOISGRMHD::By ILLINOISGRMHD::Bz ILLINOISGRMHD::failure_checker SMALLBPOYNET::smallb2 SMALLBPOYNET::minus_one_minus_u_0 INFO (CarpetIOHDF5): I/O Method 'IOHDF5_0D' registered: 0D AMR output of grid variables to HDF5 files INFO (CarpetIOHDF5): I/O Method 'IOHDF5_1D' registered: 1D AMR output of grid variables to HDF5 files INFO (CarpetIOHDF5): I/O Method 'IOHDF5_2D' registered: 2D AMR output of grid variables to HDF5 files INFO (CarpetIOHDF5): Periodic 2D AMR output requested for: ADMBASE::gxx ADMBASE::gxy ADMBASE::gxz ADMBASE::gyy ADMBASE::gyz ADMBASE::gzz ADMBASE::kxx ADMBASE::kxy ADMBASE::kxz ADMBASE::kyy ADMBASE::kyz ADMBASE::kzz ADMBASE::alp ADMBASE::betax ADMBASE::betay ADMBASE::betaz ILLINOISGRMHD::rho_b ILLINOISGRMHD::P ILLINOISGRMHD::vx ILLINOISGRMHD::vy ILLINOISGRMHD::vz ILLINOISGRMHD::Bx ILLINOISGRMHD::By ILLINOISGRMHD::Bz ML_ADMCONSTRAINTS::H INFO (CarpetIOHDF5): I/O Method 'IOHDF5_3D' registered: 3D AMR output of grid variables to HDF5 files INFO (CarpetIOScalar): Periodic scalar output requested for: ADMBASE::alp ILLINOISGRMHD::rho_b ILLINOISGRMHD::P ILLINOISGRMHD::vx ILLINOISGRMHD::vy ILLINOISGRMHD::vz ILLINOISGRMHD::Bx ILLINOISGRMHD::By ILLINOISGRMHD::Bz ML_BSSN::H ML_BSSN::M1 ML_BSSN::M2 ML_BSSN::M3 SMALLBPOYNET::smallbx SMALLBPOYNET::smallby SMALLBPOYNET::smallbz SMALLBPOYNET::smallb2 SMALLBPOYNET::minus_one_minus_u_0 INFO (CarpetRegrid2): Initialising position of centre 0 to [-15.2,0,0] INFO (CarpetRegrid2): Initialising position of centre 1 to [15.2,0,0] WARNING level 1 from host nid02984 process 0 while executing schedule bin CCTK_WRAGH, routine ML_BSSN_Helper::ML_BSSN_ParamCompat in thorn ML_BSSN_Helper, file /global/u1/c/chlin/Cactus/arrangements/McLachlan/ML_BSSN_Helper/src/Parameters.c:15: -> Forcing ML_BSSN::initial_boundary_condition="extrapolate-gammas" because ML_BSSN::my_initial_boundary_condition="extrapolate-gammas" WARNING level 1 from host nid02984 process 0 while executing schedule bin CCTK_WRAGH, routine ML_BSSN_Helper::ML_BSSN_ParamCompat in thorn ML_BSSN_Helper, file /global/u1/c/chlin/Cactus/arrangements/McLachlan/ML_BSSN_Helper/src/Parameters.c:28: -> Forcing ML_BSSN::rhs_boundary_condition="NewRad" because ML_BSSN::my_rhs_boundary_condition="NewRad" WARNING level 1 from host nid02984 process 0 while executing schedule bin CCTK_WRAGH, routine ML_BSSN_Helper::ML_BSSN_ParamCompat in thorn ML_BSSN_Helper, file /global/u1/c/chlin/Cactus/arrangements/McLachlan/ML_BSSN_Helper/src/Parameters.c:57: -> Forcing ML_BSSN::evolveA=0 because ML_BSSN::LapseACoeff=0.0 WARNING level 1 from host nid02984 process 0 while executing schedule bin CCTK_WRAGH, routine ML_BSSN_Helper::ML_BSSN_ParamCompat in thorn ML_BSSN_Helper, file /global/u1/c/chlin/Cactus/arrangements/McLachlan/ML_BSSN_Helper/src/Parameters.c:77: -> Forcing ML_BSSN::evolveB=0 because ML_BSSN::ShiftBCoeff=0.0 WARNING level 1 from host nid02984 process 0 while executing schedule bin CCTK_WRAGH, routine ML_BSSN_Helper::ML_BSSN_ParamCompat in thorn ML_BSSN_Helper, file /global/u1/c/chlin/Cactus/arrangements/McLachlan/ML_BSSN_Helper/src/Parameters.c:116: -> Forcing ML_BSSN::advectLapse=1 because ML_BSSN::LapseAdvectionCoeff=1.0 WARNING level 1 from host nid02984 process 0 while executing schedule bin CCTK_WRAGH, routine ML_BSSN_Helper::ML_BSSN_ParamCompat in thorn ML_BSSN_Helper, file /global/u1/c/chlin/Cactus/arrangements/McLachlan/ML_BSSN_Helper/src/Parameters.c:134: -> Forcing ML_BSSN::advectShift=1 because ML_BSSN::ShiftAdvectionCoeff=1.0 INFO (MoL): Using Runge-Kutta 4 as the time integrator. INFO (TerminationTrigger): Reminding you every 60 minutes about remaining walltime INFO (MoL): The maximum number of evolved variables is 769. 34 are registered. INFO (MoL): The maximum number of slow evolved variables is 769. 0 are registered. INFO (MoL): The maximum number of constrained variables is 769. 20 are registered. INFO (MoL): The maximum number of SandR variables is 769. 0 are registered. INFO (MoL): The maximum number of evolved array variables is 769. 0 are registered. INFO (MoL): The maximum number of constrained array variables is 769. 0 are registered. INFO (MoL): The maximum number of SandR array variables is 769. 0 are registered. INFO (MoL): The maximum size of any array variables is 0. WARNING level 1 from host nid02984 process 0 while executing schedule bin CCTK_PARAMCHECK, routine ML_BSSN_Helper::ML_BSSN_ParamCheck in thorn ML_BSSN_Helper, file /global/u1/c/chlin/Cactus/arrangements/McLachlan/ML_BSSN_Helper/src/Parameters.c:157: -> Parameter ML_BSSN::my_initial_boundary_condition is outdated; please update the parameter file. Do not use this parameter, and set up initial boundary conditions as usual. WARNING level 1 from host nid02984 process 0 while executing schedule bin CCTK_PARAMCHECK, routine ML_BSSN_Helper::ML_BSSN_ParamCheck in thorn ML_BSSN_Helper, file /global/u1/c/chlin/Cactus/arrangements/McLachlan/ML_BSSN_Helper/src/Parameters.c:164: -> Parameter ML_BSSN::my_rhs_boundary_condition is outdated; please update the parameter file. Do not use this parameter, and set up RHS boundary conditions as usual. WARNING level 1 from host nid02984 process 0 while executing schedule bin CCTK_PARAMCHECK, routine ML_BSSN_Helper::ML_BSSN_ParamCheck in thorn ML_BSSN_Helper, file /global/u1/c/chlin/Cactus/arrangements/McLachlan/ML_BSSN_Helper/src/Parameters.c:184: -> Parameter ML_BSSN::LapseACoeff is outdated; please update the parameter file. Instead of using this parameter, you should set ML_BSSN::evolveA. WARNING level 1 from host nid02984 process 0 while executing schedule bin CCTK_PARAMCHECK, routine ML_BSSN_Helper::ML_BSSN_ParamCheck in thorn ML_BSSN_Helper, file /global/u1/c/chlin/Cactus/arrangements/McLachlan/ML_BSSN_Helper/src/Parameters.c:191: -> Parameter ML_BSSN::ShiftBCoeff is outdated; please update the parameter file. Instead of using this parameter, you should set ML_BSSN::evolveB. WARNING level 1 from host nid02984 process 0 while executing schedule bin CCTK_PARAMCHECK, routine ML_BSSN_Helper::ML_BSSN_ParamCheck in thorn ML_BSSN_Helper, file /global/u1/c/chlin/Cactus/arrangements/McLachlan/ML_BSSN_Helper/src/Parameters.c:198: -> Parameter ML_BSSN::LapseAdvectionCoeff is outdated; please update the parameter file. Instead of using this parameter, you should set ML_BSSN::advectLapse. WARNING level 1 from host nid02984 process 0 while executing schedule bin CCTK_PARAMCHECK, routine ML_BSSN_Helper::ML_BSSN_ParamCheck in thorn ML_BSSN_Helper, file /global/u1/c/chlin/Cactus/arrangements/McLachlan/ML_BSSN_Helper/src/Parameters.c:204: -> Parameter ML_BSSN::ShiftAdvectionCoeff is outdated; please update the parameter file. Instead of using this parameter, you should set ML_BSSN::advectShift. INFO (Vectors): Testing vectorisation... [errors may result in segfaults] INFO (Vectors): 2300/2300 tests passed INFO (CarpetRegrid2): Centre 1 is at position [-15.2,0,0] with 7 levels INFO (CarpetRegrid2): Centre 2 is at position [15.2,0,0] with 7 levels INFO (CarpetRegrid2): Regridding level 0 map 0 at iteration 0 time 0 INFO (CarpetRegrid2): Regridding levels 1 and up INFO (CarpetRegrid2): Enforcing grid structure properties, iteration 0 INFO (CarpetRegrid2): Enforcing grid structure properties, iteration 0 INFO (CarpetRegrid2): Enforcing grid structure properties, iteration 0 INFO (CarpetRegrid2): Enforcing grid structure properties, iteration 0 INFO (CarpetRegrid2): Enforcing grid structure properties, iteration 0 INFO (CarpetRegrid2): Enforcing grid structure properties, iteration 0 INFO (Carpet): Grid structure (superregions, grid points): [0][0][0] exterior: [0,0,0] : [64,64,64] ([65,65,65] + PADDING) 274625 [1][0][0] exterior: [23,25,25] : [105,103,103] ([83,79,79] + PADDING) 518003 [2][0][0] exterior: [80,83,83] : [176,173,173] ([97,91,91] + PADDING) 803257 [3][0][0] exterior: [198,205,205] : [314,307,307] ([117,103,103] + PADDING) 1241253 [4][0][0] exterior: [453,467,467] : [571,557,557] ([119,91,91] + PADDING) 985439 [5][0][0] exterior: [940,967,967] : [1108,1081,1081] ([169,115,115] + PADDING) 2235025 [6][0][0] exterior: [1937,1991,1991] : [2159,2105,2105] ([223,115,115] + PADDING) 2949175 INFO (Carpet): Grid structure (superregions, coordinates): [0][0][0] exterior: [-576.000000000000000,-576.000000000000000,-576.000000000000000] : [576.000000000000000,576.000000000000000,576.000000000000000] : [18.000000000000000,18.000000000000000,18.000000000000000] [1][0][0] exterior: [-369.000000000000000,-351.000000000000000,-351.000000000000000] : [369.000000000000000,351.000000000000000,351.000000000000000] : [9.000000000000000,9.000000000000000,9.000000000000000] [2][0][0] exterior: [-216.000000000000000,-202.500000000000000,-202.500000000000000] : [216.000000000000000,202.500000000000000,202.500000000000000] : [4.500000000000000,4.500000000000000,4.500000000000000] [3][0][0] exterior: [-130.500000000000000,-114.750000000000000,-114.750000000000000] : [130.500000000000000,114.750000000000000,114.750000000000000] : [2.250000000000000,2.250000000000000,2.250000000000000] [4][0][0] exterior: [-66.375000000000000,-50.625000000000000,-50.625000000000000] : [66.375000000000000,50.625000000000000,50.625000000000000] : [1.125000000000000,1.125000000000000,1.125000000000000] [5][0][0] exterior: [-47.250000000000000,-32.062500000000000,-32.062500000000000] : [47.250000000000000,32.062500000000000,32.062500000000000] : [0.562500000000000,0.562500000000000,0.562500000000000] [6][0][0] exterior: [-31.218750000000000,-16.031250000000000,-16.031250000000000] : [31.218750000000000,16.031250000000000,16.031250000000000] : [0.281250000000000,0.281250000000000,0.281250000000000] INFO (Carpet): There are not enough time levels for the desired temporal prolongation order in the grid function group "HYDROBASE::RHO". With Carpet::prolongation_order_time=2, you need at least 3 time levels. INFO (Carpet): There are not enough time levels for the desired temporal prolongation order in the grid function group "HYDROBASE::PRESS". With Carpet::prolongation_order_time=2, you need at least 3 time levels. INFO (Carpet): There are not enough time levels for the desired temporal prolongation order in the grid function group "HYDROBASE::EPS". With Carpet::prolongation_order_time=2, you need at least 3 time levels. INFO (Carpet): There are not enough time levels for the desired temporal prolongation order in the grid function group "HYDROBASE::VEL". With Carpet::prolongation_order_time=2, you need at least 3 time levels. INFO (Carpet): There are not enough time levels for the desired temporal prolongation order in the grid function group "HYDROBASE::W_LORENTZ". With Carpet::prolongation_order_time=2, you need at least 3 time levels. INFO (Carpet): There are not enough time levels for the desired temporal prolongation order in the grid function group "HYDROBASE::BVEC". With Carpet::prolongation_order_time=2, you need at least 3 time levels. INFO (Carpet): There are not enough time levels for the desired temporal prolongation order in the grid function group "HYDROBASE::AVEC". With Carpet::prolongation_order_time=2, you need at least 3 time levels. INFO (Carpet): There are not enough time levels for the desired temporal prolongation order in the grid function group "HYDROBASE::APHI". With Carpet::prolongation_order_time=2, you need at least 3 time levels. INFO (Carpet): Global grid structure statistics: INFO (Carpet): GF: rhs: 1452k active, 2271k owned (+56%), 2680k total (+18%), 64 steps/time INFO (Carpet): GF: vars: 370, pts: 2024M active, 3307M owned (+63%), 4021M total (+22%), 1.0 comp/proc INFO (Carpet): GA: vars: 463, pts: 0M active, 0M total (+-0%) INFO (Carpet): Total required memory: 32.085 GByte (for GAs and currently active GFs) INFO (Carpet): Load balance: min avg max sdv max/avg-1 INFO (Carpet): Level 0: 37M 38M 39M 1M owned 2% INFO (Carpet): Level 1: 95M 96M 97M 1M owned 1% INFO (Carpet): Level 2: 147M 149M 150M 2M owned 1% INFO (Carpet): Level 3: 228M 230M 232M 2M owned 1% INFO (Carpet): Level 4: 181M 182M 184M 2M owned 1% INFO (Carpet): Level 5: 411M 413M 416M 2M owned 1% INFO (Carpet): Level 6: 543M 546M 548M 2M owned 0% INFO (CartGrid3D): Grid Spacings: INFO (CartGrid3D): dx=>1.8000000e+01 dy=>1.8000000e+01 dz=>1.8000000e+01 INFO (CartGrid3D): Computational Coordinates: INFO (CartGrid3D): x=>[-576.000,576.000] y=>[-576.000,576.000] z=>[-576.000,576.000] INFO (CartGrid3D): Indices of Physical Coordinates: INFO (CartGrid3D): x=>[0,64] y=>[0,64] z=>[0,64] INFO (AHFinderDirect): setting up AHFinderDirect data structures INFO (AHFinderDirect): to search for 1 horizon on 2 processors INFO (AHFinderDirect): constructing full sphere patch system INFO (AHFinderDirect): at origin (0,0,0) INFO (AHFinderDirect): with 18 angular zones per right angle INFO (AHFinderDirect): ==> 2166 nominal, 3174 ghosted angular grid points INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (SphericalSurface): SphericalSurface[0]: finest refinement-level that fully contains this surface: sf_maxreflevel = -1 INFO (SphericalSurface): SphericalSurface[0]: finest refinement-level that overlaps with this surface : sf_minreflevel = -1 INFO (SphericalSurface): SphericalSurface[0]: sf_delta_theta = 0.089760, sf_delta_phi = 0.087266 INFO (SphericalSurface): SphericalSurface[1]: finest refinement-level that fully contains this surface: sf_maxreflevel = -1 INFO (SphericalSurface): SphericalSurface[1]: finest refinement-level that overlaps with this surface : sf_minreflevel = -1 INFO (SphericalSurface): SphericalSurface[1]: sf_delta_theta = 0.089760, sf_delta_phi = 0.087266 INFO (SphericalSurface): SphericalSurface[2]: finest refinement-level that fully contains this surface: sf_maxreflevel = -1 INFO (SphericalSurface): SphericalSurface[2]: finest refinement-level that overlaps with this surface : sf_minreflevel = -1 INFO (SphericalSurface): SphericalSurface[2]: sf_delta_theta = 0.089760, sf_delta_phi = 0.087266 INFO (Time): Timestep set to 4.5 (courant_static) INFO (Meudon_Bin_NS): Setting up LORENE Bin_NS initial data INFO (Meudon_Bin_NS): Setting up coordinates INFO (Meudon_Bin_NS): Meudon_Bin_NS will use the 2D_Polytrope equation of state. INFO (Meudon_Bin_NS): Meudon_Bin_NS will use the 1 eos handle INFO (Meudon_Bin_NS): Reading from file "/global/u1/c/chlin/lorene_data/test_gr_09_063.d" Binary system read in file : ---------------------------- Binary system ============= Orbital angular velocity : 842.9922849147580 rad/s Coordinate separation between the two stellar centers : 80.0000000000000 km Absolute coordinate X of the rotation axis : 0.0000000000000 km Star 1 : ====== Relativistic star ----------------- Number of domains occupied by the star : 1 Equation of state : Star 1 EOS EOS of class Eos_poly (relativistic polytrope) : Adiabatic index gamma : 2.0000000000000 Pressure coefficient kappa : 0.1000000000000 rho_nuc c^2 / n_nuc^gamma Mean particle mass : 1.0000000000000 m_B Relativistic chemical potential at zero pressure : 1.0000000000000 m_B c^2 Central enthalpy : 0.0900000000000 c^2 Central proper baryon density : 0.4708714185261 x 0.1 fm^-3 Central proper energy density : 0.4930434078045 rho_nuc c^2 Central pressure : 0.0221719892785 rho_nuc c^2 Regularization index of the gravitational potential : k_div = 0 Central lapse N : 0.8213808112110 Central value of A^2 : 1.4514630408533 Coordinate equatorial radius (phi=0) a1 = 31.5924006423005 km Coordinate equatorial radius (phi=pi/2) a2 = 27.7774975562862 km Coordinate equatorial radius (phi=pi): 30.3543794804130 km Coordinate polar radius a3 = 27.7919576947017 km Axis ratio a2/a1 = 0.8792461791933 a3/a1 = 0.8797038885829 Baryon mass : 1.6707909842769 M_sol Gravitational mass : 1.5865459517723 M_sol Star in a binary system ----------------------- irrotational configuration Absolute abscidia of the stellar center: -34.3158126713210 km Absolute abscidia of the barycenter of the baryon density : -34.1807445379738 km d_tilde : 2.1758549531469 Orientation with respect to the absolute frame : 0.0000000000000 rad Central value of gam_euler : 1.0036741742355 Central u_euler (U^X, U^Y, U^Z) [c] : -0.0000000000000 -0.0709573993086 0.0000000000000 Central d_psi (X, Y, Z) [c] : -0.0000000000000 -0.1131052894769 0.0000000000000 Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.0000000000000 0.0003935153316 0.0000000000000 Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : 0.3718331383477 0.4943138902066 0.0012024646722 Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.3718331383477 -0.2470525597791 -0.0012024646722 Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : 0.0435878798883 -0.0425606873520 -0.0006667894195 Central value of loggam : 0.0000001114528 Central value of log(N) auto, comp : -0.1735698739165 -0.0231985648935 Central value of beta=log(AN) auto, comp : -0.0100589249522 -0.0004234931416 Central value of shift (N^X, N^Y, N^Z) [c] : -0.0000000000000 -0.0378905962741 0.0000000000000 ... shift_auto part of it [c] : -0.0000000000000 -0.0489796785016 0.0000000000000 ... shift_comp part of it [c] : 0.0000000000000 0.0110890822274 0.0000000000000 ... w_shift (NB: components in the star Cartesian frame) [c] : -0.0000000000000 -0.0534744284212 0.0000000000000 Central value of khi_shift [km c] : -0.0000000000000 Central value of (B^X, B^Y, B^Z)/N [c] : -0.0000000000000 0.0713465773081 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : -0.0000362989755 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_comp) [km^{-1}] : -0.0002981883450 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : -0.0000112315610 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_comp) [km^{-1}] : -0.0000163565138 0.0000000000000 0.0000000000000 Central A^2 K^{ij} [c/km] : A^2 K^{xx} auto, comp : -0.0000000000000 0.0000000000000 A^2 K^{xy} auto, comp : -0.0000046453614 -0.0000650115950 A^2 K^{xz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yy} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{zz} auto, comp : 0.0000000000000 0.0000000000000 Central A^2 K_{ij} K^{ij} [c^2/km^2] : A^2 K_{ij} K^{ij} auto, comp : 0.0000000000574 0.0000000008767 Star 2 : ====== Relativistic star ----------------- Number of domains occupied by the star : 1 Equation of state : Star 1 EOS EOS of class Eos_poly (relativistic polytrope) : Adiabatic index gamma : 2.0000000000000 Pressure coefficient kappa : 0.1000000000000 rho_nuc c^2 / n_nuc^gamma Mean particle mass : 1.0000000000000 m_B Relativistic chemical potential at zero pressure : 1.0000000000000 m_B c^2 Central enthalpy : 0.0630000000000 c^2 Central proper baryon density : 0.3251341961565 x 0.1 fm^-3 Central proper energy density : 0.3357054207076 rho_nuc c^2 Central pressure : 0.0105712245510 rho_nuc c^2 Regularization index of the gravitational potential : k_div = 0 Central lapse N : 0.8560863131981 Central value of A^2 : 1.3462211527558 Coordinate equatorial radius (phi=0) a1 = 36.8227739723403 km Coordinate equatorial radius (phi=pi/2) a2 = 28.1647579751040 km Coordinate equatorial radius (phi=pi): 33.4745976028041 km Coordinate polar radius a3 = 28.0515599671080 km Axis ratio a2/a1 = 0.7648733361658 a3/a1 = 0.7617992057899 Baryon mass : 1.2830422471997 M_sol Gravitational mass : 1.2341346226771 M_sol Star in a binary system ----------------------- irrotational configuration Absolute abscidia of the stellar center: 45.6841873286790 km Absolute abscidia of the barycenter of the baryon density : 45.2667839099569 km d_tilde : 2.5042244486690 Orientation with respect to the absolute frame : 3.1415926535898 rad Central value of gam_euler : 1.0072303464165 Central u_euler (U^X, U^Y, U^Z) [c] : -0.0000000000000 0.1030839742961 0.0000000000000 Central d_psi (X, Y, Z) [c] : -0.0000000000000 0.1488664791184 0.0000000000000 Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.0000000000000 -0.0013151835884 0.0000000000000 Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : 0.3604065576231 0.2005959662083 0.0039290006645 Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.3604065576231 -0.5244769501283 -0.0039290006645 Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : -0.0435970713608 0.0424962994979 -0.0014660267039 Central value of loggam : 0.0000011471727 Central value of log(N) auto, comp : -0.1258082259951 -0.0295758487532 Central value of beta=log(AN) auto, comp : -0.0063507225771 -0.0003825914672 Central value of shift (N^X, N^Y, N^Z) [c] : -0.0000000000000 0.0391100446608 0.0000000000000 ... shift_auto part of it [c] : -0.0000000000000 0.0499120259801 0.0000000000000 ... shift_comp part of it [c] : 0.0000000000000 -0.0108019813193 0.0000000000000 ... w_shift (NB: components in the star Cartesian frame) [c] : 0.0000000000000 -0.0546332508576 0.0000000000000 Central value of khi_shift [km c] : -0.0000000000000 Central value of (B^X, B^Y, B^Z)/N [c] : -0.0000000000000 -0.1043705700727 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : 0.0000507123452 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_comp) [km^{-1}] : 0.0003740463855 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : 0.0000113135293 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_comp) [km^{-1}] : 0.0000156642372 0.0000000000000 0.0000000000000 Central A^2 K^{ij} [c/km] : A^2 K^{xx} auto, comp : -0.0000000000000 0.0000000000000 A^2 K^{xy} auto, comp : -0.0000046934496 -0.0000635848483 A^2 K^{xz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yy} auto, comp : 0.0000000000000 -0.0000000000000 A^2 K^{yz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{zz} auto, comp : -0.0000000000000 0.0000000000000 Central A^2 K_{ij} K^{ij} [c^2/km^2] : A^2 K_{ij} K^{ij} auto, comp : 0.0000000000497 0.0000000008035 Summary : ------- Separation d : 80.0000000000000 km Separation d_G : 79.4475284479308 km Omega : 842.9922849147580 rad/s Baryon mass of star 1 : 1.6707909842769 M_sol Baryon mass of star 2 : 1.2830422471997 M_sol ADM mass of the system : 2.8404508262857 M_sol Total angular momentum : 10.0839081191459 G M_sol^2 / c Radius of star 1 (x_comp) : 31.5924006423005 km Radius of star 1 (y) : 27.7774975562862 km Radius of star 1 (z) : 27.7919576947017 km Radius of star 1 (x_opp) : 30.3543794804130 km Radius of star 2 (x_comp) : 36.8227739723403 km Radius of star 2 (y) : 28.1647579751040 km Radius of star 2 (z) : 28.0515599671080 km Radius of star 2 (x_opp) : 33.4745976028041 km INFO (Meudon_Bin_NS): omega [rad/s]: 842.992 INFO (Meudon_Bin_NS): dist [km]: 80 INFO (Meudon_Bin_NS): dist_mass [km]: 79.4475 INFO (Meudon_Bin_NS): mass1_b [M_sun]: 1.67079 INFO (Meudon_Bin_NS): mass2_b [M_sun]: 1.28304 INFO (Meudon_Bin_NS): mass_ADM [M_sun]: 2.84045 INFO (Meudon_Bin_NS): L_tot [G M_sun^2/c]: 10.0839 INFO (Meudon_Bin_NS): rad1_x_comp [km]: 31.5924 INFO (Meudon_Bin_NS): rad1_y [km]: 27.7775 INFO (Meudon_Bin_NS): rad1_z [km]: 27.792 INFO (Meudon_Bin_NS): rad1_x_opp [km]: 30.3544 INFO (Meudon_Bin_NS): rad2_x_comp [km]: 36.8228 INFO (Meudon_Bin_NS): rad2_y [km]: 28.1648 INFO (Meudon_Bin_NS): rad2_z [km]: 28.0516 INFO (Meudon_Bin_NS): rad2_x_opp [km]: 33.4746 INFO (Meudon_Bin_NS): K [ET unit]: 372.101387884248 INFO (Meudon_Bin_NS): Filling in Cactus grid points INFO (Meudon_Bin_NS): Done. INFO (Seed_Magnetic_Fields_BNS): Seeding magnetic fields INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 0 NumPts= 152100 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 2.410e-16, ErrDenom: 5.571e-04 | 2.00 iters/gridpt INFO (IllinoisGRMHD): C2P: Lev: 0 NumPts= 152100 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 2.990e-16, ErrDenom: 5.571e-04 | 2.00 iters/gridpt INFO (VolumeIntegrals_GRMHD): Outputting to directory 'volume_integration/' INFO (particle_tracerET): Outputting to directory 'ABEoutput_directory-nsns_test/' INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 0 NumPts= 152100 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 1.557e-16, ErrDenom: 5.571e-04 | 2.00 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (Time): Timestep set to 4.5 (courant_static) INFO (Meudon_Bin_NS): Setting up LORENE Bin_NS initial data INFO (Meudon_Bin_NS): Setting up coordinates INFO (Meudon_Bin_NS): Meudon_Bin_NS will use the 2D_Polytrope equation of state. INFO (Meudon_Bin_NS): Meudon_Bin_NS will use the 1 eos handle INFO (Meudon_Bin_NS): Reading from file "/global/u1/c/chlin/lorene_data/test_gr_09_063.d" Binary system read in file : ---------------------------- Binary system ============= Orbital angular velocity : 842.9922849147580 rad/s Coordinate separation between the two stellar centers : 80.0000000000000 km Absolute coordinate X of the rotation axis : 0.0000000000000 km Star 1 : ====== Relativistic star ----------------- Number of domains occupied by the star : 1 Equation of state : Star 1 EOS EOS of class Eos_poly (relativistic polytrope) : Adiabatic index gamma : 2.0000000000000 Pressure coefficient kappa : 0.1000000000000 rho_nuc c^2 / n_nuc^gamma Mean particle mass : 1.0000000000000 m_B Relativistic chemical potential at zero pressure : 1.0000000000000 m_B c^2 Central enthalpy : 0.0900000000000 c^2 Central proper baryon density : 0.4708714185261 x 0.1 fm^-3 Central proper energy density : 0.4930434078045 rho_nuc c^2 Central pressure : 0.0221719892785 rho_nuc c^2 Regularization index of the gravitational potential : k_div = 0 Central lapse N : 0.8213808112110 Central value of A^2 : 1.4514630408533 Coordinate equatorial radius (phi=0) a1 = 31.5924006423005 km Coordinate equatorial radius (phi=pi/2) a2 = 27.7774975562862 km Coordinate equatorial radius (phi=pi): 30.3543794804130 km Coordinate polar radius a3 = 27.7919576947017 km Axis ratio a2/a1 = 0.8792461791933 a3/a1 = 0.8797038885829 Baryon mass : 1.6707909842769 M_sol Gravitational mass : 1.5865459517723 M_sol Star in a binary system ----------------------- irrotational configuration Absolute abscidia of the stellar center: -34.3158126713210 km Absolute abscidia of the barycenter of the baryon density : -34.1807445379738 km d_tilde : 2.1758549531469 Orientation with respect to the absolute frame : 0.0000000000000 rad Central value of gam_euler : 1.0036741742355 Central u_euler (U^X, U^Y, U^Z) [c] : -0.0000000000000 -0.0709573993086 0.0000000000000 Central d_psi (X, Y, Z) [c] : -0.0000000000000 -0.1131052894769 0.0000000000000 Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.0000000000000 0.0003935153316 0.0000000000000 Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : 0.3718331383477 0.4943138902066 0.0012024646722 Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.3718331383477 -0.2470525597791 -0.0012024646722 Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : 0.0435878798883 -0.0425606873520 -0.0006667894195 Central value of loggam : 0.0000001114528 Central value of log(N) auto, comp : -0.1735698739165 -0.0231985648935 Central value of beta=log(AN) auto, comp : -0.0100589249522 -0.0004234931416 Central value of shift (N^X, N^Y, N^Z) [c] : -0.0000000000000 -0.0378905962741 0.0000000000000 ... shift_auto part of it [c] : -0.0000000000000 -0.0489796785016 0.0000000000000 ... shift_comp part of it [c] : 0.0000000000000 0.0110890822274 0.0000000000000 ... w_shift (NB: components in the star Cartesian frame) [c] : -0.0000000000000 -0.0534744284212 0.0000000000000 Central value of khi_shift [km c] : -0.0000000000000 Central value of (B^X, B^Y, B^Z)/N [c] : -0.0000000000000 0.0713465773081 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : -0.0000362989755 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_comp) [km^{-1}] : -0.0002981883450 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : -0.0000112315610 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_comp) [km^{-1}] : -0.0000163565138 0.0000000000000 0.0000000000000 Central A^2 K^{ij} [c/km] : A^2 K^{xx} auto, comp : -0.0000000000000 0.0000000000000 A^2 K^{xy} auto, comp : -0.0000046453614 -0.0000650115950 A^2 K^{xz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yy} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{zz} auto, comp : 0.0000000000000 0.0000000000000 Central A^2 K_{ij} K^{ij} [c^2/km^2] : A^2 K_{ij} K^{ij} auto, comp : 0.0000000000574 0.0000000008767 Star 2 : ====== Relativistic star ----------------- Number of domains occupied by the star : 1 Equation of state : Star 1 EOS EOS of class Eos_poly (relativistic polytrope) : Adiabatic index gamma : 2.0000000000000 Pressure coefficient kappa : 0.1000000000000 rho_nuc c^2 / n_nuc^gamma Mean particle mass : 1.0000000000000 m_B Relativistic chemical potential at zero pressure : 1.0000000000000 m_B c^2 Central enthalpy : 0.0630000000000 c^2 Central proper baryon density : 0.3251341961565 x 0.1 fm^-3 Central proper energy density : 0.3357054207076 rho_nuc c^2 Central pressure : 0.0105712245510 rho_nuc c^2 Regularization index of the gravitational potential : k_div = 0 Central lapse N : 0.8560863131981 Central value of A^2 : 1.3462211527558 Coordinate equatorial radius (phi=0) a1 = 36.8227739723403 km Coordinate equatorial radius (phi=pi/2) a2 = 28.1647579751040 km Coordinate equatorial radius (phi=pi): 33.4745976028041 km Coordinate polar radius a3 = 28.0515599671080 km Axis ratio a2/a1 = 0.7648733361658 a3/a1 = 0.7617992057899 Baryon mass : 1.2830422471997 M_sol Gravitational mass : 1.2341346226771 M_sol Star in a binary system ----------------------- irrotational configuration Absolute abscidia of the stellar center: 45.6841873286790 km Absolute abscidia of the barycenter of the baryon density : 45.2667839099569 km d_tilde : 2.5042244486690 Orientation with respect to the absolute frame : 3.1415926535898 rad Central value of gam_euler : 1.0072303464165 Central u_euler (U^X, U^Y, U^Z) [c] : -0.0000000000000 0.1030839742961 0.0000000000000 Central d_psi (X, Y, Z) [c] : -0.0000000000000 0.1488664791184 0.0000000000000 Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.0000000000000 -0.0013151835884 0.0000000000000 Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : 0.3604065576231 0.2005959662083 0.0039290006645 Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.3604065576231 -0.5244769501283 -0.0039290006645 Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : -0.0435970713608 0.0424962994979 -0.0014660267039 Central value of loggam : 0.0000011471727 Central value of log(N) auto, comp : -0.1258082259951 -0.0295758487532 Central value of beta=log(AN) auto, comp : -0.0063507225771 -0.0003825914672 Central value of shift (N^X, N^Y, N^Z) [c] : -0.0000000000000 0.0391100446608 0.0000000000000 ... shift_auto part of it [c] : -0.0000000000000 0.0499120259801 0.0000000000000 ... shift_comp part of it [c] : 0.0000000000000 -0.0108019813193 0.0000000000000 ... w_shift (NB: components in the star Cartesian frame) [c] : 0.0000000000000 -0.0546332508576 0.0000000000000 Central value of khi_shift [km c] : -0.0000000000000 Central value of (B^X, B^Y, B^Z)/N [c] : -0.0000000000000 -0.1043705700727 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : 0.0000507123452 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_comp) [km^{-1}] : 0.0003740463855 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : 0.0000113135293 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_comp) [km^{-1}] : 0.0000156642372 0.0000000000000 0.0000000000000 Central A^2 K^{ij} [c/km] : A^2 K^{xx} auto, comp : -0.0000000000000 0.0000000000000 A^2 K^{xy} auto, comp : -0.0000046934496 -0.0000635848483 A^2 K^{xz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yy} auto, comp : 0.0000000000000 -0.0000000000000 A^2 K^{yz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{zz} auto, comp : -0.0000000000000 0.0000000000000 Central A^2 K_{ij} K^{ij} [c^2/km^2] : A^2 K_{ij} K^{ij} auto, comp : 0.0000000000497 0.0000000008035 Summary : ------- Separation d : 80.0000000000000 km Separation d_G : 79.4475284479308 km Omega : 842.9922849147580 rad/s Baryon mass of star 1 : 1.6707909842769 M_sol Baryon mass of star 2 : 1.2830422471997 M_sol ADM mass of the system : 2.8404508262857 M_sol Total angular momentum : 10.0839081191459 G M_sol^2 / c Radius of star 1 (x_comp) : 31.5924006423005 km Radius of star 1 (y) : 27.7774975562862 km Radius of star 1 (z) : 27.7919576947017 km Radius of star 1 (x_opp) : 30.3543794804130 km Radius of star 2 (x_comp) : 36.8227739723403 km Radius of star 2 (y) : 28.1647579751040 km Radius of star 2 (z) : 28.0515599671080 km Radius of star 2 (x_opp) : 33.4745976028041 km INFO (Meudon_Bin_NS): omega [rad/s]: 842.992 INFO (Meudon_Bin_NS): dist [km]: 80 INFO (Meudon_Bin_NS): dist_mass [km]: 79.4475 INFO (Meudon_Bin_NS): mass1_b [M_sun]: 1.67079 INFO (Meudon_Bin_NS): mass2_b [M_sun]: 1.28304 INFO (Meudon_Bin_NS): mass_ADM [M_sun]: 2.84045 INFO (Meudon_Bin_NS): L_tot [G M_sun^2/c]: 10.0839 INFO (Meudon_Bin_NS): rad1_x_comp [km]: 31.5924 INFO (Meudon_Bin_NS): rad1_y [km]: 27.7775 INFO (Meudon_Bin_NS): rad1_z [km]: 27.792 INFO (Meudon_Bin_NS): rad1_x_opp [km]: 30.3544 INFO (Meudon_Bin_NS): rad2_x_comp [km]: 36.8228 INFO (Meudon_Bin_NS): rad2_y [km]: 28.1648 INFO (Meudon_Bin_NS): rad2_z [km]: 28.0516 INFO (Meudon_Bin_NS): rad2_x_opp [km]: 33.4746 INFO (Meudon_Bin_NS): K [ET unit]: 372.101387884248 INFO (Meudon_Bin_NS): Filling in Cactus grid points INFO (Meudon_Bin_NS): Done. INFO (Seed_Magnetic_Fields_BNS): Seeding magnetic fields INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 1 NumPts= 346800 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 3.095e-15, ErrDenom: 3.667e-03 | 2.00 iters/gridpt INFO (IllinoisGRMHD): C2P: Lev: 1 NumPts= 346800 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 3.039e-16, ErrDenom: 3.667e-03 | 2.00 iters/gridpt INFO (VolumeIntegrals_GRMHD): Outputting to directory 'volume_integration/' INFO (particle_tracerET): Outputting to directory 'ABEoutput_directory-nsns_test/' INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 1 NumPts= 346800 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 8.463e-15, ErrDenom: 3.667e-03 | 2.00 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (Time): Timestep set to 2.25 (courant_static) INFO (Meudon_Bin_NS): Setting up LORENE Bin_NS initial data INFO (Meudon_Bin_NS): Setting up coordinates INFO (Meudon_Bin_NS): Meudon_Bin_NS will use the 2D_Polytrope equation of state. INFO (Meudon_Bin_NS): Meudon_Bin_NS will use the 1 eos handle INFO (Meudon_Bin_NS): Reading from file "/global/u1/c/chlin/lorene_data/test_gr_09_063.d" Binary system read in file : ---------------------------- Binary system ============= Orbital angular velocity : 842.9922849147580 rad/s Coordinate separation between the two stellar centers : 80.0000000000000 km Absolute coordinate X of the rotation axis : 0.0000000000000 km Star 1 : ====== Relativistic star ----------------- Number of domains occupied by the star : 1 Equation of state : Star 1 EOS EOS of class Eos_poly (relativistic polytrope) : Adiabatic index gamma : 2.0000000000000 Pressure coefficient kappa : 0.1000000000000 rho_nuc c^2 / n_nuc^gamma Mean particle mass : 1.0000000000000 m_B Relativistic chemical potential at zero pressure : 1.0000000000000 m_B c^2 Central enthalpy : 0.0900000000000 c^2 Central proper baryon density : 0.4708714185261 x 0.1 fm^-3 Central proper energy density : 0.4930434078045 rho_nuc c^2 Central pressure : 0.0221719892785 rho_nuc c^2 Regularization index of the gravitational potential : k_div = 0 Central lapse N : 0.8213808112110 Central value of A^2 : 1.4514630408533 Coordinate equatorial radius (phi=0) a1 = 31.5924006423005 km Coordinate equatorial radius (phi=pi/2) a2 = 27.7774975562862 km Coordinate equatorial radius (phi=pi): 30.3543794804130 km Coordinate polar radius a3 = 27.7919576947017 km Axis ratio a2/a1 = 0.8792461791933 a3/a1 = 0.8797038885829 Baryon mass : 1.6707909842769 M_sol Gravitational mass : 1.5865459517723 M_sol Star in a binary system ----------------------- irrotational configuration Absolute abscidia of the stellar center: -34.3158126713210 km Absolute abscidia of the barycenter of the baryon density : -34.1807445379738 km d_tilde : 2.1758549531469 Orientation with respect to the absolute frame : 0.0000000000000 rad Central value of gam_euler : 1.0036741742355 Central u_euler (U^X, U^Y, U^Z) [c] : -0.0000000000000 -0.0709573993086 0.0000000000000 Central d_psi (X, Y, Z) [c] : -0.0000000000000 -0.1131052894769 0.0000000000000 Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.0000000000000 0.0003935153316 0.0000000000000 Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : 0.3718331383477 0.4943138902066 0.0012024646722 Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.3718331383477 -0.2470525597791 -0.0012024646722 Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : 0.0435878798883 -0.0425606873520 -0.0006667894195 Central value of loggam : 0.0000001114528 Central value of log(N) auto, comp : -0.1735698739165 -0.0231985648935 Central value of beta=log(AN) auto, comp : -0.0100589249522 -0.0004234931416 Central value of shift (N^X, N^Y, N^Z) [c] : -0.0000000000000 -0.0378905962741 0.0000000000000 ... shift_auto part of it [c] : -0.0000000000000 -0.0489796785016 0.0000000000000 ... shift_comp part of it [c] : 0.0000000000000 0.0110890822274 0.0000000000000 ... w_shift (NB: components in the star Cartesian frame) [c] : -0.0000000000000 -0.0534744284212 0.0000000000000 Central value of khi_shift [km c] : -0.0000000000000 Central value of (B^X, B^Y, B^Z)/N [c] : -0.0000000000000 0.0713465773081 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : -0.0000362989755 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_comp) [km^{-1}] : -0.0002981883450 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : -0.0000112315610 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_comp) [km^{-1}] : -0.0000163565138 0.0000000000000 0.0000000000000 Central A^2 K^{ij} [c/km] : A^2 K^{xx} auto, comp : -0.0000000000000 0.0000000000000 A^2 K^{xy} auto, comp : -0.0000046453614 -0.0000650115950 A^2 K^{xz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yy} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{zz} auto, comp : 0.0000000000000 0.0000000000000 Central A^2 K_{ij} K^{ij} [c^2/km^2] : A^2 K_{ij} K^{ij} auto, comp : 0.0000000000574 0.0000000008767 Star 2 : ====== Relativistic star ----------------- Number of domains occupied by the star : 1 Equation of state : Star 1 EOS EOS of class Eos_poly (relativistic polytrope) : Adiabatic index gamma : 2.0000000000000 Pressure coefficient kappa : 0.1000000000000 rho_nuc c^2 / n_nuc^gamma Mean particle mass : 1.0000000000000 m_B Relativistic chemical potential at zero pressure : 1.0000000000000 m_B c^2 Central enthalpy : 0.0630000000000 c^2 Central proper baryon density : 0.3251341961565 x 0.1 fm^-3 Central proper energy density : 0.3357054207076 rho_nuc c^2 Central pressure : 0.0105712245510 rho_nuc c^2 Regularization index of the gravitational potential : k_div = 0 Central lapse N : 0.8560863131981 Central value of A^2 : 1.3462211527558 Coordinate equatorial radius (phi=0) a1 = 36.8227739723403 km Coordinate equatorial radius (phi=pi/2) a2 = 28.1647579751040 km Coordinate equatorial radius (phi=pi): 33.4745976028041 km Coordinate polar radius a3 = 28.0515599671080 km Axis ratio a2/a1 = 0.7648733361658 a3/a1 = 0.7617992057899 Baryon mass : 1.2830422471997 M_sol Gravitational mass : 1.2341346226771 M_sol Star in a binary system ----------------------- irrotational configuration Absolute abscidia of the stellar center: 45.6841873286790 km Absolute abscidia of the barycenter of the baryon density : 45.2667839099569 km d_tilde : 2.5042244486690 Orientation with respect to the absolute frame : 3.1415926535898 rad Central value of gam_euler : 1.0072303464165 Central u_euler (U^X, U^Y, U^Z) [c] : -0.0000000000000 0.1030839742961 0.0000000000000 Central d_psi (X, Y, Z) [c] : -0.0000000000000 0.1488664791184 0.0000000000000 Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.0000000000000 -0.0013151835884 0.0000000000000 Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : 0.3604065576231 0.2005959662083 0.0039290006645 Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.3604065576231 -0.5244769501283 -0.0039290006645 Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : -0.0435970713608 0.0424962994979 -0.0014660267039 Central value of loggam : 0.0000011471727 Central value of log(N) auto, comp : -0.1258082259951 -0.0295758487532 Central value of beta=log(AN) auto, comp : -0.0063507225771 -0.0003825914672 Central value of shift (N^X, N^Y, N^Z) [c] : -0.0000000000000 0.0391100446608 0.0000000000000 ... shift_auto part of it [c] : -0.0000000000000 0.0499120259801 0.0000000000000 ... shift_comp part of it [c] : 0.0000000000000 -0.0108019813193 0.0000000000000 ... w_shift (NB: components in the star Cartesian frame) [c] : 0.0000000000000 -0.0546332508576 0.0000000000000 Central value of khi_shift [km c] : -0.0000000000000 Central value of (B^X, B^Y, B^Z)/N [c] : -0.0000000000000 -0.1043705700727 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : 0.0000507123452 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_comp) [km^{-1}] : 0.0003740463855 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : 0.0000113135293 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_comp) [km^{-1}] : 0.0000156642372 0.0000000000000 0.0000000000000 Central A^2 K^{ij} [c/km] : A^2 K^{xx} auto, comp : -0.0000000000000 0.0000000000000 A^2 K^{xy} auto, comp : -0.0000046934496 -0.0000635848483 A^2 K^{xz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yy} auto, comp : 0.0000000000000 -0.0000000000000 A^2 K^{yz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{zz} auto, comp : -0.0000000000000 0.0000000000000 Central A^2 K_{ij} K^{ij} [c^2/km^2] : A^2 K_{ij} K^{ij} auto, comp : 0.0000000000497 0.0000000008035 Summary : ------- Separation d : 80.0000000000000 km Separation d_G : 79.4475284479308 km Omega : 842.9922849147580 rad/s Baryon mass of star 1 : 1.6707909842769 M_sol Baryon mass of star 2 : 1.2830422471997 M_sol ADM mass of the system : 2.8404508262857 M_sol Total angular momentum : 10.0839081191459 G M_sol^2 / c Radius of star 1 (x_comp) : 31.5924006423005 km Radius of star 1 (y) : 27.7774975562862 km Radius of star 1 (z) : 27.7919576947017 km Radius of star 1 (x_opp) : 30.3543794804130 km Radius of star 2 (x_comp) : 36.8227739723403 km Radius of star 2 (y) : 28.1647579751040 km Radius of star 2 (z) : 28.0515599671080 km Radius of star 2 (x_opp) : 33.4745976028041 km INFO (Meudon_Bin_NS): omega [rad/s]: 842.992 INFO (Meudon_Bin_NS): dist [km]: 80 INFO (Meudon_Bin_NS): dist_mass [km]: 79.4475 INFO (Meudon_Bin_NS): mass1_b [M_sun]: 1.67079 INFO (Meudon_Bin_NS): mass2_b [M_sun]: 1.28304 INFO (Meudon_Bin_NS): mass_ADM [M_sun]: 2.84045 INFO (Meudon_Bin_NS): L_tot [G M_sun^2/c]: 10.0839 INFO (Meudon_Bin_NS): rad1_x_comp [km]: 31.5924 INFO (Meudon_Bin_NS): rad1_y [km]: 27.7775 INFO (Meudon_Bin_NS): rad1_z [km]: 27.792 INFO (Meudon_Bin_NS): rad1_x_opp [km]: 30.3544 INFO (Meudon_Bin_NS): rad2_x_comp [km]: 36.8228 INFO (Meudon_Bin_NS): rad2_y [km]: 28.1648 INFO (Meudon_Bin_NS): rad2_z [km]: 28.0516 INFO (Meudon_Bin_NS): rad2_x_opp [km]: 33.4746 INFO (Meudon_Bin_NS): K [ET unit]: 372.101387884248 INFO (Meudon_Bin_NS): Filling in Cactus grid points INFO (Meudon_Bin_NS): Done. INFO (Seed_Magnetic_Fields_BNS): Seeding magnetic fields INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 2 NumPts= 517495 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 9.520e-16, ErrDenom: 2.105e-02 | 2.00 iters/gridpt INFO (IllinoisGRMHD): C2P: Lev: 2 NumPts= 517495 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 2.357e-16, ErrDenom: 2.105e-02 | 2.00 iters/gridpt INFO (VolumeIntegrals_GRMHD): Outputting to directory 'volume_integration/' INFO (particle_tracerET): Outputting to directory 'ABEoutput_directory-nsns_test/' INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 2 NumPts= 517495 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 2.202e-15, ErrDenom: 2.105e-02 | 2.00 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (Time): Timestep set to 1.125 (courant_static) INFO (Meudon_Bin_NS): Setting up LORENE Bin_NS initial data INFO (Meudon_Bin_NS): Setting up coordinates INFO (Meudon_Bin_NS): Meudon_Bin_NS will use the 2D_Polytrope equation of state. INFO (Meudon_Bin_NS): Meudon_Bin_NS will use the 1 eos handle INFO (Meudon_Bin_NS): Reading from file "/global/u1/c/chlin/lorene_data/test_gr_09_063.d" Binary system read in file : ---------------------------- Binary system ============= Orbital angular velocity : 842.9922849147580 rad/s Coordinate separation between the two stellar centers : 80.0000000000000 km Absolute coordinate X of the rotation axis : 0.0000000000000 km Star 1 : ====== Relativistic star ----------------- Number of domains occupied by the star : 1 Equation of state : Star 1 EOS EOS of class Eos_poly (relativistic polytrope) : Adiabatic index gamma : 2.0000000000000 Pressure coefficient kappa : 0.1000000000000 rho_nuc c^2 / n_nuc^gamma Mean particle mass : 1.0000000000000 m_B Relativistic chemical potential at zero pressure : 1.0000000000000 m_B c^2 Central enthalpy : 0.0900000000000 c^2 Central proper baryon density : 0.4708714185261 x 0.1 fm^-3 Central proper energy density : 0.4930434078045 rho_nuc c^2 Central pressure : 0.0221719892785 rho_nuc c^2 Regularization index of the gravitational potential : k_div = 0 Central lapse N : 0.8213808112110 Central value of A^2 : 1.4514630408533 Coordinate equatorial radius (phi=0) a1 = 31.5924006423005 km Coordinate equatorial radius (phi=pi/2) a2 = 27.7774975562862 km Coordinate equatorial radius (phi=pi): 30.3543794804130 km Coordinate polar radius a3 = 27.7919576947017 km Axis ratio a2/a1 = 0.8792461791933 a3/a1 = 0.8797038885829 Baryon mass : 1.6707909842769 M_sol Gravitational mass : 1.5865459517723 M_sol Star in a binary system ----------------------- irrotational configuration Absolute abscidia of the stellar center: -34.3158126713210 km Absolute abscidia of the barycenter of the baryon density : -34.1807445379738 km d_tilde : 2.1758549531469 Orientation with respect to the absolute frame : 0.0000000000000 rad Central value of gam_euler : 1.0036741742355 Central u_euler (U^X, U^Y, U^Z) [c] : -0.0000000000000 -0.0709573993086 0.0000000000000 Central d_psi (X, Y, Z) [c] : -0.0000000000000 -0.1131052894769 0.0000000000000 Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.0000000000000 0.0003935153316 0.0000000000000 Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : 0.3718331383477 0.4943138902066 0.0012024646722 Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.3718331383477 -0.2470525597791 -0.0012024646722 Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : 0.0435878798883 -0.0425606873520 -0.0006667894195 Central value of loggam : 0.0000001114528 Central value of log(N) auto, comp : -0.1735698739165 -0.0231985648935 Central value of beta=log(AN) auto, comp : -0.0100589249522 -0.0004234931416 Central value of shift (N^X, N^Y, N^Z) [c] : -0.0000000000000 -0.0378905962741 0.0000000000000 ... shift_auto part of it [c] : -0.0000000000000 -0.0489796785016 0.0000000000000 ... shift_comp part of it [c] : 0.0000000000000 0.0110890822274 0.0000000000000 ... w_shift (NB: components in the star Cartesian frame) [c] : -0.0000000000000 -0.0534744284212 0.0000000000000 Central value of khi_shift [km c] : -0.0000000000000 Central value of (B^X, B^Y, B^Z)/N [c] : -0.0000000000000 0.0713465773081 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : -0.0000362989755 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_comp) [km^{-1}] : -0.0002981883450 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : -0.0000112315610 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_comp) [km^{-1}] : -0.0000163565138 0.0000000000000 0.0000000000000 Central A^2 K^{ij} [c/km] : A^2 K^{xx} auto, comp : -0.0000000000000 0.0000000000000 A^2 K^{xy} auto, comp : -0.0000046453614 -0.0000650115950 A^2 K^{xz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yy} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{zz} auto, comp : 0.0000000000000 0.0000000000000 Central A^2 K_{ij} K^{ij} [c^2/km^2] : A^2 K_{ij} K^{ij} auto, comp : 0.0000000000574 0.0000000008767 Star 2 : ====== Relativistic star ----------------- Number of domains occupied by the star : 1 Equation of state : Star 1 EOS EOS of class Eos_poly (relativistic polytrope) : Adiabatic index gamma : 2.0000000000000 Pressure coefficient kappa : 0.1000000000000 rho_nuc c^2 / n_nuc^gamma Mean particle mass : 1.0000000000000 m_B Relativistic chemical potential at zero pressure : 1.0000000000000 m_B c^2 Central enthalpy : 0.0630000000000 c^2 Central proper baryon density : 0.3251341961565 x 0.1 fm^-3 Central proper energy density : 0.3357054207076 rho_nuc c^2 Central pressure : 0.0105712245510 rho_nuc c^2 Regularization index of the gravitational potential : k_div = 0 Central lapse N : 0.8560863131981 Central value of A^2 : 1.3462211527558 Coordinate equatorial radius (phi=0) a1 = 36.8227739723403 km Coordinate equatorial radius (phi=pi/2) a2 = 28.1647579751040 km Coordinate equatorial radius (phi=pi): 33.4745976028041 km Coordinate polar radius a3 = 28.0515599671080 km Axis ratio a2/a1 = 0.7648733361658 a3/a1 = 0.7617992057899 Baryon mass : 1.2830422471997 M_sol Gravitational mass : 1.2341346226771 M_sol Star in a binary system ----------------------- irrotational configuration Absolute abscidia of the stellar center: 45.6841873286790 km Absolute abscidia of the barycenter of the baryon density : 45.2667839099569 km d_tilde : 2.5042244486690 Orientation with respect to the absolute frame : 3.1415926535898 rad Central value of gam_euler : 1.0072303464165 Central u_euler (U^X, U^Y, U^Z) [c] : -0.0000000000000 0.1030839742961 0.0000000000000 Central d_psi (X, Y, Z) [c] : -0.0000000000000 0.1488664791184 0.0000000000000 Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.0000000000000 -0.0013151835884 0.0000000000000 Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : 0.3604065576231 0.2005959662083 0.0039290006645 Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.3604065576231 -0.5244769501283 -0.0039290006645 Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : -0.0435970713608 0.0424962994979 -0.0014660267039 Central value of loggam : 0.0000011471727 Central value of log(N) auto, comp : -0.1258082259951 -0.0295758487532 Central value of beta=log(AN) auto, comp : -0.0063507225771 -0.0003825914672 Central value of shift (N^X, N^Y, N^Z) [c] : -0.0000000000000 0.0391100446608 0.0000000000000 ... shift_auto part of it [c] : -0.0000000000000 0.0499120259801 0.0000000000000 ... shift_comp part of it [c] : 0.0000000000000 -0.0108019813193 0.0000000000000 ... w_shift (NB: components in the star Cartesian frame) [c] : 0.0000000000000 -0.0546332508576 0.0000000000000 Central value of khi_shift [km c] : -0.0000000000000 Central value of (B^X, B^Y, B^Z)/N [c] : -0.0000000000000 -0.1043705700727 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : 0.0000507123452 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_comp) [km^{-1}] : 0.0003740463855 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : 0.0000113135293 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_comp) [km^{-1}] : 0.0000156642372 0.0000000000000 0.0000000000000 Central A^2 K^{ij} [c/km] : A^2 K^{xx} auto, comp : -0.0000000000000 0.0000000000000 A^2 K^{xy} auto, comp : -0.0000046934496 -0.0000635848483 A^2 K^{xz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yy} auto, comp : 0.0000000000000 -0.0000000000000 A^2 K^{yz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{zz} auto, comp : -0.0000000000000 0.0000000000000 Central A^2 K_{ij} K^{ij} [c^2/km^2] : A^2 K_{ij} K^{ij} auto, comp : 0.0000000000497 0.0000000008035 Summary : ------- Separation d : 80.0000000000000 km Separation d_G : 79.4475284479308 km Omega : 842.9922849147580 rad/s Baryon mass of star 1 : 1.6707909842769 M_sol Baryon mass of star 2 : 1.2830422471997 M_sol ADM mass of the system : 2.8404508262857 M_sol Total angular momentum : 10.0839081191459 G M_sol^2 / c Radius of star 1 (x_comp) : 31.5924006423005 km Radius of star 1 (y) : 27.7774975562862 km Radius of star 1 (z) : 27.7919576947017 km Radius of star 1 (x_opp) : 30.3543794804130 km Radius of star 2 (x_comp) : 36.8227739723403 km Radius of star 2 (y) : 28.1647579751040 km Radius of star 2 (z) : 28.0515599671080 km Radius of star 2 (x_opp) : 33.4745976028041 km INFO (Meudon_Bin_NS): omega [rad/s]: 842.992 INFO (Meudon_Bin_NS): dist [km]: 80 INFO (Meudon_Bin_NS): dist_mass [km]: 79.4475 INFO (Meudon_Bin_NS): mass1_b [M_sun]: 1.67079 INFO (Meudon_Bin_NS): mass2_b [M_sun]: 1.28304 INFO (Meudon_Bin_NS): mass_ADM [M_sun]: 2.84045 INFO (Meudon_Bin_NS): L_tot [G M_sun^2/c]: 10.0839 INFO (Meudon_Bin_NS): rad1_x_comp [km]: 31.5924 INFO (Meudon_Bin_NS): rad1_y [km]: 27.7775 INFO (Meudon_Bin_NS): rad1_z [km]: 27.792 INFO (Meudon_Bin_NS): rad1_x_opp [km]: 30.3544 INFO (Meudon_Bin_NS): rad2_x_comp [km]: 36.8228 INFO (Meudon_Bin_NS): rad2_y [km]: 28.1648 INFO (Meudon_Bin_NS): rad2_z [km]: 28.0516 INFO (Meudon_Bin_NS): rad2_x_opp [km]: 33.4746 INFO (Meudon_Bin_NS): K [ET unit]: 372.101387884248 INFO (Meudon_Bin_NS): Filling in Cactus grid points INFO (Meudon_Bin_NS): Done. INFO (Seed_Magnetic_Fields_BNS): Seeding magnetic fields INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 3 NumPts= 772265 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 4.245e-16, ErrDenom: 1.642e-01 | 2.00 iters/gridpt INFO (IllinoisGRMHD): C2P: Lev: 3 NumPts= 772265 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 2.270e-16, ErrDenom: 1.642e-01 | 2.00 iters/gridpt INFO (VolumeIntegrals_GRMHD): Outputting to directory 'volume_integration/' INFO (particle_tracerET): Outputting to directory 'ABEoutput_directory-nsns_test/' INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 3 NumPts= 772265 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 5.998e-16, ErrDenom: 1.642e-01 | 2.00 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (Time): Timestep set to 0.5625 (courant_static) INFO (Meudon_Bin_NS): Setting up LORENE Bin_NS initial data INFO (Meudon_Bin_NS): Setting up coordinates INFO (Meudon_Bin_NS): Meudon_Bin_NS will use the 2D_Polytrope equation of state. INFO (Meudon_Bin_NS): Meudon_Bin_NS will use the 1 eos handle INFO (Meudon_Bin_NS): Reading from file "/global/u1/c/chlin/lorene_data/test_gr_09_063.d" Binary system read in file : ---------------------------- Binary system ============= Orbital angular velocity : 842.9922849147580 rad/s Coordinate separation between the two stellar centers : 80.0000000000000 km Absolute coordinate X of the rotation axis : 0.0000000000000 km Star 1 : ====== Relativistic star ----------------- Number of domains occupied by the star : 1 Equation of state : Star 1 EOS EOS of class Eos_poly (relativistic polytrope) : Adiabatic index gamma : 2.0000000000000 Pressure coefficient kappa : 0.1000000000000 rho_nuc c^2 / n_nuc^gamma Mean particle mass : 1.0000000000000 m_B Relativistic chemical potential at zero pressure : 1.0000000000000 m_B c^2 Central enthalpy : 0.0900000000000 c^2 Central proper baryon density : 0.4708714185261 x 0.1 fm^-3 Central proper energy density : 0.4930434078045 rho_nuc c^2 Central pressure : 0.0221719892785 rho_nuc c^2 Regularization index of the gravitational potential : k_div = 0 Central lapse N : 0.8213808112110 Central value of A^2 : 1.4514630408533 Coordinate equatorial radius (phi=0) a1 = 31.5924006423005 km Coordinate equatorial radius (phi=pi/2) a2 = 27.7774975562862 km Coordinate equatorial radius (phi=pi): 30.3543794804130 km Coordinate polar radius a3 = 27.7919576947017 km Axis ratio a2/a1 = 0.8792461791933 a3/a1 = 0.8797038885829 Baryon mass : 1.6707909842769 M_sol Gravitational mass : 1.5865459517723 M_sol Star in a binary system ----------------------- irrotational configuration Absolute abscidia of the stellar center: -34.3158126713210 km Absolute abscidia of the barycenter of the baryon density : -34.1807445379738 km d_tilde : 2.1758549531469 Orientation with respect to the absolute frame : 0.0000000000000 rad Central value of gam_euler : 1.0036741742355 Central u_euler (U^X, U^Y, U^Z) [c] : -0.0000000000000 -0.0709573993086 0.0000000000000 Central d_psi (X, Y, Z) [c] : -0.0000000000000 -0.1131052894769 0.0000000000000 Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.0000000000000 0.0003935153316 0.0000000000000 Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : 0.3718331383477 0.4943138902066 0.0012024646722 Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.3718331383477 -0.2470525597791 -0.0012024646722 Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : 0.0435878798883 -0.0425606873520 -0.0006667894195 Central value of loggam : 0.0000001114528 Central value of log(N) auto, comp : -0.1735698739165 -0.0231985648935 Central value of beta=log(AN) auto, comp : -0.0100589249522 -0.0004234931416 Central value of shift (N^X, N^Y, N^Z) [c] : -0.0000000000000 -0.0378905962741 0.0000000000000 ... shift_auto part of it [c] : -0.0000000000000 -0.0489796785016 0.0000000000000 ... shift_comp part of it [c] : 0.0000000000000 0.0110890822274 0.0000000000000 ... w_shift (NB: components in the star Cartesian frame) [c] : -0.0000000000000 -0.0534744284212 0.0000000000000 Central value of khi_shift [km c] : -0.0000000000000 Central value of (B^X, B^Y, B^Z)/N [c] : -0.0000000000000 0.0713465773081 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : -0.0000362989755 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_comp) [km^{-1}] : -0.0002981883450 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : -0.0000112315610 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_comp) [km^{-1}] : -0.0000163565138 0.0000000000000 0.0000000000000 Central A^2 K^{ij} [c/km] : A^2 K^{xx} auto, comp : -0.0000000000000 0.0000000000000 A^2 K^{xy} auto, comp : -0.0000046453614 -0.0000650115950 A^2 K^{xz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yy} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{zz} auto, comp : 0.0000000000000 0.0000000000000 Central A^2 K_{ij} K^{ij} [c^2/km^2] : A^2 K_{ij} K^{ij} auto, comp : 0.0000000000574 0.0000000008767 Star 2 : ====== Relativistic star ----------------- Number of domains occupied by the star : 1 Equation of state : Star 1 EOS EOS of class Eos_poly (relativistic polytrope) : Adiabatic index gamma : 2.0000000000000 Pressure coefficient kappa : 0.1000000000000 rho_nuc c^2 / n_nuc^gamma Mean particle mass : 1.0000000000000 m_B Relativistic chemical potential at zero pressure : 1.0000000000000 m_B c^2 Central enthalpy : 0.0630000000000 c^2 Central proper baryon density : 0.3251341961565 x 0.1 fm^-3 Central proper energy density : 0.3357054207076 rho_nuc c^2 Central pressure : 0.0105712245510 rho_nuc c^2 Regularization index of the gravitational potential : k_div = 0 Central lapse N : 0.8560863131981 Central value of A^2 : 1.3462211527558 Coordinate equatorial radius (phi=0) a1 = 36.8227739723403 km Coordinate equatorial radius (phi=pi/2) a2 = 28.1647579751040 km Coordinate equatorial radius (phi=pi): 33.4745976028041 km Coordinate polar radius a3 = 28.0515599671080 km Axis ratio a2/a1 = 0.7648733361658 a3/a1 = 0.7617992057899 Baryon mass : 1.2830422471997 M_sol Gravitational mass : 1.2341346226771 M_sol Star in a binary system ----------------------- irrotational configuration Absolute abscidia of the stellar center: 45.6841873286790 km Absolute abscidia of the barycenter of the baryon density : 45.2667839099569 km d_tilde : 2.5042244486690 Orientation with respect to the absolute frame : 3.1415926535898 rad Central value of gam_euler : 1.0072303464165 Central u_euler (U^X, U^Y, U^Z) [c] : -0.0000000000000 0.1030839742961 0.0000000000000 Central d_psi (X, Y, Z) [c] : -0.0000000000000 0.1488664791184 0.0000000000000 Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.0000000000000 -0.0013151835884 0.0000000000000 Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : 0.3604065576231 0.2005959662083 0.0039290006645 Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.3604065576231 -0.5244769501283 -0.0039290006645 Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : -0.0435970713608 0.0424962994979 -0.0014660267039 Central value of loggam : 0.0000011471727 Central value of log(N) auto, comp : -0.1258082259951 -0.0295758487532 Central value of beta=log(AN) auto, comp : -0.0063507225771 -0.0003825914672 Central value of shift (N^X, N^Y, N^Z) [c] : -0.0000000000000 0.0391100446608 0.0000000000000 ... shift_auto part of it [c] : -0.0000000000000 0.0499120259801 0.0000000000000 ... shift_comp part of it [c] : 0.0000000000000 -0.0108019813193 0.0000000000000 ... w_shift (NB: components in the star Cartesian frame) [c] : 0.0000000000000 -0.0546332508576 0.0000000000000 Central value of khi_shift [km c] : -0.0000000000000 Central value of (B^X, B^Y, B^Z)/N [c] : -0.0000000000000 -0.1043705700727 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : 0.0000507123452 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_comp) [km^{-1}] : 0.0003740463855 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : 0.0000113135293 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_comp) [km^{-1}] : 0.0000156642372 0.0000000000000 0.0000000000000 Central A^2 K^{ij} [c/km] : A^2 K^{xx} auto, comp : -0.0000000000000 0.0000000000000 A^2 K^{xy} auto, comp : -0.0000046934496 -0.0000635848483 A^2 K^{xz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yy} auto, comp : 0.0000000000000 -0.0000000000000 A^2 K^{yz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{zz} auto, comp : -0.0000000000000 0.0000000000000 Central A^2 K_{ij} K^{ij} [c^2/km^2] : A^2 K_{ij} K^{ij} auto, comp : 0.0000000000497 0.0000000008035 Summary : ------- Separation d : 80.0000000000000 km Separation d_G : 79.4475284479308 km Omega : 842.9922849147580 rad/s Baryon mass of star 1 : 1.6707909842769 M_sol Baryon mass of star 2 : 1.2830422471997 M_sol ADM mass of the system : 2.8404508262857 M_sol Total angular momentum : 10.0839081191459 G M_sol^2 / c Radius of star 1 (x_comp) : 31.5924006423005 km Radius of star 1 (y) : 27.7774975562862 km Radius of star 1 (z) : 27.7919576947017 km Radius of star 1 (x_opp) : 30.3543794804130 km Radius of star 2 (x_comp) : 36.8227739723403 km Radius of star 2 (y) : 28.1647579751040 km Radius of star 2 (z) : 28.0515599671080 km Radius of star 2 (x_opp) : 33.4745976028041 km INFO (Meudon_Bin_NS): omega [rad/s]: 842.992 INFO (Meudon_Bin_NS): dist [km]: 80 INFO (Meudon_Bin_NS): dist_mass [km]: 79.4475 INFO (Meudon_Bin_NS): mass1_b [M_sun]: 1.67079 INFO (Meudon_Bin_NS): mass2_b [M_sun]: 1.28304 INFO (Meudon_Bin_NS): mass_ADM [M_sun]: 2.84045 INFO (Meudon_Bin_NS): L_tot [G M_sun^2/c]: 10.0839 INFO (Meudon_Bin_NS): rad1_x_comp [km]: 31.5924 INFO (Meudon_Bin_NS): rad1_y [km]: 27.7775 INFO (Meudon_Bin_NS): rad1_z [km]: 27.792 INFO (Meudon_Bin_NS): rad1_x_opp [km]: 30.3544 INFO (Meudon_Bin_NS): rad2_x_comp [km]: 36.8228 INFO (Meudon_Bin_NS): rad2_y [km]: 28.1648 INFO (Meudon_Bin_NS): rad2_z [km]: 28.0516 INFO (Meudon_Bin_NS): rad2_x_opp [km]: 33.4746 INFO (Meudon_Bin_NS): K [ET unit]: 372.101387884248 INFO (Meudon_Bin_NS): Filling in Cactus grid points INFO (Meudon_Bin_NS): Done. INFO (Seed_Magnetic_Fields_BNS): Seeding magnetic fields INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 4 NumPts= 620994 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 5.032e-16, ErrDenom: 1.313e+00 | 2.03 iters/gridpt INFO (IllinoisGRMHD): C2P: Lev: 4 NumPts= 620994 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 2.435e-16, ErrDenom: 1.313e+00 | 2.03 iters/gridpt INFO (VolumeIntegrals_GRMHD): Outputting to directory 'volume_integration/' INFO (particle_tracerET): Outputting to directory 'ABEoutput_directory-nsns_test/' INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 4 NumPts= 620994 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 7.979e-16, ErrDenom: 1.313e+00 | 2.03 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (Time): Timestep set to 0.28125 (courant_static) INFO (Meudon_Bin_NS): Setting up LORENE Bin_NS initial data INFO (Meudon_Bin_NS): Setting up coordinates INFO (Meudon_Bin_NS): Meudon_Bin_NS will use the 2D_Polytrope equation of state. INFO (Meudon_Bin_NS): Meudon_Bin_NS will use the 1 eos handle INFO (Meudon_Bin_NS): Reading from file "/global/u1/c/chlin/lorene_data/test_gr_09_063.d" Binary system read in file : ---------------------------- Binary system ============= Orbital angular velocity : 842.9922849147580 rad/s Coordinate separation between the two stellar centers : 80.0000000000000 km Absolute coordinate X of the rotation axis : 0.0000000000000 km Star 1 : ====== Relativistic star ----------------- Number of domains occupied by the star : 1 Equation of state : Star 1 EOS EOS of class Eos_poly (relativistic polytrope) : Adiabatic index gamma : 2.0000000000000 Pressure coefficient kappa : 0.1000000000000 rho_nuc c^2 / n_nuc^gamma Mean particle mass : 1.0000000000000 m_B Relativistic chemical potential at zero pressure : 1.0000000000000 m_B c^2 Central enthalpy : 0.0900000000000 c^2 Central proper baryon density : 0.4708714185261 x 0.1 fm^-3 Central proper energy density : 0.4930434078045 rho_nuc c^2 Central pressure : 0.0221719892785 rho_nuc c^2 Regularization index of the gravitational potential : k_div = 0 Central lapse N : 0.8213808112110 Central value of A^2 : 1.4514630408533 Coordinate equatorial radius (phi=0) a1 = 31.5924006423005 km Coordinate equatorial radius (phi=pi/2) a2 = 27.7774975562862 km Coordinate equatorial radius (phi=pi): 30.3543794804130 km Coordinate polar radius a3 = 27.7919576947017 km Axis ratio a2/a1 = 0.8792461791933 a3/a1 = 0.8797038885829 Baryon mass : 1.6707909842769 M_sol Gravitational mass : 1.5865459517723 M_sol Star in a binary system ----------------------- irrotational configuration Absolute abscidia of the stellar center: -34.3158126713210 km Absolute abscidia of the barycenter of the baryon density : -34.1807445379738 km d_tilde : 2.1758549531469 Orientation with respect to the absolute frame : 0.0000000000000 rad Central value of gam_euler : 1.0036741742355 Central u_euler (U^X, U^Y, U^Z) [c] : -0.0000000000000 -0.0709573993086 0.0000000000000 Central d_psi (X, Y, Z) [c] : -0.0000000000000 -0.1131052894769 0.0000000000000 Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.0000000000000 0.0003935153316 0.0000000000000 Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : 0.3718331383477 0.4943138902066 0.0012024646722 Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.3718331383477 -0.2470525597791 -0.0012024646722 Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : 0.0435878798883 -0.0425606873520 -0.0006667894195 Central value of loggam : 0.0000001114528 Central value of log(N) auto, comp : -0.1735698739165 -0.0231985648935 Central value of beta=log(AN) auto, comp : -0.0100589249522 -0.0004234931416 Central value of shift (N^X, N^Y, N^Z) [c] : -0.0000000000000 -0.0378905962741 0.0000000000000 ... shift_auto part of it [c] : -0.0000000000000 -0.0489796785016 0.0000000000000 ... shift_comp part of it [c] : 0.0000000000000 0.0110890822274 0.0000000000000 ... w_shift (NB: components in the star Cartesian frame) [c] : -0.0000000000000 -0.0534744284212 0.0000000000000 Central value of khi_shift [km c] : -0.0000000000000 Central value of (B^X, B^Y, B^Z)/N [c] : -0.0000000000000 0.0713465773081 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : -0.0000362989755 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_comp) [km^{-1}] : -0.0002981883450 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : -0.0000112315610 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_comp) [km^{-1}] : -0.0000163565138 0.0000000000000 0.0000000000000 Central A^2 K^{ij} [c/km] : A^2 K^{xx} auto, comp : -0.0000000000000 0.0000000000000 A^2 K^{xy} auto, comp : -0.0000046453614 -0.0000650115950 A^2 K^{xz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yy} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{zz} auto, comp : 0.0000000000000 0.0000000000000 Central A^2 K_{ij} K^{ij} [c^2/km^2] : A^2 K_{ij} K^{ij} auto, comp : 0.0000000000574 0.0000000008767 Star 2 : ====== Relativistic star ----------------- Number of domains occupied by the star : 1 Equation of state : Star 1 EOS EOS of class Eos_poly (relativistic polytrope) : Adiabatic index gamma : 2.0000000000000 Pressure coefficient kappa : 0.1000000000000 rho_nuc c^2 / n_nuc^gamma Mean particle mass : 1.0000000000000 m_B Relativistic chemical potential at zero pressure : 1.0000000000000 m_B c^2 Central enthalpy : 0.0630000000000 c^2 Central proper baryon density : 0.3251341961565 x 0.1 fm^-3 Central proper energy density : 0.3357054207076 rho_nuc c^2 Central pressure : 0.0105712245510 rho_nuc c^2 Regularization index of the gravitational potential : k_div = 0 Central lapse N : 0.8560863131981 Central value of A^2 : 1.3462211527558 Coordinate equatorial radius (phi=0) a1 = 36.8227739723403 km Coordinate equatorial radius (phi=pi/2) a2 = 28.1647579751040 km Coordinate equatorial radius (phi=pi): 33.4745976028041 km Coordinate polar radius a3 = 28.0515599671080 km Axis ratio a2/a1 = 0.7648733361658 a3/a1 = 0.7617992057899 Baryon mass : 1.2830422471997 M_sol Gravitational mass : 1.2341346226771 M_sol Star in a binary system ----------------------- irrotational configuration Absolute abscidia of the stellar center: 45.6841873286790 km Absolute abscidia of the barycenter of the baryon density : 45.2667839099569 km d_tilde : 2.5042244486690 Orientation with respect to the absolute frame : 3.1415926535898 rad Central value of gam_euler : 1.0072303464165 Central u_euler (U^X, U^Y, U^Z) [c] : -0.0000000000000 0.1030839742961 0.0000000000000 Central d_psi (X, Y, Z) [c] : -0.0000000000000 0.1488664791184 0.0000000000000 Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.0000000000000 -0.0013151835884 0.0000000000000 Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : 0.3604065576231 0.2005959662083 0.0039290006645 Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.3604065576231 -0.5244769501283 -0.0039290006645 Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : -0.0435970713608 0.0424962994979 -0.0014660267039 Central value of loggam : 0.0000011471727 Central value of log(N) auto, comp : -0.1258082259951 -0.0295758487532 Central value of beta=log(AN) auto, comp : -0.0063507225771 -0.0003825914672 Central value of shift (N^X, N^Y, N^Z) [c] : -0.0000000000000 0.0391100446608 0.0000000000000 ... shift_auto part of it [c] : -0.0000000000000 0.0499120259801 0.0000000000000 ... shift_comp part of it [c] : 0.0000000000000 -0.0108019813193 0.0000000000000 ... w_shift (NB: components in the star Cartesian frame) [c] : 0.0000000000000 -0.0546332508576 0.0000000000000 Central value of khi_shift [km c] : -0.0000000000000 Central value of (B^X, B^Y, B^Z)/N [c] : -0.0000000000000 -0.1043705700727 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : 0.0000507123452 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_comp) [km^{-1}] : 0.0003740463855 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : 0.0000113135293 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_comp) [km^{-1}] : 0.0000156642372 0.0000000000000 0.0000000000000 Central A^2 K^{ij} [c/km] : A^2 K^{xx} auto, comp : -0.0000000000000 0.0000000000000 A^2 K^{xy} auto, comp : -0.0000046934496 -0.0000635848483 A^2 K^{xz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yy} auto, comp : 0.0000000000000 -0.0000000000000 A^2 K^{yz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{zz} auto, comp : -0.0000000000000 0.0000000000000 Central A^2 K_{ij} K^{ij} [c^2/km^2] : A^2 K_{ij} K^{ij} auto, comp : 0.0000000000497 0.0000000008035 Summary : ------- Separation d : 80.0000000000000 km Separation d_G : 79.4475284479308 km Omega : 842.9922849147580 rad/s Baryon mass of star 1 : 1.6707909842769 M_sol Baryon mass of star 2 : 1.2830422471997 M_sol ADM mass of the system : 2.8404508262857 M_sol Total angular momentum : 10.0839081191459 G M_sol^2 / c Radius of star 1 (x_comp) : 31.5924006423005 km Radius of star 1 (y) : 27.7774975562862 km Radius of star 1 (z) : 27.7919576947017 km Radius of star 1 (x_opp) : 30.3543794804130 km Radius of star 2 (x_comp) : 36.8227739723403 km Radius of star 2 (y) : 28.1647579751040 km Radius of star 2 (z) : 28.0515599671080 km Radius of star 2 (x_opp) : 33.4745976028041 km INFO (Meudon_Bin_NS): omega [rad/s]: 842.992 INFO (Meudon_Bin_NS): dist [km]: 80 INFO (Meudon_Bin_NS): dist_mass [km]: 79.4475 INFO (Meudon_Bin_NS): mass1_b [M_sun]: 1.67079 INFO (Meudon_Bin_NS): mass2_b [M_sun]: 1.28304 INFO (Meudon_Bin_NS): mass_ADM [M_sun]: 2.84045 INFO (Meudon_Bin_NS): L_tot [G M_sun^2/c]: 10.0839 INFO (Meudon_Bin_NS): rad1_x_comp [km]: 31.5924 INFO (Meudon_Bin_NS): rad1_y [km]: 27.7775 INFO (Meudon_Bin_NS): rad1_z [km]: 27.792 INFO (Meudon_Bin_NS): rad1_x_opp [km]: 30.3544 INFO (Meudon_Bin_NS): rad2_x_comp [km]: 36.8228 INFO (Meudon_Bin_NS): rad2_y [km]: 28.1648 INFO (Meudon_Bin_NS): rad2_z [km]: 28.0516 INFO (Meudon_Bin_NS): rad2_x_opp [km]: 33.4746 INFO (Meudon_Bin_NS): K [ET unit]: 372.101387884248 INFO (Meudon_Bin_NS): Filling in Cactus grid points INFO (Meudon_Bin_NS): Done. INFO (Seed_Magnetic_Fields_BNS): Seeding magnetic fields INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 5 NumPts= 1332331 | Fixes: Font= 0 VL= 0 rho*= 1077 | Failures: 0 InHoriz= 0 / 0 | Error: 1.908e-05, ErrDenom: 1.050e+01 | 2.13 iters/gridpt INFO (IllinoisGRMHD): C2P: Lev: 5 NumPts= 1332331 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 2.407e-16, ErrDenom: 1.050e+01 | 2.13 iters/gridpt INFO (VolumeIntegrals_GRMHD): Outputting to directory 'volume_integration/' INFO (particle_tracerET): Outputting to directory 'ABEoutput_directory-nsns_test/' INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 5 NumPts= 1332331 | Fixes: Font= 0 VL= 0 rho*= 1077 | Failures: 0 InHoriz= 0 / 0 | Error: 1.908e-05, ErrDenom: 1.050e+01 | 2.13 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (Time): Timestep set to 0.140625 (courant_static) INFO (Meudon_Bin_NS): Setting up LORENE Bin_NS initial data INFO (Meudon_Bin_NS): Setting up coordinates INFO (Meudon_Bin_NS): Meudon_Bin_NS will use the 2D_Polytrope equation of state. INFO (Meudon_Bin_NS): Meudon_Bin_NS will use the 1 eos handle INFO (Meudon_Bin_NS): Reading from file "/global/u1/c/chlin/lorene_data/test_gr_09_063.d" Binary system read in file : ---------------------------- Binary system ============= Orbital angular velocity : 842.9922849147580 rad/s Coordinate separation between the two stellar centers : 80.0000000000000 km Absolute coordinate X of the rotation axis : 0.0000000000000 km Star 1 : ====== Relativistic star ----------------- Number of domains occupied by the star : 1 Equation of state : Star 1 EOS EOS of class Eos_poly (relativistic polytrope) : Adiabatic index gamma : 2.0000000000000 Pressure coefficient kappa : 0.1000000000000 rho_nuc c^2 / n_nuc^gamma Mean particle mass : 1.0000000000000 m_B Relativistic chemical potential at zero pressure : 1.0000000000000 m_B c^2 Central enthalpy : 0.0900000000000 c^2 Central proper baryon density : 0.4708714185261 x 0.1 fm^-3 Central proper energy density : 0.4930434078045 rho_nuc c^2 Central pressure : 0.0221719892785 rho_nuc c^2 Regularization index of the gravitational potential : k_div = 0 Central lapse N : 0.8213808112110 Central value of A^2 : 1.4514630408533 Coordinate equatorial radius (phi=0) a1 = 31.5924006423005 km Coordinate equatorial radius (phi=pi/2) a2 = 27.7774975562862 km Coordinate equatorial radius (phi=pi): 30.3543794804130 km Coordinate polar radius a3 = 27.7919576947017 km Axis ratio a2/a1 = 0.8792461791933 a3/a1 = 0.8797038885829 Baryon mass : 1.6707909842769 M_sol Gravitational mass : 1.5865459517723 M_sol Star in a binary system ----------------------- irrotational configuration Absolute abscidia of the stellar center: -34.3158126713210 km Absolute abscidia of the barycenter of the baryon density : -34.1807445379738 km d_tilde : 2.1758549531469 Orientation with respect to the absolute frame : 0.0000000000000 rad Central value of gam_euler : 1.0036741742355 Central u_euler (U^X, U^Y, U^Z) [c] : -0.0000000000000 -0.0709573993086 0.0000000000000 Central d_psi (X, Y, Z) [c] : -0.0000000000000 -0.1131052894769 0.0000000000000 Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.0000000000000 0.0003935153316 0.0000000000000 Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : 0.3718331383477 0.4943138902066 0.0012024646722 Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.3718331383477 -0.2470525597791 -0.0012024646722 Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : 0.0435878798883 -0.0425606873520 -0.0006667894195 Central value of loggam : 0.0000001114528 Central value of log(N) auto, comp : -0.1735698739165 -0.0231985648935 Central value of beta=log(AN) auto, comp : -0.0100589249522 -0.0004234931416 Central value of shift (N^X, N^Y, N^Z) [c] : -0.0000000000000 -0.0378905962741 0.0000000000000 ... shift_auto part of it [c] : -0.0000000000000 -0.0489796785016 0.0000000000000 ... shift_comp part of it [c] : 0.0000000000000 0.0110890822274 0.0000000000000 ... w_shift (NB: components in the star Cartesian frame) [c] : -0.0000000000000 -0.0534744284212 0.0000000000000 Central value of khi_shift [km c] : -0.0000000000000 Central value of (B^X, B^Y, B^Z)/N [c] : -0.0000000000000 0.0713465773081 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : -0.0000362989755 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_comp) [km^{-1}] : -0.0002981883450 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : -0.0000112315610 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_comp) [km^{-1}] : -0.0000163565138 0.0000000000000 0.0000000000000 Central A^2 K^{ij} [c/km] : A^2 K^{xx} auto, comp : -0.0000000000000 0.0000000000000 A^2 K^{xy} auto, comp : -0.0000046453614 -0.0000650115950 A^2 K^{xz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yy} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{zz} auto, comp : 0.0000000000000 0.0000000000000 Central A^2 K_{ij} K^{ij} [c^2/km^2] : A^2 K_{ij} K^{ij} auto, comp : 0.0000000000574 0.0000000008767 Star 2 : ====== Relativistic star ----------------- Number of domains occupied by the star : 1 Equation of state : Star 1 EOS EOS of class Eos_poly (relativistic polytrope) : Adiabatic index gamma : 2.0000000000000 Pressure coefficient kappa : 0.1000000000000 rho_nuc c^2 / n_nuc^gamma Mean particle mass : 1.0000000000000 m_B Relativistic chemical potential at zero pressure : 1.0000000000000 m_B c^2 Central enthalpy : 0.0630000000000 c^2 Central proper baryon density : 0.3251341961565 x 0.1 fm^-3 Central proper energy density : 0.3357054207076 rho_nuc c^2 Central pressure : 0.0105712245510 rho_nuc c^2 Regularization index of the gravitational potential : k_div = 0 Central lapse N : 0.8560863131981 Central value of A^2 : 1.3462211527558 Coordinate equatorial radius (phi=0) a1 = 36.8227739723403 km Coordinate equatorial radius (phi=pi/2) a2 = 28.1647579751040 km Coordinate equatorial radius (phi=pi): 33.4745976028041 km Coordinate polar radius a3 = 28.0515599671080 km Axis ratio a2/a1 = 0.7648733361658 a3/a1 = 0.7617992057899 Baryon mass : 1.2830422471997 M_sol Gravitational mass : 1.2341346226771 M_sol Star in a binary system ----------------------- irrotational configuration Absolute abscidia of the stellar center: 45.6841873286790 km Absolute abscidia of the barycenter of the baryon density : 45.2667839099569 km d_tilde : 2.5042244486690 Orientation with respect to the absolute frame : 3.1415926535898 rad Central value of gam_euler : 1.0072303464165 Central u_euler (U^X, U^Y, U^Z) [c] : -0.0000000000000 0.1030839742961 0.0000000000000 Central d_psi (X, Y, Z) [c] : -0.0000000000000 0.1488664791184 0.0000000000000 Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.0000000000000 -0.0013151835884 0.0000000000000 Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : 0.3604065576231 0.2005959662083 0.0039290006645 Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.3604065576231 -0.5244769501283 -0.0039290006645 Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : -0.0435970713608 0.0424962994979 -0.0014660267039 Central value of loggam : 0.0000011471727 Central value of log(N) auto, comp : -0.1258082259951 -0.0295758487532 Central value of beta=log(AN) auto, comp : -0.0063507225771 -0.0003825914672 Central value of shift (N^X, N^Y, N^Z) [c] : -0.0000000000000 0.0391100446608 0.0000000000000 ... shift_auto part of it [c] : -0.0000000000000 0.0499120259801 0.0000000000000 ... shift_comp part of it [c] : 0.0000000000000 -0.0108019813193 0.0000000000000 ... w_shift (NB: components in the star Cartesian frame) [c] : 0.0000000000000 -0.0546332508576 0.0000000000000 Central value of khi_shift [km c] : -0.0000000000000 Central value of (B^X, B^Y, B^Z)/N [c] : -0.0000000000000 -0.1043705700727 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : 0.0000507123452 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(logn_comp) [km^{-1}] : 0.0003740463855 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : 0.0000113135293 0.0000000000000 0.0000000000000 Central (d/dX,d/dY,d/dZ)(beta_comp) [km^{-1}] : 0.0000156642372 0.0000000000000 0.0000000000000 Central A^2 K^{ij} [c/km] : A^2 K^{xx} auto, comp : -0.0000000000000 0.0000000000000 A^2 K^{xy} auto, comp : -0.0000046934496 -0.0000635848483 A^2 K^{xz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{yy} auto, comp : 0.0000000000000 -0.0000000000000 A^2 K^{yz} auto, comp : 0.0000000000000 0.0000000000000 A^2 K^{zz} auto, comp : -0.0000000000000 0.0000000000000 Central A^2 K_{ij} K^{ij} [c^2/km^2] : A^2 K_{ij} K^{ij} auto, comp : 0.0000000000497 0.0000000008035 Summary : ------- Separation d : 80.0000000000000 km Separation d_G : 79.4475284479308 km Omega : 842.9922849147580 rad/s Baryon mass of star 1 : 1.6707909842769 M_sol Baryon mass of star 2 : 1.2830422471997 M_sol ADM mass of the system : 2.8404508262857 M_sol Total angular momentum : 10.0839081191459 G M_sol^2 / c Radius of star 1 (x_comp) : 31.5924006423005 km Radius of star 1 (y) : 27.7774975562862 km Radius of star 1 (z) : 27.7919576947017 km Radius of star 1 (x_opp) : 30.3543794804130 km Radius of star 2 (x_comp) : 36.8227739723403 km Radius of star 2 (y) : 28.1647579751040 km Radius of star 2 (z) : 28.0515599671080 km Radius of star 2 (x_opp) : 33.4745976028041 km INFO (Meudon_Bin_NS): omega [rad/s]: 842.992 INFO (Meudon_Bin_NS): dist [km]: 80 INFO (Meudon_Bin_NS): dist_mass [km]: 79.4475 INFO (Meudon_Bin_NS): mass1_b [M_sun]: 1.67079 INFO (Meudon_Bin_NS): mass2_b [M_sun]: 1.28304 INFO (Meudon_Bin_NS): mass_ADM [M_sun]: 2.84045 INFO (Meudon_Bin_NS): L_tot [G M_sun^2/c]: 10.0839 INFO (Meudon_Bin_NS): rad1_x_comp [km]: 31.5924 INFO (Meudon_Bin_NS): rad1_y [km]: 27.7775 INFO (Meudon_Bin_NS): rad1_z [km]: 27.792 INFO (Meudon_Bin_NS): rad1_x_opp [km]: 30.3544 INFO (Meudon_Bin_NS): rad2_x_comp [km]: 36.8228 INFO (Meudon_Bin_NS): rad2_y [km]: 28.1648 INFO (Meudon_Bin_NS): rad2_z [km]: 28.0516 INFO (Meudon_Bin_NS): rad2_x_opp [km]: 33.4746 INFO (Meudon_Bin_NS): K [ET unit]: 372.101387884248 INFO (Meudon_Bin_NS): Filling in Cactus grid points INFO (Meudon_Bin_NS): Done. INFO (Seed_Magnetic_Fields_BNS): Seeding magnetic fields INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 6 NumPts= 1727638 | Fixes: Font= 0 VL= 0 rho*= 25462 | Failures: 0 InHoriz= 0 / 0 | Error: 3.116e-05, ErrDenom: 7.426e+01 | 2.62 iters/gridpt INFO (IllinoisGRMHD): C2P: Lev: 6 NumPts= 1727638 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 2.468e-16, ErrDenom: 7.427e+01 | 2.63 iters/gridpt INFO (VolumeIntegrals_GRMHD): Outputting to directory 'volume_integration/' INFO (particle_tracerET): Outputting to directory 'ABEoutput_directory-nsns_test/' INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 6 NumPts= 1727638 | Fixes: Font= 0 VL= 0 rho*= 25462 | Failures: 0 InHoriz= 0 / 0 | Error: 3.116e-05, ErrDenom: 7.426e+01 | 2.62 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 0 NumPts= 152100 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 3.172e-16, ErrDenom: 5.571e-04 | 2.00 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 1 NumPts= 346800 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 8.529e-15, ErrDenom: 3.667e-03 | 2.00 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 2 NumPts= 517495 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 2.194e-15, ErrDenom: 2.105e-02 | 2.00 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 3 NumPts= 772265 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 5.566e-16, ErrDenom: 1.642e-01 | 2.00 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 4 NumPts= 620994 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 7.661e-16, ErrDenom: 1.313e+00 | 2.03 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 5 NumPts= 1332331 | Fixes: Font= 0 VL= 0 rho*= 1077 | Failures: 0 InHoriz= 0 / 0 | Error: 1.908e-05, ErrDenom: 1.050e+01 | 2.13 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 6 NumPts= 1727638 | Fixes: Font= 0 VL= 0 rho*= 25462 | Failures: 0 InHoriz= 0 / 0 | Error: 3.116e-05, ErrDenom: 7.426e+01 | 2.62 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 0 NumPts= 152100 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 1.752e-16, ErrDenom: 5.571e-04 | 2.00 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 1 NumPts= 346800 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 8.429e-15, ErrDenom: 3.667e-03 | 2.00 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 2 NumPts= 517495 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 2.160e-15, ErrDenom: 2.105e-02 | 2.00 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 3 NumPts= 772265 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 5.370e-16, ErrDenom: 1.642e-01 | 2.00 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 4 NumPts= 620994 | Fixes: Font= 0 VL= 0 rho*= 0 | Failures: 0 InHoriz= 0 / 0 | Error: 7.467e-16, ErrDenom: 1.313e+00 | 2.03 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 5 NumPts= 1332331 | Fixes: Font= 0 VL= 0 rho*= 1077 | Failures: 0 InHoriz= 0 / 0 | Error: 1.908e-05, ErrDenom: 1.050e+01 | 2.13 iters/gridpt INFO (IllinoisGRMHD): Setting Symmetry = none... at iteration = 0 INFO (IllinoisGRMHD): C2P: Lev: 6 NumPts= 1727638 | Fixes: Font= 0 VL= 0 rho*= 25462 | Failures: 0 InHoriz= 0 / 0 | Error: 3.116e-05, ErrDenom: 7.426e+01 | 2.62 iters/gridpt INFO (TerminationTrigger): Remaining wallclock time for your job is 59.0333 minutes. Triggering termination... particle_tracerET: Iteration #0: Need to specify 4981 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.80e-03. Need ~ 262 more iterations. particle_tracerET: Iteration #1: Need to specify 4966 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.40e-03. Need ~ 292 more iterations. particle_tracerET: Iteration #2: Need to specify 4943 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.80e-03. Need ~ 260 more iterations. particle_tracerET: Iteration #3: Need to specify 4925 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.75e-03. Need ~ 262 more iterations. particle_tracerET: Iteration #4: Need to specify 4904 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 255 more iterations. particle_tracerET: Iteration #5: Need to specify 4887 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.77e-03. Need ~ 259 more iterations. particle_tracerET: Iteration #6: Need to specify 4873 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.63e-03. Need ~ 268 more iterations. particle_tracerET: Iteration #7: Need to specify 4852 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.70e-03. Need ~ 262 more iterations. particle_tracerET: Iteration #8: Need to specify 4842 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.51e-03. Need ~ 275 more iterations. particle_tracerET: Iteration #9: Need to specify 4821 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.58e-03. Need ~ 269 more iterations. particle_tracerET: Iteration #10: Need to specify 4804 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.56e-03. Need ~ 269 more iterations. particle_tracerET: Iteration #11: Need to specify 4777 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.72e-03. Need ~ 257 more iterations. particle_tracerET: Iteration #12: Need to specify 4751 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 248 more iterations. particle_tracerET: Iteration #13: Need to specify 4730 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 245 more iterations. particle_tracerET: Iteration #14: Need to specify 4717 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.77e-03. Need ~ 250 more iterations. particle_tracerET: Iteration #15: Need to specify 4695 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.81e-03. Need ~ 246 more iterations. particle_tracerET: Iteration #16: Need to specify 4679 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.78e-03. Need ~ 247 more iterations. particle_tracerET: Iteration #17: Need to specify 4665 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.72e-03. Need ~ 250 more iterations. particle_tracerET: Iteration #18: Need to specify 4649 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.69e-03. Need ~ 251 more iterations. particle_tracerET: Iteration #19: Need to specify 4630 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.70e-03. Need ~ 250 more iterations. particle_tracerET: Iteration #20: Need to specify 4607 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.74e-03. Need ~ 246 more iterations. particle_tracerET: Iteration #21: Need to specify 4582 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.80e-03. Need ~ 241 more iterations. particle_tracerET: Iteration #22: Need to specify 4563 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.80e-03. Need ~ 240 more iterations. particle_tracerET: Iteration #23: Need to specify 4539 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 236 more iterations. particle_tracerET: Iteration #24: Need to specify 4520 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 235 more iterations. particle_tracerET: Iteration #25: Need to specify 4503 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 235 more iterations. particle_tracerET: Iteration #26: Need to specify 4489 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 237 more iterations. particle_tracerET: Iteration #27: Need to specify 4469 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 235 more iterations. particle_tracerET: Iteration #28: Need to specify 4452 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.78e-03. Need ~ 235 more iterations. particle_tracerET: Iteration #29: Need to specify 4431 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 233 more iterations. particle_tracerET: Iteration #30: Need to specify 4412 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 232 more iterations. particle_tracerET: Iteration #31: Need to specify 4391 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.81e-03. Need ~ 230 more iterations. particle_tracerET: Iteration #32: Need to specify 4370 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 228 more iterations. particle_tracerET: Iteration #33: Need to specify 4352 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.81e-03. Need ~ 228 more iterations. particle_tracerET: Iteration #34: Need to specify 4338 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.78e-03. Need ~ 229 more iterations. particle_tracerET: Iteration #35: Need to specify 4324 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.76e-03. Need ~ 230 more iterations. particle_tracerET: Iteration #36: Need to specify 4304 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.76e-03. Need ~ 228 more iterations. particle_tracerET: Iteration #37: Need to specify 4282 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.78e-03. Need ~ 226 more iterations. particle_tracerET: Iteration #38: Need to specify 4265 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.77e-03. Need ~ 226 more iterations. particle_tracerET: Iteration #39: Need to specify 4242 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 223 more iterations. particle_tracerET: Iteration #40: Need to specify 4230 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.76e-03. Need ~ 225 more iterations. particle_tracerET: Iteration #41: Need to specify 4214 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.74e-03. Need ~ 225 more iterations. particle_tracerET: Iteration #42: Need to specify 4199 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.73e-03. Need ~ 225 more iterations. particle_tracerET: Iteration #43: Need to specify 4182 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.72e-03. Need ~ 224 more iterations. particle_tracerET: Iteration #44: Need to specify 4168 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.70e-03. Need ~ 225 more iterations. particle_tracerET: Iteration #45: Need to specify 4151 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.69e-03. Need ~ 224 more iterations. particle_tracerET: Iteration #46: Need to specify 4129 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.71e-03. Need ~ 222 more iterations. particle_tracerET: Iteration #47: Need to specify 4104 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.73e-03. Need ~ 219 more iterations. particle_tracerET: Iteration #48: Need to specify 4093 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.70e-03. Need ~ 221 more iterations. particle_tracerET: Iteration #49: Need to specify 4074 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.70e-03. Need ~ 219 more iterations. particle_tracerET: Iteration #50: Need to specify 4049 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.73e-03. Need ~ 217 more iterations. particle_tracerET: Iteration #51: Need to specify 4027 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.74e-03. Need ~ 215 more iterations. particle_tracerET: Iteration #52: Need to specify 4004 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.76e-03. Need ~ 213 more iterations. particle_tracerET: Iteration #53: Need to specify 3987 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.75e-03. Need ~ 212 more iterations. particle_tracerET: Iteration #54: Need to specify 3966 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.76e-03. Need ~ 210 more iterations. particle_tracerET: Iteration #55: Need to specify 3947 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.76e-03. Need ~ 209 more iterations. particle_tracerET: Iteration #56: Need to specify 3931 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.75e-03. Need ~ 209 more iterations. particle_tracerET: Iteration #57: Need to specify 3909 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.76e-03. Need ~ 207 more iterations. particle_tracerET: Iteration #58: Need to specify 3889 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.77e-03. Need ~ 206 more iterations. particle_tracerET: Iteration #59: Need to specify 3874 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.75e-03. Need ~ 206 more iterations. particle_tracerET: Iteration #60: Need to specify 3843 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 202 more iterations. particle_tracerET: Iteration #61: Need to specify 3825 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 201 more iterations. particle_tracerET: Iteration #62: Need to specify 3809 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.78e-03. Need ~ 201 more iterations. particle_tracerET: Iteration #63: Need to specify 3789 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.78e-03. Need ~ 200 more iterations. particle_tracerET: Iteration #64: Need to specify 3769 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 199 more iterations. particle_tracerET: Iteration #65: Need to specify 3750 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 198 more iterations. particle_tracerET: Iteration #66: Need to specify 3731 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 196 more iterations. particle_tracerET: Iteration #67: Need to specify 3714 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.78e-03. Need ~ 196 more iterations. particle_tracerET: Iteration #68: Need to specify 3693 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 194 more iterations. particle_tracerET: Iteration #69: Need to specify 3668 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.81e-03. Need ~ 192 more iterations. particle_tracerET: Iteration #70: Need to specify 3647 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.81e-03. Need ~ 191 more iterations. particle_tracerET: Iteration #71: Need to specify 3626 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 190 more iterations. particle_tracerET: Iteration #72: Need to specify 3610 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.81e-03. Need ~ 189 more iterations. particle_tracerET: Iteration #73: Need to specify 3594 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.80e-03. Need ~ 189 more iterations. particle_tracerET: Iteration #74: Need to specify 3578 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 188 more iterations. particle_tracerET: Iteration #75: Need to specify 3552 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.81e-03. Need ~ 186 more iterations. particle_tracerET: Iteration #76: Need to specify 3529 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 184 more iterations. particle_tracerET: Iteration #77: Need to specify 3507 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 183 more iterations. particle_tracerET: Iteration #78: Need to specify 3486 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 181 more iterations. particle_tracerET: Iteration #79: Need to specify 3466 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 180 more iterations. particle_tracerET: Iteration #80: Need to specify 3447 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 179 more iterations. particle_tracerET: Iteration #81: Need to specify 3428 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 178 more iterations. particle_tracerET: Iteration #82: Need to specify 3415 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 178 more iterations. particle_tracerET: Iteration #83: Need to specify 3399 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.81e-03. Need ~ 178 more iterations. particle_tracerET: Iteration #84: Need to specify 3377 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 176 more iterations. particle_tracerET: Iteration #85: Need to specify 3356 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 175 more iterations. particle_tracerET: Iteration #86: Need to specify 3333 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 173 more iterations. particle_tracerET: Iteration #87: Need to specify 3308 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 172 more iterations. particle_tracerET: Iteration #88: Need to specify 3298 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 172 more iterations. particle_tracerET: Iteration #89: Need to specify 3280 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 171 more iterations. particle_tracerET: Iteration #90: Need to specify 3263 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 170 more iterations. particle_tracerET: Iteration #91: Need to specify 3246 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.81e-03. Need ~ 170 more iterations. particle_tracerET: Iteration #92: Need to specify 3229 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.81e-03. Need ~ 169 more iterations. particle_tracerET: Iteration #93: Need to specify 3211 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.81e-03. Need ~ 168 more iterations. particle_tracerET: Iteration #94: Need to specify 3192 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.81e-03. Need ~ 167 more iterations. particle_tracerET: Iteration #95: Need to specify 3177 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.80e-03. Need ~ 167 more iterations. particle_tracerET: Iteration #96: Need to specify 3161 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 166 more iterations. particle_tracerET: Iteration #97: Need to specify 3143 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 165 more iterations. particle_tracerET: Iteration #98: Need to specify 3118 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.80e-03. Need ~ 164 more iterations. particle_tracerET: Iteration #99: Need to specify 3105 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 163 more iterations. particle_tracerET: Iteration #100: Need to specify 3081 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.80e-03. Need ~ 162 more iterations. particle_tracerET: Iteration #101: Need to specify 3055 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.81e-03. Need ~ 160 more iterations. particle_tracerET: Iteration #102: Need to specify 3035 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 159 more iterations. particle_tracerET: Iteration #103: Need to specify 3016 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 158 more iterations. particle_tracerET: Iteration #104: Need to specify 3001 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.81e-03. Need ~ 157 more iterations. particle_tracerET: Iteration #105: Need to specify 2981 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.81e-03. Need ~ 156 more iterations. particle_tracerET: Iteration #106: Need to specify 2967 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.80e-03. Need ~ 156 more iterations. particle_tracerET: Iteration #107: Need to specify 2950 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.80e-03. Need ~ 155 more iterations. particle_tracerET: Iteration #108: Need to specify 2925 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.81e-03. Need ~ 153 more iterations. particle_tracerET: Iteration #109: Need to specify 2909 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.80e-03. Need ~ 153 more iterations. particle_tracerET: Iteration #110: Need to specify 2895 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 152 more iterations. particle_tracerET: Iteration #111: Need to specify 2878 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 151 more iterations. particle_tracerET: Iteration #112: Need to specify 2855 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.80e-03. Need ~ 150 more iterations. particle_tracerET: Iteration #113: Need to specify 2845 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.78e-03. Need ~ 150 more iterations. particle_tracerET: Iteration #114: Need to specify 2824 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.78e-03. Need ~ 149 more iterations. particle_tracerET: Iteration #115: Need to specify 2810 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.78e-03. Need ~ 148 more iterations. particle_tracerET: Iteration #116: Need to specify 2791 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.78e-03. Need ~ 147 more iterations. particle_tracerET: Iteration #117: Need to specify 2770 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.78e-03. Need ~ 146 more iterations. particle_tracerET: Iteration #118: Need to specify 2737 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.80e-03. Need ~ 143 more iterations. particle_tracerET: Iteration #119: Need to specify 2719 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.80e-03. Need ~ 143 more iterations. particle_tracerET: Iteration #120: Need to specify 2707 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.79e-03. Need ~ 142 more iterations. particle_tracerET: Iteration #121: Need to specify 2685 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.80e-03. Need ~ 141 more iterations. particle_tracerET: Iteration #122: Need to specify 2664 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.80e-03. Need ~ 140 more iterations. particle_tracerET: Iteration #123: Need to specify 2643 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.80e-03. Need ~ 139 more iterations. particle_tracerET: Iteration #124: Need to specify 2612 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 136 more iterations. particle_tracerET: Iteration #125: Need to specify 2591 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 135 more iterations. particle_tracerET: Iteration #126: Need to specify 2567 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 133 more iterations. particle_tracerET: Iteration #127: Need to specify 2550 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 133 more iterations. particle_tracerET: Iteration #128: Need to specify 2530 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 132 more iterations. particle_tracerET: Iteration #129: Need to specify 2512 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 131 more iterations. particle_tracerET: Iteration #130: Need to specify 2485 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 129 more iterations. particle_tracerET: Iteration #131: Need to specify 2464 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 128 more iterations. particle_tracerET: Iteration #132: Need to specify 2442 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 126 more iterations. particle_tracerET: Iteration #133: Need to specify 2425 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 126 more iterations. particle_tracerET: Iteration #134: Need to specify 2407 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 125 more iterations. particle_tracerET: Iteration #135: Need to specify 2390 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 124 more iterations. particle_tracerET: Iteration #136: Need to specify 2377 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 124 more iterations. particle_tracerET: Iteration #137: Need to specify 2355 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 122 more iterations. particle_tracerET: Iteration #138: Need to specify 2338 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 122 more iterations. particle_tracerET: Iteration #139: Need to specify 2324 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 121 more iterations. particle_tracerET: Iteration #140: Need to specify 2310 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 121 more iterations. particle_tracerET: Iteration #141: Need to specify 2284 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 119 more iterations. particle_tracerET: Iteration #142: Need to specify 2252 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 117 more iterations. particle_tracerET: Iteration #143: Need to specify 2234 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 116 more iterations. particle_tracerET: Iteration #144: Need to specify 2214 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 115 more iterations. particle_tracerET: Iteration #145: Need to specify 2188 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 113 more iterations. particle_tracerET: Iteration #146: Need to specify 2167 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 112 more iterations. particle_tracerET: Iteration #147: Need to specify 2155 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 112 more iterations. particle_tracerET: Iteration #148: Need to specify 2129 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 110 more iterations. particle_tracerET: Iteration #149: Need to specify 2110 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 109 more iterations. particle_tracerET: Iteration #150: Need to specify 2093 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 108 more iterations. particle_tracerET: Iteration #151: Need to specify 2082 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 108 more iterations. particle_tracerET: Iteration #152: Need to specify 2050 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 106 more iterations. particle_tracerET: Iteration #153: Need to specify 2035 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 105 more iterations. particle_tracerET: Iteration #154: Need to specify 2012 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 104 more iterations. particle_tracerET: Iteration #155: Need to specify 1996 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 103 more iterations. particle_tracerET: Iteration #156: Need to specify 1982 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 103 more iterations. particle_tracerET: Iteration #157: Need to specify 1977 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 103 more iterations. particle_tracerET: Iteration #158: Need to specify 1954 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 101 more iterations. particle_tracerET: Iteration #159: Need to specify 1941 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 101 more iterations. particle_tracerET: Iteration #160: Need to specify 1923 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 100 more iterations. particle_tracerET: Iteration #161: Need to specify 1897 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 99 more iterations. particle_tracerET: Iteration #162: Need to specify 1872 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 97 more iterations. particle_tracerET: Iteration #163: Need to specify 1852 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 96 more iterations. particle_tracerET: Iteration #164: Need to specify 1828 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 95 more iterations. particle_tracerET: Iteration #165: Need to specify 1818 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 94 more iterations. particle_tracerET: Iteration #166: Need to specify 1800 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 93 more iterations. particle_tracerET: Iteration #167: Need to specify 1791 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.82e-03. Need ~ 93 more iterations. particle_tracerET: Iteration #168: Need to specify 1767 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 92 more iterations. particle_tracerET: Iteration #169: Need to specify 1745 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 91 more iterations. particle_tracerET: Iteration #170: Need to specify 1729 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 90 more iterations. particle_tracerET: Iteration #171: Need to specify 1705 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 89 more iterations. particle_tracerET: Iteration #172: Need to specify 1687 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 88 more iterations. particle_tracerET: Iteration #173: Need to specify 1671 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 87 more iterations. particle_tracerET: Iteration #174: Need to specify 1642 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 85 more iterations. particle_tracerET: Iteration #175: Need to specify 1622 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 84 more iterations. particle_tracerET: Iteration #176: Need to specify 1599 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 83 more iterations. particle_tracerET: Iteration #177: Need to specify 1568 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 81 more iterations. particle_tracerET: Iteration #178: Need to specify 1543 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 79 more iterations. particle_tracerET: Iteration #179: Need to specify 1529 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 79 more iterations. particle_tracerET: Iteration #180: Need to specify 1513 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 78 more iterations. particle_tracerET: Iteration #181: Need to specify 1494 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 77 more iterations. particle_tracerET: Iteration #182: Need to specify 1477 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 76 more iterations. particle_tracerET: Iteration #183: Need to specify 1458 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 75 more iterations. particle_tracerET: Iteration #184: Need to specify 1438 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 74 more iterations. particle_tracerET: Iteration #185: Need to specify 1417 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 73 more iterations. particle_tracerET: Iteration #186: Need to specify 1396 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 72 more iterations. particle_tracerET: Iteration #187: Need to specify 1371 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 71 more iterations. particle_tracerET: Iteration #188: Need to specify 1351 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 69 more iterations. particle_tracerET: Iteration #189: Need to specify 1329 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 68 more iterations. particle_tracerET: Iteration #190: Need to specify 1308 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.87e-03. Need ~ 67 more iterations. particle_tracerET: Iteration #191: Need to specify 1287 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.87e-03. Need ~ 66 more iterations. particle_tracerET: Iteration #192: Need to specify 1267 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.87e-03. Need ~ 65 more iterations. particle_tracerET: Iteration #193: Need to specify 1249 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.87e-03. Need ~ 64 more iterations. particle_tracerET: Iteration #194: Need to specify 1237 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 64 more iterations. particle_tracerET: Iteration #195: Need to specify 1219 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 63 more iterations. particle_tracerET: Iteration #196: Need to specify 1194 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 61 more iterations. particle_tracerET: Iteration #197: Need to specify 1175 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 60 more iterations. particle_tracerET: Iteration #198: Need to specify 1153 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.87e-03. Need ~ 59 more iterations. particle_tracerET: Iteration #199: Need to specify 1136 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 58 more iterations. particle_tracerET: Iteration #200: Need to specify 1118 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 57 more iterations. particle_tracerET: Iteration #201: Need to specify 1094 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.87e-03. Need ~ 56 more iterations. particle_tracerET: Iteration #202: Need to specify 1075 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.87e-03. Need ~ 55 more iterations. particle_tracerET: Iteration #203: Need to specify 1059 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 54 more iterations. particle_tracerET: Iteration #204: Need to specify 1037 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.87e-03. Need ~ 53 more iterations. particle_tracerET: Iteration #205: Need to specify 1009 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.87e-03. Need ~ 52 more iterations. particle_tracerET: Iteration #206: Need to specify 996 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.87e-03. Need ~ 51 more iterations. particle_tracerET: Iteration #207: Need to specify 966 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.88e-03. Need ~ 49 more iterations. particle_tracerET: Iteration #208: Need to specify 949 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.88e-03. Need ~ 48 more iterations. particle_tracerET: Iteration #209: Need to specify 934 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.87e-03. Need ~ 48 more iterations. particle_tracerET: Iteration #210: Need to specify 911 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.88e-03. Need ~ 47 more iterations. particle_tracerET: Iteration #211: Need to specify 899 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.87e-03. Need ~ 46 more iterations. particle_tracerET: Iteration #212: Need to specify 886 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 45 more iterations. particle_tracerET: Iteration #213: Need to specify 873 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 45 more iterations. particle_tracerET: Iteration #214: Need to specify 852 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 44 more iterations. particle_tracerET: Iteration #215: Need to specify 830 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 42 more iterations. particle_tracerET: Iteration #216: Need to specify 815 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 42 more iterations. particle_tracerET: Iteration #217: Need to specify 803 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 41 more iterations. particle_tracerET: Iteration #218: Need to specify 792 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 41 more iterations. particle_tracerET: Iteration #219: Need to specify 771 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 40 more iterations. particle_tracerET: Iteration #220: Need to specify 753 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 39 more iterations. particle_tracerET: Iteration #221: Need to specify 733 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 38 more iterations. particle_tracerET: Iteration #222: Need to specify 712 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 37 more iterations. particle_tracerET: Iteration #223: Need to specify 689 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 35 more iterations. particle_tracerET: Iteration #224: Need to specify 668 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 34 more iterations. particle_tracerET: Iteration #225: Need to specify 651 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 33 more iterations. particle_tracerET: Iteration #226: Need to specify 630 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 32 more iterations. particle_tracerET: Iteration #227: Need to specify 613 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 31 more iterations. particle_tracerET: Iteration #228: Need to specify 589 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 30 more iterations. particle_tracerET: Iteration #229: Need to specify 568 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 29 more iterations. particle_tracerET: Iteration #230: Need to specify 548 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 28 more iterations. particle_tracerET: Iteration #231: Need to specify 531 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 27 more iterations. particle_tracerET: Iteration #232: Need to specify 516 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 26 more iterations. particle_tracerET: Iteration #233: Need to specify 492 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 25 more iterations. particle_tracerET: Iteration #234: Need to specify 480 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 24 more iterations. particle_tracerET: Iteration #235: Need to specify 456 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 23 more iterations. particle_tracerET: Iteration #236: Need to specify 432 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 22 more iterations. particle_tracerET: Iteration #237: Need to specify 417 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 21 more iterations. particle_tracerET: Iteration #238: Need to specify 390 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.86e-03. Need ~ 20 more iterations. particle_tracerET: Iteration #239: Need to specify 375 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 19 more iterations. particle_tracerET: Iteration #240: Need to specify 356 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 18 more iterations. particle_tracerET: Iteration #241: Need to specify 339 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 17 more iterations. particle_tracerET: Iteration #242: Need to specify 326 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 16 more iterations. particle_tracerET: Iteration #243: Need to specify 312 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 16 more iterations. particle_tracerET: Iteration #244: Need to specify 292 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 15 more iterations. particle_tracerET: Iteration #245: Need to specify 277 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 14 more iterations. particle_tracerET: Iteration #246: Need to specify 264 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 13 more iterations. particle_tracerET: Iteration #247: Need to specify 240 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 12 more iterations. particle_tracerET: Iteration #248: Need to specify 226 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.83e-03. Need ~ 11 more iterations. particle_tracerET: Iteration #249: Need to specify 201 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 10 more iterations. particle_tracerET: Iteration #250: Need to specify 181 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 9 more iterations. particle_tracerET: Iteration #251: Need to specify 157 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 8 more iterations. particle_tracerET: Iteration #252: Need to specify 140 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 7 more iterations. particle_tracerET: Iteration #253: Need to specify 121 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 6 more iterations. particle_tracerET: Iteration #254: Need to specify 101 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 5 more iterations. particle_tracerET: Iteration #255: Need to specify 80 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 4 more iterations. particle_tracerET: Iteration #256: Need to specify 63 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.84e-03. Need ~ 3 more iterations. particle_tracerET: Iteration #257: Need to specify 31 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 1 more iterations. particle_tracerET: Iteration #258: Need to specify 18 more particle locations. Central condensation parameter = 1.000000e+00. Success rate = 3.85e-03. Need ~ 0 more iterations. SHOULD BE ALL DONE! Inside particle_tracerET: I think CCTK_DELTA_TIME=1.406250e-01 1HEY PART TRACE: -1.239893e-05 2.151692e-01 1.232224e-04 INFO (TerminationTrigger): Remaining wallclock time for your job is 58.65 minutes. Triggering termination... INFO (TerminationTrigger): Remaining wallclock time for your job is 58.5167 minutes. Triggering termination... INFO (TerminationTrigger): Remaining wallclock time for your job is 58.35 minutes. Triggering termination... INFO (TerminationTrigger): Remaining wallclock time for your job is 58.1167 minutes. Triggering termination... INFO (TerminationTrigger): Remaining wallclock time for your job is 57.9167 minutes. Triggering termination... VolumeIntegrals_GRMHD: Iter 1, num_reductions=1, Integ. quantity=one, sphere moves/tracks AMR centre=-1/-1 | INSIDE center x,y,z=0.000000e+00,0.000000e+00,0.000000e+00 ; r=0.000000e+00 | OUTSIDE center x,y,z=0.000000e+00,0.000000e+00,0.000000e+00 ; r=0.000000e+00 VolumeIntegrals_GRMHD: Iter 2, num_reductions=4, Integ. quantity=centerofmass, sphere moves/tracks AMR centre=0/-1 | INSIDE center x,y,z=-1.520000e+01,0.000000e+00,0.000000e+00 ; r=1.350000e+01 | OUTSIDE center x,y,z=-1.520000e+01,0.000000e+00,0.000000e+00 ; r=0.000000e+00 VolumeIntegrals_GRMHD: AMR centre #2 tracks Integral 0: (x,y,z)=(-1.858315e+01,1.139904e-17,-1.053221e-15) [norm=8.817717e-01]. Prev centre @ (-1.520000e+01,0.000000e+00,0.000000e+00). VolumeIntegrals_GRMHD: Iter 3, num_reductions=1, Integ. quantity=one, sphere moves/tracks AMR centre=0/-1 | INSIDE center x,y,z=-1.520000e+01,0.000000e+00,0.000000e+00 ; r=1.350000e+01 | OUTSIDE center x,y,z=-1.520000e+01,0.000000e+00,0.000000e+00 ; r=0.000000e+00 VolumeIntegrals_GRMHD: Iter 4, num_reductions=4, Integ. quantity=centerofmass, sphere moves/tracks AMR centre=1/-1 | INSIDE center x,y,z=1.520000e+01,0.000000e+00,0.000000e+00 ; r=1.350000e+01 | OUTSIDE center x,y,z=1.520000e+01,0.000000e+00,0.000000e+00 ; r=0.000000e+00 VolumeIntegrals_GRMHD: AMR centre #4 tracks Integral 1: (x,y,z)=(2.079096e+01,9.475074e-16,-1.124290e-15) [norm=3.163531e-01]. Prev centre @ (1.520000e+01,0.000000e+00,0.000000e+00). VolumeIntegrals_GRMHD: Iter 5, num_reductions=1, Integ. quantity=one, sphere moves/tracks AMR centre=1/-1 | INSIDE center x,y,z=1.520000e+01,0.000000e+00,0.000000e+00 ; r=1.350000e+01 | OUTSIDE center x,y,z=1.520000e+01,0.000000e+00,0.000000e+00 ; r=0.000000e+00 VolumeIntegrals_GRMHD: Iter 6, num_reductions=1, Integ. quantity=restmass, sphere moves/tracks AMR centre=-1/-1 | INSIDE center x,y,z=0.000000e+00,0.000000e+00,0.000000e+00 ; r=0.000000e+00 | OUTSIDE center x,y,z=0.000000e+00,0.000000e+00,0.000000e+00 ; r=0.000000e+00 INFO (TerminationTrigger): Remaining wallclock time for your job is 57.5667 minutes. Triggering termination... ---------------------------------------------------------------------------------------------- Iteration Time | *me_per_hour | *RMHD::rho_b | *ure_checker | ML_BSSN::phi | ML_BSSN::trK | | maximum | maximum | maximum | maximum ---------------------------------------------------------------------------------------------- 0 0.000 | 0.0000000 | 0.0001265 | 10.0000000 | 0.0932081 | 2.005774e-18 -------------------------------------------------------------------------------- Done. Stopping: Thu Jul 19 11:27:37 PDT 2018 Done. Thu Jul 19 11:27:37 PDT 2018 Simfactory Done at date: 0