Add epsilon operator with rotation wrappers

operators-epsilon.toml

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