Issue 1:

Missing symmetrization in the conformal Ricci tensor inside the constraints/analysis block.

Where the bug is

In `fun_bssn_cons`, the code uses

    + Delta[uc] * Gammat[la, lb, lc]

but in the BSSN Ricci decomposition the corresponding term is the symmetrized object

    Delta^c * Gammat_(ab)c

which means

    (1/2) * Delta^--

Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2916/include-cottonmouth