Vector dPc

\(\mathcal{P}_k=\operatorname{span}\left\{\prod_{i=1}^dx_i^{p_i}\middle|\sum_{i=1}^dp_i\leqslant k\right\}\)

Before trying this example, you must install Basix.UFL:

pip3 install git+https://github.com/FEniCS/basix.git git+https://github.com/FEniCS/ufl.git

This element can then be created with the following lines of Python:

import basix
import basix.ufl

# Create vector dPc order 1 on a quadrilateral
element = basix.ufl.element(basix.ElementFamily.DPC, basix.CellType.quadrilateral, 1, dpc_variant=basix.DPCVariant.simplex_equispaced, discontinuous=True, shape=(2,))

# Create vector dPc order 2 on a quadrilateral
element = basix.ufl.element(basix.ElementFamily.DPC, basix.CellType.quadrilateral, 2, dpc_variant=basix.DPCVariant.simplex_equispaced, discontinuous=True, shape=(2,))

# Create vector dPc order 3 on a quadrilateral
element = basix.ufl.element(basix.ElementFamily.DPC, basix.CellType.quadrilateral, 3, dpc_variant=basix.DPCVariant.simplex_equispaced, discontinuous=True, shape=(2,))

# Create vector dPc order 1 on a hexahedron
element = basix.ufl.element(basix.ElementFamily.DPC, basix.CellType.hexahedron, 1, dpc_variant=basix.DPCVariant.simplex_equispaced, discontinuous=True, shape=(3,))

# Create vector dPc order 2 on a hexahedron
element = basix.ufl.element(basix.ElementFamily.DPC, basix.CellType.hexahedron, 2, dpc_variant=basix.DPCVariant.simplex_equispaced, discontinuous=True, shape=(3,))

Before trying this example, you must install Symfem:

pip3 install symfem

This element can then be created with the following lines of Python:

import symfem

# Create vector dPc order 1 on a quadrilateral
element = symfem.create_element("quadrilateral", "vector dPc", 1)

# Create vector dPc order 2 on a quadrilateral
element = symfem.create_element("quadrilateral", "vector dPc", 2)

# Create vector dPc order 3 on a quadrilateral
element = symfem.create_element("quadrilateral", "vector dPc", 3)

# Create vector dPc order 1 on a hexahedron
element = symfem.create_element("hexahedron", "vector dPc", 1)

# Create vector dPc order 2 on a hexahedron
element = symfem.create_element("hexahedron", "vector dPc", 2)

Before trying this example, you must install (legacy) UFL:

pip3 install fenics-ufl-legacy

This element can then be created with the following lines of Python:

import ufl_legacy

# Create vector dPc order 1 on a quadrilateral
element = ufl_legacy.VectorElement("DPC", "quadrilateral", 1)

# Create vector dPc order 2 on a quadrilateral
element = ufl_legacy.VectorElement("DPC", "quadrilateral", 2)

# Create vector dPc order 3 on a quadrilateral
element = ufl_legacy.VectorElement("DPC", "quadrilateral", 3)

# Create vector dPc order 1 on a hexahedron
element = ufl_legacy.VectorElement("DPC", "hexahedron", 1)

# Create vector dPc order 2 on a hexahedron
element = ufl_legacy.VectorElement("DPC", "hexahedron", 2)