计算经验

1. 选用 Global/Local Lax-Friedrichs flux 时，可以把常数乘上1.1。原因在于计算时可能出现一些overshoot导致不是单调的flux。
2. 用DG计算时，在初始时刻会计算精确解的 $L^2$ projection，积分应该选用尽量精确的quadrature rule计算，这样才更接近真实的DG；在最后计算精确解和数值解的误差时，积分应该选用均匀布点计算，原因为了防止某些点处出现superconvergence现象影响观察到真实的收敛阶。
3. CFL数取 $\frac{1}{2k+1}$，其中 $k$ 为有限元空间的多项式次数。另外传播速度 $c$ 应该在每一步都重新计算，来调节时间步长 $dt$。