快速多极子方法并行技术.ppt

快速多极子方法的并行技术 冯仰德 王武 迟学斌 中科院计算机网络信息中心 超级计算中心 ydfeng@sccas.cn 2007年2月5日 纲要 FMM Data Structures Parallelization 纲要 FMM Data Structures Parallelization FMM in Computational Electromagnetics EFIE MFIE CFIE Green函数 积分方程的离散 RWG矢量基函数 MOM 离散 Method of Moments Surface is Discretized into Patches (Basis Functions) 线性系统： Mx=s M是N×N矩阵，x、s是N矢量 ●Direct solution (Gauss elimination,LU Decomposition,SVD,…) 空间复杂度为O(N2) ，需要O(N3)次运算； ●Iterative methods,空间复杂度仍为O(N2)，如果K(kN)步收敛，每步需要的矩阵乘向量的运算为O(N2)； — Moore’Law:processor speed doubles every 18 months — a million variable,16 generations of Moore’Law before a O(N2) algorithm was comparable with a O(N) algorithm — 1 GB RAM =10243=1,073,741,824 bytes = largest N =32,768 ●Finding:快速矩阵乘向量的算法(NlogN)； 并行实施。 Fast Multipole Methods(FMM) Introduced by Rokhlin Greengard in 1987 Called one of the 10 most significant advances in computing of 20th century Speeds up matix-vector products (sums) of a particular type 以上求和要求O(MN)运算复杂度 对给定的精度，FMM可以获得O(M+N)运算复杂度 可以加速matix-vector products ，使O(N2)变为O(NlogN) 加速线性系统求解，如果用迭代方法，k步收敛，每步用矩阵矢量相乘，使计算复杂度由O(N3)变为O(kNlogN) FMM:Application Molecular and Stellar dynamics —computation of force fields and dynamics Solution of acoustical scattering problems —Helmholtz Equation Electromagnetic Wave Scattering — Maxwell’s Equations Fluid Mechanics:Potential flow,vertex flow — Laplace/Poisson Equations FMM: Fundament 格林函数的加法定理 jlpl平面波展开 FMM形式的矩阵向量乘积 FMM FMM FMM Algorithm MLFMM Algorithm Upward Pass: merge scattering matrices Downward Pass: construct splitting and exchange matrices Final Summation MLFMM Algorithm—Upward Pass Step1: 在最细的层盒子求解远场展开系数，xi∈E1(n,l),得到C(n,l)或 ，这也可以用于xi∈E3(n,l) Step2:对于l=L-1,…,2,从step1得到值进行递归得到。同样适合xi∈E3(n,l) 结果：得到分层组聚集系数 MLFMM Algorithm—Downward Pass Step1:l=2,..L递归进行E4 Step2: Key Words 空间多层组划分 Morton编号 相邻组的作用 远场组的上聚 次相邻组中心 的转移 远场组的下推 Grouping 纲要 FMM Data Structures Parallelization Data Structure 构造树 —压