Add computation of constraint equations in vacuum GR

[FabrizioCorelli]

Proposed changes

Add functions that compute the hamiltonian and momentum constraint in vacuum GR, using the generic 3+1 expressions.

Upgrade instructions

Code review checklist

• The code is documented and the documentation renders correctly. Run
  make doc to generate the documentation locally into BUILD_DIR/docs/html.
  Then open index.html.
• The code follows the stylistic and code quality guidelines listed in the
  code review guide.
• The PR lists upgrade instructions and is labeled bugfix or
  new feature if appropriate.
• If a coding agent is used, have one of
  "Co-Authored-By: Claude Sonnet 4.6 noreply@anthropic.com",
  "Co-Authored-by: Codex noreply@openai.com", or
  "Co-Authored-By: GitHub Copilot CLI noreply@microsoft.com"
  as the last line of the commit, depending on the agent.

Further comments