General Relativity and Quantum Cosmology
[Submitted on 7 Aug 1995 (v1), last revised 16 May 1996 (this version, v4)]
Title:Finding Apparent Horizons in Numerical Relativity
View PDFAbstract: This paper presents a detailed discussion of the ``Newton's method'' algorithm for finding apparent horizons in 3+1 numerical relativity.
We describe a method for computing the Jacobian matrix of the finite differenced $H(h)$ function by symbolically differentiating the finite difference equations, giving the Jacobian elements directly in terms of the finite difference molecule coefficients used in computing $H(h)$. Assuming the finite differencing scheme commutes with linearization, we show how the Jacobian elements may be computed by first linearizing the continuum $H(h)$ equations, then finite differencing the linearized (continuum) equations.
We find this symbolic differentiation method of computing the $H(h)$ Jacobian to be {\em much} more efficient than the usual numerical perturbation method, and also much easier to implement than is commonly thought.
When solving the discrete $H(h) = 0$ equations, we find that Newton's method generally converges very rapidly. However, if the initial guess for the horizon position contains significant high-spatial-frequency error components, Newton's method has a small (poor) radius of convergence. This is {\em not} an artifact of insufficient resolution in the finite difference grid; rather, it appears to be caused by a strong nonlinearity in the continuum $H(h)$ function for high-spatial-frequency error components in $h$. Robust variants of Newton's method can boost the radius of convergence by O(1) factors, but the underlying nonlinearity remains, and appears to worsen rapidly with increasing initial-guess-error spatial frequency.
Using 4th~order finite differencing, we find typical accuracies for computed horizon positions in the $10^{-5}$ range for $\Delta\theta = \frac{\pi/2}{50}$.
Submission history
From: Jonathan Thornburg [view email][v1] Mon, 7 Aug 1995 09:47:50 UTC (1 KB) (withdrawn)
[v2] Tue, 12 Dec 1995 06:55:01 UTC (1 KB) (withdrawn)
[v3] Wed, 13 Dec 1995 08:49:59 UTC (1 KB) (withdrawn)
[v4] Thu, 16 May 1996 22:48:49 UTC (64 KB)
References & Citations
Bibliographic and Citation Tools
Bibliographic Explorer (What is the Explorer?)
Connected Papers (What is Connected Papers?)
Litmaps (What is Litmaps?)
scite Smart Citations (What are Smart Citations?)
Code, Data and Media Associated with this Article
alphaXiv (What is alphaXiv?)
CatalyzeX Code Finder for Papers (What is CatalyzeX?)
DagsHub (What is DagsHub?)
Gotit.pub (What is GotitPub?)
Hugging Face (What is Huggingface?)
Papers with Code (What is Papers with Code?)
ScienceCast (What is ScienceCast?)
Demos
Recommenders and Search Tools
Influence Flower (What are Influence Flowers?)
CORE Recommender (What is CORE?)
IArxiv Recommender
(What is IArxiv?)
arXivLabs: experimental projects with community collaborators
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.