Function.assign on a subset fails with a Constant, a Real Function or abs in the expression

Abierto
  • ¿Algo no está bien?
#5,643 0 comentarios 0 reacciones 1 asignado Ver en GitHub

Los mantenedores suelen responder en 1 día

@connorjward ya está trabajando en esto.

Desde el 13/11/2026.

Evaluación

Este issue todavía no se ha evaluado.

Descripción

On current main (ed5fcad68), Function.assign with a subset fails when the expression contains a Constant or a Function in a Real space, and when it contains abs. The same expressions work without a subset.

from firedrake import *

mesh = UnitSquareMesh(2, 2)
V = FunctionSpace(mesh, "CG", 1)
R = FunctionSpace(mesh, "R", 0)
f, g, h = Function(V), Function(V), Function(V)
c = Constant(2.0)
r = Function(R).assign(2.0)
subset = DirichletBC(V, 0, 1).node_set

for name, expr in [("c*g", c*g), ("g + c", g + c), ("r*g", r*g), ("g + r", g + r), ("abs(g - h)", abs(g - h))]:
    f.assign(expr)
    try:
        f.assign(expr, subset=subset)
        print(name, "ok")
    except Exception as e:
        print(name, type(e).__name__)

This prints

c*g IncompatibleAxisTargetException
g + c IncompatibleAxisTargetException
r*g IncompatibleAxisTargetException
g + r IncompatibleAxisTargetException
abs(g - h) TypeError

The IncompatibleAxisTargetException is raised in pyop3/visitors/concretize.py (line 78) during concretize. The TypeError is No handler defined for pyop3.expr.base.Abs from pyop3/visitors/compress_indirections.py.

Connor, you said on #5621 that this is a bug, so I am assigning you as you suggested.

Lenguaje dominante
Python
Estrellas
684
Forks
204
Merge medio
4 d 14 h
PR fusionados (30 d)
46

Preparar el entorno

  • Sin Dockerfile ni archivo de Docker Compose
  • Tiene una plantilla de pull request
  • Sin guía de contribución

Primeros pasos

  1. Lee el issue completo y luego la guía de contribución del proyecto.
  2. Comenta en el issue que vas a ocuparte — evita que dos personas hagan lo mismo.
  3. Haz un fork del repositorio y trabaja en una rama.
  4. Abre un pull request que haga referencia al número del issue.

Más de firedrakeproject/firedrake

Todos los issues de firedrakeproject/firedrake

Issues similares

Más issues de Python

Recibe los nuevos issues en tu correo

Un resumen breve de issues de GitHub para principiantes.