Stability of the iterated Crank-Nicholson method in numerical relativity
Stability of the iterated Crank-Nicholson method in numerical relativity
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 …