How to get the `fem.Form` of the derivative of a `fem.functionspace` with `basix.ufl.quadrature_element`?

Errors
2

Can be reproduced with a much simpler code:

import basix.ufl
import dolfinx
import ufl
from mpi4py import MPI

mesh = dolfinx.mesh.create_unit_square(MPI.COMM_WORLD, 10, 10)
V = dolfinx.fem.functionspace(mesh, ("Lagrange", 1))
u, v = dolfinx.fem.Function(V), ufl.TestFunction(V)
q_el = basix.ufl.quadrature_element(mesh.topology.cell_name(),value_shape=(),
                                    scheme="default", degree=1)
W = dolfinx.fem.functionspace(mesh, q_el)
c = dolfinx.fem.Function(W)
a = c * ufl.dot(u, v) * ufl.dx(domain=mesh,
                               metadata={"quadrature_degree": q_el.degree})
dolfinx.fem.form(a)

Issue reported at: Basix Quadrature element not complete · Issue #725 · FEniCS/basix · GitHub

2 Likes