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.

Efficient Chebyshev Spectral Method for Solving Linear Elliptic Pdes Using Quasi-Inverse Technique

Efficient Spectral-Galerkin Method I. Direct Solvers of Second- and Fourth-Order Equations Using Legendre Polynomials

Spectral Monic Chebyshev Approximation for Higher Order Differential Equations

A Meshless Method for Some Inverse Problems Associated with the Helmholtz Equation

Efficient Spectral and Spectral Element Methods for Eigenvalue Problems of Schrödinger Equations with an Inverse Square Potential

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

An Efficient and Accurate Numerical Method for the Spectral Fractional Laplacian Equation

A spectral element solution of the Poisson equation with shifted boundary polynomial corrections: influence of the surrogate to true boundary mapping and an asymptotically preserving Robin formulation

Efficient spectral–Galerkin methods for fractional partial differential equations with variable coefficients

Efficient and Accurate Spectral Method Using Generalized Jacobi Functions for Solving Riesz Fractional Differential Equations

A robust spectral treatment of a class of initial value problems using modified Chebyshev polynomials

A Highly Efficient Spectral-Galerkin Method Based on Tensor Product for Fourth-Order Steklov Equation with Boundary Eigenvalue

Chebyshev polynomial derivative-based spectral tau approach for solving high-order differential equations

Solving time-dependent PDEs with the ultraspherical spectral method

Efficient Spectral Methods for Quasi-Equilibrium Closure Approximations of Symmetric Problems on Unit Circle and Sphere

On spectral numerical method for variable-order partial differential equations

A new mixed finite-element method for $H^2$ elliptic problems

EFFICIENT SCALING AND MOVING TECHNIQUES FOR SPECTRAL METHODS IN UNBOUNDED DOMAINS

An Efficient High Order Heterogeneous Multiscale Method for Elliptic Problems.

Distinct Numerical Solutions for Elliptic Cross-Interface Problems Using Finite Element and Finite Difference Methods

Efficient Exponential Integrator Finite Element Method for Semilinear Parabolic Equations