FFC fails with conditionals which are functions of a dolfin Function

Using conditionals which are functions of a dolfin Function no longer work. In this example specifically:

kappa = conditional(u > 1.0, u, 1.0)

The following code fails unless a uflacs representation is chosen:

from dolfin import *

import matplotlib.pyplot as plt

# Set uflacs representation, and the code runs
# parameters['form_compiler']['representation'] = 'uflacs'

mesh = UnitSquareMesh(32, 32)
V = FunctionSpace(mesh, 'Lagrange', 1)

u, v = Function(V), TestFunction(V)
du = TrialFunction(V)

kappa = conditional(u > 1.0, u, 1.0)
residual = dot(kappa*grad(u), grad(v))*dx
residual -= Constant(1.0)*v*dx
J = derivative(residual, u, du)

def on_bdry(x, on):
    return on
bcs = [DirichletBC(V, Constant(1.0), on_bdry)]

problem = NonlinearVariationalProblem(residual, u, bcs=bcs, J=J)
solver = NonlinearVariationalSolver(problem)                                                                         
solver.solve()

plot(u, backend='matplotlib')
plt.show()

The error reported:

File "lib/python2.7/site-packages/FFC-1.7.0.dev0-py2.7.egg/ffc/quadrature/quadraturetransformer.py", line 338, in conditional
    "True value of Condtional should only be one function: " + repr(true))
Exception: True value of Condtional should only be one function: {((1, 'k', 3, 3),): 'FE0[0][k]'}

Comments (5)