Borked H(div) trace tabulation when asking for derivatives
miklos1 opened this issue · comments
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.]])}
Thanks for the issue. I'll have a close look at this tomorrow.
This is being addressed on Bitbucket: FIAT PR #31.