Segmentation fault when assembling forms with high-order derivatives on hexahedral elements
Issue #1072
new
Some forms involving high-order derivatives result in segmentation faults during assembly on hexahedral meshes. The following MWE is the simplest example I could come up with to reproduce the behavior (in all versions supporting hexahedral elements):
from dolfin import *
mesh = UnitCubeMesh.create(1,1,1,CellType.Type.hexahedron)
V = FunctionSpace(mesh,"Lagrange",3)
u = TrialFunction(V)
v = TestFunction(V)
A = (u.dx(0).dx(0).dx(0))*(v.dx(0).dx(0).dx(0))*dx
assemble(A)