MATLAB求解拉普拉斯代码二维边界积分算子的高阶局部校正梯形规则
这是论文随附的“zeta
校正正交”的
MATLAB
代码:
B.
Wu
和
PG
Martinsson,Zeta
校正:构建奇异积分算子的校正梯形规则的新方法。
(2020
年,arxiv)
它包含修正的梯形规则,这些规则对于平滑闭合曲线上的拉普拉斯、亥姆霍兹和斯托克斯层势是稳定且高阶准确的。
作者:吴博伟,2020/7,更新于2021/1
还包含从
Alex
Barnett
的包修改的支持功能
另见表面积分的相关代码。
依赖的注意事项:
最多
42
个订单的
Zeta
更正使用预先计算的权重。
(它们对应于kapur_rokhlin_sep_log.m函数的输入k
<=
21
)对于更高阶的更正,需要来自
Symbolic
Math
Toolbox
的vpa函数。
(此存储库中的测试文件不需要
Symbolic
Toolbox。)
例子
以下是使用高阶
zeta
正交求解斯托克斯和亥姆霍兹外边界值问题的测试文件的一些示例输出。
收敛性根据边界上使用的点数显示。
主要测试文件说明:
test_sto2d_bie.m
2021-11-20 20:31:58
196KB
系统开源
1