Preconditioning for a Class of Spectral Differentiation Matrices

  • Authors:

  • Weiming Cao;Ronald D. Haynes;Manfred R. Trummer

  • Affiliations:

  • Department of Applied Mathematics, University of Texas at San Antonio, San Antonio, USA 78249;Department of Mathematics, Simon Fraser University, Burnaby, Canada V5A 1S6;Department of Mathematics and Centre for Scientific Computing, Simon Fraser University, Burnaby, Canada V5A 1S6

  • Venue:
  • Journal of Scientific Computing
  • Year:
  • 2005

Quantified Score

Hi-index 0.00

Visualization

Abstract

We propose an efficient preconditioning technique for the numerical solution of first-order partial differential equations (PDEs). This study has been motivated by the computation of an invariant torus of a system of ordinary differential equations. We find the torus by discretizing a nonlinear first-order PDE with a full two-dimensional Fourier spectral method and by applying Newton's method. This leads to large nonsymmetric linear algebraic systems. The sparsity pattern of these systems makes the use of direct solvers prohibitively expensive. Commonly used iterative methods, e.g., GMRes, BiCGStab and CGNR (Conjugate Gradient applied to the normal equations), are quite slow to converge. Our preconditioner is derived from the solution of a PDE with constant coefficients; it has a fast implementation based on the Fast Fourier Transform (FFT). It effectively increases the clustering of the spectrum, and speeds up convergence significantly. We demonstrate the performance of the preconditioner in a number of linear PDEs and the nonlinear PDE arising from the Van der Pol oscillator