diff --git a/test/heatequation/heatequation_dg.mini b/test/heatequation/heatequation_dg.mini
index 47559b75a86880bddb1134be30a31dfa4494b8e0..cb5b3ccbf593b27121fe25385ce2aff9620654cf 100644
--- a/test/heatequation/heatequation_dg.mini
+++ b/test/heatequation/heatequation_dg.mini
@@ -14,7 +14,7 @@ extension = vtu
 [formcompiler]
 explicit_time_stepping = 0, 1 | expand scheme
 exact_solution_expression = g
-compare_l2errorsquared = 1e-5
+compare_l2errorsquared = 1e-7
 
 # Disable explicit tests for now
 {__exec_suffix} == explicit | exclude
diff --git a/test/heatequation/heatequation_dg.ufl b/test/heatequation/heatequation_dg.ufl
index 80114e5461b1be6fa898d8a2f28b72d1df5d2e16..21496b69afe311614ec53312c1be5eb032e03133 100644
--- a/test/heatequation/heatequation_dg.ufl
+++ b/test/heatequation/heatequation_dg.ufl
@@ -9,18 +9,21 @@ v = TestFunction(V)
 
 n = FacetNormal(cell)('+')
 
+# penalty factor
 gamma = 1.0
+
+# SSIP: -1.0, IIPG: 0.0, NIPG: 1.0
 theta = 1.0
 
 poisson = inner(grad(u), grad(v))*dx \
   + inner(n, avg(grad(u)))*jump(v)*dS \
   + gamma*jump(u)*jump(v)*dS \
-  + theta*jump(u)*inner(avg(grad(v)), n)*dS \
+  - theta*jump(u)*inner(avg(grad(v)), n)*dS \
   - inner(n, grad(u))*v*ds \
   + gamma*u*v*ds \
-  - theta*u*inner(grad(v), n)*ds \
+  + theta*u*inner(grad(v), n)*ds \
   - f*v*dx \
-  + theta*g*inner(grad(v), n)*ds \
+  - theta*g*inner(grad(v), n)*ds \
   - gamma*g*v*ds
 
 mass = (u*v)*dx