Hierarchical interpolative factorization preconditioner for parabolic equations
From MaRDI portal
Abstract: This note proposes an efficient preconditioner for solving linear and semi-linear parabolic equations. With the Crank-Nicholson time stepping method, the algebraic system of equations at each time step is solved with the conjugate gradient method, preconditioned with hierarchical interpolative factorization. Stiffness matrices arising in the discretization of parabolic equations typically have large condition numbers, and therefore preconditioning becomes essential, especially for large time steps. We propose to use the hierarchical interpolative factorization as the preconditioning for the conjugate gradient iteration. Computed only once, the hierarchical interpolative factorization offers an efficient and accurate approximate inverse of the linear system. As a result, the preconditioned conjugate gradient iteration converges in a small number of iterations. Compared to other classical exact and approximate factorizations such as Cholesky or incomplete Cholesky, the hierarchical interpolative factorization can be computed in linear time and the application of its inverse has linear complexity. Numerical experiments demonstrate the performance of the method and the reduction of conjugate gradient iterations.
Recommendations
- Recursively preconditioned hierarchical interpolative factorization for elliptic partial differential equations
- Multilevel space-time block diagonal preconditioners for parabolic problems
- Efficient preconditioning of variational time discretization methods for parabolic partial differential equations
- scientific article; zbMATH DE number 219982
- Reuse of standard preconditioners for higher-order time discretizations of parabolic PDEs
Cites work
- A fast algorithm for particle simulations
- A fast compact exponential time differencing method for semilinear parabolic equations with Neumann boundary conditions
- A fast Galerkin method for parabolic space-time boundary integral equations
- A fast method for solving the heat equation by layer potentials
- An adaptive fast Gauss transform in two dimensions
- An Algebraic Sparsified Nested Dissection Algorithm Using Low-Rank Approximations
- An Incomplete Factorization Technique for Positive Definite Linear Systems
- Estimating Extremal Eigenvalues and Condition Numbers of Matrices
- Estimating the Largest Eigenvalue by the Power and Lanczos Algorithms with a Random Start
- Explicit unconditionally stable methods for the heat equation via potential theory
- Fast explicit integration factor methods for semilinear parabolic equations
- Fast high-order compact exponential time differencing Runge-Kutta methods for second-order semilinear parabolic equations
- Fast high-order integral equation methods for solving boundary value problems of two dimensional heat equation in complex geometry
- Hierarchical interpolative factorization for elliptic operators: differential equations
- Hierarchical interpolative factorization for elliptic operators: integral equations
- Multidimensional fast Gauss transforms by Chebyshev expansions
- On the Compression of Low Rank Matrices
- Recursively preconditioned hierarchical interpolative factorization for elliptic partial differential equations
- Semi-discrete time-dependent fourth-order problems on an interval: error estimate
- Superfast Multifrontal Method for Large Structured Linear Systems of Equations
- The chemical basis of morphogenesis
- The Fast Gauss Transform
Cited in
(5)- Recursively preconditioned hierarchical interpolative factorization for elliptic partial differential equations
- Hierarchical orthogonal factorization: sparse least squares problems
- A Hierarchical Preconditioner for Wave Problems in Quasilinear Complexity
- Hierarchical orthogonal factorization: sparse square matrices
- Second‐order accurate hierarchical approximate factorizations for solving sparse linear systems
This page was built for publication: Hierarchical interpolative factorization preconditioner for parabolic equations
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2219805)
