Incomplete factorization is a widely used preconditioning technique for Krylov subspace methods for solving large-scale sparse linear systems. Its multilevel variants, such as ILUPACK, are more robust for many symmetric or unsymmetric linear systems than the traditional, single-level incomplete LU (or ILU) techniques. However, the previous multilevel ILU techniques still lacked robustness and efficiency for some large-scale saddle-point problems, which often arise from systems of partial differential equations (PDEs). We introduce HILUCSI, or Hierarchical Incomplete LU-Crout with Scalability-oriented and Inverse-based dropping. As a multilevel preconditioner, HILUCSI statically and dynamically permutes individual rows and columns to the next level for deferred factorization. Unlike ILUPACK, HILUCSI applies symmetric preprocessing techniques at the top levels but always uses unsymmetric preprocessing and unsymmetric factorization at the coarser levels. The deferring combined with mixed preprocessing enabled a unified treatment for nearly or partially symmetric systems, and simplified the implementation by avoiding mixed $1\times 1$ and $2\times 2$ pivots for symmetric indefinite systems. We show that this combination improves robustness for indefinite systems without compromising efficiency. Furthermore, to enable superior efficiency for large-scale systems with millions or more unknowns, HILUCSI introduces a scalability-oriented dropping in conjunction with a variant of inverse-based dropping. We demonstrate the effectiveness of HILUCSI for dozens of benchmark problems, including those from the 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.
翻译:不完全的因子化是用于解决大规模分散线性系统的Krylov亚空间方法的一种广泛使用的先决条件技术。 它的多层次变体,如 ILUPACK, 与传统、 单级不完全LU( 或 ILUU) 技术相比,对于许多对称或不对称线性线性系统来说,比传统、 单级不完全LU( 或 ILUU) 的技术更加强大。 然而, 以前的多层次 ILU技术对于一些大型的垫点处理技术仍然缺乏强力和效率, 这些问题经常来自部分差异方程式( PDE) 系统。 我们引入了HILUSI, 或高等级的LUTR( 即高压性 LUF), 推迟和混合前一级LULO( 即低压性能性能), 也使得高或部分的IMLIL( 高压性能) 和高压性压性硬性能系统( ) 显示高或低度性硬性硬性硬性硬性硬性硬性硬性硬性硬性硬化的硬性硬性硬性系统 。