Models steepening of high-amplitude sound into shockwaves. Uses sixth-order compact scheme with periodic tridiagonal solver to solve the Euler equation in 1D. Includes filtering via compact scheme but no flux limiting, so it tends to develop oscillations unless lots of filtering is used.
© 2016 Ricky Reusser. WTFPL.