Bugs in quad scheme of quadrature element

Issue #93 new
Jan Blechta created an issue

The following code should pass:

from ufl import *

element = FiniteElement("Quadrature", "tetrahedron", degree=2, quad_scheme='default')
assert element.quadrature_scheme() == 'default'
for e in element.sub_elements():
    assert e.quadrature_scheme == 'default'

element = VectorElement("Quadrature", "tetrahedron", degree=2, quad_scheme='default')
assert element.quadrature_scheme() == 'default'
for e in element.sub_elements():
    assert e.quadrature_scheme == 'default'

element = TensorElement("Quadrature", "tetrahedron", degree=2, quad_scheme='default')
assert element.quadrature_scheme() == 'default'
for e in element.sub_elements():
    assert e.quadrature_scheme == 'default'

The bugs result in problems like https://bitbucket.org/fenics-project/ffc/issues/138.

Comments (2)

  1. Log in to comment