Vector-valued conditionals fail in graph rebuild (uflacs)

Issue #105 resolved
Martin Sandve Alnæs created an issue
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)

  1. Log in to comment