Tensorium
Loading...
Searching...
No Matches
tensorium_RG::ConstraintSolver< T > Class Template Reference

#include <BSSNConstraintsSolver.hpp>

Collaboration diagram for tensorium_RG::ConstraintSolver< T >:

Static Public Member Functions

static tensorium::Tensor< T, 3 > solveLichnerowicz (const tensorium::Tensor< T, 5 > &Atilde, const tensorium::Tensor< T, 5 > &g_tilde_inv, T dx, T dy, T dz, int max_iter=1000, T tol=T(1e-8))
 Solve the Lichnerowicz equation to enforce the Hamiltonian constraint in BSSN.

Member Function Documentation

◆ solveLichnerowicz()

template<typename T>
tensorium::Tensor< T, 3 > tensorium_RG::ConstraintSolver< T >::solveLichnerowicz ( const tensorium::Tensor< T, 5 > & Atilde,
const tensorium::Tensor< T, 5 > & g_tilde_inv,
T dx,
T dy,
T dz,
int max_iter = 1000,
T tol = T(1e-8) )
inlinestatic

Solve the Lichnerowicz equation to enforce the Hamiltonian constraint in BSSN.

Parameters
AtildeThe trace-free conformal extrinsic curvature tensor Ã_{ij} at each grid point. Shape: [NX][NY][NZ][3][3]
g_tilde_invThe inverse conformal metric γ̃^{ij} at each grid point. Shape: [NX][NY][NZ][3][3]
dx,dy,dzGrid spacings
max_iterMaximum number of iterations
tolConvergence threshold
Returns
psi The conformal factor ψ at each grid point. Shape: [NX][NY][NZ]

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

Here is the call graph for this function:
The documentation for this class was generated from the following file: