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.

Eigenvalue Analysis and Applications of the Legendre Dual-Petrov-Galerkin Methods for Initial Value Problems.

An Ultraweak-Local Discontinuous Galerkin Method for PDEs with High Order Spatial Derivatives.

Error estimates of a space–time Legendre spectral method for solving the Korteweg–de Vries equation

Efficient Legendre-Laguerre spectral element methods for problems on unbounded domains with diagonalization technique

A Legendre spectral method based on a hybrid format and its error estimation for fourth-order eigenvalue problems

$hp$-version $C^1$-continuous Petrov-Galerkin method for nonlinear second-order initial value problems with application to wave equations

A Local Discontinuous Galerkin Method for Nonlinear Parabolic SPDEs

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

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

A spectral Galerkin exponential Euler time-stepping scheme for parabolic SPDEs on two-dimensional domains with a C2-boundary

Numerical Analysis for Sturm–Liouville Problems with Nonlocal Generalized Boundary Conditions

Analysis of Local Discontinuous Galerkin Methods with Implicit-Explicit Time Marching for Linearized KdV Equations

An absolutely stable discontinuous Galerkin method for the indefinite time-harmonic Maxwell equations with large wave number

Enhancing accuracy with an adaptive discretization for the non-local integro-partial differential equations involving initial time singularities

A staggered discontinuous Galerkin method for a class of nonlinear elliptic equations

Analysis of Discontinuous Galerkin Methods with Upwind-Biased Fluxes for One Dimensional Linear Hyperbolic Equations with Degenerate Variable Coefficients

A high-order local discontinuous Galerkin method for the $p$-Laplace equation

Exponentially accurate spectral and spectral element methods for fractional ODEs

Discontinuous Galerkin Methods for the Helmholtz Equation with Large Wave Number

Error estimates and post-processing of local discontinuous Galerkin method for Schrödinger equations.

A priori and a posteriori error estimates of a DG-CG method for the wave equation in second order formulation