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.

Poisson Preconditioning for Self-Adjoint Elliptic Problems

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

Preconditioners for Elliptic Problems Via Non-Uniform Meshes.

Optimal operator preconditioning for pseudodifferential boundary problems

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

A Neural-preconditioned Poisson Solver for Mixed Dirichlet and Neumann Boundary Conditions

Optimal preconditioners for Nitsche-XFEM discretizations of interface problems

A Novel Preconditioner Based On Csl Operator For Solving The Helmholtz Equations

An efficient neural-network and finite-difference hybrid method for elliptic interface problems with applications

Multilevel Preconditioner with Stable Coarse Grid Corrections for the Helmholtz Equation.

Two-level hybrid Schwarz Preconditioners for The Helmholtz Equation with high wave number

Elliptic Preconditioner for Accelerating the Self-Consistent Field Iteration in Kohn-Sham Density Functional Theory.

A CONSTRAINT PRECONDITIONER FOR SOLVING SYMMETRIC POSITIVE DEFINITE SYSTEMS AND APPLICATION TO THE HELMHOLTZ EQUATIONS AND POISSON EQUATIONS

Convergence and Quasi-Optimality of Adaptive FEM with Inhomogeneous Dirichlet Data

Preconditioners for Adaptive Integral Method Implementation

Preconditioned NonSymmetric/Symmetric Discontinuous Galerkin Method for Elliptic Problem with Reconstructed Discontinuous Approximation

Two-Level Schwarz Methods for a Discontinuous Galerkin Approximation of Elliptic Problems with Jump Coefficients

Preconditioners for Higher Order Finite Element Discretizations of H(div)-elliptic Problem

Positive stable preconditioners for symmetric indefinite linear systems arising from Helmholtz equations

Nonoverlapping Domain Decomposition Methods with A Simple Coarse Space for Elliptic Problems

Domain Decomposition Interface Preconditioners For 4th-Order Elliptic Problems