Use anisotropic quadrature order in test and small improvements

python/dune/perftool/sumfact/tabulation.py

@@ -208,7 +208,7 @@ def quadrature_points_per_direction():
         q = (q,) * world_dimension()
 
     # Quadrature points in per direction
-    nb_qp = tuple([(lambda x: x // 2 + 1)(order) for order in q])
+    nb_qp = tuple(order // 2 +1 for order in q)
 
     return nb_qp
 
@@ -229,7 +229,7 @@ def local_quadrature_points_per_direction():
 
 def polynomial_degree():
     form = get_global_context_value("formdata").preprocessed_form
-    degree = form.coefficients()[0].ufl_element()._degree
+    degree = form.coefficients()[0].ufl_element().degree()
     if isinstance(degree, int):
         degree = (degree,) * world_dimension()
     return degree

test/sumfact/poisson/poisson_2d.mini

@@ -18,7 +18,7 @@ numerical_jacobian = 1, 0 | expand num
 compare_l2errorsquared = 4e-5, 4e-9 | expand deg
 sumfact = 1
 vectorize_grads = 1, 0 | expand grad
-# quadrature_order = 2, 4
+quadrature_order = 2, 4
 
 [formcompiler.ufl_variants]
 degree = 1, 2 | expand deg