LFP方程的AP方案 该项目提供了用于paper()的代码。 预计算分数拉普拉斯矩阵 可以使用以下代码来计算不同参数的计算矩阵; 我们还提供了本文中使用的所有预先计算的矩阵。 s = 0.5(alpha = 1) M_pre_pl_alp1.m s!= 0.5且0 <s <1 M_pre_pl_new.m 计算分数拉普拉斯算子 佛的例子 图1,2 compute_fl(128,3,0.4,300,2) 图3 FL_Error_vs_N(0.4,3,300,2) 空间同质的情况 对于图〜4、5、6,以下示例约为s = 0.4 LF_hom(128,0.8,1,3,300,1) 空间不均匀的情况 当\ epsilon = 1时,将AP方案与IMEX方案进行比较 对于图〜7,8,以下示例约为s = 0.8,即高斯初始条件 LFP_AP_NHE_ssn(100,128,pi,3,0.8,