matlab的欧拉方法代码高效PIDE-2D
作者:Gujji
Reddy,Alan
Seitenfuss,Debora
Medeiros,Luca
Meacci,MiltonAssunção和Michael
Vynnycky
为二维空间上非结构化网格上的抛物线积分微分方程的数值解提供了简短的MATLAB实现。
三角形上的分段线性有限元空间用于空间离散化,而时间离散化则基于后向欧拉法和Crank-Nicolson方法。
选择正交规则以离散化Volterra积分项,以便与时间步长方案一致。
此外,在组装过程中使用矢量化技术介绍了该代码的有效版本,并进行了比较研究。
数值例子证明了该代码的灵活性。
有六个zip文件,其中包含必要的MATLAB文件:
BE-LRR-unvectorized.zip
(后向Euler,左矩形规则,未向量化);
BE-LRR-vectorized.zip
(向后Euler,左矩形规则,矢量化);
BE-RRR-unvectorized.zip
(后向Euler,直角矩形规则,未向量化);
BE-RRR-vectorized.zip
(向后Euler,直角矩形规则
2022-09-09 09:23:40
23KB
系统开源
1