diff --git a/test/sumfact/poisson/CMakeLists.txt b/test/sumfact/poisson/CMakeLists.txt
index 824f2721f163a46d6a76ac9c256c95b3f9885c4c..a09ff1c4c47cf6c36fe598b4f6a06ced2c321603 100644
--- a/test/sumfact/poisson/CMakeLists.txt
+++ b/test/sumfact/poisson/CMakeLists.txt
@@ -34,6 +34,24 @@ dune_add_formcompiler_system_test(UFLFILE poisson_dg_3d.ufl
                                   INIFILE poisson_dg_3d.mini
                                   )
 
+
+#===============================
+# Poisson on 'unstructured grid'
+#===============================
+dune_add_formcompiler_system_test(UFLFILE poisson_2d.ufl
+                                  BASENAME sumfact_poisson_2d_unstructured
+                                  INIFILE poisson_2d_unstructured.mini
+                                  )
+
+#==================================
+# Poisson DG on 'unstructured grid'
+#==================================
+dune_add_formcompiler_system_test(UFLFILE poisson_dg_2d.ufl
+                                  BASENAME sumfact_poisson_dg_2d_unstructured
+                                  INIFILE poisson_dg_2d_unstructured.mini
+                                  )
+
+
 #=============================================
 # Poisson DG using FastDGGridOperator in 2D/3D
 #=============================================
diff --git a/test/sumfact/poisson/poisson_2d_unstructured.mini b/test/sumfact/poisson/poisson_2d_unstructured.mini
new file mode 100644
index 0000000000000000000000000000000000000000..b2df475e90ee5b11b2f511d2b84df347bc31c58e
--- /dev/null
+++ b/test/sumfact/poisson/poisson_2d_unstructured.mini
@@ -0,0 +1,33 @@
+__name = sumfact_poisson_2d_unstructured_{__exec_suffix}
+__exec_suffix = {deg_suffix}_{diff_suffix}_{quadvec_suffix}_{gradvec_suffix}
+
+diff_suffix = numdiff, symdiff | expand num
+quadvec_suffix = quadvec, nonquadvec | expand quad
+gradvec_suffix = gradvec, nongradvec | expand grad
+deg_suffix = deg{formcompiler.ufl_variants.degree}
+
+lowerleft = 0.0 0.0
+upperright = 1.0 1.0
+elements = 16 16
+elementType = quadrilateral
+
+[wrapper.vtkcompare]
+name = {__name}
+extension = vtu
+
+[formcompiler]
+compare_l2errorsquared = 5e-5, 5e-7 | expand deg
+grid_unstructured = 1
+
+[formcompiler.r]
+numerical_jacobian = 1, 0 | expand num
+sumfact = 1
+vectorization_quadloop = 1, 0 | expand quad
+vectorization_strategy = explicit, none | expand grad
+
+# palpo TODO
+# diagonal_transformation_matrix = 1
+# constant_transformation_matrix = 1
+
+[formcompiler.ufl_variants]
+degree = 1, 2 | expand deg
diff --git a/test/sumfact/poisson/poisson_dg_2d_unstructured.mini b/test/sumfact/poisson/poisson_dg_2d_unstructured.mini
new file mode 100644
index 0000000000000000000000000000000000000000..f5c76c1352b7a269598a6289d754720a40cb411a
--- /dev/null
+++ b/test/sumfact/poisson/poisson_dg_2d_unstructured.mini
@@ -0,0 +1,32 @@
+__name = sumfact_poisson_dg_2d_unstructured_{__exec_suffix}
+__exec_suffix = {deg_suffix}_{diff_suffix}_{quadvec_suffix}_{gradvec_suffix}
+
+diff_suffix = numdiff, symdiff | expand num
+quadvec_suffix = quadvec, nonquadvec | expand quad
+gradvec_suffix = gradvec, nongradvec | expand grad
+deg_suffix = deg{formcompiler.ufl_variants.degree}
+
+lowerleft = 0.0 0.0
+upperright = 1.0 1.0
+elements = 16 16
+elementType = quadrilateral
+
+[wrapper.vtkcompare]
+name = {__name}
+extension = vtu
+
+[formcompiler]
+compare_l2errorsquared = 5e-5, 5e-7 | expand deg
+grid_unstructured = 1
+
+[formcompiler.r]
+numerical_jacobian = 1, 0 | expand num
+sumfact = 1
+vectorization_quadloop = 1, 0 | expand quad
+vectorization_strategy = explicit, none | expand grad
+# palpo TODO
+diagonal_transformation_matrix = 1
+constant_transformation_matrix = 1
+
+[formcompiler.ufl_variants]
+degree = 1, 2 | expand deg