HILUCSI: Simple, Robust, and Fast Multilevel ILU for Large-Scale Saddle-Point Problems from PDEs

11/22/2019
by   Qiao Chen, et al.
0

Incomplete factorization is a widely used preconditioning technique for Krylov subspace methods for solving large-scale sparse linear systems. Its multilevel variants, such as those in ILUPACK and ARMS, have been shown to be more robust for many symmetric or unsymmetric linear systems than the traditional, single-level incomplete LU (or ILU) techniques. However, multilevel ILU still lacked robustness and efficiency for some large-scale saddle-point problems, which often arise from systems of partial differential equations (PDEs). In this work, we introduce HILUCSI, or Hierarchical Incomplete LU-Crout with Scalability-oriented and Inverse-based dropping, which is specifically designed to take advantage of some special features of such systems. HILUCSI differs from the state-of-the-art ILU techniques in two main aspects. First, HILUCSI leverages the near or partial symmetry of the underlying problems and the inherent block structures of multilevel ILU to improve robustness and to simplify the treatment of indefinite systems. Second, HILUCSI introduces a scalability-oriented dropping in conjunction with a variant of inverse-based dropping to improve the efficiency for large-scale problems from PDEs. We demonstrate the effectiveness of HILUCSI for a number of benchmark problems, including those from mixed formulation of the Poisson equation, Stokes equations, and Navier-Stokes equations. We also compare its performance with ILUPACK, the supernodal ILUTP in SuperLU, and multithreaded direct solvers in PARDISO and MUMPS.

READ FULL TEXT
research
11/22/2019

HILUCSI: Simple, Robust, and Fast Multilevel ILU with Mixed Symmetric and Unsymmetric Processing

Incomplete factorization is a widely used preconditioning technique for ...
research
11/14/2020

Robust and Efficient Multilevel-ILU Preconditioned Newton-GMRES for Incompressible Navier-Stokes

We introduce a new preconditioned Newton-GMRES method for solving the no...
research
11/11/2020

A staggered-grid multilevel incomplete LU for steady incompressible flows

Algorithms for studying transitions and instabilities in incompressible ...
research
06/18/2021

HIFIR: Hybrid Incomplete Factorization with Iterative Refinement for Preconditioning Ill-conditioned and Singular Systems

We introduce a software package called HIFIR for preconditioning sparse,...
research
04/25/2016

Extreme-scale Multigrid Components within PETSc

Elliptic partial differential equations (PDEs) frequently arise in conti...
research
08/01/2023

A Dual-space Multilevel Kernel-splitting Framework for Discrete and Continuous Convolution

We introduce a new class of multilevel, adaptive, dual-space methods for...
research
04/01/2019

Multigrid Solvers in Reconfigurable Hardware

The problem of finding the solution of Partial Differential Equations (P...

Please sign up or login with your details

Forgot password? Click here to reset