稀疏近似逆与多层块ILU+预条件技术论文.pdf

1 稀疏近似逆与多层块 ILU 预条件技术 a),b) a) b) 谷同祥 迟学斌 刘兴平 a)中科院计算机网络信息中心超级计算中心，北京 349 信箱，100080 b)北京应用物理与计算数学研究所，北京 8009 信箱，100088 E-mail: gtx@jupiter.cnc.ac.cn 摘 要 本文设计了一种求解一般稀疏矩阵的健壮且有效的可并行化预条件子，这种预条件 子涉及在多层块 ILU 预条件子(BILUM) 中使用稀疏近似逆技术。所得的预条件子保持了 BILUM 的健壮性。当使用大块(子区域)来形成块无关集时，提供了控制填充数量的一种新 方法。它比标准的 BILUM 预条件子有两点优势：控制稀疏性的能力和增强了并行性。数值 例子显示了新预条件子的有效性和效率。 关健词 稀疏矩阵, 不完全 LU 分解, 多层块 ILU 预条件子, 稀疏近似逆, Krylov 子空间 方法。 1. 引言 目前，迭代求解技术流行的选择是预条件Krylov 子空间方法，其中一些专门用于求解 非对称(无结构)线性方程组，见如[20]。预条件迭代法的收敛速度通常由预条件子的质量来 决定，有了高质量的预条件子，Krylov 子空间加速就不再是关键[12, 22]，这使得人们由着 重寻找有效的迭代解法器转向构造健壮的预条件子[12] 。并行和分布式处理机的出现与流行 使得寻求适合高性能计算机的可并行化预条件子变得越来越重要。 构造“一般目的”的预条件子需要权衡其有效性、一般性、并行性和健壮性。一般地， 高精度不完全 LU(ILU)类预条件子(具有大量填充)是更健壮的，但它比其低精度形式花费更 多的构造时间。高精度预条件子由于增加了行和列之间的耦合而并行性更差。 对可并行化的预条件子，近年来出现了几种稀疏近似逆技术[3, 8, 9, 10, 13, 14]，这些预 条件子可提供较大的并行性且对某类问题被证明是有效的。然而，对更多的一般稀疏矩阵， 它们的健壮性并不强，未达到可与传统串行预条件子(如 ILUT[18])相媲美的水平。另一方 面，ILUT 直接实现中并行性是有限的，特别是当需要用高精度分解时。 许多作者[10, 13, 27]已经注意到稀疏近似逆技术的一些缺陷和限制。其中之一是由于 内在的局部耦合性使它们对小规模问题更有效[10, 27]。这至少对基于 Frobenius 范数极小 化的稀疏近似逆技术是正确的。因此，这些作者建议在块形式整体预条件子中，如块 SSOR 和基于 Schur 补技术的预条件子中，将稀疏近似逆用做局部预条件子。 多层块 ILU 预条件子(BILUM)[22, 24]推广了具有多层结构的逐次无关集排序的概念， 具有较好的并行度，对某类问题可以导致一个几乎与问题大小无关的收敛速度。BILUM 中无 关块的求解使用精确逆或基于奇异值分解的正规逆(regularized-inverse)，这些策略对小 的块是有效的，然而，求逆的花费随着块大小的增加而急剧增加。对具有大块的基于区域分 解的实现，预条件子的这些构造花费是不能容忍的。另一个相关的问题是较高的存储花费， 因为大块不再是稠密的而其逆却是稠密的。 1 国家重点基础研究专项（G1999032805）和国家高技术研究发展计划（2001AA111043）资助项目 1 我们考察在 BILUM 的构造中用稀疏近似逆技术求大块的逆。对最终的约化系统也构造 一个稀疏近似逆并用作对最终层上近似解的一个解法器。事实上，几乎所有基于 Schur 补技 术的多层方法都使用了某种稀疏近似逆技术，其中一些是非常简单的。对求解某些特殊结构 矩阵的某些多层方法，尽管基于高阶多项式的一些稀疏近似逆是最优的，但它们的实际使用 和有效实现成问题。因此，为了有效性特别是在并行计算机上，它们通常被次最优的稀疏近 似逆所代替。另外，稀疏近似逆也已与多重网格方