gh Namespace Reference

Items related to evolving the first-order generalized harmonic system. More...

Namespaces
namespace	BoundaryConditions
	Boundary conditions for the generalized harmonic system.
namespace	BoundaryCorrections
	Boundary corrections/numerical fluxes.
namespace	ConstraintDamping
	Holds classes implementing DampingFunction (functions $R^n\to R$).
namespace	gauges
	Gauge conditions for generalized harmonic evolution systems.
namespace	OptionTags
	Input option tags for the generalized harmonic evolution system.
namespace	Tags
	Tags for the generalized harmonic formulation of Einstein equations.

Classes
struct	CharacteristicFieldsCompute
struct	CharacteristicSpeedsCompute
struct	CharacteristicSpeedsOnStrahlkorper
struct	CharacteristicSpeedsOnStrahlkorperCompute
struct	ComputeNormalDotFluxes
	Set the normal dot the flux to zero since the generalized harmonic system has no fluxes and they're currently still needed for the evolution scheme. More...
struct	EvolvedFieldsFromCharacteristicFieldsCompute
class	NumericInitialData
	Numeric initial data loaded from volume data files. More...
struct	System
struct	TimeDerivative
	Compute the RHS of the Generalized Harmonic formulation of Einstein's equations. More...

Functions
template<size_t Dim>
void	initial_gh_variables_from_adm (gsl::not_null< tnsr::aa< DataVector, Dim > * > spacetime_metric, gsl::not_null< tnsr::aa< DataVector, Dim > * > pi, gsl::not_null< tnsr::iaa< DataVector, Dim > * > phi, const tnsr::ii< DataVector, Dim > &spatial_metric, const Scalar< DataVector > &lapse, const tnsr::I< DataVector, Dim > &shift, const tnsr::ii< DataVector, Dim > &extrinsic_curvature, const Mesh< Dim > &mesh, const InverseJacobian< DataVector, Dim, Frame::ElementLogical, Frame::Inertial > &inv_jacobian)
	Compute initial GH variables from ADM variables. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::ijj< DataType, SpatialDim, Frame >	covariant_deriv_of_extrinsic_curvature (const tnsr::ii< DataType, SpatialDim, Frame > &extrinsic_curvature, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_unit_normal_vector, const tnsr::Ijj< DataType, SpatialDim, Frame > &spatial_christoffel_second_kind, const tnsr::AA< DataType, SpatialDim, Frame > &inverse_spacetime_metric, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::iaa< DataType, SpatialDim, Frame > &d_pi, const tnsr::ijaa< DataType, SpatialDim, Frame > &d_phi)
	Computes the covariant derivative of extrinsic curvature from generalized harmonic variables and the spacetime normal vector. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	covariant_deriv_of_extrinsic_curvature (gsl::not_null< tnsr::ijj< DataType, SpatialDim, Frame > * > d_extrinsic_curvature, const tnsr::ii< DataType, SpatialDim, Frame > &extrinsic_curvature, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_unit_normal_vector, const tnsr::Ijj< DataType, SpatialDim, Frame > &spatial_christoffel_second_kind, const tnsr::AA< DataType, SpatialDim, Frame > &inverse_spacetime_metric, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::iaa< DataType, SpatialDim, Frame > &d_pi, const tnsr::ijaa< DataType, SpatialDim, Frame > &d_phi)
template<typename DataType , size_t SpatialDim, typename Frame >
void	spacetime_derivative_of_spacetime_metric (gsl::not_null< tnsr::abb< DataType, SpatialDim, Frame > * > da_spacetime_metric, const Scalar< DataType > &lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::aa< DataType, SpatialDim, Frame > &pi, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes the spacetime derivative of the spacetime metric, ${\partial }_a{g}_{bc}$. More...
template<size_t Dim, typename Frame >
std::array< DataVector, 4 >	characteristic_speeds (const Scalar< DataVector > &gamma_1, const Scalar< DataVector > &lapse, const tnsr::I< DataVector, Dim, Frame > &shift, const tnsr::i< DataVector, Dim, Frame > &unit_normal_one_form, const std::optional< tnsr::I< DataVector, Dim, Frame > > &mesh_velocity)
	Compute the characteristic speeds for the generalized harmonic system. More...
template<size_t Dim, typename Frame >
void	characteristic_speeds (gsl::not_null< std::array< DataVector, 4 > * > char_speeds, const Scalar< DataVector > &gamma_1, const Scalar< DataVector > &lapse, const tnsr::I< DataVector, Dim, Frame > &shift, const tnsr::i< DataVector, Dim, Frame > &unit_normal_one_form, const std::optional< tnsr::I< DataVector, Dim, Frame > > &mesh_velocity)
	Compute the characteristic speeds for the generalized harmonic system. More...
template<size_t Dim, typename Frame >
Tags::CharacteristicFields< DataVector, Dim, Frame >::type	characteristic_fields (const Scalar< DataVector > &gamma_2, const tnsr::II< DataVector, Dim, Frame > &inverse_spatial_metric, const tnsr::aa< DataVector, Dim, Frame > &spacetime_metric, const tnsr::aa< DataVector, Dim, Frame > &pi, const tnsr::iaa< DataVector, Dim, Frame > &phi, const tnsr::i< DataVector, Dim, Frame > &unit_normal_one_form)
	Computes characteristic fields from evolved fields. More...
template<size_t Dim, typename Frame >
void	characteristic_fields (gsl::not_null< typename Tags::CharacteristicFields< DataVector, Dim, Frame >::type * > char_fields, const Scalar< DataVector > &gamma_2, const tnsr::II< DataVector, Dim, Frame > &inverse_spatial_metric, const tnsr::aa< DataVector, Dim, Frame > &spacetime_metric, const tnsr::aa< DataVector, Dim, Frame > &pi, const tnsr::iaa< DataVector, Dim, Frame > &phi, const tnsr::i< DataVector, Dim, Frame > &unit_normal_one_form)
	Computes characteristic fields from evolved fields. More...
template<size_t Dim, typename Frame >
Tags::EvolvedFieldsFromCharacteristicFields< DataVector, Dim, Frame >::type	evolved_fields_from_characteristic_fields (const Scalar< DataVector > &gamma_2, const tnsr::aa< DataVector, Dim, Frame > &u_psi, const tnsr::iaa< DataVector, Dim, Frame > &u_zero, const tnsr::aa< DataVector, Dim, Frame > &u_plus, const tnsr::aa< DataVector, Dim, Frame > &u_minus, const tnsr::i< DataVector, Dim, Frame > &unit_normal_one_form)
	For expressions used here to compute evolved fields from characteristic ones, see CharacteristicFieldsCompute.
template<size_t Dim, typename Frame >
void	evolved_fields_from_characteristic_fields (gsl::not_null< typename Tags::EvolvedFieldsFromCharacteristicFields< DataVector, Dim, Frame >::type * > evolved_fields, const Scalar< DataVector > &gamma_2, const tnsr::aa< DataVector, Dim, Frame > &u_psi, const tnsr::iaa< DataVector, Dim, Frame > &u_zero, const tnsr::aa< DataVector, Dim, Frame > &u_plus, const tnsr::aa< DataVector, Dim, Frame > &u_minus, const tnsr::i< DataVector, Dim, Frame > &unit_normal_one_form)
	For expressions used here to compute evolved fields from characteristic ones, see CharacteristicFieldsCompute.
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::iaa< DataType, SpatialDim, Frame >	three_index_constraint (const tnsr::iaa< DataType, SpatialDim, Frame > &d_spacetime_metric, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes the generalized-harmonic 3-index constraint. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	three_index_constraint (gsl::not_null< tnsr::iaa< DataType, SpatialDim, Frame > * > constraint, const tnsr::iaa< DataType, SpatialDim, Frame > &d_spacetime_metric, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes the generalized-harmonic 3-index constraint. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::a< DataType, SpatialDim, Frame >	gauge_constraint (const tnsr::a< DataType, SpatialDim, Frame > &gauge_function, const tnsr::a< DataType, SpatialDim, Frame > &spacetime_normal_one_form, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_normal_vector, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const tnsr::AA< DataType, SpatialDim, Frame > &inverse_spacetime_metric, const tnsr::aa< DataType, SpatialDim, Frame > &pi, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes the generalized-harmonic gauge constraint. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	gauge_constraint (gsl::not_null< tnsr::a< DataType, SpatialDim, Frame > * > constraint, const tnsr::a< DataType, SpatialDim, Frame > &gauge_function, const tnsr::a< DataType, SpatialDim, Frame > &spacetime_normal_one_form, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_normal_vector, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const tnsr::AA< DataType, SpatialDim, Frame > &inverse_spacetime_metric, const tnsr::aa< DataType, SpatialDim, Frame > &pi, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes the generalized-harmonic gauge constraint. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::ia< DataType, SpatialDim, Frame >	two_index_constraint (const tnsr::ab< DataType, SpatialDim, Frame > &spacetime_d_gauge_function, const tnsr::a< DataType, SpatialDim, Frame > &spacetime_normal_one_form, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_normal_vector, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const tnsr::AA< DataType, SpatialDim, Frame > &inverse_spacetime_metric, const tnsr::aa< DataType, SpatialDim, Frame > &pi, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::iaa< DataType, SpatialDim, Frame > &d_pi, const tnsr::ijaa< DataType, SpatialDim, Frame > &d_phi, const Scalar< DataType > &gamma2, const tnsr::iaa< DataType, SpatialDim, Frame > &three_index_constraint)
	Computes the generalized-harmonic 2-index constraint. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	two_index_constraint (gsl::not_null< tnsr::ia< DataType, SpatialDim, Frame > * > constraint, const tnsr::ab< DataType, SpatialDim, Frame > &spacetime_d_gauge_function, const tnsr::a< DataType, SpatialDim, Frame > &spacetime_normal_one_form, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_normal_vector, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const tnsr::AA< DataType, SpatialDim, Frame > &inverse_spacetime_metric, const tnsr::aa< DataType, SpatialDim, Frame > &pi, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::iaa< DataType, SpatialDim, Frame > &d_pi, const tnsr::ijaa< DataType, SpatialDim, Frame > &d_phi, const Scalar< DataType > &gamma2, const tnsr::iaa< DataType, SpatialDim, Frame > &three_index_constraint)
	Computes the generalized-harmonic 2-index constraint. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::iaa< DataType, SpatialDim, Frame >	four_index_constraint (const tnsr::ijaa< DataType, SpatialDim, Frame > &d_phi)
	Computes the generalized-harmonic 4-index constraint. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	four_index_constraint (gsl::not_null< tnsr::iaa< DataType, SpatialDim, Frame > * > constraint, const tnsr::ijaa< DataType, SpatialDim, Frame > &d_phi)
	Computes the generalized-harmonic 4-index constraint. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::a< DataType, SpatialDim, Frame >	f_constraint (const tnsr::a< DataType, SpatialDim, Frame > &gauge_function, const tnsr::ab< DataType, SpatialDim, Frame > &spacetime_d_gauge_function, const tnsr::a< DataType, SpatialDim, Frame > &spacetime_normal_one_form, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_normal_vector, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const tnsr::AA< DataType, SpatialDim, Frame > &inverse_spacetime_metric, const tnsr::aa< DataType, SpatialDim, Frame > &pi, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::iaa< DataType, SpatialDim, Frame > &d_pi, const tnsr::ijaa< DataType, SpatialDim, Frame > &d_phi, const Scalar< DataType > &gamma2, const tnsr::iaa< DataType, SpatialDim, Frame > &three_index_constraint)
	Computes the generalized-harmonic F constraint. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	f_constraint (gsl::not_null< tnsr::a< DataType, SpatialDim, Frame > * > constraint, const tnsr::a< DataType, SpatialDim, Frame > &gauge_function, const tnsr::ab< DataType, SpatialDim, Frame > &spacetime_d_gauge_function, const tnsr::a< DataType, SpatialDim, Frame > &spacetime_normal_one_form, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_normal_vector, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const tnsr::AA< DataType, SpatialDim, Frame > &inverse_spacetime_metric, const tnsr::aa< DataType, SpatialDim, Frame > &pi, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::iaa< DataType, SpatialDim, Frame > &d_pi, const tnsr::ijaa< DataType, SpatialDim, Frame > &d_phi, const Scalar< DataType > &gamma2, const tnsr::iaa< DataType, SpatialDim, Frame > &three_index_constraint)
	Computes the generalized-harmonic F constraint. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::a< DataType, SpatialDim, Frame >	f_constraint (const tnsr::a< DataType, SpatialDim, Frame > &gauge_function, const tnsr::ab< DataType, SpatialDim, Frame > &spacetime_d_gauge_function, const tnsr::a< DataType, SpatialDim, Frame > &spacetime_normal_one_form, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_normal_vector, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const tnsr::AA< DataType, SpatialDim, Frame > &inverse_spacetime_metric, const tnsr::aa< DataType, SpatialDim, Frame > &pi, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::iaa< DataType, SpatialDim, Frame > &d_pi, const tnsr::ijaa< DataType, SpatialDim, Frame > &d_phi, const Scalar< DataType > &gamma2, const tnsr::iaa< DataType, SpatialDim, Frame > &three_index_constraint, const tnsr::aa< DataType, SpatialDim, Frame > &trace_reversed_stress_energy)
	Computes the generalized-harmonic F constraint. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	f_constraint (gsl::not_null< tnsr::a< DataType, SpatialDim, Frame > * > constraint, const tnsr::a< DataType, SpatialDim, Frame > &gauge_function, const tnsr::ab< DataType, SpatialDim, Frame > &spacetime_d_gauge_function, const tnsr::a< DataType, SpatialDim, Frame > &spacetime_normal_one_form, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_normal_vector, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const tnsr::AA< DataType, SpatialDim, Frame > &inverse_spacetime_metric, const tnsr::aa< DataType, SpatialDim, Frame > &pi, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::iaa< DataType, SpatialDim, Frame > &d_pi, const tnsr::ijaa< DataType, SpatialDim, Frame > &d_phi, const Scalar< DataType > &gamma2, const tnsr::iaa< DataType, SpatialDim, Frame > &three_index_constraint, const tnsr::aa< DataType, SpatialDim, Frame > &trace_reversed_stress_energy)
	Computes the generalized-harmonic F constraint. More...
template<typename DataType , size_t SpatialDim, typename Frame >
Scalar< DataType >	constraint_energy (const tnsr::a< DataType, SpatialDim, Frame > &gauge_constraint, const tnsr::a< DataType, SpatialDim, Frame > &f_constraint, const tnsr::ia< DataType, SpatialDim, Frame > &two_index_constraint, const tnsr::iaa< DataType, SpatialDim, Frame > &three_index_constraint, const tnsr::iaa< DataType, SpatialDim, Frame > &four_index_constraint, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const Scalar< DataType > &spatial_metric_determinant, double gauge_constraint_multiplier=1.0, double two_index_constraint_multiplier=1.0, double three_index_constraint_multiplier=1.0, double four_index_constraint_multiplier=1.0)
	Computes the generalized-harmonic (unnormalized) constraint energy. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	constraint_energy (gsl::not_null< Scalar< DataType > * > energy, const tnsr::a< DataType, SpatialDim, Frame > &gauge_constraint, const tnsr::a< DataType, SpatialDim, Frame > &f_constraint, const tnsr::ia< DataType, SpatialDim, Frame > &two_index_constraint, const tnsr::iaa< DataType, SpatialDim, Frame > &three_index_constraint, const tnsr::iaa< DataType, SpatialDim, Frame > &four_index_constraint, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const Scalar< DataType > &spatial_metric_determinant, double gauge_constraint_multiplier=1.0, double two_index_constraint_multiplier=1.0, double three_index_constraint_multiplier=1.0, double four_index_constraint_multiplier=1.0)
	Computes the generalized-harmonic (unnormalized) constraint energy. More...
template<typename DataType , size_t SpatialDim, typename Frame >
Scalar< DataType >	constraint_energy_normalization (const tnsr::iaa< DataType, SpatialDim, Frame > &d_spacetime_metric, const tnsr::iaa< DataType, SpatialDim, Frame > &d_pi, const tnsr::ijaa< DataType, SpatialDim, Frame > &d_phi, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const Scalar< DataType > &sqrt_spatial_metric_determinant, double dimensional_constant)
	Computes the generalized-harmonic normalized constraint energy. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	constraint_energy_normalization (gsl::not_null< Scalar< DataType > * > energy_norm, const tnsr::iaa< DataType, SpatialDim, Frame > &d_spacetime_metric, const tnsr::iaa< DataType, SpatialDim, Frame > &d_pi, const tnsr::ijaa< DataType, SpatialDim, Frame > &d_phi, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const Scalar< DataType > &sqrt_spatial_metric_determinant, double dimensional_constant)
	Computes the generalized-harmonic normalized constraint energy. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	christoffel_second_kind (const gsl::not_null< tnsr::Ijj< DataType, SpatialDim, Frame > * > christoffel, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::II< DataType, SpatialDim, Frame > &inv_metric)
	Computes spatial Christoffel symbol of the 2nd kind from the the generalized harmonic spatial derivative variable and the inverse spatial metric. More...
template<typename DataType , size_t SpatialDim, typename Frame >
auto	christoffel_second_kind (const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::II< DataType, SpatialDim, Frame > &inv_metric) -> tnsr::Ijj< DataType, SpatialDim, Frame >
	Computes spatial Christoffel symbol of the 2nd kind from the the generalized harmonic spatial derivative variable and the inverse spatial metric. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::a< DataType, SpatialDim, Frame >	trace_christoffel (const tnsr::a< DataType, SpatialDim, Frame > &spacetime_normal_one_form, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_normal_vector, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const tnsr::AA< DataType, SpatialDim, Frame > &inverse_spacetime_metric, const tnsr::aa< DataType, SpatialDim, Frame > &pi, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Compute ${\mathrm{\Gamma }}_a$ from the generalized harmonic evolved variables. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	trace_christoffel (gsl::not_null< tnsr::a< DataType, SpatialDim, Frame > * > trace, const tnsr::a< DataType, SpatialDim, Frame > &spacetime_normal_one_form, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_normal_vector, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const tnsr::AA< DataType, SpatialDim, Frame > &inverse_spacetime_metric, const tnsr::aa< DataType, SpatialDim, Frame > &pi, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Compute ${\mathrm{\Gamma }}_a$ from the generalized harmonic evolved variables. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	deriv_spatial_metric (gsl::not_null< tnsr::ijj< DataType, SpatialDim, Frame > * > d_spatial_metric, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes spatial derivatives of the spatial metric from the generalized harmonic spatial derivative variable. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::ijj< DataType, SpatialDim, Frame >	deriv_spatial_metric (const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes spatial derivatives of the spatial metric from the generalized harmonic spatial derivative variable. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	extrinsic_curvature (gsl::not_null< tnsr::ii< DataType, SpatialDim, Frame > * > ex_curv, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_normal_vector, const tnsr::aa< DataType, SpatialDim, Frame > &pi, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes extrinsic curvature from generalized harmonic variables and the spacetime normal vector. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::ii< DataType, SpatialDim, Frame >	extrinsic_curvature (const tnsr::A< DataType, SpatialDim, Frame > &spacetime_normal_vector, const tnsr::aa< DataType, SpatialDim, Frame > &pi, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes extrinsic curvature from generalized harmonic variables and the spacetime normal vector. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	gauge_source (gsl::not_null< tnsr::a< DataType, SpatialDim, Frame > * > gauge_source_h, const Scalar< DataType > &lapse, const Scalar< DataType > &dt_lapse, const tnsr::i< DataType, SpatialDim, Frame > &deriv_lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::I< DataType, SpatialDim, Frame > &dt_shift, const tnsr::iJ< DataType, SpatialDim, Frame > &deriv_shift, const tnsr::ii< DataType, SpatialDim, Frame > &spatial_metric, const Scalar< DataType > &trace_extrinsic_curvature, const tnsr::i< DataType, SpatialDim, Frame > &trace_christoffel_last_indices)
	Computes generalized harmonic gauge source function. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::a< DataType, SpatialDim, Frame >	gauge_source (const Scalar< DataType > &lapse, const Scalar< DataType > &dt_lapse, const tnsr::i< DataType, SpatialDim, Frame > &deriv_lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::I< DataType, SpatialDim, Frame > &dt_shift, const tnsr::iJ< DataType, SpatialDim, Frame > &deriv_shift, const tnsr::ii< DataType, SpatialDim, Frame > &spatial_metric, const Scalar< DataType > &trace_extrinsic_curvature, const tnsr::i< DataType, SpatialDim, Frame > &trace_christoffel_last_indices)
	Computes generalized harmonic gauge source function. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	phi (gsl::not_null< tnsr::iaa< DataType, SpatialDim, Frame > * > phi, const Scalar< DataType > &lapse, const tnsr::i< DataType, SpatialDim, Frame > &deriv_lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::iJ< DataType, SpatialDim, Frame > &deriv_shift, const tnsr::ii< DataType, SpatialDim, Frame > &spatial_metric, const tnsr::ijj< DataType, SpatialDim, Frame > &deriv_spatial_metric)
	Computes the auxiliary variable ${\mathrm{\Phi }}_{iab}$ used by the generalized harmonic formulation of Einstein's equations. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::iaa< DataType, SpatialDim, Frame >	phi (const Scalar< DataType > &lapse, const tnsr::i< DataType, SpatialDim, Frame > &deriv_lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::iJ< DataType, SpatialDim, Frame > &deriv_shift, const tnsr::ii< DataType, SpatialDim, Frame > &spatial_metric, const tnsr::ijj< DataType, SpatialDim, Frame > &deriv_spatial_metric)
	Computes the auxiliary variable ${\mathrm{\Phi }}_{iab}$ used by the generalized harmonic formulation of Einstein's equations. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	pi (gsl::not_null< tnsr::aa< DataType, SpatialDim, Frame > * > pi, const Scalar< DataType > &lapse, const Scalar< DataType > &dt_lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::I< DataType, SpatialDim, Frame > &dt_shift, const tnsr::ii< DataType, SpatialDim, Frame > &spatial_metric, const tnsr::ii< DataType, SpatialDim, Frame > &dt_spatial_metric, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes the conjugate momentum ${\mathrm{\Pi }}_{ab}$ of the spacetime metric $g_{ab}$. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::aa< DataType, SpatialDim, Frame >	pi (const Scalar< DataType > &lapse, const Scalar< DataType > &dt_lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::I< DataType, SpatialDim, Frame > &dt_shift, const tnsr::ii< DataType, SpatialDim, Frame > &spatial_metric, const tnsr::ii< DataType, SpatialDim, Frame > &dt_spatial_metric, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes the conjugate momentum ${\mathrm{\Pi }}_{ab}$ of the spacetime metric $g_{ab}$. More...
template<typename DataType , size_t VolumeDim, typename Frame >
void	spatial_ricci_tensor (gsl::not_null< tnsr::ii< DataType, VolumeDim, Frame > * > ricci, const tnsr::iaa< DataType, VolumeDim, Frame > &phi, const tnsr::ijaa< DataType, VolumeDim, Frame > &deriv_phi, const tnsr::II< DataType, VolumeDim, Frame > &inverse_spatial_metric)
	Compute spatial Ricci tensor using evolved variables and their first derivatives. More...
template<typename DataType , size_t VolumeDim, typename Frame >
tnsr::ii< DataType, VolumeDim, Frame >	spatial_ricci_tensor (const tnsr::iaa< DataType, VolumeDim, Frame > &phi, const tnsr::ijaa< DataType, VolumeDim, Frame > &deriv_phi, const tnsr::II< DataType, VolumeDim, Frame > &inverse_spatial_metric)
	Compute spatial Ricci tensor using evolved variables and their first derivatives. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	second_time_deriv_of_spacetime_metric (gsl::not_null< tnsr::aa< DataType, SpatialDim, Frame > * > d2t2_spacetime_metric, const Scalar< DataType > &lapse, const Scalar< DataType > &dt_lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::I< DataType, SpatialDim, Frame > &dt_shift, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::iaa< DataType, SpatialDim, Frame > &dt_phi, const tnsr::aa< DataType, SpatialDim, Frame > &pi, const tnsr::aa< DataType, SpatialDim, Frame > &dt_pi)
	Computes the second time derivative of the spacetime metric from the generalized harmonic variables, lapse, shift, and the spacetime unit normal 1-form. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::aa< DataType, SpatialDim, Frame >	second_time_deriv_of_spacetime_metric (const Scalar< DataType > &lapse, const Scalar< DataType > &dt_lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::I< DataType, SpatialDim, Frame > &dt_shift, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::iaa< DataType, SpatialDim, Frame > &dt_phi, const tnsr::aa< DataType, SpatialDim, Frame > &pi, const tnsr::aa< DataType, SpatialDim, Frame > &dt_pi)
	Computes the second time derivative of the spacetime metric from the generalized harmonic variables, lapse, shift, and the spacetime unit normal 1-form. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	spacetime_deriv_of_det_spatial_metric (gsl::not_null< tnsr::a< DataType, SpatialDim, Frame > * > d4_det_spatial_metric, const Scalar< DataType > &sqrt_det_spatial_metric, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const tnsr::ii< DataType, SpatialDim, Frame > &dt_spatial_metric, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes spacetime derivatives of the determinant of spatial metric, using the generalized harmonic variables, spatial metric, and its time derivative. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::a< DataType, SpatialDim, Frame >	spacetime_deriv_of_det_spatial_metric (const Scalar< DataType > &sqrt_det_spatial_metric, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const tnsr::ii< DataType, SpatialDim, Frame > &dt_spatial_metric, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes spacetime derivatives of the determinant of spatial metric, using the generalized harmonic variables, spatial metric, and its time derivative. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	spacetime_deriv_of_norm_of_shift (gsl::not_null< tnsr::a< DataType, SpatialDim, Frame > * > d4_norm_of_shift, const Scalar< DataType > &lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::ii< DataType, SpatialDim, Frame > &spatial_metric, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const tnsr::AA< DataType, SpatialDim, Frame > &inverse_spacetime_metric, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_unit_normal, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::aa< DataType, SpatialDim, Frame > &pi)
	Computes spacetime derivatives of the norm of the shift vector. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::a< DataType, SpatialDim, Frame >	spacetime_deriv_of_norm_of_shift (const Scalar< DataType > &lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::ii< DataType, SpatialDim, Frame > &spatial_metric, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const tnsr::AA< DataType, SpatialDim, Frame > &inverse_spacetime_metric, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_unit_normal, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::aa< DataType, SpatialDim, Frame > &pi)
	Computes spacetime derivatives of the norm of the shift vector. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	spatial_deriv_of_lapse (gsl::not_null< tnsr::i< DataType, SpatialDim, Frame > * > deriv_lapse, const Scalar< DataType > &lapse, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_unit_normal, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes spatial derivatives of lapse ( $\alpha$) from the generalized harmonic variables and spacetime unit normal 1-form. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::i< DataType, SpatialDim, Frame >	spatial_deriv_of_lapse (const Scalar< DataType > &lapse, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_unit_normal, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes spatial derivatives of lapse ( $\alpha$) from the generalized harmonic variables and spacetime unit normal 1-form. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	spatial_deriv_of_shift (gsl::not_null< tnsr::iJ< DataType, SpatialDim, Frame > * > deriv_shift, const Scalar< DataType > &lapse, const tnsr::AA< DataType, SpatialDim, Frame > &inverse_spacetime_metric, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_unit_normal, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes spatial derivatives of the shift vector from the generalized harmonic and geometric variables. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::iJ< DataType, SpatialDim, Frame >	spatial_deriv_of_shift (const Scalar< DataType > &lapse, const tnsr::AA< DataType, SpatialDim, Frame > &inverse_spacetime_metric, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_unit_normal, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes spatial derivatives of the shift vector from the generalized harmonic and geometric variables. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	time_derivative_of_spacetime_metric (gsl::not_null< tnsr::aa< DataType, SpatialDim, Frame > * > dt_spacetime_metric, const Scalar< DataType > &lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::aa< DataType, SpatialDim, Frame > &pi, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes the time derivative of the spacetime metric from the generalized harmonic quantities ${\mathrm{\Pi }}_{ab}$, ${\mathrm{\Phi }}_{iab}$, and the lapse $\alpha$ and shift ${\beta }^i$. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::aa< DataType, SpatialDim, Frame >	time_derivative_of_spacetime_metric (const Scalar< DataType > &lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::aa< DataType, SpatialDim, Frame > &pi, const tnsr::iaa< DataType, SpatialDim, Frame > &phi)
	Computes the time derivative of the spacetime metric from the generalized harmonic quantities ${\mathrm{\Pi }}_{ab}$, ${\mathrm{\Phi }}_{iab}$, and the lapse $\alpha$ and shift ${\beta }^i$. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	time_deriv_of_lapse (gsl::not_null< Scalar< DataType > * > dt_lapse, const Scalar< DataType > &lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_unit_normal, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::aa< DataType, SpatialDim, Frame > &pi)
	Computes time derivative of lapse ( $\alpha$) from the generalized harmonic variables, lapse, shift and the spacetime unit normal 1-form. More...
template<typename DataType , size_t SpatialDim, typename Frame >
Scalar< DataType >	time_deriv_of_lapse (const Scalar< DataType > &lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_unit_normal, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::aa< DataType, SpatialDim, Frame > &pi)
	Computes time derivative of lapse ( $\alpha$) from the generalized harmonic variables, lapse, shift and the spacetime unit normal 1-form. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	time_deriv_of_lower_shift (gsl::not_null< tnsr::i< DataType, SpatialDim, Frame > * > dt_lower_shift, const Scalar< DataType > &lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::ii< DataType, SpatialDim, Frame > &spatial_metric, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_unit_normal, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::aa< DataType, SpatialDim, Frame > &pi)
	Computes time derivative of index lowered shift from generalized harmonic variables, spatial metric and its time derivative. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::i< DataType, SpatialDim, Frame >	time_deriv_of_lower_shift (const Scalar< DataType > &lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::ii< DataType, SpatialDim, Frame > &spatial_metric, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_unit_normal, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::aa< DataType, SpatialDim, Frame > &pi)
	Computes time derivative of index lowered shift from generalized harmonic variables, spatial metric and its time derivative. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	time_deriv_of_shift (gsl::not_null< tnsr::I< DataType, SpatialDim, Frame > * > dt_shift, const Scalar< DataType > &lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_unit_normal, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::aa< DataType, SpatialDim, Frame > &pi)
	Computes time derivative of the shift vector from the generalized harmonic and geometric variables. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::I< DataType, SpatialDim, Frame >	time_deriv_of_shift (const Scalar< DataType > &lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::II< DataType, SpatialDim, Frame > &inverse_spatial_metric, const tnsr::A< DataType, SpatialDim, Frame > &spacetime_unit_normal, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::aa< DataType, SpatialDim, Frame > &pi)
	Computes time derivative of the shift vector from the generalized harmonic and geometric variables. More...
template<typename DataType , size_t SpatialDim, typename Frame >
void	time_deriv_of_spatial_metric (gsl::not_null< tnsr::ii< DataType, SpatialDim, Frame > * > dt_spatial_metric, const Scalar< DataType > &lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::aa< DataType, SpatialDim, Frame > &pi)
	Computes time derivative of the spatial metric. More...
template<typename DataType , size_t SpatialDim, typename Frame >
tnsr::ii< DataType, SpatialDim, Frame >	time_deriv_of_spatial_metric (const Scalar< DataType > &lapse, const tnsr::I< DataType, SpatialDim, Frame > &shift, const tnsr::iaa< DataType, SpatialDim, Frame > &phi, const tnsr::aa< DataType, SpatialDim, Frame > &pi)
	Computes time derivative of the spatial metric. More...

Detailed Description

Items related to evolving the first-order generalized harmonic system.

Function Documentation

characteristic_fields() [1/2]

template<size_t Dim, typename Frame >

Tags::CharacteristicFields< DataVector, Dim, Frame >::type gh::characteristic_fields	(	const Scalar< DataVector > &	gamma_2,
		const tnsr::II< DataVector, Dim, Frame > &	inverse_spatial_metric,
		const tnsr::aa< DataVector, Dim, Frame > &	spacetime_metric,
		const tnsr::aa< DataVector, Dim, Frame > &	pi,
		const tnsr::iaa< DataVector, Dim, Frame > &	phi,
		const tnsr::i< DataVector, Dim, Frame > &	unit_normal_one_form
	)

Computes characteristic fields from evolved fields.

CharacteristicFieldsCompute and EvolvedFieldsFromCharacteristicFieldsCompute convert between characteristic and evolved fields for the generalized harmonic system.

CharacteristicFieldsCompute computes characteristic fields as described in "A New Generalized Harmonic Evolution System" by Lindblom et. al [127] . Their names used here differ from this paper:

$$\begin{array}{rlr}\mathrm{S}\mathrm{p}\mathrm{E}\mathrm{C}\mathrm{T}\mathrm{R}\mathrm{E}& & \mathrm{L}\mathrm{i}\mathrm{n}\mathrm{d}\mathrm{b}\mathrm{l}\mathrm{o}\mathrm{m}\\ u_{ab}^{\psi }& & u_{ab}^{\hat{0}}\\ u_{iab}^0& & u_{iab}^{\hat{2}}\\ u_{ab}^{\pm }& & u_{ab}^{\hat{1}\pm }\end{array}$$

The characteristic fields $u$ are given in terms of the evolved fields by Eq.(32) - (34) of [127], respectively:

$$\begin{array}{rl}{u}_{ab}^{\psi }=& {\psi }_{ab}\\ u_{iab}^0=& ({\delta }_i^k-n_i{n}^k){\mathrm{\Phi }}_{kab}:=P_i^k{\mathrm{\Phi }}_{kab}\\ u_{ab}^{\pm }=& {\mathrm{\Pi }}_{ab}\pm n^i{\mathrm{\Phi }}_{iab}-{\gamma }_2{\psi }_{ab}\end{array}$$

where ${\psi }_{ab}$ is the spacetime metric, ${\mathrm{\Pi }}_{ab}$ and ${\mathrm{\Phi }}_{iab}$ are evolved generalized harmonic fields introduced by first derivatives of ${\psi }_{ab}$, ${\gamma }_2$ is a constraint damping parameter, and $n_k$ is the unit normal to the surface.

EvolvedFieldsFromCharacteristicFieldsCompute computes evolved fields $w$ in terms of the characteristic fields. This uses the inverse of above relations:

$$\begin{array}{rl}{\psi }_{ab}=& u_{ab}^{\psi },\\ {\mathrm{\Pi }}_{ab}=& \frac{1}{2}(u_{ab}^{+}+u_{ab}^{-})+{\gamma }_2{u}_{ab}^{\psi },\\ {\mathrm{\Phi }}_{iab}=& \frac{1}{2}(u_{ab}^{+}-u_{ab}^{-})n_i+u_{iab}^0.\end{array}$$

The corresponding characteristic speeds $v$ are computed by CharacteristicSpeedsCompute .

characteristic_fields() [2/2]

template<size_t Dim, typename Frame >

void gh::characteristic_fields	(	gsl::not_null< typename Tags::CharacteristicFields< DataVector, Dim, Frame >::type * >	char_fields,
		const Scalar< DataVector > &	gamma_2,
		const tnsr::II< DataVector, Dim, Frame > &	inverse_spatial_metric,
		const tnsr::aa< DataVector, Dim, Frame > &	spacetime_metric,
		const tnsr::aa< DataVector, Dim, Frame > &	pi,
		const tnsr::iaa< DataVector, Dim, Frame > &	phi,
		const tnsr::i< DataVector, Dim, Frame > &	unit_normal_one_form
	)

Computes characteristic fields from evolved fields.

CharacteristicFieldsCompute and EvolvedFieldsFromCharacteristicFieldsCompute convert between characteristic and evolved fields for the generalized harmonic system.

CharacteristicFieldsCompute computes characteristic fields as described in "A New Generalized Harmonic Evolution System" by Lindblom et. al [127] . Their names used here differ from this paper:

$$\begin{array}{rlr}\mathrm{S}\mathrm{p}\mathrm{E}\mathrm{C}\mathrm{T}\mathrm{R}\mathrm{E}& & \mathrm{L}\mathrm{i}\mathrm{n}\mathrm{d}\mathrm{b}\mathrm{l}\mathrm{o}\mathrm{m}\\ u_{ab}^{\psi }& & u_{ab}^{\hat{0}}\\ u_{iab}^0& & u_{iab}^{\hat{2}}\\ u_{ab}^{\pm }& & u_{ab}^{\hat{1}\pm }\end{array}$$

The characteristic fields $u$ are given in terms of the evolved fields by Eq.(32) - (34) of [127], respectively:

$$\begin{array}{rl}{u}_{ab}^{\psi }=& {\psi }_{ab}\\ u_{iab}^0=& ({\delta }_i^k-n_i{n}^k){\mathrm{\Phi }}_{kab}:=P_i^k{\mathrm{\Phi }}_{kab}\\ u_{ab}^{\pm }=& {\mathrm{\Pi }}_{ab}\pm n^i{\mathrm{\Phi }}_{iab}-{\gamma }_2{\psi }_{ab}\end{array}$$

where ${\psi }_{ab}$ is the spacetime metric, ${\mathrm{\Pi }}_{ab}$ and ${\mathrm{\Phi }}_{iab}$ are evolved generalized harmonic fields introduced by first derivatives of ${\psi }_{ab}$, ${\gamma }_2$ is a constraint damping parameter, and $n_k$ is the unit normal to the surface.

EvolvedFieldsFromCharacteristicFieldsCompute computes evolved fields $w$ in terms of the characteristic fields. This uses the inverse of above relations:

$$\begin{array}{rl}{\psi }_{ab}=& u_{ab}^{\psi },\\ {\mathrm{\Pi }}_{ab}=& \frac{1}{2}(u_{ab}^{+}+u_{ab}^{-})+{\gamma }_2{u}_{ab}^{\psi },\\ {\mathrm{\Phi }}_{iab}=& \frac{1}{2}(u_{ab}^{+}-u_{ab}^{-})n_i+u_{iab}^0.\end{array}$$

The corresponding characteristic speeds $v$ are computed by CharacteristicSpeedsCompute .

characteristic_speeds() [1/2]

template<size_t Dim, typename Frame >

std::array< DataVector, 4 > gh::characteristic_speeds	(	const Scalar< DataVector > &	gamma_1,
		const Scalar< DataVector > &	lapse,
		const tnsr::I< DataVector, Dim, Frame > &	shift,
		const tnsr::i< DataVector, Dim, Frame > &	unit_normal_one_form,
		const std::optional< tnsr::I< DataVector, Dim, Frame > > &	mesh_velocity
	)

Compute the characteristic speeds for the generalized harmonic system.

Computes the speeds as described in "A New Generalized Harmonic Evolution System" by Lindblom et. al [127] [see text following Eq.(34)]. The characteristic fields' names used here differ from this paper:

$$\begin{array}{rlr}\mathrm{S}\mathrm{p}\mathrm{E}\mathrm{C}\mathrm{T}\mathrm{R}\mathrm{E}& & \mathrm{L}\mathrm{i}\mathrm{n}\mathrm{d}\mathrm{b}\mathrm{l}\mathrm{o}\mathrm{m}\\ u_{ab}^{\psi }& & u_{ab}^{\hat{0}}\\ u_{iab}^0& & u_{iab}^{\hat{2}}\\ u_{ab}^{\pm }& & u_{ab}^{\hat{1}\pm }\end{array}$$

The corresponding characteristic speeds $v$ are given in the text between Eq.(34) and Eq.(35) of [127] :

$$\begin{array}{rl}{v}_{\psi }=& -(1+{\gamma }_1)n_k{\beta }^k-{\gamma }_1{n}_k{v}_g^k\\ v_0=& -n_k{\beta }^k\\ v_{\pm }=& -n_k{\beta }^k\pm \alpha \end{array}$$

where $\alpha ,{\beta }^k$ are the lapse and shift respectively, ${\gamma }_1$ is a constraint damping parameter, $n_k$ is the unit normal to the surface, and $v_g^k$ is the (optional) mesh velocity.

characteristic_speeds() [2/2]

template<size_t Dim, typename Frame >

void gh::characteristic_speeds	(	gsl::not_null< std::array< DataVector, 4 > * >	char_speeds,
		const Scalar< DataVector > &	gamma_1,
		const Scalar< DataVector > &	lapse,
		const tnsr::I< DataVector, Dim, Frame > &	shift,
		const tnsr::i< DataVector, Dim, Frame > &	unit_normal_one_form,
		const std::optional< tnsr::I< DataVector, Dim, Frame > > &	mesh_velocity
	)

Compute the characteristic speeds for the generalized harmonic system.

Computes the speeds as described in "A New Generalized Harmonic Evolution System" by Lindblom et. al [127] [see text following Eq.(34)]. The characteristic fields' names used here differ from this paper:

$$\begin{array}{rlr}\mathrm{S}\mathrm{p}\mathrm{E}\mathrm{C}\mathrm{T}\mathrm{R}\mathrm{E}& & \mathrm{L}\mathrm{i}\mathrm{n}\mathrm{d}\mathrm{b}\mathrm{l}\mathrm{o}\mathrm{m}\\ u_{ab}^{\psi }& & u_{ab}^{\hat{0}}\\ u_{iab}^0& & u_{iab}^{\hat{2}}\\ u_{ab}^{\pm }& & u_{ab}^{\hat{1}\pm }\end{array}$$

The corresponding characteristic speeds $v$ are given in the text between Eq.(34) and Eq.(35) of [127] :

$$\begin{array}{rl}{v}_{\psi }=& -(1+{\gamma }_1)n_k{\beta }^k-{\gamma }_1{n}_k{v}_g^k\\ v_0=& -n_k{\beta }^k\\ v_{\pm }=& -n_k{\beta }^k\pm \alpha \end{array}$$

where $\alpha ,{\beta }^k$ are the lapse and shift respectively, ${\gamma }_1$ is a constraint damping parameter, $n_k$ is the unit normal to the surface, and $v_g^k$ is the (optional) mesh velocity.

constraint_energy() [1/2]

template<typename DataType , size_t SpatialDim, typename Frame >

Scalar< DataType > gh::constraint_energy	(	const tnsr::a< DataType, SpatialDim, Frame > &	gauge_constraint,
		const tnsr::a< DataType, SpatialDim, Frame > &	f_constraint,
		const tnsr::ia< DataType, SpatialDim, Frame > &	two_index_constraint,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	three_index_constraint,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	four_index_constraint,
		const tnsr::II< DataType, SpatialDim, Frame > &	inverse_spatial_metric,
		const Scalar< DataType > &	spatial_metric_determinant,
		double	gauge_constraint_multiplier = 1.0,
		double	two_index_constraint_multiplier = 1.0,
		double	three_index_constraint_multiplier = 1.0,
		double	four_index_constraint_multiplier = 1.0
	)

Computes the generalized-harmonic (unnormalized) constraint energy.

Details

Computes the generalized-harmonic unnormalized constraint energy [Eq. (53) of [127] with $m^{ab}={\delta }^{ab}$ and with each term in the sum scaled by an arbitrary coefficient],

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ \text{(1)}& & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}+K_4{C}_{ikab}{C}_{jlab}{g}^{ij}{g}^{kl}.\end{array}$$

Here $C_a$ is the gauge constraint, $F_a$ is the f constraint, $C_{ia}$ is the two-index constraint, $C_{iab}$ is the three-index constraint, $C_{ikab}$ is the four-index constraint, $g^{ij}$ is the inverse spatial metric, and $K_1$, $K_2$, $K_3$, and $K_4$ are constant multipliers for each term that each default to a value of 1.0. Note that in this equation, spacetime indices $a,b$ are raised and lowered with the Kronecker delta.

Also note that the argument four_index_constraint is a rank-3 tensor. This is because gh::four_index_constraint() takes advantage of the antisymmetry of the four-index constraint's first two indices to only compute and return the independent components of $C_{ijab}$, which can be written as

$$\begin{array}{}\text{(2)}& D_{iab}\equiv \frac{1}{2}{ϵ}_i{}^{jk}{C}_{jkab},\end{array}$$

where ${ϵ}_{ijk}$ is the flat-space Levi-Civita tensor, whose inidces are raised and lowered with the Kronecker delta. The result is

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(3)}& & +& 2K_{4}g{D}_{iab}{D}_{jab}{g}^{ij},\end{array}$$

where $g$ is the determinant of the spatial metric.

To derive this expression for the constraint energy implemented here, Eq.~(53) of [127] is

$$\begin{array}{rcl}{S}_{AB}d{c}^{A}d{c}^B& =& m^{ab}[d{F}_{a}d{F}_b+g^{ij}(d{C}_{ia}d{C}_{jb}+g^{kl}{m}^{cd}d{C}_{ikac}d{C}_{jlbd})\\ \text{(4)}& & +& {\mathrm{\Lambda }}^2(d{C}_{a}d{C}_b+g^{ij}{m}^{cd}d{C}_{iac}d{C}_{jbd})].\end{array}$$

Replace $d{c}^A\to c^A$ to get

$$\begin{array}{rcl}E& =& m^{ab}[F_a{F}_b+g^{ij}(C_{ia}{C}_{jb}+g^{kl}{m}^{cd}{C}_{ikac}{C}_{jlbd})\\ & +& {\mathrm{\Lambda }}^2(C_a{C}_b+g^{ij}{m}^{cd}{C}_{iac}{C}_{jbd})]\\ & =& m^{ab}{F}_a{F}_b+m^{ab}{g}^{ij}{C}_{ia}{C}_{jb}+m^{ab}{g}^{ij}{g}^{kl}{m}^{cd}{C}_{ikac}{C}_{jlbd}\\ \text{(5)}& & +& m^{ab}{\mathrm{\Lambda }}^2{C}_a{C}_b+m^{ab}{\mathrm{\Lambda }}^2{g}^{ij}{m}^{cd}{C}_{iac}{C}_{jbd}.\end{array}$$

Here $m^{ab}$ is an arbitrary positive-definite matrix, and $\mathrm{\Lambda }$ is an arbitrary real scalar. Choose $m^{ab}={\delta }^{ab}$ but allow an arbitrary coefficient to be placed in front of each term. Then, absorb ${\mathrm{\Lambda }}^2$ into one of these coefficients, to get

$$\begin{array}{rcl}E& =& K_F{\delta }^{ab}{F}_a{F}_b+K_2{\delta }^{ab}{g}^{ij}{C}_{ia}{C}_{jb}+K_4{\delta }^{ab}{g}^{ij}{g}^{kl}{\delta }^{cd}{C}_{ikac}{C}_{jlbd}\\ \text{(6)}& & +& K_1{\delta }^{ab}{C}_a{C}_b+K_3{\delta }^{ab}{g}^{ij}{\delta }^{cd}{C}_{iac}{C}_{jbd}.\end{array}$$

Adopting a Euclidean norm for the constraint space (i.e., choosing to raise and lower spacetime indices with Kronecker deltas) gives

$$\begin{array}{rcl}E& =& K_F{F}_a{F}_a+K_2{g}^{ij}{C}_{ia}{C}_{ja}+K_4{g}^{ij}{g}^{kl}{C}_{ikac}{C}_{jlac}\\ \text{(7)}& & +& K_1{C}_a{C}_a+K_3{g}^{ij}{C}_{iac}{C}_{jac}.\end{array}$$

The two-index constraint and f constraint can be viewed as the time and space components of a combined spacetime constraint. So next choose $K_F=K_2$, giving

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ \text{(8)}& & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}+K_4{C}_{ikab}{C}_{jlab}{g}^{ij}{g}^{kl}.\end{array}$$

Note that $C_{ikab}$ is antisymmetric on the first two indices. Next, replace the four-index constraint $C_{ijab}$ with $D_{iab}$, which contains the independent components of $C_{ijab}$. Specifically,

$$\begin{array}{}\text{(9)}& D_{iab}\equiv \frac{1}{2}{ϵ}_i{}^{jk}{C}_{jkab}.\end{array}$$

The inverse relationship is

$$\begin{array}{}\text{(10)}& C_{jkab}={ϵ}^i{}_{jk}{D}_{iab},\end{array}$$

where ${ϵ}_{ijk}$ is the flat-space Levi-Civita tensor, whose indices are raised and lowered with the Kronecker delta. Inserting this relation to replace $C_{jkab}$ with $D_{iab}$ gives

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(11)}& & +& K_4{D}_{mab}{D}_{nab}{ϵ}^m{}_{ik}{ϵ}^n{}_{jl}{g}^{ij}{g}^{kl}.\end{array}$$

There's a subtle point here: $g^{ij}$ is the inverse spatial metric, which is not necessarily flat. But ${ϵ}^i{}_{jk}$ is the flat space Levi-Civita tensor. In order to raise and lower indices of the Levi-Civita tensor with the inverse spatial metrics, put in the appropriate factors of $\sqrt{g}$, where $g$ is the metric determinant, to make the curved-space Levi-Civita tensor compatible with $g^{ij}$. Let ${\epsilon }^{ijk}$ represent the curved space Levi-Civita tensor compatible with $g^{ij}$:

$$\begin{array}{}\text{(12)}& {\epsilon }^{mik}=g^{-1/2}{ϵ}^{mik}\text{(13)}& {\epsilon }_{mik}=g^{1/2}{ϵ}_{mik}.\end{array}$$

Then we can write the constraint energy as

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(14)}& & +& K_4{D}_{mab}{D}_{nab}g{g}^{-1/2}{ϵ}^m{}_{ik}{g}^{-1/2}{ϵ}^n{}_{jl}{g}^{ij}{g}^{kl}.\end{array}$$

The factors of $g^{-1/2}$ make the Levi-Civita tensor compatible with $g^{ij}$. Swapping which summed indices are raised and which are lowered gives

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(15)}& & +& K_4{D}_{mab}{D}_{nab}g{g}^{-1/2}{ϵ}^{mik}{g}^{-1/2}{ϵ}^{njl}{g}_{ij}{g}_{kl},\end{array}$$

or

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(16)}& & +& K_4{D}_{mab}{D}_{nab}g{\epsilon }^{mik}{\epsilon }^{njl}{g}_{ij}{g}_{kl},\end{array}$$

or, reversing up and down repeated indices again,

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(17)}& & +& K_4{D}_{mab}{D}_{nab}g{\epsilon }^m{}_{ik}{\epsilon }^n{}_{jl}{g}^{ij}{g}^{kl}.\end{array}$$

The metric raises and lowers the indices of ${\epsilon }^{ijk}$, so this can be written as

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(18)}& & +& K_{4}g{D}_{mab}{D}^n{}_{ab}{\epsilon }^{mjl}{\epsilon }_{njl}.\end{array}$$

Now, in flat space (Eq. (1.23) of [194]),

$$\begin{array}{}\text{(19)}& {ϵ}^{mjl}{ϵ}_{njl}={\delta }_{nj}^{mj}={\delta }_n^m{\delta }_j^j-{\delta }_j^m{\delta }_n^j=2{\delta }_n^m.\end{array}$$

But this holds for curved space as well: multiply the left hand side by $1=g^{1/2}{g}^{-1/2}$ to get

$$\begin{array}{}\text{(20)}& g^{-1/2}{ϵ}^{mjl}{g}^{1/2}{ϵ}_{njl}={\epsilon }^{mjl}{\epsilon }_{njl}={\delta }_{nj}^{mj}={\delta }_n^m{\delta }_j^j-{\delta }_j^m{\delta }_n^j=2{\delta }_n^m.\end{array}$$

So the constraint energy is

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(21)}& & +& 2K_4{D}_{mab}{D}^n{}_{ab}{\delta }_n^m.\end{array}$$

Simplifying gives the formula implemented here:

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(22)}& & +& 2K_{4}g{D}_{iab}{D}_{jab}{g}^{ij}.\end{array}$$

constraint_energy() [2/2]

template<typename DataType , size_t SpatialDim, typename Frame >

void gh::constraint_energy	(	gsl::not_null< Scalar< DataType > * >	energy,
		const tnsr::a< DataType, SpatialDim, Frame > &	gauge_constraint,
		const tnsr::a< DataType, SpatialDim, Frame > &	f_constraint,
		const tnsr::ia< DataType, SpatialDim, Frame > &	two_index_constraint,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	three_index_constraint,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	four_index_constraint,
		const tnsr::II< DataType, SpatialDim, Frame > &	inverse_spatial_metric,
		const Scalar< DataType > &	spatial_metric_determinant,
		double	gauge_constraint_multiplier = 1.0,
		double	two_index_constraint_multiplier = 1.0,
		double	three_index_constraint_multiplier = 1.0,
		double	four_index_constraint_multiplier = 1.0
	)

Computes the generalized-harmonic (unnormalized) constraint energy.

Details

Computes the generalized-harmonic unnormalized constraint energy [Eq. (53) of [127] with $m^{ab}={\delta }^{ab}$ and with each term in the sum scaled by an arbitrary coefficient],

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ \text{(23)}& & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}+K_4{C}_{ikab}{C}_{jlab}{g}^{ij}{g}^{kl}.\end{array}$$

Here $C_a$ is the gauge constraint, $F_a$ is the f constraint, $C_{ia}$ is the two-index constraint, $C_{iab}$ is the three-index constraint, $C_{ikab}$ is the four-index constraint, $g^{ij}$ is the inverse spatial metric, and $K_1$, $K_2$, $K_3$, and $K_4$ are constant multipliers for each term that each default to a value of 1.0. Note that in this equation, spacetime indices $a,b$ are raised and lowered with the Kronecker delta.

Also note that the argument four_index_constraint is a rank-3 tensor. This is because gh::four_index_constraint() takes advantage of the antisymmetry of the four-index constraint's first two indices to only compute and return the independent components of $C_{ijab}$, which can be written as

$$\begin{array}{}\text{(24)}& D_{iab}\equiv \frac{1}{2}{ϵ}_i{}^{jk}{C}_{jkab},\end{array}$$

where ${ϵ}_{ijk}$ is the flat-space Levi-Civita tensor, whose inidces are raised and lowered with the Kronecker delta. The result is

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(25)}& & +& 2K_{4}g{D}_{iab}{D}_{jab}{g}^{ij},\end{array}$$

where $g$ is the determinant of the spatial metric.

To derive this expression for the constraint energy implemented here, Eq.~(53) of [127] is

$$\begin{array}{rcl}{S}_{AB}d{c}^{A}d{c}^B& =& m^{ab}[d{F}_{a}d{F}_b+g^{ij}(d{C}_{ia}d{C}_{jb}+g^{kl}{m}^{cd}d{C}_{ikac}d{C}_{jlbd})\\ \text{(26)}& & +& {\mathrm{\Lambda }}^2(d{C}_{a}d{C}_b+g^{ij}{m}^{cd}d{C}_{iac}d{C}_{jbd})].\end{array}$$

Replace $d{c}^A\to c^A$ to get

$$\begin{array}{rcl}E& =& m^{ab}[F_a{F}_b+g^{ij}(C_{ia}{C}_{jb}+g^{kl}{m}^{cd}{C}_{ikac}{C}_{jlbd})\\ & +& {\mathrm{\Lambda }}^2(C_a{C}_b+g^{ij}{m}^{cd}{C}_{iac}{C}_{jbd})]\\ & =& m^{ab}{F}_a{F}_b+m^{ab}{g}^{ij}{C}_{ia}{C}_{jb}+m^{ab}{g}^{ij}{g}^{kl}{m}^{cd}{C}_{ikac}{C}_{jlbd}\\ \text{(27)}& & +& m^{ab}{\mathrm{\Lambda }}^2{C}_a{C}_b+m^{ab}{\mathrm{\Lambda }}^2{g}^{ij}{m}^{cd}{C}_{iac}{C}_{jbd}.\end{array}$$

Here $m^{ab}$ is an arbitrary positive-definite matrix, and $\mathrm{\Lambda }$ is an arbitrary real scalar. Choose $m^{ab}={\delta }^{ab}$ but allow an arbitrary coefficient to be placed in front of each term. Then, absorb ${\mathrm{\Lambda }}^2$ into one of these coefficients, to get

$$\begin{array}{rcl}E& =& K_F{\delta }^{ab}{F}_a{F}_b+K_2{\delta }^{ab}{g}^{ij}{C}_{ia}{C}_{jb}+K_4{\delta }^{ab}{g}^{ij}{g}^{kl}{\delta }^{cd}{C}_{ikac}{C}_{jlbd}\\ \text{(28)}& & +& K_1{\delta }^{ab}{C}_a{C}_b+K_3{\delta }^{ab}{g}^{ij}{\delta }^{cd}{C}_{iac}{C}_{jbd}.\end{array}$$

Adopting a Euclidean norm for the constraint space (i.e., choosing to raise and lower spacetime indices with Kronecker deltas) gives

$$\begin{array}{rcl}E& =& K_F{F}_a{F}_a+K_2{g}^{ij}{C}_{ia}{C}_{ja}+K_4{g}^{ij}{g}^{kl}{C}_{ikac}{C}_{jlac}\\ \text{(29)}& & +& K_1{C}_a{C}_a+K_3{g}^{ij}{C}_{iac}{C}_{jac}.\end{array}$$

The two-index constraint and f constraint can be viewed as the time and space components of a combined spacetime constraint. So next choose $K_F=K_2$, giving

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ \text{(30)}& & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}+K_4{C}_{ikab}{C}_{jlab}{g}^{ij}{g}^{kl}.\end{array}$$

Note that $C_{ikab}$ is antisymmetric on the first two indices. Next, replace the four-index constraint $C_{ijab}$ with $D_{iab}$, which contains the independent components of $C_{ijab}$. Specifically,

$$\begin{array}{}\text{(31)}& D_{iab}\equiv \frac{1}{2}{ϵ}_i{}^{jk}{C}_{jkab}.\end{array}$$

The inverse relationship is

$$\begin{array}{}\text{(32)}& C_{jkab}={ϵ}^i{}_{jk}{D}_{iab},\end{array}$$

where ${ϵ}_{ijk}$ is the flat-space Levi-Civita tensor, whose indices are raised and lowered with the Kronecker delta. Inserting this relation to replace $C_{jkab}$ with $D_{iab}$ gives

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(33)}& & +& K_4{D}_{mab}{D}_{nab}{ϵ}^m{}_{ik}{ϵ}^n{}_{jl}{g}^{ij}{g}^{kl}.\end{array}$$

There's a subtle point here: $g^{ij}$ is the inverse spatial metric, which is not necessarily flat. But ${ϵ}^i{}_{jk}$ is the flat space Levi-Civita tensor. In order to raise and lower indices of the Levi-Civita tensor with the inverse spatial metrics, put in the appropriate factors of $\sqrt{g}$, where $g$ is the metric determinant, to make the curved-space Levi-Civita tensor compatible with $g^{ij}$. Let ${\epsilon }^{ijk}$ represent the curved space Levi-Civita tensor compatible with $g^{ij}$:

$$\begin{array}{}\text{(34)}& {\epsilon }^{mik}=g^{-1/2}{ϵ}^{mik}\text{(35)}& {\epsilon }_{mik}=g^{1/2}{ϵ}_{mik}.\end{array}$$

Then we can write the constraint energy as

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(36)}& & +& K_4{D}_{mab}{D}_{nab}g{g}^{-1/2}{ϵ}^m{}_{ik}{g}^{-1/2}{ϵ}^n{}_{jl}{g}^{ij}{g}^{kl}.\end{array}$$

The factors of $g^{-1/2}$ make the Levi-Civita tensor compatible with $g^{ij}$. Swapping which summed indices are raised and which are lowered gives

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(37)}& & +& K_4{D}_{mab}{D}_{nab}g{g}^{-1/2}{ϵ}^{mik}{g}^{-1/2}{ϵ}^{njl}{g}_{ij}{g}_{kl},\end{array}$$

or

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(38)}& & +& K_4{D}_{mab}{D}_{nab}g{\epsilon }^{mik}{\epsilon }^{njl}{g}_{ij}{g}_{kl},\end{array}$$

or, reversing up and down repeated indices again,

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(39)}& & +& K_4{D}_{mab}{D}_{nab}g{\epsilon }^m{}_{ik}{\epsilon }^n{}_{jl}{g}^{ij}{g}^{kl}.\end{array}$$

The metric raises and lowers the indices of ${\epsilon }^{ijk}$, so this can be written as

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(40)}& & +& K_{4}g{D}_{mab}{D}^n{}_{ab}{\epsilon }^{mjl}{\epsilon }_{njl}.\end{array}$$

Now, in flat space (Eq. (1.23) of [194]),

$$\begin{array}{}\text{(41)}& {ϵ}^{mjl}{ϵ}_{njl}={\delta }_{nj}^{mj}={\delta }_n^m{\delta }_j^j-{\delta }_j^m{\delta }_n^j=2{\delta }_n^m.\end{array}$$

But this holds for curved space as well: multiply the left hand side by $1=g^{1/2}{g}^{-1/2}$ to get

$$\begin{array}{}\text{(42)}& g^{-1/2}{ϵ}^{mjl}{g}^{1/2}{ϵ}_{njl}={\epsilon }^{mjl}{\epsilon }_{njl}={\delta }_{nj}^{mj}={\delta }_n^m{\delta }_j^j-{\delta }_j^m{\delta }_n^j=2{\delta }_n^m.\end{array}$$

So the constraint energy is

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(43)}& & +& 2K_4{D}_{mab}{D}^n{}_{ab}{\delta }_n^m.\end{array}$$

Simplifying gives the formula implemented here:

$$\begin{array}{rcl}E& =& K_1{C}_a{C}_a+K_2(F_a{F}_a+C_{ia}{C}_{ja}{g}^{ij})\\ & +& K_3{C}_{iab}{C}_{jab}{g}^{ij}\\ \text{(44)}& & +& 2K_{4}g{D}_{iab}{D}_{jab}{g}^{ij}.\end{array}$$

constraint_energy_normalization() [1/2]

template<typename DataType , size_t SpatialDim, typename Frame >

Scalar< DataType > gh::constraint_energy_normalization	(	const tnsr::iaa< DataType, SpatialDim, Frame > &	d_spacetime_metric,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	d_pi,
		const tnsr::ijaa< DataType, SpatialDim, Frame > &	d_phi,
		const tnsr::II< DataType, SpatialDim, Frame > &	inverse_spatial_metric,
		const Scalar< DataType > &	sqrt_spatial_metric_determinant,
		double	dimensional_constant
	)

Computes the generalized-harmonic normalized constraint energy.

Details

Computes the generalized-harmonic normalized constraint energy integrand [Eq. (70) of [127] with $m^{ab}={\delta }^{ab}$],

$$\begin{array}{}\text{(45)}& ||E||={\gamma }^{ij}{\delta }^{ab}{\delta }^{cd}({\mathrm{\Lambda }}^2{\partial }_i{g}_{ac}{\partial }_j{g}_{bd}+{\partial }_i{\mathrm{\Pi }}_{ac}{\partial }_j{\mathrm{\Pi }}_{bd}+{\gamma }^{kl}{\partial }_i{\mathrm{\Phi }}_{kac}{\partial }_j{\mathrm{\Phi }}_{lbd}){\gamma }^{1/2}\end{array}$$

Here ${\gamma }^{ij}$ is the inverse spatial metric, and ${\mathrm{\Lambda }}^2$ is the squared dimensional constant. Note that in this equation, spacetime indices $a,b$ are raised and lowered with the Kronecker delta.

constraint_energy_normalization() [2/2]

template<typename DataType , size_t SpatialDim, typename Frame >

void gh::constraint_energy_normalization	(	gsl::not_null< Scalar< DataType > * >	energy_norm,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	d_spacetime_metric,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	d_pi,
		const tnsr::ijaa< DataType, SpatialDim, Frame > &	d_phi,
		const tnsr::II< DataType, SpatialDim, Frame > &	inverse_spatial_metric,
		const Scalar< DataType > &	sqrt_spatial_metric_determinant,
		double	dimensional_constant
	)

Computes the generalized-harmonic normalized constraint energy.

Details

Computes the generalized-harmonic normalized constraint energy integrand [Eq. (70) of [127] with $m^{ab}={\delta }^{ab}$],

$$\begin{array}{}\text{(46)}& ||E||={\gamma }^{ij}{\delta }^{ab}{\delta }^{cd}({\mathrm{\Lambda }}^2{\partial }_i{g}_{ac}{\partial }_j{g}_{bd}+{\partial }_i{\mathrm{\Pi }}_{ac}{\partial }_j{\mathrm{\Pi }}_{bd}+{\gamma }^{kl}{\partial }_i{\mathrm{\Phi }}_{kac}{\partial }_j{\mathrm{\Phi }}_{lbd}){\gamma }^{1/2}\end{array}$$

Here ${\gamma }^{ij}$ is the inverse spatial metric, and ${\mathrm{\Lambda }}^2$ is the squared dimensional constant. Note that in this equation, spacetime indices $a,b$ are raised and lowered with the Kronecker delta.

f_constraint() [1/4]

template<typename DataType , size_t SpatialDim, typename Frame >

tnsr::a< DataType, SpatialDim, Frame > gh::f_constraint	(	const tnsr::a< DataType, SpatialDim, Frame > &	gauge_function,
		const tnsr::ab< DataType, SpatialDim, Frame > &	spacetime_d_gauge_function,
		const tnsr::a< DataType, SpatialDim, Frame > &	spacetime_normal_one_form,
		const tnsr::A< DataType, SpatialDim, Frame > &	spacetime_normal_vector,
		const tnsr::II< DataType, SpatialDim, Frame > &	inverse_spatial_metric,
		const tnsr::AA< DataType, SpatialDim, Frame > &	inverse_spacetime_metric,
		const tnsr::aa< DataType, SpatialDim, Frame > &	pi,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	phi,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	d_pi,
		const tnsr::ijaa< DataType, SpatialDim, Frame > &	d_phi,
		const Scalar< DataType > &	gamma2,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	three_index_constraint
	)

Computes the generalized-harmonic F constraint.

Details

Computes the generalized-harmonic F constraint [Eq. (43) of [127]],

$$\begin{array}{rcl}{\mathcal{F}}_a& \equiv & \frac{1}{2}{g}_a^i{\psi }^{bc}{\partial }_i{\mathrm{\Pi }}_{bc}-g^{ij}{\partial }_i{\mathrm{\Pi }}_{ja}-g^{ij}{t}^b{\partial }_i{\mathrm{\Phi }}_{jba}+\frac{1}{2}{t}_a{\psi }^{bc}{g}^{ij}{\partial }_i{\mathrm{\Phi }}_{jbc}\\ & & +t_a{g}^{ij}{\partial }_i{H}_j+g_a^i{\mathrm{\Phi }}_{ijb}{g}^{jk}{\mathrm{\Phi }}_{kcd}{\psi }^{bd}{t}^c-\frac{1}{2}{g}_a^i{\mathrm{\Phi }}_{ijb}{g}^{jk}{\mathrm{\Phi }}_{kcd}{\psi }^{cd}{t}^b\\ & & -g_a^i{t}^b{\partial }_i{H}_b+g^{ij}{\mathrm{\Phi }}_{icd}{\mathrm{\Phi }}_{jba}{\psi }^{bc}{t}^d-\frac{1}{2}{t}_a{g}^{ij}{g}^{mn}{\mathrm{\Phi }}_{imc}{\mathrm{\Phi }}_{njd}{\psi }^{cd}\\ & & -\frac{1}{4}{t}_a{g}^{ij}{\mathrm{\Phi }}_{icd}{\mathrm{\Phi }}_{jbe}{\psi }^{cb}{\psi }^{de}+\frac{1}{4}{t}_a{\mathrm{\Pi }}_{cd}{\mathrm{\Pi }}_{be}{\psi }^{cb}{\psi }^{de}-g^{ij}{H}_i{\mathrm{\Pi }}_{ja}\\ & & -t^b{g}^{ij}{\mathrm{\Pi }}_{bi}{\mathrm{\Pi }}_{ja}-\frac{1}{4}{g}_a^i{\mathrm{\Phi }}_{icd}{t}^c{t}^d{\mathrm{\Pi }}_{be}{\psi }^{be}+\frac{1}{2}{t}_a{\mathrm{\Pi }}_{cd}{\mathrm{\Pi }}_{be}{\psi }^{ce}{t}^d{t}^b\\ & & +g_a^i{\mathrm{\Phi }}_{icd}{\mathrm{\Pi }}_{be}{t}^c{t}^b{\psi }^{de}-g^{ij}{\mathrm{\Phi }}_{iba}{t}^b{\mathrm{\Pi }}_{je}{t}^e-\frac{1}{2}{g}^{ij}{\mathrm{\Phi }}_{icd}{t}^c{t}^d{\mathrm{\Pi }}_{ja}\\ & & -g^{ij}{H}_i{\mathrm{\Phi }}_{jba}{t}^b+g_a^i{\mathrm{\Phi }}_{icd}{H}_b{\psi }^{bc}{t}^d+{\gamma }_2(g^{id}{\mathcal{C}}_{ida}-\frac{1}{2}{g}_a^i{\psi }^{cd}{\mathcal{C}}_{icd})\\ & & +\frac{1}{2}{t}_a{\mathrm{\Pi }}_{cd}{\psi }^{cd}{H}_b{t}^b-t_a{g}^{ij}{\mathrm{\Phi }}_{ijc}{H}_d{\psi }^{cd}+\frac{1}{2}{t}_a{g}^{ij}{H}_i{\mathrm{\Phi }}_{jcd}{\psi }^{cd}\\ \text{(47)}& & & -16\pi t^a{T}_{ab}\end{array}$$

where $H_a$ is the gauge function, ${\psi }_{ab}$ is the spacetime metric, ${\mathrm{\Pi }}_{ab}=-t^c{\partial }_c{\psi }_{ab}$, and ${\mathrm{\Phi }}_{iab}={\partial }_i{\psi }_{ab}$; $t^a$ is the timelike unit normal vector to the spatial slice, $g^{ij}$ is the inverse spatial metric, $g_c^b={\delta }_c^b+t^b{t}_c$, and $T_{ab}$ is the stress-energy tensor if nonzero (if using the overload with no stress-energy tensor provided, the stress energy term is omitted).

To justify the stress-energy contribution to the F constraint, note that the stress-energy tensor appears in the dynamics of the Generalized Harmonic system only through ${\partial }_t{\mathrm{\Pi }}_{ab}$. That dependence arises from (using $\dots$ to indicate collections of terms that are known to be independent of the stress-energy tensor):

$${\mathcal{F}}_a=\dots {\alpha }^{-1}({\partial }_t{\mathcal{C}}_a),$$

where

$${\mathcal{C}}_a=H_a+g^{ij}{\mathrm{\Phi }}_{ija}+t^b{\mathrm{\Pi }}_{ba}-\frac{1}{2}{g}_a{}^i{\psi }^{bc}{\mathrm{\Phi }}_{ibc}-\frac{1}{2}{t}_a{\psi }^{bc}{\mathrm{\Pi }}_{bc}.$$

.

Therefore, the Stress-energy contribution can be calculated from the trace-reversed contribution appearing in grmhd::GhValenciaDivClean::add_stress_energy_term_to_dt_pi – the trace reversal in that function and the trace-reversal that appears explicitly in ${\mathcal{C}}_a$ cancel.

f_constraint() [2/4]

template<typename DataType , size_t SpatialDim, typename Frame >

tnsr::a< DataType, SpatialDim, Frame > gh::f_constraint	(	const tnsr::a< DataType, SpatialDim, Frame > &	gauge_function,
		const tnsr::ab< DataType, SpatialDim, Frame > &	spacetime_d_gauge_function,
		const tnsr::a< DataType, SpatialDim, Frame > &	spacetime_normal_one_form,
		const tnsr::A< DataType, SpatialDim, Frame > &	spacetime_normal_vector,
		const tnsr::II< DataType, SpatialDim, Frame > &	inverse_spatial_metric,
		const tnsr::AA< DataType, SpatialDim, Frame > &	inverse_spacetime_metric,
		const tnsr::aa< DataType, SpatialDim, Frame > &	pi,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	phi,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	d_pi,
		const tnsr::ijaa< DataType, SpatialDim, Frame > &	d_phi,
		const Scalar< DataType > &	gamma2,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	three_index_constraint,
		const tnsr::aa< DataType, SpatialDim, Frame > &	trace_reversed_stress_energy
	)

Computes the generalized-harmonic F constraint.

Details

Computes the generalized-harmonic F constraint [Eq. (43) of [127]],

$$\begin{array}{rcl}{\mathcal{F}}_a& \equiv & \frac{1}{2}{g}_a^i{\psi }^{bc}{\partial }_i{\mathrm{\Pi }}_{bc}-g^{ij}{\partial }_i{\mathrm{\Pi }}_{ja}-g^{ij}{t}^b{\partial }_i{\mathrm{\Phi }}_{jba}+\frac{1}{2}{t}_a{\psi }^{bc}{g}^{ij}{\partial }_i{\mathrm{\Phi }}_{jbc}\\ & & +t_a{g}^{ij}{\partial }_i{H}_j+g_a^i{\mathrm{\Phi }}_{ijb}{g}^{jk}{\mathrm{\Phi }}_{kcd}{\psi }^{bd}{t}^c-\frac{1}{2}{g}_a^i{\mathrm{\Phi }}_{ijb}{g}^{jk}{\mathrm{\Phi }}_{kcd}{\psi }^{cd}{t}^b\\ & & -g_a^i{t}^b{\partial }_i{H}_b+g^{ij}{\mathrm{\Phi }}_{icd}{\mathrm{\Phi }}_{jba}{\psi }^{bc}{t}^d-\frac{1}{2}{t}_a{g}^{ij}{g}^{mn}{\mathrm{\Phi }}_{imc}{\mathrm{\Phi }}_{njd}{\psi }^{cd}\\ & & -\frac{1}{4}{t}_a{g}^{ij}{\mathrm{\Phi }}_{icd}{\mathrm{\Phi }}_{jbe}{\psi }^{cb}{\psi }^{de}+\frac{1}{4}{t}_a{\mathrm{\Pi }}_{cd}{\mathrm{\Pi }}_{be}{\psi }^{cb}{\psi }^{de}-g^{ij}{H}_i{\mathrm{\Pi }}_{ja}\\ & & -t^b{g}^{ij}{\mathrm{\Pi }}_{bi}{\mathrm{\Pi }}_{ja}-\frac{1}{4}{g}_a^i{\mathrm{\Phi }}_{icd}{t}^c{t}^d{\mathrm{\Pi }}_{be}{\psi }^{be}+\frac{1}{2}{t}_a{\mathrm{\Pi }}_{cd}{\mathrm{\Pi }}_{be}{\psi }^{ce}{t}^d{t}^b\\ & & +g_a^i{\mathrm{\Phi }}_{icd}{\mathrm{\Pi }}_{be}{t}^c{t}^b{\psi }^{de}-g^{ij}{\mathrm{\Phi }}_{iba}{t}^b{\mathrm{\Pi }}_{je}{t}^e-\frac{1}{2}{g}^{ij}{\mathrm{\Phi }}_{icd}{t}^c{t}^d{\mathrm{\Pi }}_{ja}\\ & & -g^{ij}{H}_i{\mathrm{\Phi }}_{jba}{t}^b+g_a^i{\mathrm{\Phi }}_{icd}{H}_b{\psi }^{bc}{t}^d+{\gamma }_2(g^{id}{\mathcal{C}}_{ida}-\frac{1}{2}{g}_a^i{\psi }^{cd}{\mathcal{C}}_{icd})\\ & & +\frac{1}{2}{t}_a{\mathrm{\Pi }}_{cd}{\psi }^{cd}{H}_b{t}^b-t_a{g}^{ij}{\mathrm{\Phi }}_{ijc}{H}_d{\psi }^{cd}+\frac{1}{2}{t}_a{g}^{ij}{H}_i{\mathrm{\Phi }}_{jcd}{\psi }^{cd}\\ \text{(48)}& & & -16\pi t^a{T}_{ab}\end{array}$$

where $H_a$ is the gauge function, ${\psi }_{ab}$ is the spacetime metric, ${\mathrm{\Pi }}_{ab}=-t^c{\partial }_c{\psi }_{ab}$, and ${\mathrm{\Phi }}_{iab}={\partial }_i{\psi }_{ab}$; $t^a$ is the timelike unit normal vector to the spatial slice, $g^{ij}$ is the inverse spatial metric, $g_c^b={\delta }_c^b+t^b{t}_c$, and $T_{ab}$ is the stress-energy tensor if nonzero (if using the overload with no stress-energy tensor provided, the stress energy term is omitted).

To justify the stress-energy contribution to the F constraint, note that the stress-energy tensor appears in the dynamics of the Generalized Harmonic system only through ${\partial }_t{\mathrm{\Pi }}_{ab}$. That dependence arises from (using $\dots$ to indicate collections of terms that are known to be independent of the stress-energy tensor):

$${\mathcal{F}}_a=\dots {\alpha }^{-1}({\partial }_t{\mathcal{C}}_a),$$

where

$${\mathcal{C}}_a=H_a+g^{ij}{\mathrm{\Phi }}_{ija}+t^b{\mathrm{\Pi }}_{ba}-\frac{1}{2}{g}_a{}^i{\psi }^{bc}{\mathrm{\Phi }}_{ibc}-\frac{1}{2}{t}_a{\psi }^{bc}{\mathrm{\Pi }}_{bc}.$$

.

Therefore, the Stress-energy contribution can be calculated from the trace-reversed contribution appearing in grmhd::GhValenciaDivClean::add_stress_energy_term_to_dt_pi – the trace reversal in that function and the trace-reversal that appears explicitly in ${\mathcal{C}}_a$ cancel.

f_constraint() [3/4]

template<typename DataType , size_t SpatialDim, typename Frame >

void gh::f_constraint	(	gsl::not_null< tnsr::a< DataType, SpatialDim, Frame > * >	constraint,
		const tnsr::a< DataType, SpatialDim, Frame > &	gauge_function,
		const tnsr::ab< DataType, SpatialDim, Frame > &	spacetime_d_gauge_function,
		const tnsr::a< DataType, SpatialDim, Frame > &	spacetime_normal_one_form,
		const tnsr::A< DataType, SpatialDim, Frame > &	spacetime_normal_vector,
		const tnsr::II< DataType, SpatialDim, Frame > &	inverse_spatial_metric,
		const tnsr::AA< DataType, SpatialDim, Frame > &	inverse_spacetime_metric,
		const tnsr::aa< DataType, SpatialDim, Frame > &	pi,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	phi,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	d_pi,
		const tnsr::ijaa< DataType, SpatialDim, Frame > &	d_phi,
		const Scalar< DataType > &	gamma2,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	three_index_constraint
	)

Computes the generalized-harmonic F constraint.

Details

Computes the generalized-harmonic F constraint [Eq. (43) of [127]],

$$\begin{array}{rcl}{\mathcal{F}}_a& \equiv & \frac{1}{2}{g}_a^i{\psi }^{bc}{\partial }_i{\mathrm{\Pi }}_{bc}-g^{ij}{\partial }_i{\mathrm{\Pi }}_{ja}-g^{ij}{t}^b{\partial }_i{\mathrm{\Phi }}_{jba}+\frac{1}{2}{t}_a{\psi }^{bc}{g}^{ij}{\partial }_i{\mathrm{\Phi }}_{jbc}\\ & & +t_a{g}^{ij}{\partial }_i{H}_j+g_a^i{\mathrm{\Phi }}_{ijb}{g}^{jk}{\mathrm{\Phi }}_{kcd}{\psi }^{bd}{t}^c-\frac{1}{2}{g}_a^i{\mathrm{\Phi }}_{ijb}{g}^{jk}{\mathrm{\Phi }}_{kcd}{\psi }^{cd}{t}^b\\ & & -g_a^i{t}^b{\partial }_i{H}_b+g^{ij}{\mathrm{\Phi }}_{icd}{\mathrm{\Phi }}_{jba}{\psi }^{bc}{t}^d-\frac{1}{2}{t}_a{g}^{ij}{g}^{mn}{\mathrm{\Phi }}_{imc}{\mathrm{\Phi }}_{njd}{\psi }^{cd}\\ & & -\frac{1}{4}{t}_a{g}^{ij}{\mathrm{\Phi }}_{icd}{\mathrm{\Phi }}_{jbe}{\psi }^{cb}{\psi }^{de}+\frac{1}{4}{t}_a{\mathrm{\Pi }}_{cd}{\mathrm{\Pi }}_{be}{\psi }^{cb}{\psi }^{de}-g^{ij}{H}_i{\mathrm{\Pi }}_{ja}\\ & & -t^b{g}^{ij}{\mathrm{\Pi }}_{bi}{\mathrm{\Pi }}_{ja}-\frac{1}{4}{g}_a^i{\mathrm{\Phi }}_{icd}{t}^c{t}^d{\mathrm{\Pi }}_{be}{\psi }^{be}+\frac{1}{2}{t}_a{\mathrm{\Pi }}_{cd}{\mathrm{\Pi }}_{be}{\psi }^{ce}{t}^d{t}^b\\ & & +g_a^i{\mathrm{\Phi }}_{icd}{\mathrm{\Pi }}_{be}{t}^c{t}^b{\psi }^{de}-g^{ij}{\mathrm{\Phi }}_{iba}{t}^b{\mathrm{\Pi }}_{je}{t}^e-\frac{1}{2}{g}^{ij}{\mathrm{\Phi }}_{icd}{t}^c{t}^d{\mathrm{\Pi }}_{ja}\\ & & -g^{ij}{H}_i{\mathrm{\Phi }}_{jba}{t}^b+g_a^i{\mathrm{\Phi }}_{icd}{H}_b{\psi }^{bc}{t}^d+{\gamma }_2(g^{id}{\mathcal{C}}_{ida}-\frac{1}{2}{g}_a^i{\psi }^{cd}{\mathcal{C}}_{icd})\\ & & +\frac{1}{2}{t}_a{\mathrm{\Pi }}_{cd}{\psi }^{cd}{H}_b{t}^b-t_a{g}^{ij}{\mathrm{\Phi }}_{ijc}{H}_d{\psi }^{cd}+\frac{1}{2}{t}_a{g}^{ij}{H}_i{\mathrm{\Phi }}_{jcd}{\psi }^{cd}\\ \text{(49)}& & & -16\pi t^a{T}_{ab}\end{array}$$

where $H_a$ is the gauge function, ${\psi }_{ab}$ is the spacetime metric, ${\mathrm{\Pi }}_{ab}=-t^c{\partial }_c{\psi }_{ab}$, and ${\mathrm{\Phi }}_{iab}={\partial }_i{\psi }_{ab}$; $t^a$ is the timelike unit normal vector to the spatial slice, $g^{ij}$ is the inverse spatial metric, $g_c^b={\delta }_c^b+t^b{t}_c$, and $T_{ab}$ is the stress-energy tensor if nonzero (if using the overload with no stress-energy tensor provided, the stress energy term is omitted).

To justify the stress-energy contribution to the F constraint, note that the stress-energy tensor appears in the dynamics of the Generalized Harmonic system only through ${\partial }_t{\mathrm{\Pi }}_{ab}$. That dependence arises from (using $\dots$ to indicate collections of terms that are known to be independent of the stress-energy tensor):

$${\mathcal{F}}_a=\dots {\alpha }^{-1}({\partial }_t{\mathcal{C}}_a),$$

where

$${\mathcal{C}}_a=H_a+g^{ij}{\mathrm{\Phi }}_{ija}+t^b{\mathrm{\Pi }}_{ba}-\frac{1}{2}{g}_a{}^i{\psi }^{bc}{\mathrm{\Phi }}_{ibc}-\frac{1}{2}{t}_a{\psi }^{bc}{\mathrm{\Pi }}_{bc}.$$

.

Therefore, the Stress-energy contribution can be calculated from the trace-reversed contribution appearing in grmhd::GhValenciaDivClean::add_stress_energy_term_to_dt_pi – the trace reversal in that function and the trace-reversal that appears explicitly in ${\mathcal{C}}_a$ cancel.

f_constraint() [4/4]

template<typename DataType , size_t SpatialDim, typename Frame >

void gh::f_constraint	(	gsl::not_null< tnsr::a< DataType, SpatialDim, Frame > * >	constraint,
		const tnsr::a< DataType, SpatialDim, Frame > &	gauge_function,
		const tnsr::ab< DataType, SpatialDim, Frame > &	spacetime_d_gauge_function,
		const tnsr::a< DataType, SpatialDim, Frame > &	spacetime_normal_one_form,
		const tnsr::A< DataType, SpatialDim, Frame > &	spacetime_normal_vector,
		const tnsr::II< DataType, SpatialDim, Frame > &	inverse_spatial_metric,
		const tnsr::AA< DataType, SpatialDim, Frame > &	inverse_spacetime_metric,
		const tnsr::aa< DataType, SpatialDim, Frame > &	pi,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	phi,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	d_pi,
		const tnsr::ijaa< DataType, SpatialDim, Frame > &	d_phi,
		const Scalar< DataType > &	gamma2,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	three_index_constraint,
		const tnsr::aa< DataType, SpatialDim, Frame > &	trace_reversed_stress_energy
	)

Computes the generalized-harmonic F constraint.

Details

Computes the generalized-harmonic F constraint [Eq. (43) of [127]],

$$\begin{array}{rcl}{\mathcal{F}}_a& \equiv & \frac{1}{2}{g}_a^i{\psi }^{bc}{\partial }_i{\mathrm{\Pi }}_{bc}-g^{ij}{\partial }_i{\mathrm{\Pi }}_{ja}-g^{ij}{t}^b{\partial }_i{\mathrm{\Phi }}_{jba}+\frac{1}{2}{t}_a{\psi }^{bc}{g}^{ij}{\partial }_i{\mathrm{\Phi }}_{jbc}\\ & & +t_a{g}^{ij}{\partial }_i{H}_j+g_a^i{\mathrm{\Phi }}_{ijb}{g}^{jk}{\mathrm{\Phi }}_{kcd}{\psi }^{bd}{t}^c-\frac{1}{2}{g}_a^i{\mathrm{\Phi }}_{ijb}{g}^{jk}{\mathrm{\Phi }}_{kcd}{\psi }^{cd}{t}^b\\ & & -g_a^i{t}^b{\partial }_i{H}_b+g^{ij}{\mathrm{\Phi }}_{icd}{\mathrm{\Phi }}_{jba}{\psi }^{bc}{t}^d-\frac{1}{2}{t}_a{g}^{ij}{g}^{mn}{\mathrm{\Phi }}_{imc}{\mathrm{\Phi }}_{njd}{\psi }^{cd}\\ & & -\frac{1}{4}{t}_a{g}^{ij}{\mathrm{\Phi }}_{icd}{\mathrm{\Phi }}_{jbe}{\psi }^{cb}{\psi }^{de}+\frac{1}{4}{t}_a{\mathrm{\Pi }}_{cd}{\mathrm{\Pi }}_{be}{\psi }^{cb}{\psi }^{de}-g^{ij}{H}_i{\mathrm{\Pi }}_{ja}\\ & & -t^b{g}^{ij}{\mathrm{\Pi }}_{bi}{\mathrm{\Pi }}_{ja}-\frac{1}{4}{g}_a^i{\mathrm{\Phi }}_{icd}{t}^c{t}^d{\mathrm{\Pi }}_{be}{\psi }^{be}+\frac{1}{2}{t}_a{\mathrm{\Pi }}_{cd}{\mathrm{\Pi }}_{be}{\psi }^{ce}{t}^d{t}^b\\ & & +g_a^i{\mathrm{\Phi }}_{icd}{\mathrm{\Pi }}_{be}{t}^c{t}^b{\psi }^{de}-g^{ij}{\mathrm{\Phi }}_{iba}{t}^b{\mathrm{\Pi }}_{je}{t}^e-\frac{1}{2}{g}^{ij}{\mathrm{\Phi }}_{icd}{t}^c{t}^d{\mathrm{\Pi }}_{ja}\\ & & -g^{ij}{H}_i{\mathrm{\Phi }}_{jba}{t}^b+g_a^i{\mathrm{\Phi }}_{icd}{H}_b{\psi }^{bc}{t}^d+{\gamma }_2(g^{id}{\mathcal{C}}_{ida}-\frac{1}{2}{g}_a^i{\psi }^{cd}{\mathcal{C}}_{icd})\\ & & +\frac{1}{2}{t}_a{\mathrm{\Pi }}_{cd}{\psi }^{cd}{H}_b{t}^b-t_a{g}^{ij}{\mathrm{\Phi }}_{ijc}{H}_d{\psi }^{cd}+\frac{1}{2}{t}_a{g}^{ij}{H}_i{\mathrm{\Phi }}_{jcd}{\psi }^{cd}\\ \text{(50)}& & & -16\pi t^a{T}_{ab}\end{array}$$

where $H_a$ is the gauge function, ${\psi }_{ab}$ is the spacetime metric, ${\mathrm{\Pi }}_{ab}=-t^c{\partial }_c{\psi }_{ab}$, and ${\mathrm{\Phi }}_{iab}={\partial }_i{\psi }_{ab}$; $t^a$ is the timelike unit normal vector to the spatial slice, $g^{ij}$ is the inverse spatial metric, $g_c^b={\delta }_c^b+t^b{t}_c$, and $T_{ab}$ is the stress-energy tensor if nonzero (if using the overload with no stress-energy tensor provided, the stress energy term is omitted).

To justify the stress-energy contribution to the F constraint, note that the stress-energy tensor appears in the dynamics of the Generalized Harmonic system only through ${\partial }_t{\mathrm{\Pi }}_{ab}$. That dependence arises from (using $\dots$ to indicate collections of terms that are known to be independent of the stress-energy tensor):

$${\mathcal{F}}_a=\dots {\alpha }^{-1}({\partial }_t{\mathcal{C}}_a),$$

where

$${\mathcal{C}}_a=H_a+g^{ij}{\mathrm{\Phi }}_{ija}+t^b{\mathrm{\Pi }}_{ba}-\frac{1}{2}{g}_a{}^i{\psi }^{bc}{\mathrm{\Phi }}_{ibc}-\frac{1}{2}{t}_a{\psi }^{bc}{\mathrm{\Pi }}_{bc}.$$

.

Therefore, the Stress-energy contribution can be calculated from the trace-reversed contribution appearing in grmhd::GhValenciaDivClean::add_stress_energy_term_to_dt_pi – the trace reversal in that function and the trace-reversal that appears explicitly in ${\mathcal{C}}_a$ cancel.

four_index_constraint() [1/2]

template<typename DataType , size_t SpatialDim, typename Frame >

tnsr::iaa< DataType, SpatialDim, Frame > gh::four_index_constraint

(

const tnsr::ijaa< DataType, SpatialDim, Frame > &

d_phi

)

Computes the generalized-harmonic 4-index constraint.

Details

Computes the independent components of the generalized-harmonic 4-index constraint. The constraint itself is given by Eq. (45) of [127],

$$\begin{array}{}\text{(51)}& C_{ijab}=2{\partial }_{[i}{\mathrm{\Phi }}_{j]ab},\end{array}$$

where ${\mathrm{\Phi }}_{iab}={\partial }_i{\psi }_{ab}$. Because the constraint is antisymmetric on the two spatial indices, here we compute and store only the independent components of $C_{ijab}$. Specifically, we compute

$$\begin{array}{}\text{(52)}& D_{iab}\equiv \frac{1}{2}{ϵ}_i{}^{jk}{C}_{jkab}={ϵ}_i{}^{jk}{\partial }_j{\mathrm{\Phi }}_{kab},\end{array}$$

where ${ϵ}_{ijk}$ is the flat-space Levi-Civita symbol, which is raised and lowered with the Kronecker delta. In terms of $D_{iab}$, the full 4-index constraint is

$$\begin{array}{}\text{(53)}& C_{jkab}={ϵ}^i{}_{jk}{D}_{iab}.\end{array}$$

four_index_constraint() [2/2]

template<typename DataType , size_t SpatialDim, typename Frame >

void gh::four_index_constraint	(	gsl::not_null< tnsr::iaa< DataType, SpatialDim, Frame > * >	constraint,
		const tnsr::ijaa< DataType, SpatialDim, Frame > &	d_phi
	)

Computes the generalized-harmonic 4-index constraint.

Details

Computes the independent components of the generalized-harmonic 4-index constraint. The constraint itself is given by Eq. (45) of [127],

$$\begin{array}{}\text{(54)}& C_{ijab}=2{\partial }_{[i}{\mathrm{\Phi }}_{j]ab},\end{array}$$

where ${\mathrm{\Phi }}_{iab}={\partial }_i{\psi }_{ab}$. Because the constraint is antisymmetric on the two spatial indices, here we compute and store only the independent components of $C_{ijab}$. Specifically, we compute

$$\begin{array}{}\text{(55)}& D_{iab}\equiv \frac{1}{2}{ϵ}_i{}^{jk}{C}_{jkab}={ϵ}_i{}^{jk}{\partial }_j{\mathrm{\Phi }}_{kab},\end{array}$$

where ${ϵ}_{ijk}$ is the flat-space Levi-Civita symbol, which is raised and lowered with the Kronecker delta. In terms of $D_{iab}$, the full 4-index constraint is

$$\begin{array}{}\text{(56)}& C_{jkab}={ϵ}^i{}_{jk}{D}_{iab}.\end{array}$$

gauge_constraint() [1/2]

template<typename DataType , size_t SpatialDim, typename Frame >

tnsr::a< DataType, SpatialDim, Frame > gh::gauge_constraint	(	const tnsr::a< DataType, SpatialDim, Frame > &	gauge_function,
		const tnsr::a< DataType, SpatialDim, Frame > &	spacetime_normal_one_form,
		const tnsr::A< DataType, SpatialDim, Frame > &	spacetime_normal_vector,
		const tnsr::II< DataType, SpatialDim, Frame > &	inverse_spatial_metric,
		const tnsr::AA< DataType, SpatialDim, Frame > &	inverse_spacetime_metric,
		const tnsr::aa< DataType, SpatialDim, Frame > &	pi,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	phi
	)

Computes the generalized-harmonic gauge constraint.

Details

Computes the generalized-harmonic gauge constraint [Eq. (40) of [127]],

$$C_a=H_a+g^{ij}{\mathrm{\Phi }}_{ija}+t^b{\mathrm{\Pi }}_{ba}-\frac{1}{2}{g}_a^i{\psi }^{bc}{\mathrm{\Phi }}_{ibc}-\frac{1}{2}{t}_a{\psi }^{bc}{\mathrm{\Pi }}_{bc},$$

where $H_a$ is the gauge function, ${\psi }_{ab}$ is the spacetime metric, ${\mathrm{\Pi }}_{ab}=-t^c{\partial }_c{\psi }_{ab}$, and ${\mathrm{\Phi }}_{iab}={\partial }_i{\psi }_{ab}$; $t^a$ is the timelike unit normal vector to the spatial slice, $g^{ij}$ is the inverse spatial metric, and $g_c^b={\delta }_c^b+t^b{t}_c$.

gauge_constraint() [2/2]

template<typename DataType , size_t SpatialDim, typename Frame >

void gh::gauge_constraint	(	gsl::not_null< tnsr::a< DataType, SpatialDim, Frame > * >	constraint,
		const tnsr::a< DataType, SpatialDim, Frame > &	gauge_function,
		const tnsr::a< DataType, SpatialDim, Frame > &	spacetime_normal_one_form,
		const tnsr::A< DataType, SpatialDim, Frame > &	spacetime_normal_vector,
		const tnsr::II< DataType, SpatialDim, Frame > &	inverse_spatial_metric,
		const tnsr::AA< DataType, SpatialDim, Frame > &	inverse_spacetime_metric,
		const tnsr::aa< DataType, SpatialDim, Frame > &	pi,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	phi
	)

Computes the generalized-harmonic gauge constraint.

Details

Computes the generalized-harmonic gauge constraint [Eq. (40) of [127]],

$$C_a=H_a+g^{ij}{\mathrm{\Phi }}_{ija}+t^b{\mathrm{\Pi }}_{ba}-\frac{1}{2}{g}_a^i{\psi }^{bc}{\mathrm{\Phi }}_{ibc}-\frac{1}{2}{t}_a{\psi }^{bc}{\mathrm{\Pi }}_{bc},$$

where $H_a$ is the gauge function, ${\psi }_{ab}$ is the spacetime metric, ${\mathrm{\Pi }}_{ab}=-t^c{\partial }_c{\psi }_{ab}$, and ${\mathrm{\Phi }}_{iab}={\partial }_i{\psi }_{ab}$; $t^a$ is the timelike unit normal vector to the spatial slice, $g^{ij}$ is the inverse spatial metric, and $g_c^b={\delta }_c^b+t^b{t}_c$.

initial_gh_variables_from_adm()

template<size_t Dim>

void gh::initial_gh_variables_from_adm	(	gsl::not_null< tnsr::aa< DataVector, Dim > * >	spacetime_metric,
		gsl::not_null< tnsr::aa< DataVector, Dim > * >	pi,
		gsl::not_null< tnsr::iaa< DataVector, Dim > * >	phi,
		const tnsr::ii< DataVector, Dim > &	spatial_metric,
		const Scalar< DataVector > &	lapse,
		const tnsr::I< DataVector, Dim > &	shift,
		const tnsr::ii< DataVector, Dim > &	extrinsic_curvature,
		const Mesh< Dim > &	mesh,
		const InverseJacobian< DataVector, Dim, Frame::ElementLogical, Frame::Inertial > &	inv_jacobian
	)

Compute initial GH variables from ADM variables.

• The spacetime metric is assembled from the spatial metric, lapse, and shift. See gr::spacetime_metric for details.
• Phi is set to the numerical derivative of the spacetime metric. This ensures that the 3-index constraint is initially satisfied.
• Pi is computed by choosing the time derivatives of lapse and shift to be zero. The gh::gauges::SetPiAndPhiFromConstraints mutator exists to override Pi later in the algorithm (it should be combined with this function).

three_index_constraint() [1/2]

template<typename DataType , size_t SpatialDim, typename Frame >

tnsr::iaa< DataType, SpatialDim, Frame > gh::three_index_constraint	(	const tnsr::iaa< DataType, SpatialDim, Frame > &	d_spacetime_metric,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	phi
	)

Computes the generalized-harmonic 3-index constraint.

Details

Computes the generalized-harmonic 3-index constraint, $C_{iab}={\partial }_i{\psi }_{ab}-{\mathrm{\Phi }}_{iab},$ which is given by Eq. (26) of [127]

three_index_constraint() [2/2]

template<typename DataType , size_t SpatialDim, typename Frame >

void gh::three_index_constraint	(	gsl::not_null< tnsr::iaa< DataType, SpatialDim, Frame > * >	constraint,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	d_spacetime_metric,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	phi
	)

Computes the generalized-harmonic 3-index constraint.

Details

Computes the generalized-harmonic 3-index constraint, $C_{iab}={\partial }_i{\psi }_{ab}-{\mathrm{\Phi }}_{iab},$ which is given by Eq. (26) of [127]

trace_christoffel() [1/2]

template<typename DataType , size_t SpatialDim, typename Frame >

tnsr::a< DataType, SpatialDim, Frame > gh::trace_christoffel	(	const tnsr::a< DataType, SpatialDim, Frame > &	spacetime_normal_one_form,
		const tnsr::A< DataType, SpatialDim, Frame > &	spacetime_normal_vector,
		const tnsr::II< DataType, SpatialDim, Frame > &	inverse_spatial_metric,
		const tnsr::AA< DataType, SpatialDim, Frame > &	inverse_spacetime_metric,
		const tnsr::aa< DataType, SpatialDim, Frame > &	pi,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	phi
	)

Compute ${\mathrm{\Gamma }}_a$ from the generalized harmonic evolved variables.

Starting from Eq. (40) of [127] we get

$$\begin{array}{rl}{\mathrm{\Gamma }}_a& ={\gamma }^{ij}{\mathrm{\Phi }}_{ija}+n^b{\mathrm{\Pi }}_{ab}-\frac{1}{2}{\gamma }^i{}_a{g}^{bc}{\mathrm{\Phi }}_{ibc}-\frac{1}{2}{n}_a{g}^{bc}{\mathrm{\Pi }}_{bc}\\ {\mathrm{\Gamma }}_a& ={\gamma }^{ij}{\mathrm{\Phi }}_{ija}+n^b{\mathrm{\Pi }}_{ab}-\frac{1}{2}{n}_a{g}^{bc}(n^i{\mathrm{\Phi }}_{ibc}+{\mathrm{\Pi }}_{bc})-\frac{1}{2}{\delta }_a^i{g}^{bc}{\mathrm{\Phi }}_{ibc}\end{array}$$

trace_christoffel() [2/2]

template<typename DataType , size_t SpatialDim, typename Frame >

void gh::trace_christoffel	(	gsl::not_null< tnsr::a< DataType, SpatialDim, Frame > * >	trace,
		const tnsr::a< DataType, SpatialDim, Frame > &	spacetime_normal_one_form,
		const tnsr::A< DataType, SpatialDim, Frame > &	spacetime_normal_vector,
		const tnsr::II< DataType, SpatialDim, Frame > &	inverse_spatial_metric,
		const tnsr::AA< DataType, SpatialDim, Frame > &	inverse_spacetime_metric,
		const tnsr::aa< DataType, SpatialDim, Frame > &	pi,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	phi
	)

Compute ${\mathrm{\Gamma }}_a$ from the generalized harmonic evolved variables.

Starting from Eq. (40) of [127] we get

$$\begin{array}{rl}{\mathrm{\Gamma }}_a& ={\gamma }^{ij}{\mathrm{\Phi }}_{ija}+n^b{\mathrm{\Pi }}_{ab}-\frac{1}{2}{\gamma }^i{}_a{g}^{bc}{\mathrm{\Phi }}_{ibc}-\frac{1}{2}{n}_a{g}^{bc}{\mathrm{\Pi }}_{bc}\\ {\mathrm{\Gamma }}_a& ={\gamma }^{ij}{\mathrm{\Phi }}_{ija}+n^b{\mathrm{\Pi }}_{ab}-\frac{1}{2}{n}_a{g}^{bc}(n^i{\mathrm{\Phi }}_{ibc}+{\mathrm{\Pi }}_{bc})-\frac{1}{2}{\delta }_a^i{g}^{bc}{\mathrm{\Phi }}_{ibc}\end{array}$$

two_index_constraint() [1/2]

template<typename DataType , size_t SpatialDim, typename Frame >

tnsr::ia< DataType, SpatialDim, Frame > gh::two_index_constraint	(	const tnsr::ab< DataType, SpatialDim, Frame > &	spacetime_d_gauge_function,
		const tnsr::a< DataType, SpatialDim, Frame > &	spacetime_normal_one_form,
		const tnsr::A< DataType, SpatialDim, Frame > &	spacetime_normal_vector,
		const tnsr::II< DataType, SpatialDim, Frame > &	inverse_spatial_metric,
		const tnsr::AA< DataType, SpatialDim, Frame > &	inverse_spacetime_metric,
		const tnsr::aa< DataType, SpatialDim, Frame > &	pi,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	phi,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	d_pi,
		const tnsr::ijaa< DataType, SpatialDim, Frame > &	d_phi,
		const Scalar< DataType > &	gamma2,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	three_index_constraint
	)

Computes the generalized-harmonic 2-index constraint.

Details

Computes the generalized-harmonic 2-index constraint [Eq. (44) of [127]],

$$\begin{array}{rcl}{C}_{ia}& \equiv & g^{jk}{\partial }_j{\mathrm{\Phi }}_{ika}-\frac{1}{2}{g}_a^j{\psi }^{cd}{\partial }_j{\mathrm{\Phi }}_{icd}+t^b{\partial }_i{\mathrm{\Pi }}_{ba}-\frac{1}{2}{t}_a{\psi }^{cd}{\partial }_i{\mathrm{\Pi }}_{cd}\\ & & +{\partial }_i{H}_a+\frac{1}{2}{g}_a^j{\mathrm{\Phi }}_{jcd}{\mathrm{\Phi }}_{ief}{\psi }^{ce}{\psi }^{df}+\frac{1}{2}{g}^{jk}{\mathrm{\Phi }}_{jcd}{\mathrm{\Phi }}_{ike}{\psi }^{cd}{t}^e{t}_a\\ & & -g^{jk}{g}^{mn}{\mathrm{\Phi }}_{jma}{\mathrm{\Phi }}_{ikn}+\frac{1}{2}{\mathrm{\Phi }}_{icd}{\mathrm{\Pi }}_{be}{t}_a({\psi }^{cb}{\psi }^{de}+\frac{1}{2}{\psi }^{be}{t}^c{t}^d)\\ \text{(57)}& & & -{\mathrm{\Phi }}_{icd}{\mathrm{\Pi }}_{ba}{t}^c({\psi }^{bd}+\frac{1}{2}{t}^b{t}^d)+\frac{1}{2}{\gamma }_2(t_a{\psi }^{cd}-2{\delta }_a^c{t}^d)C_{icd}.\end{array}$$

where $H_a$ is the gauge function, ${\psi }_{ab}$ is the spacetime metric, ${\mathrm{\Pi }}_{ab}=-t^c{\partial }_c{\psi }_{ab}$, and ${\mathrm{\Phi }}_{iab}={\partial }_i{\psi }_{ab}$; $t^a$ is the timelike unit normal vector to the spatial slice, $g^{ij}$ is the inverse spatial metric, and $g_c^b={\delta }_c^b+t^b{t}_c$.

two_index_constraint() [2/2]

template<typename DataType , size_t SpatialDim, typename Frame >

void gh::two_index_constraint	(	gsl::not_null< tnsr::ia< DataType, SpatialDim, Frame > * >	constraint,
		const tnsr::ab< DataType, SpatialDim, Frame > &	spacetime_d_gauge_function,
		const tnsr::a< DataType, SpatialDim, Frame > &	spacetime_normal_one_form,
		const tnsr::A< DataType, SpatialDim, Frame > &	spacetime_normal_vector,
		const tnsr::II< DataType, SpatialDim, Frame > &	inverse_spatial_metric,
		const tnsr::AA< DataType, SpatialDim, Frame > &	inverse_spacetime_metric,
		const tnsr::aa< DataType, SpatialDim, Frame > &	pi,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	phi,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	d_pi,
		const tnsr::ijaa< DataType, SpatialDim, Frame > &	d_phi,
		const Scalar< DataType > &	gamma2,
		const tnsr::iaa< DataType, SpatialDim, Frame > &	three_index_constraint
	)

Computes the generalized-harmonic 2-index constraint.

Details

Computes the generalized-harmonic 2-index constraint [Eq. (44) of [127]],

$$\begin{array}{rcl}{C}_{ia}& \equiv & g^{jk}{\partial }_j{\mathrm{\Phi }}_{ika}-\frac{1}{2}{g}_a^j{\psi }^{cd}{\partial }_j{\mathrm{\Phi }}_{icd}+t^b{\partial }_i{\mathrm{\Pi }}_{ba}-\frac{1}{2}{t}_a{\psi }^{cd}{\partial }_i{\mathrm{\Pi }}_{cd}\\ & & +{\partial }_i{H}_a+\frac{1}{2}{g}_a^j{\mathrm{\Phi }}_{jcd}{\mathrm{\Phi }}_{ief}{\psi }^{ce}{\psi }^{df}+\frac{1}{2}{g}^{jk}{\mathrm{\Phi }}_{jcd}{\mathrm{\Phi }}_{ike}{\psi }^{cd}{t}^e{t}_a\\ & & -g^{jk}{g}^{mn}{\mathrm{\Phi }}_{jma}{\mathrm{\Phi }}_{ikn}+\frac{1}{2}{\mathrm{\Phi }}_{icd}{\mathrm{\Pi }}_{be}{t}_a({\psi }^{cb}{\psi }^{de}+\frac{1}{2}{\psi }^{be}{t}^c{t}^d)\\ \text{(58)}& & & -{\mathrm{\Phi }}_{icd}{\mathrm{\Pi }}_{ba}{t}^c({\psi }^{bd}+\frac{1}{2}{t}^b{t}^d)+\frac{1}{2}{\gamma }_2(t_a{\psi }^{cd}-2{\delta }_a^c{t}^d)C_{icd}.\end{array}$$

where $H_a$ is the gauge function, ${\psi }_{ab}$ is the spacetime metric, ${\mathrm{\Pi }}_{ab}=-t^c{\partial }_c{\psi }_{ab}$, and ${\mathrm{\Phi }}_{iab}={\partial }_i{\psi }_{ab}$; $t^a$ is the timelike unit normal vector to the spatial slice, $g^{ij}$ is the inverse spatial metric, and $g_c^b={\delta }_c^b+t^b{t}_c$.