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.

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

Parallel-in-Time Solver for the All-at-Once Runge–Kutta Discretization

Parallel-in-Time Solver for the All-at-Once Runge--Kutta Discretization

Preconditioning of spectral methods via Birkhoff interpolation

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

A preconditioned MINRES method for optimal control of wave equations and its asymptotic spectral distribution theory

Preconditioning Low Rank Generalized Minimal Residual Method (GMRES) for Implicit Discretizations of Matrix Differential Equations

A robust and time-parallel preconditioner for parabolic reconstruction problems using Isogeometric Analysis

Fine-grained Analysis and Faster Algorithms for Iteratively Solving Linear Systems

Algebraic Spectral Analysis of the DSSR Preconditioner.

Spectral Analysis and Preconditioned Iterative Solvers for Large Structured Linear Systems

Faster Linear Systems and Matrix Norm Approximation via Multi-level Sketched Preconditioning

Optimal and Low-Memory Near-Optimal Preconditioning of Fully Implicit Runge-Kutta Schemes for Parabolic PDEs

Parallelization in time by diagonalization

Analysis of any order Runge-Kutta Spectral Volume Schemes for 1D Hyperbolic Equations

Adaptive rational Krylov methods for exponential Runge--Kutta integrators

Preconditioning techniques for generalized Sylvester matrix equations

Adaptive Rational Krylov Methods for Exponential Runge–Kutta Integrators

A -Preconditioner for Space Fractional Diffusion Equation with Non-separable Variable Coefficients

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

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