Implement poisson_neumann_symdiff!

and split mini files to avoid reuse of old vtu files!

python/dune/perftool/pdelab/localoperator.py

@@ -213,7 +213,8 @@ class AssemblyMethod(ClassMember):
         from cgen import LiteralLines, Block
         content = signature
         content.append('{')
-        content.extend('  ' + l for l in generate_code(kernel)[0].split('\n'))
+        if kernel is not None:
+            content.extend('  ' + l for l in generate_code(kernel)[0].split('\n'))
         content.append('}')
         ClassMember.__init__(self, content)
 
@@ -314,6 +315,14 @@ def generate_localoperator_kernels(formdata, namedata):
                         kernel = generate_kernel(integral)
                     operator_kernels[(integral.integral_type(), 'jacobian')] = kernel
 
+                # Generate dummy functions for those kernels, that vanished in the differentiation process
+                # We *could* solve this problem by using lambda_* terms but we do not really want that, so
+                # we use empty jacobian assembly methods instead
+                alpha_measures = set(i.integral_type() for i in form.integrals())
+                jacobian_measures = set(i.integral_type() for i in jacform.integrals())
+                for it in alpha_measures - jacobian_measures:
+                    operator_kernels[(it, 'jacobian')] = None
+
         # TODO: JacobianApply for matrix-free computations.
 
     # Return the set of generated kernels

test/poisson/CMakeLists.txt

@@ -5,7 +5,7 @@ add_generated_executable(UFLFILE poisson.ufl
                          )
 
 dune_add_system_test(TARGET poisson_numdiff
-                     INIFILE poisson.mini
+                     INIFILE poisson_numdiff.mini
                      SCRIPT dune_vtkcompare.py
                      )
 
@@ -19,7 +19,7 @@ add_generated_executable(UFLFILE poisson.ufl
                          )
 
 dune_add_system_test(TARGET poisson_symdiff
-                     INIFILE poisson.mini
+                     INIFILE poisson_symdiff.mini
                      SCRIPT dune_vtkcompare.py
                      )
 
@@ -30,6 +30,15 @@ add_generated_executable(UFLFILE poisson_neumann.ufl
                          )
 
 dune_add_system_test(TARGET poisson_neumann_numdiff
-                     INIFILE poisson.mini
+                     INIFILE poisson_neumann_numdiff.mini
                      SCRIPT dune_vtkcompare.py
                      )
+
+# 4. Poisson Test Case: source f, mixed neumann/dirichlet boundary + symbolic differentiation
+add_generated_executable(UFLFILE poisson_neumann.ufl
+                         TARGET poisson_neumann_symdiff
+                         )
+
+dune_add_system_test(TARGET poisson_neumann_symdiff
+                     INIFILE poisson_neumann_symdiff.mini
+                     )

test/poisson/poisson_neumann_numdiff.mini

+__name = poisson_neumann_numdiff
+
+lowerleft = 0.0 0.0
+upperright = 1.0 1.0
+elements = 32 32
+elementType = simplical
+
+[wrapper.vtkcompare]
+name = {__name}
+reference = poisson_ref
+extension = vtu

test/poisson/poisson_neumann_symdiff.mini

+__name = poisson_neumann_symdiff
+
+lowerleft = 0.0 0.0
+upperright = 1.0 1.0
+elements = 32 32
+elementType = simplical
+
+[wrapper.vtkcompare]
+name = {__name}
+reference = poisson_ref
+extension = vtu

test/poisson/poisson.mini → test/poisson/poisson_numdiff.mini

-__name = laplace
+__name = poisson_numdiff
 
 lowerleft = 0.0 0.0
 upperright = 1.0 1.0
@@ -6,6 +6,6 @@ elements = 32 32
 elementType = simplical
 
 [wrapper.vtkcompare]
-name = poisson
+name = {__name}
 reference = poisson_ref
 extension = vtu

test/poisson/poisson_symdiff.mini

+__name = poisson_symdiff
+
+lowerleft = 0.0 0.0
+upperright = 1.0 1.0
+elements = 32 32
+elementType = simplical
+
+[wrapper.vtkcompare]
+name = {__name}
+reference = poisson_ref
+extension = vtu