Tensorium
|
#include <BSSNConstraintsSolver.hpp>
Static Public Member Functions | |
static tensorium::Tensor< T, 3 > | solveLichnerowicz (const tensorium::Tensor< T, 5 > Ã, 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. | |
|
inlinestatic |
Solve the Lichnerowicz equation to enforce the Hamiltonian constraint in BSSN.
Atilde | The trace-free conformal extrinsic curvature tensor Ã_{ij} at each grid point. Shape: [NX][NY][NZ][3][3] |
g_tilde_inv | The inverse conformal metric γ̃^{ij} at each grid point. Shape: [NX][NY][NZ][3][3] |
dx,dy,dz | Grid spacings |
max_iter | Maximum number of iterations |
tol | Convergence threshold |
References psi, and tensorium::Tensor< K, Rank >::shape().
Referenced by tensorium_RG::BSSN< T >::init_BSSN().