快速多极方法的概念性介绍.pdfVIP

中国科学 数学 年 第 卷 第 期 快速多极方法的概念性介绍 应乐兴 收稿日期 接受日期 美国 基金及美国国家科学基金 批准号 资助项目 摘要 本文对快速多极方法作一概念性的非严格的描述, 该方法用来评估具有源分布的卷曲核函数. 我们分别考虑了非振动核情形以及振动核情形. 对于非振动核情形, 我们概述了由 Greengard 和 Rokhlin 提出的经典快速多极方法的主要思想. 对于振动核情形, 我们介绍了最近由Engquist 和Ying 所发展的定向快速多极方法. 关键词 快速多极方法 非振动核 振动核 多极方法 主题分类 引言 许多计算物理问题都需要计算所谓的 N 体问题: 给定一个具有 N 个目标点的集合X d 、一 个具有 N 个源点的集合 Y d 、核函数G(x, y)、以及一个源点的权重集合f (y) : y Y , 我们对 每个x X 计算由下式定义的势能函数u(x), u(x) = ∑ G(x, y)f (y). y ∈Y 这类问题出现在很多领域中. 在天体物理学中, X = Y 代表N 个恒星的位置, G(x, y) = 1/x y 是万有引力势能. 在静电力学中, Y 和 X 分别代表 N 个电荷和 N 个探测器的位置, 核 G(x, y) = 1/x y 是Coulumb 势能. 在调和时间域的波散射的数值解时, X = Y 通常是一个离散点的集合, 且 G(x, y) = exp(iω x y )/x y , 其中ω 是波场的角频率. 直接计算u(x) 需要O(N ) 个步骤, 当N 很大时, 计算代价非常高. 快速多极方法 (FMM) 是在 过去的二十年间发展起来的, 该方法在计算具有可调整精度的近似解时仅需O(N logα N ) 个步骤即可, 其中α = 0, 1. 快速多极方法主要有两类, 一类是关于非振动核的, 另一类是关于振动核的. 本文中, 我 们对每类中的一种方法作一概念性的非严格的介绍, 即分别介绍关于非振动Coulumb 势能 的经 典FMM 算法以及关于振动Helmholtz 核 的定向FMM 算法. 非振动核情形 令 X = Y = P 是由拟一致性地分布于单位方框 Ω = [0, 1] 中的 N 个点构成的集合, 并令 G(x, y) = 1/x y (参见图 1). 我们的目标是快速计算所有成对的相互作用. 英文版发表信息 应乐兴 快速多极方法的概念性介绍 图 在单位方框 中拟一致分布的 个点 B A O(n )2 f (y) u(x)=uA O(n) O(n) O(1) cB cA f =∑ f (y) u =G(c ,c )f B y B A A B B 图 一