Brezzi–Douglas–Marini

\(\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:

pip3 install git+https://github.com/fenics/basix.git

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

import basix

# Create Brezzi-Douglas-Marini order 1 on a triangle
element = basix.create_element(basix.ElementFamily.BDM, basix.CellType.triangle, 1)

# Create Brezzi-Douglas-Marini order 2 on a triangle
element = basix.create_element(basix.ElementFamily.BDM, basix.CellType.triangle, 2)

# Create Brezzi-Douglas-Marini order 1 on a tetrahedron
element = basix.create_element(basix.ElementFamily.BDM, basix.CellType.tetrahedron, 1)

# Create Brezzi-Douglas-Marini order 2 on a tetrahedron
element = basix.create_element(basix.ElementFamily.BDM, basix.CellType.tetrahedron, 2)

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 Brezzi-Douglas-Marini order 1 on a triangle
element = symfem.create_element("triangle", "N2div", 1)

# Create Brezzi-Douglas-Marini order 2 on a triangle
element = symfem.create_element("triangle", "N2div", 2)

# Create Brezzi-Douglas-Marini order 1 on a tetrahedron
element = symfem.create_element("tetrahedron", "N2div", 1)

# Create Brezzi-Douglas-Marini order 2 on a tetrahedron
element = symfem.create_element("tetrahedron", "N2div", 2)

Before trying this example, you must install UFL:

pip3 install UFL

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

import ufl

# Create Brezzi-Douglas-Marini order 1 on a triangle
element = ufl.FiniteElement("BDM", "triangle", 1)

# Create Brezzi-Douglas-Marini order 2 on a triangle
element = ufl.FiniteElement("BDM", "triangle", 2)

# Create Brezzi-Douglas-Marini order 1 on a tetrahedron
element = ufl.FiniteElement("BDM", "tetrahedron", 1)

# Create Brezzi-Douglas-Marini order 2 on a tetrahedron
element = ufl.FiniteElement("BDM", "tetrahedron", 2)