Abstract:The Legendre spectral Galerkin method of self-adjoint second order elliptic equations usually results in a linear system with a dense and ill-conditioned coefficient matrix. In this paper, the linear system is solved by a preconditioned conjugate gradient (PCG) method where the preconditioner M is constructed by approximating the variable coefficients with a (T+1)-term Legendre series in each direction to a desired accuracy. A feature of the proposed PCG method is that the iteration step increases slightly with the size of the resulting matrix when reaching a certain approximation accuracy. The efficiency of the method lies in that the system with the preconditioner M is approximately solved by a one-step method based on the incomplete LU factorization technique with no fill-in, denoted by ILU(0). The ILU(0) factorization of M is an element of R(N-1)dx(N-1)d can be computed using O(T-2d N-d) operations, and the number of nonzeros in the factorization factors is of O(T-d N-d), d = 1, 2, 3. A conclusion of the algorithm is to fast solve the resulting system from the Legendre Galerkin spectral method for Poisson equations with Dirichlet boundary conditions, which has a complexity of O(N-d). To further speed up the PCG method, an algorithm is developed for fast matrix-vector multiplications by the resulting matrix of Legendre-Galerkin spectral discretization, without the need to explicitly form it. The complexity of the fast matrix-vector multiplications is of O(N-d (log(2) N)(2)). In view that T is independent of N in one dimension and is set to be of order O(log(2) N) in two and three dimensions, the PCG method has a O(N-d (log(2) N)(2d)) quasi-optimal complexity for a d dimensional domain with (N-1)(d) unknows, d = 1, 2, 3. In addition, a fast direct solver for the three-dimensional Poisson equation is developed, which is of O(N-3(log(2) N)(2)) and improves the existing results on the computational complexity. Numerical examples are given to demonstrate the efficiency of proposed preconditioners and the algorithm for fast matrix-vector multiplications.

Preconditioning of spectral methods via Birkhoff interpolation

Spectral Analysis of Preconditioned Matrices Arising from Stage-Parallel Implicit Runge–Kutta Methods of Arbitrarily High Order

Preconditioning fractional spectral collocation

Estimating the Condition Number for Multivariate Interpolation Problems

Preconditioned Chebyshev BiCG for parameterized linear systems

Preconditioner Design via Bregman Divergences

Preconditioned Legendre Spectral Galerkin Methods for the Non-separable Elliptic Equation

Spectral analysis for preconditioning of multi-dimensional Riesz fractional diffusion equations

Preconditioned discontinuous Galerkin method and convection-diffusion-reaction problems with guaranteed bounds to resulting spectra

Circulant Preconditioned Iterations for Fractional Diffusion Equations Based on Hermitian and Skew-Hermitian Splittings

Preconditioners based on Voronoi quantizers of random variable coefficients for stochastic elliptic partial differential equations

Preconditioned discontinuous Galerkin method and convection‐diffusion‐reaction problems with guaranteed bounds to resulting spectra

Chebyshev-Type Methods and Preconditioning Techniques

Preconditioner Design via the Bregman Divergence

Spectral stability and perturbation results for kernel differentiation matrices on the sphere

A Bernoulli-barycentric rational matrix collocation method with preconditioning for a class of evolutionary PDEs

Coarse Spaces Based on Higher-Order Interpolation for Schwarz Preconditioners for Helmholtz Problems

A Note on Parallel Preconditioning for the All-at-Once Solution of Riesz Fractional Diffusion Equations

BCCB preconditioners for solving linear systems from delay differential equations

Strang-type Preconditioners for Solving Fractional Diffusion Equations by Boundary Value Methods

Preconditioners for Adaptive Integral Method Implementation