How to solve "double free or corruption (!prev)" error in mixed nonlinear problem?

As far as Im aware, real elements have not been implemented in DOLFINx yet, i.e.

import ufl
from petsc4py.PETSc import ScalarType
from mpi4py import MPI
from dolfinx import fem, mesh

# Create mesh
domain = mesh.create_rectangle(MPI.COMM_WORLD, [[- 0.5, - 0.5], [0.5, 0.5]], [20, 20],
                               mesh.CellType.quadrilateral)

# Function spaces and functions
Q = ufl.FiniteElement("Real", domain.ufl_cell(), 0) # Single Lagrange multiplier
W = fem.FunctionSpace(domain, Q)

will throw an error when trying to generate this function space