第七章并行三维位势多极边界元法.docVIP

第七章 并行多极三维位势边界元法 7.1 基于FMM的Krylov空间的广义极小残值法[GMRES(m)] GMRES(m)算法是一种求解大型非对称线性方程组的Krylov子空间投影法，1986年由Y.Saad和M.H.Schltz提出。该法基于Krylov向量的完全正交化，因而在Krylov子空间上只求出一个近似解或者拟最小残值解。由于所需计算量和存储量较少等方面的优势，目前广泛应用于各种工程应用领域。 在工程问题的数值解模拟中，无论是边界元法还是有限元法都将形成离散方程矩阵，二者不同之处在于矩阵系数及其性质上。有限元法形成的刚度系数是稀疏矩阵，边界元法形成的影响系数矩阵是稠密矩阵。 7.1. 1 KrylovGMRES(m)算法的基本理论 线性方程组的求解方法有很多种。传统边界元法求解线性方程组，一般采用高斯消去法或是由高斯消去法派生的相关方法，解方程所需要的计算次数与方程秩的成比例。 广义极小残值算法是在Arnoldi算法的基础上提出的稠密系数矩阵线性方程组的解法，计算次数只与成比例。 对于大型线性代数方程组 取为任意一向量，令，则式化为 式中，。对于固定整数，令左右空间分别为 在子空间中找到式的近似解，使得 (7-4