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( 高压性能) 和高压性压性硬性能系统( ) 显示高或低度性硬性硬性硬性硬性硬性硬性硬性硬性硬性硬化的硬性硬性硬性系统 。

0
下载
关闭预览

相关内容

【图与几何深度学习】Graph and geometric deep learning,49页ppt
专知会员服务
25+阅读 · 2021年4月2日
Stabilizing Transformers for Reinforcement Learning
专知会员服务
58+阅读 · 2019年10月17日
【哈佛大学商学院课程Fall 2019】机器学习可解释性
专知会员服务
103+阅读 · 2019年10月9日
MIT新书《强化学习与最优控制》
专知会员服务
275+阅读 · 2019年10月9日
计算机类 | PLDI 2020等国际会议信息6条
Call4Papers
3+阅读 · 2019年7月8日
计算机 | 入门级EI会议ICVRIS 2019诚邀稿件
Call4Papers
10+阅读 · 2019年6月24日
深度卷积神经网络中的降采样
极市平台
12+阅读 · 2019年5月24日
Transferring Knowledge across Learning Processes
CreateAMind
27+阅读 · 2019年5月18日
ICLR2019最佳论文出炉
专知
12+阅读 · 2019年5月6日
CCF A类 | 顶级会议RTSS 2019诚邀稿件
Call4Papers
10+阅读 · 2019年4月17日
Unsupervised Learning via Meta-Learning
CreateAMind
42+阅读 · 2019年1月3日
Ray RLlib: Scalable 降龙十八掌
CreateAMind
9+阅读 · 2018年12月28日
Hierarchical Imitation - Reinforcement Learning
CreateAMind
19+阅读 · 2018年5月25日
【学习】Hierarchical Softmax
机器学习研究会
4+阅读 · 2017年8月6日
VIP会员
相关资讯
计算机类 | PLDI 2020等国际会议信息6条
Call4Papers
3+阅读 · 2019年7月8日
计算机 | 入门级EI会议ICVRIS 2019诚邀稿件
Call4Papers
10+阅读 · 2019年6月24日
深度卷积神经网络中的降采样
极市平台
12+阅读 · 2019年5月24日
Transferring Knowledge across Learning Processes
CreateAMind
27+阅读 · 2019年5月18日
ICLR2019最佳论文出炉
专知
12+阅读 · 2019年5月6日
CCF A类 | 顶级会议RTSS 2019诚邀稿件
Call4Papers
10+阅读 · 2019年4月17日
Unsupervised Learning via Meta-Learning
CreateAMind
42+阅读 · 2019年1月3日
Ray RLlib: Scalable 降龙十八掌
CreateAMind
9+阅读 · 2018年12月28日
Hierarchical Imitation - Reinforcement Learning
CreateAMind
19+阅读 · 2018年5月25日
【学习】Hierarchical Softmax
机器学习研究会
4+阅读 · 2017年8月6日
Top
微信扫码咨询专知VIP会员