adds fma mapper method to `SubstitutionMapper`

CMakeLists.txt

@@ -37,8 +37,8 @@ add_subdirectory(cmake/modules)
 add_subdirectory(python)
 
 add_subdirectory(test)
-add_subdirectory(bin)
-add_subdirectory(applications)
+#add_subdirectory(bin)
+#add_subdirectory(applications)
 
 # finalize the dune project, e.g. generating config.h etc.
 finalize_dune_project(GENERATE_CONFIG_H_CMAKE)

python/dune/codegen/loopy/symbolic.py

@@ -122,6 +122,7 @@ lp.type_inference.TypeInferenceMapper.map_vectorized_sumfact_kernel = needs_reso
 
 # FusedMultiplyAdd node
 lp.symbolic.IdentityMapper.map_fused_multiply_add = identity_map_fused_multiply_add
+lp.symbolic.SubstitutionMapper.map_fused_multiply_add = lp.symbolic.SubstitutionMapper.map_variable
 lp.symbolic.WalkMapper.map_fused_multiply_add = walk_map_fused_multiply_add
 lp.symbolic.StringifyMapper.map_fused_multiply_add = stringify_map_fused_multiply_add
 lp.symbolic.DependencyMapper.map_fused_multiply_add = dependency_map_fused_multiply_add

test/blockstructured/poisson/poisson_matrix_free.mini

@@ -13,6 +13,7 @@ compare_l2errorsquared = 1e-7
 
 [formcompiler.r]
 matrix_free = 1
+generate_jacobians = 0
 blockstructured = 1
 number_of_blocks = 4
 geometry_mixins = blockstructured_equidistant

test/poisson/CMakeLists.txt

@@ -81,6 +81,13 @@ dune_add_formcompiler_system_test(UFLFILE poisson_dg_quadrilateral.ufl
         INIFILE poisson_dg_axiparallel.mini
         )
 
+
+# 10. Poisson Test Case with a full permeability tensor
+dune_add_formcompiler_system_test(UFLFILE poisson_tensor_inv.ufl
+        BASENAME poisson_tensor_inv
+        INIFILE poisson_tensor_inv.mini
+        )
+
 # the reference vtk file
 add_executable(poisson_dg_ref reference_main.cc)
 set_target_properties(poisson_dg_ref PROPERTIES EXCLUDE_FROM_ALL 1)

test/poisson/poisson_dg_quadrilateral.ufl

-cell = "quadrilateral"
+generate_jacobian_kernelscell = "quadrilateral"
 degree = 1
 dim = 2