dolfin-x support for CG elements on mesh skeleton

LEoPart-project / leopart-x

Implementation of Lagrangian-Eulerian Particle tracking with DOLFINx

GNU Lesser General Public License v3.0

8 stars 1 forks source link

dolfin-x support for CG elements on mesh skeleton #5

Open nate-sime opened 4 years ago

nate-sime commented 4 years ago

Consider the following:

from mpi4py import MPI
from dolfinx import *
import ufl

mesh = UnitSquareMesh(MPI.COMM_WORLD, 1, 1)
V = FunctionSpace(mesh, ufl.FiniteElement("CG", mesh.ufl_cell(), 1)["facet"])

which throws the error:

RuntimeError: Cannot handle this element type: <<CG1 on a triangle>>|_{facet}

EDG and EDG-HDG methods rely on CG elements defined on the mesh skeleton.

jmmaljaars commented 4 years ago

Is this bug already posted on dolfin-x issue tracker? Also, maybe @chrisrichardson knows more about the status of this issue