GRHayLHD

Samuel Cupp
Leonardo Rosa Werneck
Terrence Pierre Jacques
Zachariah Etienne

\( \)Date\( \)

Abstract

Provides a general relativistic hydrodynamic (GRHD) evolution code using GRHayLvia the GRHayLibthorn. This is a trimmed-down version of the IllinoisGRMHDthorn where magnetic fields are set to zero.

1 Introduction

This thorn provides a GRHD evolution code built on the General Relativistic Hydrodynamics Library (GRHayL), which is included in the Einstein Toolkit via the GRHayLib thorn. The core library inherits or adapts code from the original IllinoisGRMHD thorn into a modular, infrastructure-agnostic library. While the IllinoisGRMHD thorn can be set to have zero magnetic fields, they will still be evolved. This is much slower than not evolving them, and GRHayLHD provides an implementation for evolving systems without magnetic fields.

This thorn provides support for hybrid and tabulated EOS, both with and without entropy. The core functions all come from GRHayLib, ensuring that any improvements or updates to the library can be easily adopted by the thorn. The conservative-to-primitive solver is selected at runtime, allowing for users to take advantage of any Con2Prim methods provided by GRHayL.

2 Parameters

Most of the behavior of this thorn is controlled by GRHayLib, the thorn that provides GRHayL functionality within the Einstein Toolkit. However, there are several parameters in the thorn which control diagnostic or debugging features. Most are self-explanatory, but some parameters deserve a more detailed explanation.

First, many analysis or diagnostic thorns use HydroBase variables. By default, GRHayLHD never copies the data back from GRHayLHD variables to HydroBase variables, which would prevent the usage of these thorns. The Convert_to_HydroBase_every parameter sets how often to copy back this data. Naturally, copying more frequently will slow down the simulation, so this should be set on the same frequency as the analysis or IO thorn that is using the HydroBase variable. We reduced the amount of data copying in GRHayLHD by using the HydroBase variables directly, which removes the need to copy them back and removes extra grid functions from memory. We still use distinct velocity variables because we use a different velocity than is defined in HydroBase.

Next, there are several parameters for debugging or other testing. The perturbing parameters give users the ability to affect either the initial data or evolution data at runtime without slowing down standard simulations. The perturb_initial_data parameter controls whether the initial data from HydroBase is perturbed, and the perturb_every_con2prim parameter controls whether the conservative variables are perturbed before every con2prim call. The magnitude of these perturbations is set by random_pert.

3 Parameters




convert_to_hydrobase_every
Scope: private INT



Description: How often to convert GRHayLHD primitive variables to HydroBase (Valencia formulation) primitive variables? Needed for some ET-based diagnostics. NOT needed for pure GRHayLHD runs.



Range Default: (none)
0:*
zero (disable) or positive (every N iterations)






matter_bc
Scope: private KEYWORD



Description: Chosen Matter boundary condition



Range Default: outflow
copy
Copy BCs
outflow
Copy BCs with additional outflow BC for velocities
frozen
Frozen boundaries






perturb_every_con2prim
Scope: private BOOLEAN



Description: Random perturbation applied every time con2prim runs



Default: no






perturb_initial_data
Scope: private BOOLEAN



Description: Random perturbation to initial data



Default: no






random_pert
Scope: private REAL



Description: Magnitude of random perturbations



Range Default: (none)
*:*
Anything goes.






random_seed
Scope: private INT



Description: Random seed for random, generally roundoff-level perturbation on initial data. Seeds srand(), and rand() is used for the RNG.



Range Default: (none)
0:99999999
Anything unsigned goes.






symmetry
Scope: private KEYWORD



Description: Currently only no symmetry supported, though work has begun in adding equatorial-symmetry support.



Range Default: none
none
no symmetry, full 3d domain






update_tmunu
Scope: private BOOLEAN



Description: Update Tmunu, for RHS of Einstein’s equations?



Default: yes






verbose
Scope: private KEYWORD



Description: Determines how much evolution information is output



Range Default: yes
no
Complete silence
yes
”Outputs information about conservative-to-prim itive fixes, failures, etc.”






eos_type
Scope: shared from GRHAYLIBKEYWORD






evolve_entropy
Scope: shared from GRHAYLIBBOOLEAN



4 Interfaces

General

Implements:

grhaylhd

Inherits:

admbase

tmunubase

hydrobase

grhaylib

Grid Variables

4.0.1 PRIVATE GROUPS





  Group Names     Variable Names   Details    




grmhd_velocities compact 0
vx description Components of primitive three velocity vî. Note that vî is defined in terms of 4-velocity as: vî = uî/u                                                   ˆ0  . Note that this definition differs from the Valencia formalism.
vy dimensions 3
vz distribution DEFAULT
group type GF
tags InterpNumTimelevels=1 prolongation=”none”
timelevels 1
variable type REAL




u0 u0 compact 0
description 4-velocity component uˆ0
dimensions 3
distribution DEFAULT
group type GF
tags InterpNumTimelevels=1 prolongation=”none” checkpoint=”no”
timelevels 1
variable type REAL




grmhd_conservatives compact 0
rho_star description Evolved hydrodynamic variables
tau dimensions 3
Stildex distribution DEFAULT
Stildey group type GF
Stildez tags prolongation=”ENO”
timelevels 3
variable type REAL




ent_star ent_star compact 0
description Evolved entropy variable
dimensions 3
distribution DEFAULT
group type GF
tags prolongation=”ENO”
timelevels 3
variable type REAL




ye_star Ye_star compact 0
description Evolved electron fraction variable
dimensions 3
distribution DEFAULT
group type GF
tags prolongation=”ENO”
timelevels 3
variable type REAL




grmhd_conservatives_rhs compact 0
rho_star_rhs description Right-hand side for grmhd_conservatives variables
tau_rhs dimensions 3
Stildex_rhs distribution DEFAULT
Stildey_rhs group type GF
Stildez_rhs tags prolongation=”none” Checkpoint=”no”
timelevels 1
variable type REAL








  Group Names     Variable Names   Details    




ent_star_rhs ent_star_rhs compact 0
description Right-hand side for ent_star
dimensions 3
distribution DEFAULT
group type GF
tags prolongation=”none” Checkpoint=”no”
timelevels 1
variable type REAL




ye_star_rhs Ye_star_rhs compact 0
description Right-hand side for Ye_star
dimensions 3
distribution DEFAULT
group type GF
tags prolongation=”none” Checkpoint=”no”
timelevels 1
variable type REAL




grmhd_flux_temps compact 0
rho_star_flux description Flux variables for computing RHSs
tau_flux dimensions 3
Stildex_flux distribution DEFAULT
Stildey_flux group type GF
Stildez_flux tags prolongation=”none” Checkpoint=”no”
timelevels 1
variable type REAL




ent_star_flux ent_star_flux compact 0
description Flux variable for ent_star
dimensions 3
distribution DEFAULT
group type GF
tags prolongation=”none” Checkpoint=”no”
timelevels 1
variable type REAL




ye_star_flux Ye_star_flux compact 0
description Flux variable for Ye_star
dimensions 3
distribution DEFAULT
group type GF
tags prolongation=”none” Checkpoint=”no”
timelevels 1
variable type REAL




failure_checker failure_checker compact 0
description Gridfunction to track conservative-to-primitives solver fixes. Beware that this gridfunction is overwritten at each RK substep.
dimensions 3
distribution DEFAULT
group type GF
tags prolongation=”none” Checkpoint=”no” InterpNumTimelevels=1
timelevels 1
variable type REAL




Uses header:

Symmetry.h

GRHayLib.h

5 Schedule

This section lists all the variables which are assigned storage by thorn GRHayLET/GRHayLHD. Storage can either last for the duration of the run (Always means that if this thorn is activated storage will be assigned, Conditional means that if this thorn is activated storage will be assigned for the duration of the run if some condition is met), or can be turned on for the duration of a schedule function.

Storage

 

Always: Conditional:
grmhd_conservatives[3] grmhd_velocities Ye_star[3] Ye_star_rhs Ye_star_flux
u0 grmhd_conservatives_rhs ent_star[3] ent_star_rhs ent_star_flux
grmhd_flux_temps failure_checker  
   

Scheduled Functions

MoL_Register (conditional)

  grhaylhd_registervars

  register evolved, rhs variables in grhaylhd for mol

 

  Language: c
  Options: meta
  Type: function

CCTK_BASEGRID (conditional)

  grhaylhd_initsymbound

  schedule symmetries

 

  Language: c
  Type: function

GRHayLHD_Con2Prim (conditional)

  grhaylhd_perturb_conservatives

  perturb conservatives before con2prim

 

  After: grhaylhd_sync_conservatives
  Before: grhaylhd_conservs_to_prims
  Type: group

GRHayLHD_Con2Prim (conditional)

  grhaylhd_conservs_to_prims

  compute primitive variables from conservatives

 

  After: grhaylhd_sync_conservatives
  Type: group

GRHayLHD_Con2Prim (conditional)

  grhaylhd_outer_boundaries

  apply selected boundary conditions to the primitive variables

 

  After: grhaylhd_conservs_to_prims
  Type: group

AddToTmunu (conditional)

  grhaylhd_compute_tmunu

  compute stress-energy tensor

 

  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    hydrobase::rho
    hydrobase::press
    hydrobase::eps
    grmhd_velocities
    u0
  Type: function
  Writes: tmunubase::stress_energy_scalar(everywhere)
    tmunubase::stress_energy_vector(everywhere)
    tmunubase::stress_energy_tensor(everywhere)

MoL_CalcRHS (conditional)

  grhaylhd_rhs

  evaluate rhss grhd equations

 

  Type: group

GRHayLHD_RHS (conditional)

  grhaylhd_evaluate_sources_rhs

  evaluate source terms in grhd rhss

 

  Type: group

GRHayLHD_RHS (conditional)

  grhaylhd_evaluate_fluxes_rhs

  evaluate flux terms in grhd rhss

 

  After: grhaylhd_evaluate_sources_rhs
  Type: group

GRHayLHD_RHS (conditional)

  convert_grhaylhd_to_hydrobase

  convert needed hydrobase variables for nrpyleakage

 

  After: grhaylhd_evaluate_fluxes_rhs
  Language: c
  Reads: admbase::lapse
    admbase::shift
    grmhd_velocities
  Type: function
  Writes: hydrobase::vel(everywhere)

CCTK_ANALYSIS (conditional)

  convert_grhaylhd_to_hydrobase

  convert grhaylhd-native variables to hydrobase

 

  After: ml_bssn_evolcalcgroup
  Before: compute_bi_b2_poyn_fluxet
    convert_to_mhd_3velocity
    particle_traceret
    volumeintegralgroup
  Language: c
  Options: global-early
    loop-local
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    grmhd_velocities
  Type: function
  Writes: hydrobase::vel(everywhere)
    hydrobase::w_lorentz(everywhere)

GRHayLHD_prims_to_conservs (conditional)

  grhaylhd_hybrid_entropy_prims_to_conservs

  entropy+hybrid version of grhaylhd_prims_to_conservs

 

  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    hydrobase::rho
    hydrobase::press
    hydrobase::eps
    hydrobase::entropy
    grmhd_velocities
  Sync: grmhd_conservatives
    ent_star
  Type: function
  Writes: u0(everywhere)
    grmhd_velocities(everywhere)
    grmhd_conservatives(everywhere)
    ent_star(everywhere)
    hydrobase::rho(everywhere)
    hydrobase::press(everywhere)
    hydrobase::eps(everywhere)
    hydrobase::entropy(everywhere)

HydroBase_Prim2ConInitial (conditional)

  grhaylhd_prim2con2prim

  convert hydrobase initial data to grhaylhd variables

 

  Type: group

GRHayLHD_conservs_to_prims (conditional)

  grhaylhd_hybrid_entropy_conservs_to_prims

  entropy+hybrid version of grhaylhd_conservs_to_prims

 

  Language: c
  Reads: grid::coordinates
    admbase::metric
    admbase::lapse
    admbase::shift
    grmhd_conservatives
    ent_star
  Type: function
  Writes: grmhd_conservatives(everywhere)
    ent_star(everywhere)
    failure_checker(everywhere)
    u0(everywhere)
    grmhd_velocities(everywhere)
    hydrobase::rho(everywhere)
    hydrobase::press(everywhere)
    hydrobase::eps(everywhere)
    hydrobase::entropy(everywhere)

GRHayLHD_outer_boundaries (conditional)

  grhaylhd_hybrid_entropy_outer_boundaries

  entropy+hybrid version of grhaylhd_outer_boundaries

 

  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    hydrobase::rho
    hydrobase::press
    hydrobase::entropy
    grmhd_velocities
  Sync: grmhd_velocities
    hydrobase::rho
    hydrobase::press
    hydrobase::eps
    hydrobase::entropy
  Type: function
  Writes: u0(everywhere)
    grmhd_conservatives(everywhere)
    ent_star(everywhere)
    grmhd_velocities(everywhere)
    hydrobase::rho(everywhere)
    hydrobase::press(everywhere)
    hydrobase::eps(everywhere)
    hydrobase::entropy(everywhere)

GRHayLHD_evaluate_sources_rhs (conditional)

  grhaylhd_hybrid_entropy_evaluate_sources_rhs

  entropy+hybrid version of grhaylhd_evaluate_sources_rhs

 

  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    admbase::curv
    hydrobase::rho
    hydrobase::press
    hydrobase::entropy
    grmhd_velocities
  Type: function
  Writes: grmhd_conservatives_rhs
    ent_star_rhs

GRHayLHD_evaluate_fluxes_rhs (conditional)

  grhaylhd_hybrid_entropy_evaluate_fluxes_rhs

  entropy+hybrid version of grhaylhd_evaluate_fluxes_rhs

 

  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    hydrobase::rho
    hydrobase::press
    hydrobase::entropy
    grmhd_velocities
    grmhd_conservatives_rhs(interior)
    ent_star_rhs(interior)
  Type: function
  Writes: grmhd_flux_temps
    grmhd_conservatives_rhs
    ent_star_flux
    ent_star_rhs

GRHayLHD_perturb_primitives (conditional)

  grhaylhd_hybrid_entropy_perturb_primitives

  entropy+hybrid version of grhaylhd_perturb_primitives

 

  Language: c
  Reads: hydrobase::rho
    hydrobase::press
    hydrobase::eps
    hydrobase::entropy
    grmhd_velocities
  Type: function
  Writes: hydrobase::rho(everywhere)
    hydrobase::press(everywhere)
    hydrobase::eps(everywhere)
    hydrobase::entropy
    grmhd_velocities(everywhere)

GRHayLHD_perturb_conservatives (conditional)

  grhaylhd_hybrid_entropy_perturb_conservatives

  entropy+hybrid version of grhaylhd_perturb_conservatives

 

  Language: c
  Reads: grmhd_conservatives
    ent_star
  Type: function
  Writes: grmhd_conservatives(everywhere)
    ent_star(everywhere)

GRHayLHD_prims_to_conservs (conditional)

  grhaylhd_hybrid_prims_to_conservs

  hybrid version of grhaylhd_prims_to_conservs

 

  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    hydrobase::rho
    hydrobase::press
    hydrobase::eps
    grmhd_velocities
  Sync: grmhd_conservatives
  Type: function
  Writes: u0(everywhere)
    grmhd_velocities(everywhere)
    grmhd_conservatives(everywhere)
    hydrobase::rho(everywhere)
    hydrobase::press(everywhere)
    hydrobase::eps(everywhere)

GRHayLHD_conservs_to_prims (conditional)

  grhaylhd_hybrid_conservs_to_prims

  hybrid version of grhaylhd_conservs_to_prims

 

  Language: c
  Reads: grid::coordinates
    admbase::metric
    admbase::lapse
    admbase::shift
    grmhd_conservatives
  Type: function
  Writes: grmhd_conservatives(everywhere)
    failure_checker(everywhere)
    u0(everywhere)
    grmhd_velocities(everywhere)
    hydrobase::rho(everywhere)
    hydrobase::press(everywhere)
    hydrobase::eps(everywhere)

GRHayLHD_outer_boundaries (conditional)

  grhaylhd_hybrid_outer_boundaries

  hybrid version of grhaylhd_outer_boundaries

 

  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    hydrobase::rho
    hydrobase::press
    grmhd_velocities
  Sync: grmhd_velocities
    hydrobase::rho
    hydrobase::press
    hydrobase::eps
  Type: function
  Writes: u0(everywhere)
    grmhd_conservatives(everywhere)
    grmhd_velocities(everywhere)
    hydrobase::rho(everywhere)
    hydrobase::press(everywhere)
    hydrobase::eps(everywhere)

GRHayLHD_evaluate_sources_rhs (conditional)

  grhaylhd_hybrid_evaluate_sources_rhs

  hybrid version of grhaylhd_evaluate_sources_rhs

 

  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    admbase::curv
    hydrobase::rho
    hydrobase::press
    grmhd_velocities
  Type: function
  Writes: grmhd_conservatives_rhs

GRHayLHD_Prim2Con2Prim (conditional)

  convert_hydrobase_to_grhaylhd

  convert hydrobase initial data to grhaylhd variables

 

  Language: c
  Reads: admbase::lapse
    admbase::shift
    hydrobase::vel
  Sync: grmhd_velocities
  Type: function
  Writes: grmhd_velocities(everywhere)

GRHayLHD_evaluate_fluxes_rhs (conditional)

  grhaylhd_hybrid_evaluate_fluxes_rhs

  hybrid version of grhaylhd_evaluate_fluxes_rhs

 

  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    hydrobase::rho
    hydrobase::press
    grmhd_velocities
    grmhd_conservatives_rhs(interior)
  Type: function
  Writes: grmhd_flux_temps
    grmhd_conservatives_rhs

GRHayLHD_perturb_primitives (conditional)

  grhaylhd_hybrid_perturb_primitives

  hybrid version of grhaylhd_perturb_primitives

 

  Language: c
  Reads: hydrobase::rho
    hydrobase::press
    hydrobase::eps
    grmhd_velocities
  Type: function
  Writes: hydrobase::rho(everywhere)
    hydrobase::press(everywhere)
    hydrobase::eps(everywhere)
    grmhd_velocities(everywhere)

GRHayLHD_perturb_conservatives (conditional)

  grhaylhd_hybrid_perturb_conservatives

  hybrid version of grhaylhd_perturb_conservatives

 

  Language: c
  Reads: grmhd_conservatives
  Type: function
  Writes: grmhd_conservatives(everywhere)

GRHayLHD_prims_to_conservs (conditional)

  grhaylhd_tabulated_entropy_prims_to_conservs

  entropy+tabulated version of grhaylhd_prims_to_conservs

 

  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    hydrobase::rho
    hydrobase::press
    hydrobase::eps
    hydrobase::entropy
    hydrobase::y_e
    hydrobase::temperature
    grmhd_velocities
  Sync: grmhd_conservatives
    ent_star
    ye_star
  Type: function
  Writes: u0(everywhere)
    grmhd_velocities(everywhere)
    grmhd_conservatives(everywhere)
    ent_star(everywhere)
    ye_star(everywhere)
    hydrobase::rho(everywhere)
    hydrobase::press(everywhere)
    hydrobase::eps(everywhere)
    hydrobase::entropy(everywhere)
    hydrobase::y_e(everywhere)
    hydrobase::temperature(everywhere)

GRHayLHD_conservs_to_prims (conditional)

  grhaylhd_tabulated_entropy_conservs_to_prims

  entropy+tabulated version of grhaylhd_conservs_to_prims

 

  Language: c
  Reads: grid::coordinates
    admbase::metric
    admbase::lapse
    admbase::shift
    grmhd_conservatives
    ent_star
    ye_star
  Type: function
  Writes: grmhd_conservatives(everywhere)
    ent_star(everywhere)
    ye_star(everywhere)
    u0(everywhere)
    grmhd_velocities(everywhere)
    failure_checker(everywhere)
    hydrobase::rho(everywhere)
    hydrobase::press(everywhere)
    hydrobase::eps(everywhere)
    hydrobase::entropy(everywhere)
    hydrobase::y_e(everywhere)
    hydrobase::temperature(everywhere)

GRHayLHD_outer_boundaries (conditional)

  grhaylhd_tabulated_entropy_outer_boundaries

  entropy+tabulated version of grhaylhd_outer_boundaries

 

  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    hydrobase::rho
    hydrobase::press
    hydrobase::entropy
    hydrobase::y_e
    hydrobase::temperature
    grmhd_velocities
  Sync: grmhd_velocities
    hydrobase::rho
    hydrobase::eps
    hydrobase::press
    hydrobase::entropy
    hydrobase::y_e
    hydrobase::temperature
  Type: function
  Writes: u0(everywhere)
    grmhd_conservatives(everywhere)
    ent_star(everywhere)
    ye_star(everywhere)
    grmhd_velocities(everywhere)
    hydrobase::rho(everywhere)
    hydrobase::press(everywhere)
    hydrobase::eps(everywhere)
    hydrobase::entropy(everywhere)
    hydrobase::y_e(everywhere)
    hydrobase::temperature(everywhere)

GRHayLHD_evaluate_sources_rhs (conditional)

  grhaylhd_tabulated_entropy_evaluate_sources_rhs

  entropy+tabulated version of grhaylhd_evaluate_sources_rhs

 

  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    admbase::curv
    hydrobase::rho
    hydrobase::press
    hydrobase::entropy
    hydrobase::y_e
    hydrobase::temperature
    grmhd_velocities
  Type: function
  Writes: grmhd_conservatives_rhs
    ent_star_rhs
    ye_star_rhs

GRHayLHD_evaluate_fluxes_rhs (conditional)

  grhaylhd_tabulated_entropy_evaluate_fluxes_rhs

  entropy+tabulated version of grhaylhd_evaluate_fluxes_rhs

 

  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    hydrobase::rho
    hydrobase::press
    hydrobase::entropy
    hydrobase::y_e
    hydrobase::temperature
    grmhd_velocities
    grmhd_conservatives_rhs(interior)
    ent_star_rhs(interior)
    ye_star_rhs(interior)
  Type: function
  Writes: grmhd_flux_temps
    ent_star_flux
    ye_star_flux
    grmhd_conservatives_rhs
    ent_star_rhs
    ye_star_rhs

GRHayLHD_perturb_primitives (conditional)

  grhaylhd_tabulated_entropy_perturb_primitives

  entropy+tabulated version of grhaylhd_perturb_primitives

 

  Language: c
  Reads: hydrobase::rho
    hydrobase::press
    hydrobase::eps
    hydrobase::entropy
    hydrobase::y_e
    hydrobase::temperature
    grmhd_velocities
  Type: function
  Writes: hydrobase::rho(everywhere)
    hydrobase::press(everywhere)
    hydrobase::eps(everywhere)
    hydrobase::entropy
    hydrobase::y_e(everywhere)
    hydrobase::temperature(everywhere)
    grmhd_velocities(everywhere)

GRHayLHD_perturb_conservatives (conditional)

  grhaylhd_tabulated_entropy_perturb_conservatives

  entropy+tabulated version of grhaylhd_perturb_conservatives

 

  Language: c
  Reads: grmhd_conservatives
    ent_star
    ye_star
  Type: function
  Writes: grmhd_conservatives(everywhere)
    ent_star(everywhere)
    ye_star(everywhere)

GRHayLHD_Prim2Con2Prim (conditional)

  grhaylhd_perturb_primitives

  perturb initial primitive data

 

  After: convert_hydrobase_to_grhaylhd
  Before: grhaylhd_prims_to_conservs
  Type: group

GRHayLHD_prims_to_conservs (conditional)

  grhaylhd_tabulated_prims_to_conservs

  tabulated version of grhaylhd_prims_to_conservs

 

  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    hydrobase::rho
    hydrobase::press
    hydrobase::eps
    hydrobase::y_e
    hydrobase::temperature
    grmhd_velocities
  Sync: grmhd_conservatives
    ye_star
  Type: function
  Writes: u0(everywhere)
    grmhd_velocities(everywhere)
    grmhd_conservatives(everywhere)
    ye_star(everywhere)
    hydrobase::rho(everywhere)
    hydrobase::press(everywhere)
    hydrobase::eps(everywhere)
    hydrobase::y_e(everywhere)
    hydrobase::temperature(everywhere)

GRHayLHD_conservs_to_prims (conditional)

  grhaylhd_tabulated_conservs_to_prims

  tabulated version of grhaylhd_conservs_to_prims

 

  Language: c
  Reads: grid::coordinates
    admbase::metric
    admbase::lapse
    admbase::shift
    grmhd_conservatives
    ye_star
  Type: function
  Writes: grmhd_conservatives(everywhere)
    ye_star(everywhere)
    u0(everywhere)
    grmhd_velocities(everywhere)
    failure_checker(everywhere)
    hydrobase::rho(everywhere)
    hydrobase::press(everywhere)
    hydrobase::eps(everywhere)
    hydrobase::y_e(everywhere)
    hydrobase::temperature(everywhere)

GRHayLHD_outer_boundaries (conditional)

  grhaylhd_tabulated_outer_boundaries

  tabulated version of grhaylhd_outer_boundaries

 

  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    hydrobase::rho
    hydrobase::press
    hydrobase::y_e
    hydrobase::temperature
    grmhd_velocities
  Sync: grmhd_velocities
    hydrobase::rho
    hydrobase::press
    hydrobase::eps
    hydrobase::y_e
    hydrobase::temperature
  Type: function
  Writes: u0(everywhere)
    grmhd_conservatives(everywhere)
    ye_star(everywhere)
    grmhd_velocities(everywhere)
    hydrobase::rho(everywhere)
    hydrobase::press(everywhere)
    hydrobase::eps(everywhere)
    hydrobase::y_e(everywhere)
    hydrobase::temperature(everywhere)

GRHayLHD_evaluate_sources_rhs (conditional)

  grhaylhd_tabulated_evaluate_sources_rhs

  tabulated version of grhaylhd_evaluate_sources_rhs

 

  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    admbase::curv
    hydrobase::rho
    hydrobase::press
    hydrobase::y_e
    hydrobase::temperature
    grmhd_velocities
  Type: function
  Writes: grmhd_conservatives_rhs
    ye_star_rhs

GRHayLHD_evaluate_fluxes_rhs (conditional)

  grhaylhd_tabulated_evaluate_fluxes_rhs

  tabulated version of grhaylhd_evaluate_fluxes_rhs

 

  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    hydrobase::rho
    hydrobase::press
    hydrobase::y_e
    hydrobase::temperature
    grmhd_velocities
    grmhd_conservatives_rhs(interior)
    ye_star_rhs(interior)
  Type: function
  Writes: grmhd_flux_temps
    grmhd_conservatives_rhs
    ye_star_flux
    ye_star_rhs

GRHayLHD_perturb_primitives (conditional)

  grhaylhd_tabulated_perturb_primitives

  tabulated version of grhaylhd_perturb_primitives

 

  Language: c
  Reads: hydrobase::rho
    hydrobase::press
    hydrobase::eps
    hydrobase::y_e
    hydrobase::temperature
    grmhd_velocities
  Type: function
  Writes: grmhd_velocities(everywhere)
    hydrobase::rho(everywhere)
    hydrobase::press(everywhere)
    hydrobase::eps(everywhere)
    hydrobase::y_e(everywhere)
    hydrobase::temperature(everywhere)

GRHayLHD_perturb_conservatives (conditional)

  grhaylhd_tabulated_perturb_conservatives

  tabulated version of grhaylhd_perturb_conservatives

 

  Language: c
  Reads: grmhd_conservatives
    ye_star
  Type: function
  Writes: grmhd_conservatives(everywhere)
    ye_star(everywhere)

GRHayLHD_Prim2Con2Prim (conditional)

  grhaylhd_prims_to_conservs

  compute conservative variables from primitives

 

  After: convert_hydrobase_to_grhaylhd
  Type: group

GRHayLHD_Prim2Con2Prim (conditional)

  grhaylhd_conservs_to_prims

  compute primitive variables from conservatives

 

  After: grhaylhd_prims_to_conservs
  Type: group

GRHayLHD_Prim2Con2Prim (conditional)

  convert_grhaylhd_to_hydrobase

  convert grhaylhd-native variables to hydrobase

 

  After: grhaylhd_conservs_to_prims
  Language: c
  Reads: admbase::metric
    admbase::lapse
    admbase::shift
    grmhd_velocities
  Type: function
  Writes: hydrobase::vel(everywhere)
    hydrobase::w_lorentz(everywhere)

HydroBase_Con2Prim (conditional)

  grhaylhd_con2prim

  compute primitive variables from conservatives

 

  Type: group

GRHayLHD_Con2Prim (conditional)

  grhaylhd_empty_function

  sync conservative variable ghost zones for c2p routine

 

  Language: c
  Sync: grmhd_conservatives
    ye_star
    ent_star
  Type: function

Aliased Functions

 

Alias Name:         Function Name:
GRHayLHD_empty_function GRHayLHD_sync_conservatives