Black hole initial data by numerical integration of the parabolic–hyperbolic form of the constraints

The parabolic–hyperbolic form of the constraints is integrated numerically. The applied numerical stencil is fourth-order accurate (in the spatial directions) while “time”-integration is made by us...