快速多级子算法（FMM）介绍和实现
本文摘要（由AI生成）：
本文用简洁的语言介绍了快速多级算法（FMM）的原理和实现，并展示了使用C++开发的FMM模块，用户可以自定义核函数等参数，适用于解决大规模工程问题。FMM通过树形结构减少计算量，提高计算效率。文章还介绍了FMM与边界元方法、矩量法结合解决大规模声场、电磁场问题的应用，并展示了使用Laplace核函数的性能比较，显示了FMM在解决大规模问题时的优势。
国内对FMM（Fast Multipole Method）的介绍都比较复杂，涉及大量的计算公式。
本文试图用最简单的语言介绍FMM的原理和实现，并介绍使用C 开发的快速多级算法模块，用户可自定义核函数，收敛标准，截断系数等参数，可用于实际工程，后续会介绍FMM与边界元方法（Boundary Element Method）/ 矩量法（Method of Moment）结合解决大规模声场，电磁场问题。
在前面简单介绍了快速多级算法 FMM。快速多级子算法能加快解决非对称满秩矩阵，扩展了矩量法，边界元等方法的应用规模，使其能解决较大规模的实际工程问题。
以下是电磁计算软件FEKO中关于多层快速多级算法（MultiLevelFMM）的介绍，是对FMM的一种改进：
https://www.feko.info/product-detail/numerical_methods/mlfmm
1. 概述
FMM算法的提出来源于多粒子系统相互作用的势场计算，比如带电粒子或者天体之间引力等。以静电场为例，空间中N个带点离子构成的系统，第i个离子所在位置的静电势 A（xi） 表示为：
其中
1. i不等于j；
2. xi是第i个粒子所在的坐标；
3. mj是第j个粒子所占的权重，与带电量呈正比；
4. rij是第i个粒子和第j个粒子之间的距离；
按照常规计算方法，对N个粒子实现求和问题，计算量达到O(N^2);
在BEM MOM 等数值计算中，一次这样的求和也就是一次矩阵和向量的乘法迭代。
FMM的实现基本思想是以树形结构为基础，通过多级展开和局部展开，把原对象进行分层分组，将N*N的关系转换为少数组对象之间的关系（如上图），从而减少计算量。该算法实现的核心是如何把每个对象归纳到一组对象中，这个主要是通过动态树结构来实现的。
计算过程如下：
1. 多级展开：
多级展开将叶子节点内（每个小方块内）所有多级展开系数累加，即可以得到该叶子节点的多级展开系数，展开节点为叶子节点的中心。
2.多级-多级转换（Level-2）
将源点和目标点系数进行转换
3.多级-局部转换（Level-2）
利用多级展开向局部展开系数的传递关系式，将该节点的“相互作用表”（紫色直线连接区域部分即为左上角点的相互作用表）中所有节点的多级展开系数传递，并累加到该节点的局部展开系数，局部展开点是该点的中心
4.局部-局部转换（Level-2）
5. 多级-局部转换（Level-3）
6. 最终求和
2. 实现
FMM基础模块，该模块采用C 开发，使用面向对象思想，定义了基础的核函数，并实现了几种简单的核函数，用户可以继承定义自己的核函数，从而应用到流体，电磁，声场等各种需要计算非对称满秩矩阵的数值计算方法中。
基本函数和数据结构：
主要函数：
CreateFMMTree //递归的方法构建节点，生成四叉树或者八叉树结构;
L2L/Local2Local //局部到局部转换;
M2L/Moment2Local //多级到局部转换;
M2M/Moment2Moment //多级到多级转换;
UpwardPass //树结构向上遍历;
DownwardPass //树结构向下遍历;
下图显示了采用Laplace核函数，一次迭代即 矩阵*向量，FMM和高斯迭代 所用的时间比较，可以看出FMM所用时间与N几乎呈线性，而且精度很高。考虑到求解线性方程组时需要计算N次，在N达到几万以后，高斯迭代法在PC上基本求不出结果。
声明：原创文章，欢迎留言与我讨论，如需转载留言
