I noticed that UFL does not know FunctionAXPY.
In my code I have a c++ function that return a dolfin.cpp.function.Function f. But I cannot do things like,
plot(f+f)
project(f+f)
because f+f is a FunctionAXPY. I could do
g = Function(f)
plot(g+g)
but I was wondering how hard it would be to convert the FunctionAXPY* type to regular Function coefficient on the fly. I suppose this would have to happen in swig.