Hi Bhavesh,
On 25 Jun 2021, at 14:29, Lawrence Mitchell <wence@gmx.li> wrote:
Hi Bhavesh,
two things (see below).
On 25 Jun 2021, at 14:12, Bhavesh Shrimali <bhavesh.shrimali@gmail.com> wrote:
Also, for definiteness, the complete code:
from firedrake import * parameters["form_compiler"]["quadrature_degree"]=5 mesh = UnitCubeMesh(5,5,5) P2 = FiniteElement("CG", "tetrahedron", 2) B4 = FiniteElement("B", "tetrahedron", 4) FB3 = FiniteElement("FacetBubble", "tetrahedron", 3) Vh = VectorFunctionSpace(mesh, P2 + B4 + FB3)
I realised one final subtlety here. You are gluing together elements like this, but this construction doesn't end up with an interpolatory element (so your non-zero boundary condition on the mesh face with marker 2 is wrong). You should make a NodalEnrichedElement(P2, B4, FB3) here. Lawrence