wrong result of LagrangeInterpolator for vector enriched element spaces

[Jaroslav Hron]

The LagrangeInterpolator seems to give wrong results for vector enriched element spaces.

from dolfin import *

mesh = RectangleMesh(0.0, 0.0, 1.0, 1.0, 4, 4, "crossed")

V = VectorFunctionSpace(mesh, "CG", 2)
B = VectorFunctionSpace(mesh, "Bubble", 3)
U=V+B

e=Expression(("1.0","0.0"))
u=interpolate(e,U)
plot(u)

u=Function(U)
LI=LagrangeInterpolator()
LI.interpolate(u,e)
plot(u)
interactive()

[Jaroslav Hron]

should go to dolfin....