|
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().