Abstract
We examine spectral and pseudospectral methods as well as waveform relaxation methods for the wave equation in one space dimension. Our goal is to study block Gauss-Jacobi waveform relaxation schemes which can be efficiently implemented in a parallel computing environment. These schemes are applied to semidiscrete systems written in terms of sparse or dense matrices. It is demonstrated that the spectral formulations lead to the implicit system of ordinary differential equations Wã′ = Sã′+g(t) w, with sparse matrices W and S which can be effectively solved by direct application of any Runge-Kutta method. We also examine waveform relaxation iterations based on splittings W = W 1-W 2 and S = S 1 + S 2 and demonstrate that these iterations are only linearly convergent on finite time windows. Waveform relaxation methods applied to the explicit system ã′ = W -1Sa′+g(t) W -1w are somewhat faster but less convenient to implement since the matrix W -1S is no longer sparse. The pseudospectral methods lead to the system Ũ′ = D̃Ũ + g(t) w with a differentiation matrix D̃ of order one and the corresponding waveform relaxation iterations are much faster than the iterations corresponding to the spectral cases (both implicit and explicit).
Original language | English (US) |
---|---|
Pages (from-to) | 1-28 |
Number of pages | 28 |
Journal | Journal of Scientific Computing |
Volume | 20 |
Issue number | 1 |
DOIs | |
State | Published - Feb 2004 |
Keywords
- Spectra and pseudospectra
- Spectral and pseudospectral methods
- Wave equation
- Wave-form relaxation iterations
ASJC Scopus subject areas
- Software
- Theoretical Computer Science
- Numerical Analysis
- General Engineering
- Computational Theory and Mathematics
- Computational Mathematics
- Applied Mathematics