BSSN Evolution

In this section, we integrate the 3+1 decomposition of Einstein's Equation in the BSSN formalism using SENR. The code is designed to be compact, human-readable, and easily extended.

Fig. 1: Convergence of the Hamiltonian constraint with increasing finite difference stencil order. A UIUC BH with M = 1 and a = 0, evolved to t/M = 68 using the advecting 1+log and Gamma-driver gauge conditions. The Hamiltonian constraint is measured at r/M = 10 in the wake of the passing gauge wave.

Code Features

Back to Top