Testing the accuracy and stability of spectral methods in numerical relativity

Creators: Boyle, Michael; Lindblom, Lee; Pfeiffer, Harald P.; Scheel, Mark A.; Kidder, Lawrence E.

Abstract

The accuracy and stability of the Caltech-Cornell pseudospectral code is evaluated using the Kidder, Scheel, and Teukolsky (KST) representation of the Einstein evolution equations. The basic "Mexico City tests" widely adopted by the numerical relativity community are adapted here for codes based on spectral methods. Exponential convergence of the spectral code is established, apparently limited only by numerical roundoff error or by truncation error in the time integration. A general expression for the growth of errors due to finite machine precision is derived, and it is shown that this limit is achieved here for the linear plane-wave test.

Additional Information

©2007 The American Physical Society. (Received 13 September 2006; revised 31 October 2006; published 5 January 2007) We thank Rob Owen, Luisa Buchman, and Olivier Sarbach for helpful conversations. This work was supported in part by a grant from the Sherman Fairchild Foundation to Caltech and Cornell; by NSF Grant Nos. PHY-0099568, PHY-0244906, DMS-0553302, PHY-0601459, and NASA Grant Nos. NAG5-12834, and NNG05GG52G at Caltech; and by NSF Grant Nos. PHY-0312072, PHY-0354631, and NASA Grant No. NNG05GG51G at Cornell.

Files

BOYprd07.pdf

Files (1.3 MB)

Name	Size	Download all
BOYprd07.pdf md5:453e383a3ca3ea81e15fd08b3f3ad68a	1.3 MB	Preview Download

Additional details

	All versions	This version
Views	0	0
Downloads	0	0
Data volume	0 Bytes	0 Bytes