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 spectral–Galerkin methods for fractional partial differential equations with variable coefficients

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

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

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

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

Exponentially accurate spectral and spectral element methods for fractional ODEs

Semi-implicit Galerkin–Legendre Spectral Schemes for Nonlinear Time-Space Fractional Diffusion–Reaction Equations with Smooth and Nonsmooth Solutions

High‐order finite difference/spectral‐Galerkin approximations for the nonlinear time–space fractional Ginzburg–Landau equation

A Fast Discrete Spectral Method for Stochastic Partial Differential Equations

The Spectral Collocation Method for Efficiently Solving PDEs with Fractional Laplacian.

Fourier pseudospectral methods for the spatial variable-order fractional wave equations

On spectral numerical method for variable-order partial differential equations

Two Efficient Generalized Laguerre Spectral Algorithms for Fractional Initial Value Problems

An efficient multigrid solver for two-dimensional spatial fractional diffusion equations with variable coefficients

The numerical solution of fractional integral equations via orthogonal polynomials in fractional powers

Finite element discretizations for variable-order fractional diffusion problems

Fourier-Gegenbauer Pseudospectral Method for Solving Time-Dependent One-Dimensional Fractional Partial Differential Equations with Variable Coefficients and Periodic Solutions

Spectral Collocation Method for Fractional Differential/Integral Equations with Generalized Fractional Operator

Efficient and Accurate Spectral Method for Solving Fractional Differential Equations on the Half Line Using Orthogonal Generalized Rational Jacobi Functions

Generalized Jacobi Functions and Their Applications to Fractional Differential Equations.

Numerical methods for fractional partial differential equations