Vector-valued conditionals fail in graph rebuild (uflacs)
Issue #105
resolved
from dolfin import *
import ufl
ufl.algorithms.apply_derivatives.CONDITIONAL_WORKAROUND = False
parameters["form_compiler"]["representation"] = "uflacs"
mesh = UnitIntervalMesh(1)
#d = 1 # works
d = 2 # fails
V = VectorFunctionSpace(mesh, "CG", 1, dim=d)
f = Function(V)
cond = conditional(f**2 < 1.0, f, f)[0]
print assemble(cond*dx)
Reported by @MarcoMorandini
Comments (4)
-
reporter -
Works for me. Thanks!
-
reporter - changed status to resolved
-
- removed milestone
Removing milestone: 1.7 (automated comment)
- Log in to comment
Fix in next. @MarcoMorandini can you test it with your application?