A PDE compiler for full-metric numerical relativity.
Abstract
This paper describes work on 3+1 numerical relativity with a full (all components nonzero) metric, such as arises in the Smarr/York minimal distortion gauge. The author presents a technique for overcoming the resulting complexity of the Einstein equations, by using a "PDE compiler" to automatically generate a finite differencing code from a high-level description of the tensor differential operators involved. The author has constructed a prototype PDE compiler, and a full-metric numerical relativity code based on it. The code is fully 3-covariant and uses minimal distortion spatial coordinates to evolve a vacuum axisymmetric spacetime containing a single black hole present in the initial data. It uses 4th order centered finite differencing on a non-staggered uniform grid.
- Publication:
-
Frontiers in Numerical Relativity
- Pub Date:
- 1989
- Bibcode:
- 1989fnr..book..370T
- Keywords:
-
- General Relativity: Numerical Methods