Cleanup heatequation tests

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

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