MUSCL-THINC-BVD and WENO5-THINC-BVD solvers

Comparison two hybrid solvers based on the THINC reconstruction and BVD algorithm.

The present implementation serves the purpose of comparing and summarizing the main hybrid algorithms reported in:[1] Deng, Xi, Bin Xie, and Feng Xiao. "Some practical versions of boundary variation diminishing (BVD) algorithm." arXiv preprint arXiv:1708.01148 (2017) and [2] Deng, Xi, et al. "High fidelity discontinuity-resolving reconstruction for compressible multiphase flows with moving interfaces." Journal of Computational Physics (2018). Where WENO5 and MUSCL methods are hybridized with the THINC reconstruction by means of the Boundary Value Diminishing (BVD) algorithm. As always, these snippets are written so that they are readable rather than fully optimized codes.