On monotone iteration and Schwarz methods for nonlinear parabolic PDEs
Journal of Computational and Applied Mathematics
High-performance modeling acoustic and elastic waves using the parallel Dichotomy Algorithm
Journal of Computational Physics
An Analysis of Equivalent Operator Preconditioning for Equation-Free Newton-Krylov Methods
SIAM Journal on Numerical Analysis
Hi-index | 0.01 |
For a restricted class of parabolic PDEs one can devise a practical numerical solver with a parallel complexity that is theoretically optimal. The method uses a multidimensional FFT to decouple the unknowns in the spatial domain into independent scalar ODEs. These are discretized to give recurrence relations in the time dimension solved by parallel cyclic reduction. This is the FFT/CR algorithm. We discuss the use of FFT/CR as a preconditioner to iteratively solve more general parabolic PDEs. This approach naturally leads to a waveform relaxation scheme. Waveform relaxation was developed as an iterative method for solving large systems of ODEs. It is the continuous-in-time analogue of stationary iterative methods for linear algebraic equations. Using the FFT/CR solver as a preconditioner preserves most of the potential for concurrency that accounts for the attractiveness of waveform relaxation with simple preconditioners like Jacobi or red-black Gauss--Seidel, while showing an important advantage: the convergence rate of the resulting iteration is independent of the mesh size used in the spatial discretization. The method can be accelerated by applying an appropriate scaling of the system before preconditioning.