fixing a bug in enforcing the algebraic constraint
Jacob found that currently algebraic constraint is not enforced when the determinant of spatial metric is less than 1. This does not seem to affect the bbh runs, in which the spatial metric is always larger than 1, as I still see convergence in the constraints.