Data Flow

  +---------------------------------+              +----------------------------+
  | Mathematica Notebook            |------------>>| Initial Data & Parameters  |
  | bimetric-ss-in-3+1 (cpp solver) |  Solve the   |  input.dat  (initial data) |
  +---------------------------------+  constraint  |  config.ini (parameters)   |
            |   Export the             equations   +----------------------------+
            |   C++ code                                         |
            V                                                    V
  +---------------------------------+  eom-std.h   +----------------------------+
  | eomEvolution.h eomMatter.h      |------------>>|       bim-solver.cpp       |
  | eomSources.h   eomConstraints.h |              +----------------------------+
  | eomGauge.h     eomLapseRatios.h |                            |
  | eomMiscEquations.h              |                            V
  +---------------------------------+              +----------------------------+
                                                   |         output.dat         |
     See also:  maximalSlice.h                     +----------------------------+
      
     Numerics:  methodOfLines.h, finiteDifferences.h,
                dataTypes.h, matrix.h, bandSol.h, cubicSpline.h  

Regularization of spherically symmetric evolution codes

• The equations are regularized according to: [1], [2], [4], [5], [6], [10], and [11].

• A nice overview of the parity conditions at r = 0 can be found in [11].

• Also, parity conditions for components of vectors and tensors can be found in Table 1 in [4].