Implement regularizations for Euler eq

Implements the Navier-Stokes and general regularization from J.-L. Guermond and B. Popov, “Viscous regularization of the euler equations and entropy principles,” SIAM Journal on Applied Mathematics, vol. 74, no. 2, pp. 284–305, 2014.

src/EulerEq/rhs.jl

@@ -217,18 +217,18 @@ function rhs!(env, P::Project, hrsc::HRSC.AbstractArtificialViscosity, ::Mesh1d,
   reg = prms.av_regularization
 
   if reg === :mono
-    rhs_mono!(cache, mesh, equation)
+    rhs_mono!(cache, mesh, equation, P)
   elseif reg === :navierstokes
-    rhs_navierstokes!(cache, mesh, equation)
+    rhs_navierstokes!(cache, mesh, equation, P)
   elseif reg === :general
-    rhs_general!(cache, mesh, equation)
+    rhs_general!(cache, mesh, equation, P)
   else
     TODO(reg)
   end
 end
 
 
-function rhs_mono!(cache, mesh, equation)
+function rhs_mono!(cache, mesh, equation, P)
 
   @unpack rho, q, E                = get_dynamic_variables(cache)
   @unpack flx_rho, flx_q, flx_E, p, eps,
@@ -273,7 +273,7 @@ function rhs_mono!(cache, mesh, equation)
 
   return
 end
-function rhs_navierstokes!(cache, mesh, equation)
+function rhs_navierstokes!(cache, mesh, equation, P)
 
   @unpack rho, q, E                        = get_dynamic_variables(cache)
   @unpack flx_rho, flx_q, flx_E, u,