Stability of the iterated Crank-Nicholson method in numerical relativity
Abstract
The iterated Crank-Nicholson method has become a popular algorithm in numerical relativity. We show that one should carry out exactly two iterations and no more. While the limit of an infinite number of iterations is the standard Crank-Nicholson method, it can in fact be worse to do more than two iterations, and it never helps. We explain how this paradoxical result arises.
- Publication:
-
Physical Review D
- Pub Date:
- April 2000
- DOI:
- arXiv:
- arXiv:gr-qc/9909026
- Bibcode:
- 2000PhRvD..61h7501T
- Keywords:
-
- 04.25.Dm;
- 04.20.-q;
- 04.70.-s;
- Numerical relativity;
- Classical general relativity;
- Physics of black holes;
- General Relativity and Quantum Cosmology
- E-Print:
- 2 pages