Namespaces
namespace	detail
namespace	bssn
namespace	fd
namespace	init

Classes
class	BSSNAtildeTensor
	Computes the trace-free conformal extrinsic curvature tensor \( \tilde{A}_{ij} \) in the BSSN formalism. More...
struct	ChiContext
class	BSSNChristoffel
	Compute the conformal Christoffel symbols \( \tilde{\Gamma}^k_{ij} \). More...
class	ConstraintSolver
class	BSSNContractedGamma
class	ExtrinsicCurvature
	Computes the extrinsic curvature tensor \( K_{ij} \) from BSSN variables. More...
class	RicciPhysicalTensor
	Computes the physical 3-Ricci tensor \( R_{ij} \) as the sum of conformal and conformal-factor contributions. More...
class	RicciConformalTensor
	Provides methods to compute the \( \chi \)-dependent part of the Ricci tensor in the BSSN formalism. More...
class	RicciTildeTensor
	Provides methods to compute the conformal Ricci tensor \( \tilde{R}_{ij} \) in the BSSN formulation. More...
struct	BSSNGrid
	Storage structure for all evolved BSSN variables on a single grid point or patch. More...
class	BSSN
	Driver class to initialize and store BSSN variables from an input spacetime metric. More...
class	TildeGamma
	Computes the contracted conformal Christoffel vector \( \tilde{\Gamma}^i \). More...
class	ChristoffelSym
	Stores and computes Christoffel symbols \( \Gamma^\lambda_{\mu\nu} \). More...
class	Metric
	A callable 4D metric class for general relativity (Minkowski, Schwarzschild, Kerr, etc.). More...
class	RicciTensor
	Computes the Ricci tensor and Ricci scalar from a 4D Riemann tensor. More...
class	RiemannTensor
	Computes the 4D Riemann curvature tensor \( R^\rho_{\ \sigma\mu\nu} \). More...
struct	AlignedDeleter
struct	GridDims
struct	Strides
struct	Field3D
struct	BoundaryPeriodic
struct	BoundaryClamp

Typedefs
template<typename T>
using	aligned_unique_ptr = std::unique_ptr<T[], AlignedDeleter<T>>

Enumerations
enum class	DiffMode { PARTIAL , COV , COV2 , SPEC }
enum	Sym6 { XX = 0 , XY = 1 , XZ = 2 , YY = 3 , YZ = 4 , ZZ = 5 }

Functions
template<typename T, typename ScalarFunc>
tensorium::Vector< T >	covariant_scalar (const tensorium::Vector< T > &X, T dx, T dy, T dz, ScalarFunc &&func, const tensorium::Tensor< T, 3 > &christoffel)
template<typename T, typename ScalarFunc>
tensorium::Tensor< T, 2 >	autodiff_scalar_second (const tensorium::Vector< T > &X, T dx, T dy, T dz, ScalarFunc &&func, const tensorium::Tensor< T, 3 > &christoffel)
template<typename T, typename ScalarFunc>
tensorium::Vector< T >	autodiff_rank0 (const tensorium::Vector< T > &X, T dx, T dy, T dz, ScalarFunc &&func, DiffMode mode, const tensorium::Tensor< T, 3 > &christoffel)
template<typename T, typename VectorFunc>
tensorium::Tensor< T, 2 >	autodiff_rank1 (const tensorium::Vector< T > &X, T dx, T dy, T dz, VectorFunc &&func, DiffMode mode, const tensorium::Tensor< T, 3 > &christoffel)
template<typename T, typename TensorFunc>
tensorium::Tensor< T, 3 >	autodiff_rank2_first (const tensorium::Vector< T > &X, T dx, T dy, T dz, TensorFunc &&func, DiffMode mode, const tensorium::Tensor< T, 3 > &christoffel)
template<typename T, typename FieldFunc>
auto	autodiff (const tensorium::Vector< T > &X, T dx, T dy, T dz, FieldFunc &&func, DiffMode mode, const tensorium::Tensor< T, 3 > &christoffel={})
template<typename T>
void	compute_partial_derivatives_3D (const tensorium::Tensor< T, 5 > &gamma_field, size_t i, size_t j, size_t k, T dx, T dy, T dz, tensorium::Tensor< T, 3 > &dgamma_out)
	Compute the 3D partial derivatives \( \partial_k \tilde{\gamma}_{ij} \) from a 5D field tensor.
template<typename T>
void	spectral_derivative_1D (tensorium::Vector< T > &field, tensorium::Vector< T > &dfield, T dx)
template<typename T>
void	spectral_partial_scalar_3D (const tensorium::Tensor< T, 3 > &scalar_field, T dx, T dy, T dz, tensorium::Tensor< T, 4 > &grad_out)
template<typename T, typename TensorFunc>
tensorium::Tensor< T, 3 >	spectral_partial_tensor2 (const tensorium::Vector< T > &X, T dx, T dy, T dz, TensorFunc &&func, size_t NX, size_t NY, size_t NZ)
template<typename T, typename ScalarFunc>
tensorium::Vector< T >	partial_scalar (const tensorium::Vector< T > &X, T dx, T dy, T dz, ScalarFunc &&func)
template<typename T, typename VectorFunc>
tensorium::Tensor< T, 2 >	partial_vector (const tensorium::Vector< T > &X, T dx, T dy, T dz, VectorFunc &&func)
template<typename T, typename TensorFunc>
void	compute_partial_derivatives_tensor2D (const tensorium::Vector< T > &X, T dx, T dy, T dz, TensorFunc &&func, tensorium::Tensor< T, 3 > &out)
template<typename T, typename TensorFunc>
void	compute_second_derivatives_tensor2D (const tensorium::Vector< T > &X, T dx, T dy, T dz, TensorFunc &&func, tensorium::Tensor< T, 4 > &out)
template<typename T, typename VectorFunc>
void	compute_partial_derivatives_vector (const tensorium::Vector< T > &X, T dx, T dy, T dz, VectorFunc &&func, tensorium::Tensor< T, 2 > &out)
template<typename T, typename ScalarFunc>
void	compute_partial_derivatives_scalar (const tensorium::Vector< T > &X, T dx, T dy, T dz, ScalarFunc &&func, tensorium::Vector< T > &out)
template<typename T>
tensorium::Tensor< T, 2 >	compute_dt_gamma_from_beta (const tensorium::Vector< T > &beta_cov, const tensorium::Tensor< T, 2 > &partial_beta_cov, const tensorium::Tensor< T, 3 > &christoffel)
template<typename T>
void	compute_partial_derivatives_vector3D (const tensorium::Tensor< T, 4 > &vec_field, size_t i, size_t j, size_t k, T dx, T dy, T dz, tensorium::Tensor< T, 2 > &dvec_out)
template<typename T, typename ScalarFunc>
void	compute_second_derivatives_scalar (const tensorium::Vector< T > &X, T dx, T dy, T dz, ScalarFunc &&func, tensorium::Tensor< T, 2 > &out)
template<typename T, typename ScalarFunc>
tensorium::Tensor< T, 2 >	covariant_scalar_second (const tensorium::Vector< T > &X, T dx, T dy, T dz, ScalarFunc &&func, const tensorium::Tensor< T, 3 > &christoffel)
template<typename T>
void	compute_second_derivatives_scalar3D (const tensorium::Tensor< T, 3 > &scalar_field, T dx, T dy, T dz, tensorium::Tensor< T, 5 > &hessian_out)
template<typename T>
tensorium::Tensor< T, 5 >	covariant_scalar_second_3D (const tensorium::Tensor< T, 3 > &chi, const tensorium::Tensor< T, 5 > &christoffel, T dx, T dy, T dz)
template<typename T, typename VectorFunc>
tensorium::Tensor< T, 2 >	covariant_vector (const tensorium::Vector< T > &X, T dx, T dy, T dz, VectorFunc &&func, const tensorium::Tensor< T, 3 > &Gamma)
template<typename T, typename TensorFunc>
tensorium::Tensor< T, 3 >	covariant_tensor2 (const tensorium::Vector< T > &X, T dx, T dy, T dz, TensorFunc &&func, const tensorium::Tensor< T, 3 > &Gamma)
template<typename T, typename TensorFunc>
tensorium::Tensor< T, 3 >	partial_tensor2 (const tensorium::Vector< T > &X, T dx, T dy, T dz, TensorFunc &&func)
template<typename T, typename TensorFunc>
tensorium::Tensor< T, 4 >	covariant_tensor2_second (const tensorium::Vector< T > &X, T dx, T dy, T dz, TensorFunc &&func, const tensorium::Tensor< T, 3 > &Gamma)
template<typename T, typename TensorFunc>
tensorium::Tensor< T, 3 >	populate_tensor3D_component (size_t i, size_t j, TensorFunc &&func, T dx, T dy, T dz, size_t NX, size_t NY, size_t NZ)
template<typename T>
tensorium::Tensor< T, 5 >	generate_conformal_metric_field (const tensorium_RG::Metric< T > &metric, size_t Nx, size_t Ny, size_t Nz, T dx, T dy, T dz)
	Generate the conformal 3-metric field \( \tilde{\gamma}_{ij}(x^i) \) on a 3D grid.
template<typename T>
void	print_tensor2 (const std::string &name, const tensorium::Tensor< T, 2 > &tensor)
template<typename T>
void	print_tensor3 (const std::string &name, const tensorium::Tensor< T, 3 > &tensor)
template<typename T>
void	print_vector (const std::string &name, const tensorium::Vector< T > &vec)
template<typename T>
tensorium::Tensor< T, 2 >	compute_dt_gamma_from_metric (const tensorium::Vector< T > &X, const tensorium_RG::Metric< T > &metric, T dt)
template<typename T>
	__attribute__ ((always_inline, hot, flatten)) inline tensorium
	Compute the inverse of a 2D metric tensor using local matrix inversion.
template<typename T>
constexpr size_t	simd_width_elems () noexcept
template<typename T>
size_t	pad_simd (size_t n) noexcept
template<typename T>
aligned_unique_ptr< T >	aligned_alloc_n (size_t n)
template<typename T>
void	halo_periodic (Field3D< T > &f, const GridDims &D)
int	sym6_index (int i, int j) noexcept
template<typename T>
void	store_sym6 (Field3D< T > *f6, size_t idx, T xx, T xy, T xz, T yy, T yz, T zz)
template<typename T>
void	load_sym6 (Field3D< T > *f6, size_t idx, T &xx, T &xy, T &xz, T &yy, T &yz, T &zz)
template<typename T>
Field3D< T >	make_field (const Strides< T > &st)
template<typename T>
void	halo_linear_extrapolate_full (Field3D< T > &f, const GridDims &D)
template<typename T>
T	sym6_get (const Field3D< T > *f6, size_t idx, int i, int j) noexcept
template<typename T>
T	sym6_inv_get (const Field3D< T > *f6, size_t idx, int i, int j) noexcept

Typedef Documentation

◆ aligned_unique_ptr

template<typename T>

using tensorium_RG::aligned_unique_ptr = std::unique_ptr<T[], AlignedDeleter<T>>

Enumeration Type Documentation

◆ DiffMode

enum class tensorium_RG::DiffMode

strong

Enumerator
PARTIAL
COV
COV2
SPEC

◆ Sym6

enum tensorium_RG::Sym6

Enumerator
XX
XY
XZ
YY
YZ
ZZ

Function Documentation

◆ attribute()

template<typename T>

tensorium_RG::__attribute__ ( (always_inline, hot, flatten) )

Compute the inverse of a 2D metric tensor using local matrix inversion.

Compute Christoffel symbols at an offset position along one direction.

Converts the tensor to a matrix, computes the inverse, and converts back.

Template Parameters

T	Scalar type

Parameters

g	Metric tensor \( g_{\mu\nu} \)

Returns: Inverse tensor \( g^{\mu\nu} \)

Used in numerical relativity to evaluate derivatives of Christoffel symbols.

Template Parameters

T	Scalar type

Parameters

X	Base coordinate
direction	Axis of offset
offset	Value to add (positive or negative)
h	Finite difference step size
g	Output metric tensor
g_inv	Output inverse metric
Gamma_out	Output Christoffel tensor
metric	Metric object

References tensorium::Tensor< K, Rank >::dimensions, and tensorium::Matrix< K, RowMajor >::inverse().

Here is the call graph for this function:

◆ aligned_alloc_n()

template<typename T>

aligned_unique_ptr< T > tensorium_RG::aligned_alloc_n ( size_t n )

inline

References TENSORIUM_ALIGN.

Referenced by tensorium_RG::bssn::detail::allocate_like(), and make_field().

Here is the caller graph for this function:

◆ autodiff()

template<typename T, typename FieldFunc>

auto tensorium_RG::autodiff	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		FieldFunc &&	func,
		DiffMode	mode,
		const tensorium::Tensor< T, 3 > &	christoffel = {} )

References X().

Referenced by tensorium_RG::BSSN< T >::init_BSSN().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ autodiff_rank0()

template<typename T, typename ScalarFunc>

tensorium::Vector< T > tensorium_RG::autodiff_rank0	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		ScalarFunc &&	func,
		DiffMode	mode,
		const tensorium::Tensor< T, 3 > &	christoffel )

References COV, covariant_scalar(), PARTIAL, partial_scalar(), and X().

Here is the call graph for this function:

◆ autodiff_rank1()

template<typename T, typename VectorFunc>

tensorium::Tensor< T, 2 > tensorium_RG::autodiff_rank1	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		VectorFunc &&	func,
		DiffMode	mode,
		const tensorium::Tensor< T, 3 > &	christoffel )

References COV, covariant_vector(), PARTIAL, partial_vector(), and X().

Here is the call graph for this function:

◆ autodiff_rank2_first()

template<typename T, typename TensorFunc>

tensorium::Tensor< T, 3 > tensorium_RG::autodiff_rank2_first	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		TensorFunc &&	func,
		DiffMode	mode,
		const tensorium::Tensor< T, 3 > &	christoffel )

References COV, covariant_tensor2(), PARTIAL, partial_tensor2(), SPEC, spectral_partial_tensor2(), and X().

Here is the call graph for this function:

◆ autodiff_scalar_second()

template<typename T, typename ScalarFunc>

tensorium::Tensor< T, 2 > tensorium_RG::autodiff_scalar_second	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		ScalarFunc &&	func,
		const tensorium::Tensor< T, 3 > &	christoffel )

References covariant_scalar_second(), and X().

Here is the call graph for this function:

◆ compute_dt_gamma_from_beta()

template<typename T>

tensorium::Tensor< T, 2 > tensorium_RG::compute_dt_gamma_from_beta	(	const tensorium::Vector< T > &	beta_cov,
		const tensorium::Tensor< T, 2 > &	partial_beta_cov,
		const tensorium::Tensor< T, 3 > &	christoffel )

Referenced by tensorium_RG::BSSN< T >::init_BSSN().

Here is the caller graph for this function:

◆ compute_dt_gamma_from_metric()

template<typename T>

tensorium::Tensor< T, 2 > tensorium_RG::compute_dt_gamma_from_metric	(	const tensorium::Vector< T > &	X,
		const tensorium_RG::Metric< T > &	metric,
		T	dt )

References tensorium_RG::Metric< T >::BSSN(), and X().

Referenced by tensorium_RG::BSSN< T >::init_BSSN().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ compute_partial_derivatives_3D()

template<typename T>

void tensorium_RG::compute_partial_derivatives_3D	(	const tensorium::Tensor< T, 5 > &	gamma_field,
		size_t	i,
		size_t	j,
		size_t	k,
		T	dx,
		T	dy,
		T	dz,
		tensorium::Tensor< T, 3 > &	dgamma_out )

inline

Compute the 3D partial derivatives \( \partial_k \tilde{\gamma}_{ij} \) from a 5D field tensor.

Uses centered 4th-order or 2nd-order finite differences for inner and boundary points.

\[\partial_k \tilde{\gamma}_{ij} = \begin{cases} \frac{-f(x+2h) + 8f(x+h) - 8f(x-h) + f(x-2h)}{12h} & \text{if inner} \\ \frac{f(x+h) - f(x-h)}{2h} & \text{if near boundary} \\ 0 & \text{otherwise} \end{cases} \]

Parameters

gamma_field	Field of type Tensor<T, 5> with dimensions [Nx, Ny, Nz, 3, 3]
i,j,k	Grid position
dx,dy,dz	Grid spacings
dgamma_out	Output tensor \( \partial_k \tilde{\gamma}_{ij} \)

References tensorium::Tensor< K, Rank >::resize(), and tensorium::Tensor< K, Rank >::shape().

Here is the call graph for this function:

◆ compute_partial_derivatives_scalar()

template<typename T, typename ScalarFunc>

void tensorium_RG::compute_partial_derivatives_scalar	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		ScalarFunc &&	func,
		tensorium::Vector< T > &	out )

References tensorium::Vector< K >::resize(), tensorium_RG::detail::shifted_copy(), tensorium_RG::detail::spatial_indices(), and X().

Referenced by tensorium_RG::BSSNContractedGamma< T >::compute(), covariant_scalar_second(), and covariant_scalar_second_3D().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ compute_partial_derivatives_tensor2D()

template<typename T, typename TensorFunc>

void tensorium_RG::compute_partial_derivatives_tensor2D	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		TensorFunc &&	func,
		tensorium::Tensor< T, 3 > &	out )

References tensorium::Tensor< K, Rank >::resize(), tensorium_RG::detail::shifted_copy(), tensorium_RG::detail::spatial_indices(), and X().

Here is the call graph for this function:

◆ compute_partial_derivatives_vector()

template<typename T, typename VectorFunc>

void tensorium_RG::compute_partial_derivatives_vector	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		VectorFunc &&	func,
		tensorium::Tensor< T, 2 > &	out )

References tensorium::Tensor< K, Rank >::resize(), tensorium_RG::detail::shifted_copy(), tensorium_RG::detail::spatial_indices(), and X().

Here is the call graph for this function:

◆ compute_partial_derivatives_vector3D()

template<typename T>

void tensorium_RG::compute_partial_derivatives_vector3D	(	const tensorium::Tensor< T, 4 > &	vec_field,
		size_t	i,
		size_t	j,
		size_t	k,
		T	dx,
		T	dy,
		T	dz,
		tensorium::Tensor< T, 2 > &	dvec_out )

inline

References tensorium::Tensor< K, Rank >::resize(), and tensorium::Tensor< K, Rank >::shape().

Here is the call graph for this function:

◆ compute_second_derivatives_scalar()

template<typename T, typename ScalarFunc>

void tensorium_RG::compute_second_derivatives_scalar	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		ScalarFunc &&	func,
		tensorium::Tensor< T, 2 > &	out )

inline

References tensorium::Tensor< K, Rank >::resize(), tensorium_RG::detail::shifted_copy(), tensorium_RG::detail::spatial_indices(), and X().

Referenced by tensorium_RG::ChiContext< T >::compute(), tensorium_RG::RicciTildeTensor< T >::compute_laplacian_term(), covariant_scalar_second(), and covariant_scalar_second_3D().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ compute_second_derivatives_scalar3D()

template<typename T>

void tensorium_RG::compute_second_derivatives_scalar3D	(	const tensorium::Tensor< T, 3 > &	scalar_field,
		T	dx,
		T	dy,
		T	dz,
		tensorium::Tensor< T, 5 > &	hessian_out )

References tensorium::Tensor< K, Rank >::resize(), and tensorium::Tensor< K, Rank >::shape().

Here is the call graph for this function:

◆ compute_second_derivatives_tensor2D()

template<typename T, typename TensorFunc>

void tensorium_RG::compute_second_derivatives_tensor2D	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		TensorFunc &&	func,
		tensorium::Tensor< T, 4 > &	out )

References tensorium_RG::detail::shifted_copy(), tensorium_RG::detail::spatial_indices(), and X().

Here is the call graph for this function:

◆ covariant_scalar()

template<typename T, typename ScalarFunc>

tensorium::Vector< T > tensorium_RG::covariant_scalar	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		ScalarFunc &&	func,
		const tensorium::Tensor< T, 3 > &	christoffel )

References partial_scalar(), and X().

Referenced by autodiff_rank0().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ covariant_scalar_second()

template<typename T, typename ScalarFunc>

tensorium::Tensor< T, 2 > tensorium_RG::covariant_scalar_second	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		ScalarFunc &&	func,
		const tensorium::Tensor< T, 3 > &	christoffel )

References compute_partial_derivatives_scalar(), compute_second_derivatives_scalar(), and X().

Referenced by autodiff_scalar_second().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ covariant_scalar_second_3D()

template<typename T>

tensorium::Tensor< T, 5 > tensorium_RG::covariant_scalar_second_3D	(	const tensorium::Tensor< T, 3 > &	chi,
		const tensorium::Tensor< T, 5 > &	christoffel,
		T	dx,
		T	dy,
		T	dz )

References compute_partial_derivatives_scalar(), compute_second_derivatives_scalar(), tensorium::Tensor< K, Rank >::shape(), and X().

Here is the call graph for this function:

◆ covariant_tensor2()

template<typename T, typename TensorFunc>

tensorium::Tensor< T, 3 > tensorium_RG::covariant_tensor2	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		TensorFunc &&	func,
		const tensorium::Tensor< T, 3 > &	Gamma )

References partial_tensor2(), and X().

Referenced by autodiff_rank2_first(), and covariant_tensor2_second().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ covariant_tensor2_second()

template<typename T, typename TensorFunc>

tensorium::Tensor< T, 4 > tensorium_RG::covariant_tensor2_second	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		TensorFunc &&	func,
		const tensorium::Tensor< T, 3 > &	Gamma )

References covariant_tensor2(), tensorium_RG::detail::shift_axis(), tensorium_RG::detail::spatial_indices(), and X().

Here is the call graph for this function:

◆ covariant_vector()

template<typename T, typename VectorFunc>

tensorium::Tensor< T, 2 > tensorium_RG::covariant_vector	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		VectorFunc &&	func,
		const tensorium::Tensor< T, 3 > &	Gamma )

References partial_vector(), and X().

Referenced by autodiff_rank1().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ generate_conformal_metric_field()

template<typename T>

tensorium::Tensor< T, 5 > tensorium_RG::generate_conformal_metric_field	(	const tensorium_RG::Metric< T > &	metric,
		size_t	Nx,
		size_t	Ny,
		size_t	Nz,
		T	dx,
		T	dy,
		T	dz )

Generate the conformal 3-metric field \( \tilde{\gamma}_{ij}(x^i) \) on a 3D grid.

For each spatial grid point \( (x, y, z) \), the spatial metric \( \gamma_{ij} \) is extracted from the given metric, the conformal factor \( \chi \) is computed, and the conformal metric is defined as:

\[\tilde{\gamma}_{ij} = \chi \, \gamma_{ij} \]

Template Parameters

T	Scalar type (typically double)

Parameters

metric	A metric object providing BSSN(X, α, β, γ) and compute_conformal_factor(γ)
Nx,Ny,Nz	Grid resolution in x, y, z directions
dx,dy,dz	Grid spacings in x, y, z directions

Returns: A 5D tensor with shape \( (N_x, N_y, N_z, 3, 3) \) containing \( \tilde{\gamma}_{ij} \) at each grid point

< Coordinates X^μ = (t, x, y, z)

< Time is fixed at t = 0

References tensorium_RG::Metric< T >::BSSN(), tensorium_RG::Metric< T >::compute_conformal_factor(), tensorium_RG::Metric< T >::compute_conformal_metric(), and X().

Here is the call graph for this function:

◆ halo_linear_extrapolate_full()

template<typename T>

void tensorium_RG::halo_linear_extrapolate_full	(	Field3D< T > &	f,
		const GridDims &	D )

inline

References tensorium_RG::Field3D< T >::idx(), and tensorium_RG::Field3D< T >::ptr().

Referenced by tensorium_RG::BoundaryClamp::apply().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ halo_periodic()

template<typename T>

void tensorium_RG::halo_periodic	(	Field3D< T > &	f,
		const GridDims &	D )

inline

References tensorium_RG::Field3D< T >::idx(), and tensorium_RG::Field3D< T >::ptr().

Referenced by tensorium_RG::BoundaryPeriodic::apply().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ load_sym6()

template<typename T>

void tensorium_RG::load_sym6	(	Field3D< T > *	f6,
		size_t	idx,
		T &	xx,
		T &	xy,
		T &	xz,
		T &	yy,
		T &	yz,
		T &	zz )

inline

References tensorium_RG::Field3D< T >::ptr(), XX, XY, XZ, YY, YZ, and ZZ.

Here is the call graph for this function:

◆ make_field()

template<typename T>

Field3D< T > tensorium_RG::make_field ( const Strides< T > & st )

inline

References aligned_alloc_n(), tensorium_RG::Field3D< T >::data, tensorium_RG::Strides< T >::nx_tot, tensorium_RG::Strides< T >::ny_tot, tensorium_RG::Strides< T >::nz_tot, and tensorium_RG::Field3D< T >::st.

Here is the call graph for this function:

◆ pad_simd()

template<typename T>

size_t tensorium_RG::pad_simd ( size_t n )

inlinenoexcept

References simd_width_elems().

Here is the call graph for this function:

◆ partial_scalar()

template<typename T, typename ScalarFunc>

tensorium::Vector< T > tensorium_RG::partial_scalar	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		ScalarFunc &&	func )

References tensorium_RG::detail::shift_axis(), tensorium_RG::detail::spatial_indices(), and X().

Referenced by autodiff_rank0(), tensorium_RG::ChiContext< T >::compute(), tensorium_RG::RicciTildeTensor< T >::compute_dGamma_term(), and covariant_scalar().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ partial_tensor2()

template<typename T, typename TensorFunc>

tensorium::Tensor< T, 3 > tensorium_RG::partial_tensor2	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		TensorFunc &&	func )

References tensorium_RG::detail::shift_axis(), tensorium_RG::detail::spatial_indices(), and X().

Referenced by autodiff_rank2_first(), and covariant_tensor2().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ partial_vector()

template<typename T, typename VectorFunc>

tensorium::Tensor< T, 2 > tensorium_RG::partial_vector	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		VectorFunc &&	func )

References tensorium_RG::detail::shift_axis(), tensorium_RG::detail::spatial_indices(), and X().

Referenced by autodiff_rank1(), and covariant_vector().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ populate_tensor3D_component()

template<typename T, typename TensorFunc>

tensorium::Tensor< T, 3 > tensorium_RG::populate_tensor3D_component	(	size_t	i,
		size_t	j,
		TensorFunc &&	func,
		T	dx,
		T	dy,
		T	dz,
		size_t	NX,
		size_t	NY,
		size_t	NZ )

References X().

Referenced by spectral_partial_tensor2().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ print_tensor2()

template<typename T>

void tensorium_RG::print_tensor2	(	const std::string &	name,
		const tensorium::Tensor< T, 2 > &	tensor )

inline

Referenced by tensorium_RG::BSSN< T >::init_BSSN().

Here is the caller graph for this function:

◆ print_tensor3()

template<typename T>

void tensorium_RG::print_tensor3	(	const std::string &	name,
		const tensorium::Tensor< T, 3 > &	tensor )

inline

Referenced by tensorium_RG::BSSN< T >::init_BSSN().

Here is the caller graph for this function:

◆ print_vector()

template<typename T>

void tensorium_RG::print_vector	(	const std::string &	name,
		const tensorium::Vector< T > &	vec )

inline

References tensorium::Vector< K >::size().

Referenced by tensorium_RG::BSSN< T >::init_BSSN().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ simd_width_elems()

template<typename T>

size_t tensorium_RG::simd_width_elems ( )

inlineconstexprnoexcept

References TENSORIUM_SIMD_WIDTH_F32, and TENSORIUM_SIMD_WIDTH_F64.

Referenced by pad_simd().

Here is the caller graph for this function:

◆ spectral_derivative_1D()

template<typename T>

void tensorium_RG::spectral_derivative_1D	(	tensorium::Vector< T > &	field,
		tensorium::Vector< T > &	dfield,
		T	dx )

References tensorium::Vector< K >::size().

Here is the call graph for this function:

◆ spectral_partial_scalar_3D()

template<typename T>

void tensorium_RG::spectral_partial_scalar_3D	(	const tensorium::Tensor< T, 3 > &	scalar_field,
		T	dx,
		T	dy,
		T	dz,
		tensorium::Tensor< T, 4 > &	grad_out )

References tensorium::Tensor< K, Rank >::resize(), and tensorium::Tensor< K, Rank >::shape().

Referenced by spectral_partial_tensor2().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ spectral_partial_tensor2()

template<typename T, typename TensorFunc>

tensorium::Tensor< T, 3 > tensorium_RG::spectral_partial_tensor2	(	const tensorium::Vector< T > &	X,
		T	dx,
		T	dy,
		T	dz,
		TensorFunc &&	func,
		size_t	NX,
		size_t	NY,
		size_t	NZ )

References populate_tensor3D_component(), spectral_partial_scalar_3D(), and X().

Referenced by autodiff_rank2_first().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ store_sym6()

template<typename T>

void tensorium_RG::store_sym6	(	Field3D< T > *	f6,
		size_t	idx,
		T	xx,
		T	xy,
		T	xz,
		T	yy,
		T	yz,
		T	zz )

inline

References tensorium_RG::Field3D< T >::ptr(), XX, XY, XZ, YY, YZ, and ZZ.

Here is the call graph for this function:

◆ sym6_get()

template<typename T>

T tensorium_RG::sym6_get	(	const Field3D< T > *	f6,
		size_t	idx,
		int	i,
		int	j )

inlinenoexcept

References XX, XY, XZ, YY, YZ, and ZZ.

Referenced by tensorium_RG::bssn::compute_bssn_constraints(), tensorium_RG::bssn::compute_constraint_monitor(), tensorium_RG::bssn::compute_invariant_stats(), tensorium_RG::bssn::compute_rhs_A_tilde(), tensorium_RG::bssn::compute_rhs_Gamma(), tensorium_RG::bssn::compute_rhs_K(), tensorium_RG::bssn::compute_ricci_bssn(), and sym6_inv_get().

Here is the caller graph for this function:

◆ sym6_index()

int tensorium_RG::sym6_index	(	int	i,
		int	j )

inlinenoexcept

References XX, XY, XZ, YY, YZ, and ZZ.

Referenced by tensorium_RG::bssn::compute_bssn_constraints(), tensorium_RG::bssn::compute_rhs_A_tilde(), tensorium_RG::bssn::compute_rhs_Gamma(), tensorium_RG::bssn::compute_rhs_gamma_tilde(), tensorium_RG::bssn::compute_rhs_K(), tensorium_RG::bssn::compute_ricci_bssn(), and tensorium_RG::bssn::detail::diff_gamma_tilde_inv_component().

Here is the caller graph for this function:

◆ sym6_inv_get()

template<typename T>

T tensorium_RG::sym6_inv_get	(	const Field3D< T > *	f6,
		size_t	idx,
		int	i,
		int	j )

inlinenoexcept

References sym6_get().

Here is the call graph for this function:

Namespaces

Classes

Typedefs

Enumerations

Functions

Typedef Documentation

◆ aligned_unique_ptr

Enumeration Type Documentation

◆ DiffMode

◆ Sym6

Function Documentation

◆ __attribute__()

◆ aligned_alloc_n()

◆ autodiff()

◆ autodiff_rank0()

◆ autodiff_rank1()

◆ autodiff_rank2_first()

◆ autodiff_scalar_second()

◆ compute_dt_gamma_from_beta()

◆ compute_dt_gamma_from_metric()

◆ compute_partial_derivatives_3D()

◆ compute_partial_derivatives_scalar()

◆ compute_partial_derivatives_tensor2D()

◆ compute_partial_derivatives_vector()

◆ compute_partial_derivatives_vector3D()

◆ compute_second_derivatives_scalar()

◆ compute_second_derivatives_scalar3D()

◆ compute_second_derivatives_tensor2D()

◆ covariant_scalar()

◆ covariant_scalar_second()

◆ covariant_scalar_second_3D()

◆ covariant_tensor2()

◆ covariant_tensor2_second()

◆ covariant_vector()

◆ generate_conformal_metric_field()

◆ halo_linear_extrapolate_full()

◆ halo_periodic()

◆ load_sym6()

◆ make_field()

◆ pad_simd()

◆ partial_scalar()

◆ partial_tensor2()

◆ partial_vector()

◆ populate_tensor3D_component()

◆ print_tensor2()

◆ print_tensor3()

◆ print_vector()

◆ simd_width_elems()

◆ spectral_derivative_1D()

◆ spectral_partial_scalar_3D()

◆ spectral_partial_tensor2()

◆ store_sym6()

◆ sym6_get()

◆ sym6_index()

◆ sym6_inv_get()

◆ attribute()