diff --git a/operators.toml b/operators.toml
index 78faa15765a20894eb22fc7f2ae0eb5bc648d73f..d79469f19d1a661e6d4117a68d61c59b18874619 100644
--- a/operators.toml
+++ b/operators.toml
@@ -14,6 +14,14 @@ quadrature    = 0
 loop-strategy = "cubes"
 optimizations = ["moveconstants", "vectorize"]
 
+[[diffusion]]
+trial-space   = "P2"
+test-space    = "P2"
+dimensions    = [2, 3]
+quadrature    = 2
+loop-strategy = "cubes"
+optimizations = ["moveconstants", "vectorize"]
+
 [[div_k_grad]]
 trial-space   = "P1"
 test-space    = "P1"
@@ -40,6 +48,25 @@ quadrature    = 2
 loop-strategy = "sawtooth"
 optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
 
+
+[[divergence]]
+trial-space   = "P2"
+test-space    = "P1"
+components = [0, 1, 2]
+dimensions    = [2, 3]
+quadrature    = 2
+loop-strategy = "cubes"
+optimizations = ["moveconstants", "vectorize"]
+
+[[gradient]]
+trial-space   = "P1"
+test-space    = "P2"
+components = [0, 1, 2]
+dimensions    = [2, 3]
+quadrature    = 2
+loop-strategy = "cubes"
+optimizations = ["moveconstants", "vectorize"]
+
 [[epsilon]]
 trial-space               = "P2"
 test-space                = "P2"
@@ -58,7 +85,7 @@ components-trial = [0, 1, 2]
 components-test  = [0, 1, 2]
 form-space-args.coefficient_function_space = "P2"
 dimensions      = [3]
-quadrature      = 2
+quadrature      = 3
 loop-strategy   = "sawtooth"
 optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
 blending        = "IcosahedralShellMap"