数学物理方程离散特征值问题的几何网格因式分解算法

来源：计算数学杂志2022年第4期北京时间：

作者：孙家昶

单位：中国科学院软件研究所并行软件与计算科学实验室, 北京 100080

摘要：本文提出求解数学物理方程大型离散特征值问题的几何网格预变换块因式分解算法(简称GPA算法).通过长期研究我们发现:结构化网格矩阵G满足幂等方程G^m=I_N,(m\ll N={m dim}(G)),故可在实数域或复数范围内进行因式分解;且G与有限元刚度矩阵A之间乘法存在互易性:A\cdot G=G\cdot A,利用G的几何不变性可把N阶大型矩阵A正交分解为m-块对角块矩阵异步并行是我们算法的计算数学基础.本文以正三角形、方形、平行六边形及正十七边形等结构化网格为例,特别是详细分析了六边形上的离散特征值异步并行算法及程序实现细节.文后附有若干2-3万阶量级离散矩阵特征值的桌面电脑数值计算例子(正三角形与方形网格,串行加速比分别为3-4倍),符合本文算法分析得出的"几何网格预处理的并行度与正多边形边数成正比"的结论.这类几何网格因式分解算法原则上可推广到三维乃至高维数学物理方程离散特征值计算问题，也可用于大型线性方程组的高效并行求解.

关键词：数理方程离散特征值,互易算子,几何块预处理子,特征值问题因式分解,异步并行算法

基金资助：国家重点研发计划高性能计算重点专项（2016YFB0200601），中国科学院软件研究所基础研究项目资助.