Skip to content

BUG: tensor product of Hermite and R fails with unsupported mapping #4952

@danshapero

Description

@danshapero

Describe the bug
Trying to create a tensor product function space with Hermite elements in one direction and constant (R space) in the other gives an error.

Steps to Reproduce
Here's a minimal failing example:

import firedrake

nx, nz = 16, 16
interval = firedrake.UnitIntervalMesh(nx)
mesh = firedrake.ExtrudedMesh(interval, nz)

x_element = firedrake.FiniteElement("Hermite", "interval", 3)
z_element = firedrake.FiniteElement("R", "interval", 0)
element = firedrake.TensorProductElement(x_element, z_element)
Q = firedrake.FunctionSpace(mesh, element)

Error message
I can give a full backtrace if it's useful but here's the final line:

  File "/home/daniel/.pyenv/versions/firedrake-2025-12-18/lib/python3.13/site-packages/finat/ufl/finiteelementbase.py", line 268, in pullback
    raise ValueError(f"Unsupported mapping: {self.mapping()}")
ValueError: Unsupported mapping: undefined

Metadata

Metadata

Assignees

No one assigned

    Labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions