Fixup

Author: pbrubeck

finat/element_factory.py

@@ -234,10 +234,6 @@ def convert_finiteelement(element, **kwargs):
     elif element.variant() is not None:
         finat_kwargs["variant"] = element.variant()
 
-    quad_scheme = element.quadrature_scheme()
-    if quad_scheme is not None:
-        finat_kwargs["quad_scheme"] = quad_scheme
-
     return make_finat_element(cell, element.degree(), **finat_kwargs), deps
 
 

finat/quadrature.py

@@ -26,6 +26,12 @@ def make_quadrature(ref_el, degree, scheme="default"):
     :arg degree: The degree of polynomial that the rule should
         integrate exactly.
     """
+
+    if hasattr(scheme, "dual_basis"):
+        Q, point_set = scheme.dual_basis
+        weights = numpy.ones((len(point_set.points),))
+        return QuadratureRule(point_set, weights, ref_el=ref_el)
+
     if ref_el.get_shape() == TENSORPRODUCT:
         try:
             degree = tuple(degree)