Borked H(div) trace tabulation when asking for derivatives

[miklos1]

import FIAT
from FIAT.reference_element import UFCTetrahedron

elem = FIAT.HDivTrace(UFCTetrahedron(), 0)
print elem.tabulate(1, [[1./3, 1./3]], entity=(2, 0))

prints:

{(1, 0, 0): TraceError('Gradient evaluations are illegal on trace elements.',),
 (0, 1, 0): TraceError('Gradient evaluations are illegal on trace elements.',),
 (0, 0, 0): TraceError('Gradient evaluations are illegal on trace elements.',),
 (0, 0, 1): array([[ 1.],
       [ 0.],
       [ 0.],
       [ 0.]])}

[thomasgibson]

Thanks for the issue. I'll have a close look at this tomorrow.

[thomasgibson]

This is being addressed on Bitbucket: FIAT PR #31.