朱莉娅-韦诺
任意阶的 WENO 重建算法的优化 Julia 实现。基于 Dumbser、Hidalgo 和 Zanotti 在非保守双曲系统的高阶时空自适应 WENO 有限体积方案(DOI 10.1016/j.cma.2013.09.022) 中的工作。
用法
u您的输入数据必须采用 shape数组的形式(nx,ny,nz,nvar)。 nvar是每个单元格中包含的变量数,并且nx,ny,nz分别是x、y、z 轴中的单元格数。如果您的网格是 1D 或 2D,请根据需要设置ny=1和/或nz=1。选择整数N,其中N+1是所需的精度顺序。WENO 重建的系数通过调用获得:
julia> weno(u,N);
该调用每次都会构造 WENO 系数矩阵M1,M2,M3,M4和振荡指标矩阵。Σ如果需要许多 WENO 重建,则预先计算这些实体的计算效率会更高:
julia> M1,M2,M3,M4 = coefficient_matrices(N);
julia> Σ = oscillation_indicator(N);
julia> chΣT = chol(Σ)';
julia> wen