Parallel Implementation of a Two-level Algebraic ILU(k)-based Domain Decomposition Preconditioner

ABSTRACT We discuss the parallel implementation of a two-level algebraic ILU(k)-based domain decomposition preconditioner using the PETSc library. We present strategies to improve performance and minimize communication among processes during setup and application phases. We compare our implementation with an off-the-shelf preconditioner in PETSc for solving linear systems arising in reservoir simulation problems, and show that for some cases our implementation performs better.

摘要：本文探讨了利用PETSc库实现的、基于ILU(k)的两层代数区域分解预条件子的并行化方案。文中提出了在预条件子构建与应用两个阶段提升性能、降低进程间通信开销的优化策略。针对油藏模拟问题中产生的线性系统求解任务，本文将所提实现方案与PETSc内置的现成预条件子开展对比实验，结果表明在部分场景下本文实现的预条件子具备更优的性能表现。