Add sumfact dg only volume tests

test/sumfact/poisson/CMakeLists.txt

@@ -3,9 +3,15 @@
 #                                   BASENAME sumfact_poisson
 #                                   INIFILE poisson.mini
 #                                   )
-#
-# # 2. Poisson Test Case: DG, f + pure dirichlet
-#dune_add_formcompiler_system_test(UFLFILE poisson_dg.ufl
-#                                  BASENAME sumfact_poisson_dg
-#                                  INIFILE poisson_dg.mini
-#                                  )
+
+# 2. Poisson Test Case: DG, f + pure dirichlet
+dune_add_formcompiler_system_test(UFLFILE poisson_dg.ufl
+                                 BASENAME sumfact_poisson_dg
+                                 INIFILE poisson_dg.mini
+                                 )
+
+# 2. Poisson Test Case: DG but only volume parts (for code developement)
+dune_add_formcompiler_system_test(UFLFILE poisson_dg_only_volume.ufl
+                                 BASENAME sumfact_poisson_dg_only_volume
+                                 INIFILE poisson_dg_only_volume.mini
+                                 )

test/sumfact/poisson/poisson_dg_only_volume.mini

+__name = sumfact_poisson_dg_only_volume_{__exec_suffix}
+
+__exec_suffix = {__num_suffix}_{__sumfact_suffix}
+__num_suffix  = numdiff, symdiff |expand num
+__sumfact_suffix = normal, sumfact | expand sumf
+
+cells = 1 1
+extension = 1. 1.
+
+[wrapper.vtkcompare]
+name = {__name}
+extension = vtu
+
+[formcompiler]
+numerical_jacobian = 1, 0 | expand num
+sumfact = 0, 1 | expand sumf

test/sumfact/poisson/poisson_dg_only_volume.ufl

+cell = "quadrilateral"
+
+f = Expression("return -2.0*x.size();", cell=cell)
+g = Expression("return x.two_norm2();", on_intersection=True, cell=cell)
+
+V = FiniteElement("DG", cell, 1)
+
+u = TrialFunction(V)
+v = TestFunction(V)
+
+n = FacetNormal(cell)('+')
+
+gamma = 1.0
+theta = 1.0
+
+r = inner(grad(u), grad(v))*dx \
+  - f*v*dx
+
+forms = [r]
\ No newline at end of file