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

Aperta
  • Qualcosa non va?
#5,643 0 commenti 0 reazioni 1 assegnatario Vedi su GitHub

I maintainer di solito rispondono entro 1 giorno

@connorjward ci sta già lavorando.

Dal 13/11/2026.

Valutazione

Questa issue non è ancora stata valutata.

Descrizione

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.

Lingua principale
Python
Stelle
684
Fork
204
Merge medio
4g 14h
PR unite (30g)
46

Preparare l'ambiente

  • Nessun Dockerfile né file Docker Compose
  • Ha un modello di pull request
  • Nessuna guida per i contributori

Come iniziare

  1. Leggi tutta la issue e poi la guida ai contributi del progetto.
  2. Commenta sulla issue per dire che te ne occupi tu — evita che due persone facciano lo stesso lavoro.
  3. Fai un fork del repository e lavora su un branch.
  4. Apri una pull request che faccia riferimento al numero della issue.

Altre issue di firedrakeproject/firedrake

Tutte le issue di firedrakeproject/firedrake

Issue simili

Altre issue su Python

Ricevi le nuove issue nella tua casella

Un breve riepilogo di issue GitHub adatte ai principianti.