Trac June 2023

trac@lists.einsteintoolkit.org

1 participants
81 discussions

#2497: IllinoisGRMHD is incompatible with setting TmunuBase::stress_energy_at_RHS = "no"
by Samuel Cupp 13 Jun '23

13 Jun '23

#2497: IllinoisGRMHD is incompatible with setting TmunuBase::stress_energy_at_RHS = "no" Reporter: Gabriele Bozzola Status: open Milestone: Version: development version Type: bug Priority: major Component: EinsteinToolkit thorn Comment (by Samuel Cupp): I’m not sure what your schedule.ccl currently looks like, but the master branch of \`IllinoisGRMHD\` only has those two functions \(which I’ll refer to as con2prim\) in AddToTmunu, not SetTmunu. Did you also add it to SetTmunu? All Leo did was change the scheduling from something that looks like ``` schedule con2prim in AddToTmunu ``` to ``` if(stress_energy_at_RHS) { schedule con2prim in AddToTmunu } else { schedule con2prim in HydroBase_Con2Prim } ``` ‌ -- Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2497/illinoisgrmhd-is-…

1 0

#2738: Caunda/Proca and Canuda/Scalar test faiilures after commits that "remove parameter eta_beta_dynamic"
by Miguel Zilhão 12 Jun '23

12 Jun '23

#2738: Caunda/Proca and Canuda/Scalar test faiilures after commits that "remove parameter eta_beta_dynamic" Reporter: Roland Haas Status: new Milestone: Version: Type: bug Priority: major Component: EinsteinToolkit thorn Comment (by Miguel Zilhão): I see, sorry, I should have read properly the description that this is happening in the scalar and proca arrangements :-\) I've removed the parameter that was giving problems in the following commits: * [https://bitbucket.org/canuda/proca/commits/d24ddf835b0552410175fab2f6f84fe1… * [https://bitbucket.org/canuda/scalar/commits/44df0f7bb4afa01b406a5876acd107a… I hope this is now fixed. ‌ -- Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2738/caunda-proca-and-…

1 0

#2738: Caunda/Proca and Canuda/Scalar test faiilures after commits that "remove parameter eta_beta_dynamic"
by Roland Haas 12 Jun '23

12 Jun '23

#2738: Caunda/Proca and Canuda/Scalar test faiilures after commits that "remove parameter eta_beta_dynamic" Reporter: Roland Haas Status: new Milestone: Version: Type: bug Priority: major Component: EinsteinToolkit thorn Comment (by Roland Haas): Yes, this is still happening. You can take a look at the tests at [https://einsteintoolkit.github.io/tests/](https://einsteintoolkit.github.io… The LeanBSSN\_Ei\_mu0.4\_c0.05.log test \(you can click on the “log” link to access the output file\) fails with ``` Major error in parameter file '/home/runner/simulations/TestJob01_temp_1/output-0000/arrangements/Proca/NPScalars_Proca/test/LeanBSSN_Ei_mu0.4_c0.05.par' line 150: Parameter 'LeanBSSNMoL::eta_beta_dynamic' not found ``` which should also fail on your system Did you maybe forget to push all changes? I just tried and I get a test failure for the LeabBSSN Ei test as well after a fresh checkout on my workstation: ``` ./GetComponents --parallel --shallow https://bitbucket.org/einsteintoolkit/manifest/raw/master/einsteintoolkit.th ``` Please note that there are apparently multiple tests in Canuda named `LeanBSSN_Ei_mu0.4_c0.05.par` and the one that fails is `./repos/Proca/NPScalars_Proca/test/LeanBSSN_Ei_mu0.4_c0.05.par`. -- Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2738/caunda-proca-and-…

1 0

#2738: Caunda/Proca and Canuda/Scalar test faiilures after commits that "remove parameter eta_beta_dynamic"
by Miguel Zilhão 12 Jun '23

12 Jun '23

#2738: Caunda/Proca and Canuda/Scalar test faiilures after commits that "remove parameter eta_beta_dynamic" Reporter: Roland Haas Status: new Milestone: Version: Type: bug Priority: major Component: EinsteinToolkit thorn Comment (by Miguel Zilhão): This was indeed a problem, but it should have been fixed with commit [https://bitbucket.org/canuda/lean\_public/commits/7457e1822e3a880efd67e11f1… On my machine all tests are currently passing… Is this still a problem with the current master branch? -- Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2738/caunda-proca-and-…

1 0

#2497: IllinoisGRMHD is incompatible with setting TmunuBase::stress_energy_at_RHS = "no"
by Vikram Manikantan 11 Jun '23

11 Jun '23

1 0

#2497: IllinoisGRMHD is incompatible with setting TmunuBase::stress_energy_at_RHS = "no"
by Vikram Manikantan 11 Jun '23

11 Jun '23

1 0

#2497: IllinoisGRMHD is incompatible with setting TmunuBase::stress_energy_at_RHS = "no"
by Vikram Manikantan 11 Jun '23

11 Jun '23

#2497: IllinoisGRMHD is incompatible with setting TmunuBase::stress_energy_at_RHS = "no" Reporter: Gabriele Bozzola Status: open Milestone: Version: development version Type: bug Priority: major Component: EinsteinToolkit thorn Comment (by Vikram Manikantan): Hi @{5bae587b96242d2e2b6110a4} ! Yes, I am attaching it here. I also want to follow up on a previous question. Currently there is a SetTmunu scheduled in schedule.ccl and in your solution there is an AddToTmunu in the if statement. I am wondering if both should be in the file, or just the AddToTmunu? ‌ disk.par --> ``` ################################################################################ # Script variables ################################################################################ # __file__ = disk.rpar # __name__ = __main__ # bh_movement_threshold = 0.124 # carpet_str = CartGrid3D::type = "coordbase" # Carpet::domain_from_coordbase = "yes" # # Driver::ghost_size = 4 # CoordBase::boundary_size_x_lower = 4 # CoordBase::boundary_size_y_lower = 4 # CoordBase::boundary_size_z_lower = 4 # CoordBase::boundary_size_x_upper = 4 # CoordBase::boundary_size_y_upper = 4 # CoordBase::boundary_size_z_upper = 4 # # ReflectionSymmetry::avoid_origin_x = no # ReflectionSymmetry::avoid_origin_y = no # ReflectionSymmetry::avoid_origin_z = no # # CoordBase::domainsize = "minmax" # CoordBase::xmax = 5120.001785714286 # CoordBase::ymax = 5120.001785714286 # CoordBase::zmax = 5120.001785714286 # CoordBase::dx = 102.4 # CoordBase::dy = 102.4 # CoordBase::dz = 102.4 # CoordBase::xmin = -5119.998214285714 # CoordBase::ymin = -5119.998214285714 # CoordBase::zmin = -5119.998214285714 # # CarpetRegrid2::num_centres = 3 # Carpet::max_refinement_levels = 14 # Carpet::time_refinement_factors = "[1,1,1,1,1,1,1,1,1,2,4,8,16,32]" # # CarpetRegrid2::num_levels_1 = 14 # CarpetRegrid2::position_x_1 = -6.5 # CarpetRegrid2::position_y_1 = 0 # CarpetRegrid2::position_z_1 = 0 # CarpetRegrid2::radius_1[1] = 2539.52 # CarpetRegrid2::radius_1[2] = 1269.76 # CarpetRegrid2::radius_1[3] = 634.88 # CarpetRegrid2::radius_1[4] = 317.44 # CarpetRegrid2::radius_1[5] = 158.72 # CarpetRegrid2::radius_1[6] = 79.36 # CarpetRegrid2::radius_1[7] = 39.68 # CarpetRegrid2::radius_1[8] = 19.84 # CarpetRegrid2::radius_1[9] = 9.92 # CarpetRegrid2::radius_1[10] = 4.96 # CarpetRegrid2::radius_1[11] = 2.48 # CarpetRegrid2::radius_1[12] = 1.24 # CarpetRegrid2::radius_1[13] = 0.62 # CarpetRegrid2::num_levels_2 = 14 # CarpetRegrid2::position_x_2 = 6.5 # CarpetRegrid2::position_y_2 = 0 # CarpetRegrid2::position_z_2 = 0 # CarpetRegrid2::radius_2[1] = 2539.52 # CarpetRegrid2::radius_2[2] = 1269.76 # CarpetRegrid2::radius_2[3] = 634.88 # CarpetRegrid2::radius_2[4] = 317.44 # CarpetRegrid2::radius_2[5] = 158.72 # CarpetRegrid2::radius_2[6] = 79.36 # CarpetRegrid2::radius_2[7] = 39.68 # CarpetRegrid2::radius_2[8] = 19.84 # CarpetRegrid2::radius_2[9] = 9.92 # CarpetRegrid2::radius_2[10] = 4.96 # CarpetRegrid2::radius_2[11] = 2.48 # CarpetRegrid2::radius_2[12] = 1.24 # CarpetRegrid2::radius_2[13] = 0.62 # CarpetRegrid2::num_levels_3 = 11 # CarpetRegrid2::position_x_3 = 0 # CarpetRegrid2::position_y_3 = 0 # CarpetRegrid2::position_z_3 = 0 # CarpetRegrid2::radius_3[1] = 2539.52 # CarpetRegrid2::radius_3[2] = 1269.76 # CarpetRegrid2::radius_3[3] = 634.88 # CarpetRegrid2::radius_3[4] = 317.44 # CarpetRegrid2::radius_3[5] = 158.72 # CarpetRegrid2::radius_3[6] = 79.36 # CarpetRegrid2::radius_3[7] = 39.68 # CarpetRegrid2::radius_3[8] = 19.84 # CarpetRegrid2::radius_3[9] = 9.92 # CarpetRegrid2::radius_3[10] = 4.96 # cfl_condition = 1.5 # cfl_num = 0.45 # cfl_num_on_disk = 0.45 # col_width = 30 # damp_lorenz = 8 # diss_str = Dissipation::order = 7 # Dissipation::vars = "LeanBSSNMoL::conf_fac # LeanBSSNMoL::hmetric # LeanBSSNMoL::hcurv # LeanBSSNMoL::trk # LeanBSSNMoL::gammat # ADMBase::lapse # ADMBase::shift # LeanBSSNMoL::conf_fac # LeanBSSNMoL::hmetric # LeanBSSNMoL::hcurv # LeanBSSNMoL::trk # LeanBSSNMoL::gammat" # Dissipation::epsdis_for_level[0] = 1.953125e-05 # Dissipation::epsdis_for_level[1] = 3.90625e-05 # Dissipation::epsdis_for_level[2] = 7.8125e-05 # Dissipation::epsdis_for_level[3] = 0.00015625 # Dissipation::epsdis_for_level[4] = 0.0003125 # Dissipation::epsdis_for_level[5] = 0.000625 # Dissipation::epsdis_for_level[6] = 0.00125 # Dissipation::epsdis_for_level[7] = 0.0025 # Dissipation::epsdis_for_level[8] = 0.005 # Dissipation::epsdis_for_level[9] = 0.01 # Dissipation::epsdis_for_level[10] = 0.01 # Dissipation::epsdis_for_level[11] = 0.01 # Dissipation::epsdis_for_level[12] = 0.01 # Dissipation::epsdis_for_level[13] = 0.01 # Dissipation::epsdis_for_level[14] = 0.01 # dx_coarse = 102.4 # dx_fine = 0.0125 # dx_fine_on_disk = 3.2 # e = 2.718281828459045 # effective_cfl = 0.0017578125 # eps_dis = 0.01 # eta_beta = 2 # evolution_order = 6 # expected_ah_radius_minus = 0.5 # expected_ah_radius_plus = 0.5 # expected_final_ah_radius = 1.0 # expected_merger_time = 2000 # extr_radius = 200.0 # final_time = 4000 # ghosts = 4 # gw_radii_str = Multipole::radius[0] = 100.0 # Multipole::radius[1] = 114.28571428571429 # Multipole::radius[2] = 128.57142857142858 # Multipole::radius[3] = 142.85714285714286 # Multipole::radius[4] = 157.14285714285714 # Multipole::radius[5] = 171.42857142857144 # Multipole::radius[6] = 185.71428571428572 # Multipole::radius[7] = 200.0 # horizon_every = 64 # i = 7 # inf = inf # int_outer_boundary = 50 # key = xp # lorenz_condition = 1.5 # max_dt_allowed = 0.1875 # max_gw_radius = 200 # min_gw_radius = 100 # nan = nan # num_gw_radii = 8 # num_ref_radii = 13 # number_of_synced_levels = 9 # out_every_2D = 256 # out_every_3D = 1024 # outer_boundary = 5120.0 # outer_boundary_factor = 2 # outflow_radii_str = Outflow::surface_index[0] = 5 # Outflow::override_radius[0] = "yes" # Outflow::radius[0] = 100.0 # Outflow::surface_index[1] = 5 # Outflow::override_radius[1] = "yes" # Outflow::radius[1] = 114.28571428571429 # Outflow::surface_index[2] = 5 # Outflow::override_radius[2] = "yes" # Outflow::radius[2] = 128.57142857142858 # Outflow::surface_index[3] = 5 # Outflow::override_radius[3] = "yes" # Outflow::radius[3] = 142.85714285714286 # Outflow::surface_index[4] = 5 # Outflow::override_radius[4] = "yes" # Outflow::radius[4] = 157.14285714285714 # Outflow::surface_index[5] = 5 # Outflow::override_radius[5] = "yes" # Outflow::radius[5] = 171.42857142857144 # Outflow::surface_index[6] = 5 # Outflow::override_radius[6] = "yes" # Outflow::radius[6] = 185.71428571428572 # Outflow::surface_index[7] = 5 # Outflow::override_radius[7] = "yes" # Outflow::radius[7] = 200.0 # par_b = 6.5 # pi = 3.141592653589793 # radius = 2539.52 # radius0 = 0.62 # radius_number = 1 # refinement_level_radius_factor = 2 # regrid_every = 256 # rl_synced_every = 32 # rlrf = 2 # skip_radii = 3 # step_gw_radius = 14.285714285714286 # tau = 6.283185307179586 # timers_every = 1024 # tp_str = ADMBase::initial_data = "twopuncturespowerlaw" # ADMBase::initial_lapse = "twopunctures-averaged" # ADMBase::initial_shift = "zero" # ADMBase::initial_dtlapse = "zero" # ADMBase::initial_dtshift = "zero" # # TwoPuncturesPowerLaw::give_bare_mass = "no" # TwoPuncturesPowerLaw::par_b = 6.5 # TwoPuncturesPowerLaw::target_m_plus = 0.5 # TwoPuncturesPowerLaw::target_m_minus = 0.5 # TwoPuncturesPowerLaw::par_m_plus = 0.5 # TwoPuncturesPowerLaw::par_m_minus = 0.5 # TwoPuncturesPowerLaw::par_P_plus[0] = -0.000413874209381327 # TwoPuncturesPowerLaw::par_P_plus[1] = 0.0807444483485998 # TwoPuncturesPowerLaw::par_P_plus[2] = 0 # TwoPuncturesPowerLaw::par_P_minus[0] = 0.000413874209381327 # TwoPuncturesPowerLaw::par_P_minus[1] = -0.0807444483485998 # TwoPuncturesPowerLaw::par_P_minus[2] = 0 # TwoPuncturesPowerLaw::par_S_plus[0] = 0.0 # TwoPuncturesPowerLaw::par_S_plus[1] = 0.0 # TwoPuncturesPowerLaw::par_S_plus[2] = 0.0 # TwoPuncturesPowerLaw::par_S_minus[0] = 0.0 # TwoPuncturesPowerLaw::par_S_minus[1] = 0.0 # TwoPuncturesPowerLaw::par_S_minus[2] = 0.0 # TwoPuncturesPowerLaw::center_offset[0] = 0.0 # TwoPuncturesPowerLaw::center_offset[1] = 0.0 # TwoPuncturesPowerLaw::center_offset[2] = 0.0 # val = 6.5 # wave_extract_every = 32 # xm = -6.5 # xp = 6.5 ################################################################################ # Par file ################################################################################ ActiveThorns = " ADMBase ADMCoupling ADMMacros AEILocalInterp AHFinderDirect Boundary Carpet CarpetIOASCII CarpetIOBasic CarpetIOHDF5 CarpetIOScalar CarpetInterp CarpetLib # CarpetMask CarpetReduce CarpetRegrid2 CarpetSlab CarpetTracker CartGrid3D # Christoffel Convert_to_HydroBase CoordBase CoordGauge CycleClock # DiskDiagnostics Dissipation EOS_Base EOS_Omni Fortran GenericFD GSL ID_converter_ILGRMHD IllinoisGRMHD HydroBase HDF5 hwloc IOUtil InitBase LeanBSSNMoL LocalInterp LoopControl MoL MPI Multipole NaNChecker NewRad NPScalars PunctureTracker QuasiLocalMeasures OpenBLAS Outflow Seed_Magnetic_Fields Slab SpaceMask SphericalSurface StaticConformal SummationByParts SymBase SystemStatistics SystemTopology TerminationTrigger TensorTypes TGRTensor Time Timers TimerReport TmunuBase TwoPuncturesPowerLaw Vectors VolumeIntegrals_GRMHD zlib ReflectionSymmetry " # # Cactus::presync_mode = "mixed-error" #################################################################### # Grid details #################################################################### Carpet::use_buffer_zones = "yes" Carpet::init_3_timelevels = "no" Carpet::init_fill_timelevels = "yes" Carpet::verbose = "no" Carpet::veryverbose = "no" Carpet::schedule_barriers = "no" Carpet::storage_verbose = "no" Grid::avoid_origin = "no" Time::dtfac = 0.0017578125 Cactus::terminate = time Cactus::cctk_final_time = 4000 MoL::ODE_Method = "rk4" MoL::MoL_Intermediate_Steps = 4 MoL::MoL_Num_Scratch_Levels = 1 Carpet::prolongation_order_space = 5 Carpet::prolongation_order_time = 2 Cactus::cctk_full_warnings = "yes" Cactus::highlight_warning_messages = "no" CarpetLib::poison_new_memory = "yes" CarpetLib::poison_value = 113 CarpetLib::output_bboxes = "no" CarpetLib::interleave_communications = "yes" CarpetLib::combine_sends = "yes" CarpetLib::support_staggered_operators = "yes" # IllinoisGRMHD requires staggered operators CarpetRegrid2::freeze_unaligned_levels = "yes" CarpetRegrid2::freeze_unaligned_parent_levels = "yes" CarpetRegrid2::verbose = "yes" CarpetRegrid2::regrid_every = 256 CartGrid3D::type = "coordbase" Carpet::domain_from_coordbase = "yes" Driver::ghost_size = 4 CoordBase::boundary_size_x_lower = 4 CoordBase::boundary_size_y_lower = 4 CoordBase::boundary_size_z_lower = 4 CoordBase::boundary_size_x_upper = 4 CoordBase::boundary_size_y_upper = 4 CoordBase::boundary_size_z_upper = 4 ReflectionSymmetry::avoid_origin_x = no ReflectionSymmetry::avoid_origin_y = no ReflectionSymmetry::avoid_origin_z = no CoordBase::domainsize = "minmax" CoordBase::xmax = 5120.001785714286 CoordBase::ymax = 5120.001785714286 CoordBase::zmax = 5120.001785714286 CoordBase::dx = 102.4 CoordBase::dy = 102.4 CoordBase::dz = 102.4 CoordBase::xmin = -5119.998214285714 CoordBase::ymin = -5119.998214285714 CoordBase::zmin = -5119.998214285714 CarpetRegrid2::num_centres = 3 Carpet::max_refinement_levels = 14 Carpet::time_refinement_factors = "[1,1,1,1,1,1,1,1,1,2,4,8,16,32]" CarpetRegrid2::num_levels_1 = 14 CarpetRegrid2::position_x_1 = -6.5 CarpetRegrid2::position_y_1 = 0 CarpetRegrid2::position_z_1 = 0 CarpetRegrid2::radius_1[1] = 2539.52 CarpetRegrid2::radius_1[2] = 1269.76 CarpetRegrid2::radius_1[3] = 634.88 CarpetRegrid2::radius_1[4] = 317.44 CarpetRegrid2::radius_1[5] = 158.72 CarpetRegrid2::radius_1[6] = 79.36 CarpetRegrid2::radius_1[7] = 39.68 CarpetRegrid2::radius_1[8] = 19.84 CarpetRegrid2::radius_1[9] = 9.92 CarpetRegrid2::radius_1[10] = 4.96 CarpetRegrid2::radius_1[11] = 2.48 CarpetRegrid2::radius_1[12] = 1.24 CarpetRegrid2::radius_1[13] = 0.62 CarpetRegrid2::num_levels_2 = 14 CarpetRegrid2::position_x_2 = 6.5 CarpetRegrid2::position_y_2 = 0 CarpetRegrid2::position_z_2 = 0 CarpetRegrid2::radius_2[1] = 2539.52 CarpetRegrid2::radius_2[2] = 1269.76 CarpetRegrid2::radius_2[3] = 634.88 CarpetRegrid2::radius_2[4] = 317.44 CarpetRegrid2::radius_2[5] = 158.72 CarpetRegrid2::radius_2[6] = 79.36 CarpetRegrid2::radius_2[7] = 39.68 CarpetRegrid2::radius_2[8] = 19.84 CarpetRegrid2::radius_2[9] = 9.92 CarpetRegrid2::radius_2[10] = 4.96 CarpetRegrid2::radius_2[11] = 2.48 CarpetRegrid2::radius_2[12] = 1.24 CarpetRegrid2::radius_2[13] = 0.62 CarpetRegrid2::num_levels_3 = 11 CarpetRegrid2::position_x_3 = 0 CarpetRegrid2::position_y_3 = 0 CarpetRegrid2::position_z_3 = 0 CarpetRegrid2::radius_3[1] = 2539.52 CarpetRegrid2::radius_3[2] = 1269.76 CarpetRegrid2::radius_3[3] = 634.88 CarpetRegrid2::radius_3[4] = 317.44 CarpetRegrid2::radius_3[5] = 158.72 CarpetRegrid2::radius_3[6] = 79.36 CarpetRegrid2::radius_3[7] = 39.68 CarpetRegrid2::radius_3[8] = 19.84 CarpetRegrid2::radius_3[9] = 9.92 CarpetRegrid2::radius_3[10] = 4.96 CarpetRegrid2::movement_threshold_1 = 0.124 CarpetRegrid2::movement_threshold_2 = 0.124 CarpetRegrid2::movement_threshold_3 = 0.124 #################################################################### # Initial Data #################################################################### InitBase::initial_data_setup_method = "init_some_levels" HydroBase::initial_Avec = "zero" HydroBase::initial_Aphi = "zero" HydroBase::initial_Bvec = "zero" HydroBase::initial_hydro = "twopuncturespowerlaw" TwoPuncturesPowerLaw::ell_inner = 5.15 TwoPuncturesPowerLaw::sam_disk = 0.0 TwoPuncturesPowerLaw::RoM_inner = 18.0 IllinoisGRMHD::gamma_th = 1.333333333333333333333333333 ID_converter_ILGRMHD::Gamma_Initial = 1.333333333333333333333333333 # rho_max = 1 IllinoisGRMHD::K_poly = 7.98742131429e-4 ID_converter_ILGRMHD::K_Initial = 7.98742131429e-4 IllinoisGRMHD::rho_b_atm = 1e-8 IllinoisGRMHD::tau_atm = 1e-7 IllinoisGRMHD::compute_extra_diagnostics = "yes" ADMBase::initial_data = "twopuncturespowerlaw" ADMBase::initial_lapse = "twopunctures-averaged" ADMBase::initial_shift = "zero" ADMBase::initial_dtlapse = "zero" ADMBase::initial_dtshift = "zero" TwoPuncturesPowerLaw::give_bare_mass = "no" TwoPuncturesPowerLaw::par_b = 6.5 TwoPuncturesPowerLaw::target_m_plus = 0.5 TwoPuncturesPowerLaw::target_m_minus = 0.5 TwoPuncturesPowerLaw::par_m_plus = 0.5 TwoPuncturesPowerLaw::par_m_minus = 0.5 TwoPuncturesPowerLaw::par_P_plus[0] = -0.000413874209381327 TwoPuncturesPowerLaw::par_P_plus[1] = 0.0807444483485998 TwoPuncturesPowerLaw::par_P_plus[2] = 0 TwoPuncturesPowerLaw::par_P_minus[0] = 0.000413874209381327 TwoPuncturesPowerLaw::par_P_minus[1] = -0.0807444483485998 TwoPuncturesPowerLaw::par_P_minus[2] = 0 TwoPuncturesPowerLaw::par_S_plus[0] = 0.0 TwoPuncturesPowerLaw::par_S_plus[1] = 0.0 TwoPuncturesPowerLaw::par_S_plus[2] = 0.0 TwoPuncturesPowerLaw::par_S_minus[0] = 0.0 T -- Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2497/illinoisgrmhd-is-…

1 0

#2741: Some thorns seem to incorrectly use CCTK_GFINDEX4D
by Roland Haas 09 Jun '23

09 Jun '23

#2741: Some thorns seem to incorrectly use CCTK_GFINDEX4D Reporter: Gabriele Bozzola Status: open Milestone: Version: Type: bug Priority: major Component: EinsteinToolkit thorn Comment (by Roland Haas): yup. I had not updated my own repo. Sorry about the confusion. `particle_tracerET` is actually a newly submitted module in the Schwarzschild release. Anyhow, assigned to @{557058:8bc23f2a-45c0-477d-8ac4-a5a16c734278} . -- Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2741/some-thorns-seem-…

1 0

#2741: Some thorns seem to incorrectly use CCTK_GFINDEX4D
by Gabriele Bozzola 09 Jun '23

09 Jun '23

#2741: Some thorns seem to incorrectly use CCTK_GFINDEX4D Reporter: Gabriele Bozzola Status: open Milestone: Version: Type: bug Priority: major Component: EinsteinToolkit thorn Comment (by Gabriele Bozzola): Yes, we were running an older version and sorry for opening a duplicated ticket. The problem is now only in `wvuthorns_diagnostics` ``` particle_tracerET/src/compute_particles_four_velocities.C: const double vU[3] = { vel[CCTK_GFINDEX4D(cctkGH,i,j,k,0)], particle_tracerET/src/compute_particles_four_velocities.C: vel[CCTK_GFINDEX4D(cctkGH,i,j,k,1)], particle_tracerET/src/compute_particles_four_velocities.C: vel[CCTK_GFINDEX4D(cctkGH,i,j,k,2)] }; particle_tracerET/src/convert_to_MHD_3velocity.C: double ETvx = vel[CCTK_GFINDEX4D(cctkGH,i,j,k,0)]; particle_tracerET/src/convert_to_MHD_3velocity.C: double ETvy = vel[CCTK_GFINDEX4D(cctkGH,i,j,k,1)]; particle_tracerET/src/convert_to_MHD_3velocity.C: double ETvz = vel[CCTK_GFINDEX4D(cctkGH,i,j,k,2)]; Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C: Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,0)] = -(y[index] + 0.5*dY)*A_b*pow(max(Pressure_at_Ax_stagger-P_cut,0.0),n_s); Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C: if(!have_two_NSs_along_x_axis) Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,1)] = (xL + 0.5*dX)*A_b*pow(max(Pressure_at_Ay_stagger-P_cut,0.0),n_s); Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C: if(r1<=r_NS1) Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,1)] = (x1 + 0.5*dX)*A_b*pow(max(Pressure_at_Ay_stagger-P_cut,0.0),n_s); Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C: if(r2<=r_NS2) Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,1)] = (x2 + 0.5*dX)*A_b*pow(max(Pressure_at_Ay_stagger-P_cut,0.0),n_s); Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C: if(r1>r_NS1 && r2>r_NS2) Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,1)] = 0.0; // No external B-field. Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C: Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,0)] = - (yL + 0.5 * dY) * (Ap1 + Ap2); Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C: Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,1)] = (x1 + 0.5 * dX) * Ap1 + (x2 + 0.5 * dX) * Ap2; Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C: Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,2)] = 0.0; Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C: Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,0)] = -y[index]*A_b*pow(max(PL-P_cut,0.0),n_s); Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C: Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,1)] = x[index]*A_b*pow(max(PL-P_cut,0.0),n_s); Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C: Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,2)] = 0.0; smallbPoynET/src/compute_bi_b2_Poyn_fluxET.C: double ETvx = vel[CCTK_GFINDEX4D(cctkGH,i,j,k,0)]; smallbPoynET/src/compute_bi_b2_Poyn_fluxET.C: double ETvy = vel[CCTK_GFINDEX4D(cctkGH,i,j,k,1)]; smallbPoynET/src/compute_bi_b2_Poyn_fluxET.C: double ETvz = vel[CCTK_GFINDEX4D(cctkGH,i,j,k,2)]; smallbPoynET/src/compute_bi_b2_Poyn_fluxET.C: double Bx_center = Bvec[CCTK_GFINDEX4D(cctkGH,i,j,k,0)]; smallbPoynET/src/compute_bi_b2_Poyn_fluxET.C: double By_center = Bvec[CCTK_GFINDEX4D(cctkGH,i,j,k,1)]; smallbPoynET/src/compute_bi_b2_Poyn_fluxET.C: double Bz_center = Bvec[CCTK_GFINDEX4D(cctkGH,i,j,k,2)]; ``` ‌ -- Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2741/some-thorns-seem-…

1 0

#2741: Some thorns seem to incorrectly use CCTK_GFINDEX4D
by Roland Haas 09 Jun '23

09 Jun '23

#2741: Some thorns seem to incorrectly use CCTK_GFINDEX4D Reporter: Gabriele Bozzola Status: open Milestone: Version: Type: bug Priority: major Component: EinsteinToolkit thorn Comment (by Roland Haas): Oh, it still is in the current version \(master, git hash [0cb4d23](https://bitbucket.org/zach_etienne/wvuthorns/commits/0cb4d23a59e51… "Baikal & BaikalVacuum codegen: Add pip install sympy==1.11 to the automated codegen instructions" of [wvuthorns](https://bitbucket.org/zach_etienne/wvuthorns)\). We \(@{557058:8bc23f2a-45c0-477d-8ac4-a5a16c734278} and I\) just missed it in the last round of removing `CCTK_GFINDEX4D` and also missed it when handling a slightly different issue in `ID_converter_ILGRMHD` that actually showed lines with `CCTK_GFINDEX4D` in the ticket description. -- Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2741/some-thorns-seem-…

1 0

← Newer
1
2
3
4
5
6
7
8
9
Older →

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Trac June 2023