陈 明_预条件技术在多层快速多极子方法中的应用.pdf

Altair 2014 技术大会论文集 预条件技术在多层快速多极子方法中的应用 Accelerating the Multilevel Fast Multipole Method with Parallel Preconditioner for Large-scale Scattering Problems 陈明 (华东电子工程研究所、合肥、230031) 摘 要:针对并行多层快速多极子方法在分析电大开放结构使用电场积分方程收敛慢的问 题，本文利用并行多层快速多极子方法的近场元素 ，构造了一种并行 SAI 预条件技术用来 加速线性方程组迭代求解收敛的速度，提高了并行多层快速多极子方法分析电大开放结构的 效率。数值算例结果表明当未知量的规模达到 100 万时，预条件技术仍然可以取得很好的 改善效果，收敛改善效果达5 ～10 倍。 关键词: 预条件技术 并行技术 矩量法 多层快速多极子方法 Abstract : A novel parallel framework is proposed for the iterative solution of the parallel multilevel fast multipole method (MLFMM). The inversion of the near-field impedance matrix is used as the preconditioning matrix to improve the convergence history of the ill-conditioned linear system formulated by electric field integral equation. the parallel technique is used to construct the preconditioning matrix. Our numerical experiments reveal that with an efficiently parallelized MLFMM and the effective parallel preconditioner, we are able to solve problems with tens of millions of unknowns in a few hours. Both the number of iteration steps and the overall simulation time can be saved greatly. Key words: preconditioning techniques ，parallel technique ，method of moment ， multilevel fast multipole method 1 概述 线性方程组的求解通常有两种方法，一种是直接解法，如高斯消去法、列主元高斯消去 - 1 - Altair 2014 技术大会论文集 法、全主元高斯消去法等。对于未知量为 的问题，直接解法需要 3 的运算操作和 N O( N ) O( N 2 ) 的内存消耗才能实现对系数矩阵的直接分解。直接解法的优点是只需要求解一次方 程就可以求解多右边向量的问题，缺点是计算复杂度太高。对于电大尺寸目标，如此高的运 算操作和内存消耗是限制直接解法应用范围的一个很重要的原因，因而直接解法通常用于求 解小未知量的问题。另一种是迭代解法，如雅可比迭代法(Jacobi) 、高斯-塞德尔迭代法 (Gauss-Seidel)、超松弛迭代法(SSOR) 、共轭梯度法(CG) 、广义最小余量法(GMRES) 、 双共轭梯度法(Bi-CG)等。迭代方法的优点是只需要知道矩阵矢量乘的运算就可以通过