参考文献
Fornberg, B. and Flyer, N., 2015.A primer on radial basis functions with applications to the geosciences. Society for Industrial and Applied Mathematics.
Bollig, E.F., Flyer, N. and Erlebacher, G., 2012. Solution to PDEs using radial basis function finite-differences (RBF-FD) on multiple GPUs. Journal of Computational Physics, 231(21), pp.7133-7151.
数值实现
Matlab 2019a
地球物理局 地震波动力学实验室 无网格组
声明:
# 系列文章优先满足个人研究需求
# 欢迎批评指正,禁止转载
目 录
石中居士:径向基函数生成的有限差分法(RBF-FD)原理及程序实现——目录zhuanlan.zhihu.com节点排序
RBF-FD DM的稀疏模式取决于物理域中节点的排序方式。例如,在球面设置的一个标准的
节点集通常没有特定的节点排序,从而导致图1(a)所示的
结构。
项的这种分散会降低计算效率,特别是对于并行RBF-FD代码中的内存访问,增加了直接稀疏求解器的成本。为了修改节点排序,并获得具有紧密对角化结构的
,使用了带宽减少算法(例如反向Cuthill-McKee(MATLAB命令:
symrcm
)和局部敏感的哈希算法)。后者将物理域分解为粗略的块;在每个块中,使用空间填充曲线(例如,栅格、希尔伯特、皮亚诺)对节点进行重新排序,以使物理空间中相邻的节点在索引空间(内存位置)中相邻。在部分(b)中可以看到对图5.3(a)的矩阵应用反向Cuthill-McKee的结果,而在部分(c)中可以应用局部敏感的哈希算法(按列进行栅格索引)的结果。
在某些情况下,节点最初被有利地排序,例如根据节点的空间紧密度排序。