1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
|
Description of files in this directory
======================================
$Header$
This file contains (typically) one file for each model, containing
a Fortran-callable subroutine to compute the 4-metric $g_{ab}$ and
inverse 4-metric $g^{ab}$:
subroutine Exact__model_name_coords(
$ x, y, z, t,
$ gdtt, gdtx, gdty, gdtz,
$ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx,
$ gutt, gutx, guty, gutz,
$ guxx, guyy, guzz, guxy, guyz, guzx,
$ psi, Tmunu_flag)
CCTK_REAL x, y, z, t
CCTK_REAL gdtt, gdtx, gdty, gdtz,
$ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx,
$ gutt, gutx, guty, gutz,
$ guxx, guyy, guzz, guxy, guyz, guzx
CCTK_REAL psi
LOGICAL Tmunu_flag
...
RETURN
END
psi is normally unused, but can be set by the subroutine to explicitly
specify the Cactus static conformal factor. At present only the
"Minkowski_gauge_wave" model does this
The subroutine should set Tmunu_flag to .true. if this model needs
to set the stress-energy tensor (i.e. if it's non-vacuum and/or it has
a cosmological constant, or .false. if it doesn't (i.e. if it's vacuum
and has no cosmological constant).
FIXME: Wouldn't it be cleaner to just compute $g_{ab}$,
and use a 4x4 matrix inversion to get the other one?
|