My intention with this code was to study the behavior of different numeric schemes for spatial derivatives. The main goal was to understand the trade off between high order accuracy and the need of dissipation when you might have high frequency solutions with shocks.

The way the code organization is very rudimentary but I would need a considerable amount of time to polish it nowadays. That being said, I recommend anyone to check it out only if there is no intentions of actually adding new features. That would probably be very unpleasant, to say the least.

The code can be downloaded here. You just need to run the GUI.m file and that’s it. The order of accuracy of the methods are

• WENO     – 4th or 2nd order (3 or 2 points/stencil)
• ENO         – 2nd or 1st order (3 or 2 points/stencil)
• Compact – 6th order
• Centered – 2nd order
• Upwind   – 1st order