Hi Community,

I am working with mixed finite element which contains: H(1), H(d) and H© spaces.

I have tried to use quad/hex meshes on solving my PDE,

from dolfin import *

from mshr import *

mesh = UnitSquareMesh.create(8, 8, CellType.Type.quadrilateral)

CGE = FiniteElement(“CG”, mesh.ufl_cell(), degreeCG)

CurlE = FiniteElement(“N1curl”, mesh.ufl_cell(), degreeCurl) # N1curl and N2curl

DivE = FiniteElement(“BDM”, mesh.ufl_cell(), degreeDiv) # RT and BDM

Z = FunctionSpace(mesh, MixedElement([CGE, CGE, CurlE, CurlE, DivE, DivE]))

but I’ve faced an error:

ufl.log.UFLException: Cellname “quadrilateral” invalid for “Nedelec 1st kind H(curl)” finite element.

So I went through Jack S.Hale’s post and I assumed probably H( d ) and H( c ) are not implemented yet.

It would be great if someone give me a hint on this case.

And if they are not implemented, Is there any other option to implement quad/hex meshes for H( c ) and H( d ).

Thank you in advance.