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.

Circulant Preconditioned Iterative Methods for Peridynamic Model Simulation

Block-circulant with circulant-block preconditioners for two-dimensional spatial fractional diffusion equations

Circulant Preconditioned Iterations for Fractional Diffusion Equations Based on Hermitian and Skew-Hermitian Splittings

Fast Iterative Solvers for Numerical Simulations of Scattering and Radiation on Thin Wires

On the Preconditioners for Fast Multipole Boundary Element Methods for 2D Multi-Domain Elastostatics

Preconditioned Nonlinear Conjugate Gradient Method for Real-time Interior-point Hyperelasticity

Circulant-based approximate inverse preconditioners for a class of fractional diffusion equations

A Parallel-in-time Method Based on Preconditioner for Biot's Model

Adaptive Complex Frequency with V-cycle GMRES for Preconditioning 3D Helmholtz Equation

Application of Diagonally Perturbed Incomplete Factorization Preconditioned Conjugate Gradient Algorithms for Edge Finite-Element Analysis of Helmholtz Equations

Preconditioners for Adaptive Integral Method Implementation

R. Chan's circulant-based approximate inverse preconditioning iterative method for solving second-order space fractional advection–dispersion equations with variable coefficients

Preconditioning techniques for a mixed Stokes/Darcy model in porous media applications

A fast computational framework for the linear peridynamic model

Preconditioned iterative methods for the convective FitzHugh–Nagumo equations

On K-Step CSCS-based Polynomial Preconditioners for Toeplitz Linear Systems with Application to Fractional Diffusion Equations

Improved Algebraic Preconditioning for Multiscale Electromagnetic Problems

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

Efficient and Scalable Methods for Dual-Porosity/Permeability Models in Fractured Reservoir Simulations

Matrix-based implementation and GPU acceleration of linearized ordinary state-based peridynamic models in MATLAB