http://xyst.i2p/DEBUG/cppcheck/144.html
Towards the ultimate conservative difference // scheme. II. Monotonicity and conservation combined in a second-order // scheme. Journal of computational physics, 14(4), 361-370. // Derived from the flux limiter phi as: s = phi_inv - (1 - phi) // van Albada limiter //auto sL = std::max(0.0, (2.0*delta1[c]*delta2[c] + muscl_eps) // /(delta1[c]*delta1[c] + delta2[c]*delta2[c] + muscl_eps)); //auto sR = std::max(0.0, (2.0*delta3[c]*delta2[c] + muscl_eps) // ...