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.