78 lines
1.8 KiB
Plaintext
78 lines
1.8 KiB
Plaintext
|
|
#define tetradtype 2
|
|
|
|
#define Cell
|
|
|
|
#define ghost_width 3
|
|
|
|
|
|
|
|
#define GAUGE 0
|
|
|
|
#define CPBC_ghost_width (ghost_width)
|
|
|
|
#define ABV 0
|
|
|
|
#define EScalar_CC 2
|
|
|
|
#if 0
|
|
|
|
define tetradtype
|
|
v:r; u: phi; w: theta
|
|
tetradtype 0
|
|
v^a = (x,y,z)
|
|
orthonormal order: v,u,w
|
|
m = (phi - i theta)/sqrt(2) following Frans, Eq.(8) of PRD 75, 124018(2007)
|
|
tetradtype 1
|
|
orthonormal order: w,u,v
|
|
m = (theta + i phi)/sqrt(2) following Sperhake, Eq.(3.2) of PRD 85, 124062(2012)
|
|
tetradtype 2
|
|
v_a = (x,y,z)
|
|
orthonormal order: v,u,w
|
|
m = (phi - i theta)/sqrt(2) following Frans, Eq.(8) of PRD 75, 124018(2007)
|
|
|
|
define Cell or Vertex
|
|
Cell center or Vertex center
|
|
|
|
define ghost_width
|
|
2nd order: 2
|
|
4th order: 3
|
|
6th order: 4
|
|
8th order: 5
|
|
|
|
define WithShell
|
|
use shell or not
|
|
|
|
define CPBC
|
|
use constraint preserving boundary condition or not
|
|
only affect Z4c
|
|
CPBC only supports WithShell
|
|
|
|
define GAUGE
|
|
0: B^i gauge
|
|
1: David puncture gauge
|
|
2: MB B^i gauge
|
|
3: RIT B^i gauge
|
|
4: MB beta gauge (beta gauge not means Eq.(3) of PRD 84, 124006)
|
|
5: RIT beta gauge (beta gauge not means Eq.(3) of PRD 84, 124006)
|
|
6: MGB1 B^i gauge
|
|
7: MGB2 B^i gauge
|
|
|
|
define CPBC_ghost_width (ghost_width)
|
|
buffer points for CPBC boundary
|
|
|
|
define ABV
|
|
0: using BSSN variable for constraint violation and psi4 calculation
|
|
1: using ADM variable for constraint violation and psi4 calculation
|
|
|
|
define EScalar_CC
|
|
Type of Potential and Scalar Distribution in F(R) Scalar-Tensor Theory
|
|
1: Case C of 1112.3928, V=0
|
|
2: shell with phi(r) = phi0 * a2^2/(1+a2^2), f(R) = R+a2*R^2 induced V
|
|
3: ground state of Schrodinger-Newton system, f(R) = R+a2*R^2 induced V
|
|
4: a2 = +oo and phi(r) = phi0 * 0.5 * ( tanh((r+r0)/sigma) - tanh((r-r0)/sigma) )
|
|
5: shell with phi(r) = phi0 * Exp(-(r-r0)**2/sigma), V = 0
|
|
|
|
#endif
|
|
|