diff --git a/test/sumfact/poisson/CMakeLists.txt b/test/sumfact/poisson/CMakeLists.txt index 07d4ab70639052991c9b9a1d72454a73822f9403..0fd21faffaaaab182f24409f79feef3d724de4f0 100644 --- a/test/sumfact/poisson/CMakeLists.txt +++ b/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 + ) diff --git a/test/sumfact/poisson/poisson_dg_only_volume.mini b/test/sumfact/poisson/poisson_dg_only_volume.mini new file mode 100644 index 0000000000000000000000000000000000000000..70d9a2cdd3b3bd0116f61ef5bb4a93cc0c9ecebe --- /dev/null +++ b/test/sumfact/poisson/poisson_dg_only_volume.mini @@ -0,0 +1,16 @@ +__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 diff --git a/test/sumfact/poisson/poisson_dg_only_volume.ufl b/test/sumfact/poisson/poisson_dg_only_volume.ufl new file mode 100644 index 0000000000000000000000000000000000000000..8c8caa9ccdacfb53d980822cbf3e956463fd1e37 --- /dev/null +++ b/test/sumfact/poisson/poisson_dg_only_volume.ufl @@ -0,0 +1,19 @@ +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