P_0 coefficients

Update hog hash, fixes the compiler issues in frozen velocity operators and also adds P_0 coefficients

operators-epsilon.toml

@@ -93,4 +93,16 @@ quadrature      = 2
 loop-strategy   = "sawtooth"
 operator-name   = "EpsilonP1Viscosity"
 optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending        = "ParametricP2Map"
\ No newline at end of file
+blending        = "ParametricP2Map"
+
+[[epsilon]]
+trial-space               = "P2Vector"
+test-space                = "P2Vector"
+form-space-args.coefficient_function_space = "P0"
+dimensions      = [3]
+quadrature      = 3
+loop-strategy   = "sawtooth"
+operator-name   = "EpsilonP0Viscosity"
+optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending        = "IcosahedralShellMap"
+