Interpolation in parallel

The following example crashes, because the interpolation points are not on the same process after refining the Mesh. An efficient and scalable approach is needed. Interpolation and evaluation in parallel should work seamlessly, provided the points are within the global mesh.

from dolfin import *
mesh = UnitSquareMesh(10,10)
Q = FunctionSpace(mesh, "CG", 1)
F = Function(Q)
F.interpolate(Expression("x[0]"))

mesh2 = refine(mesh, redistribute=True)
Q2 = FunctionSpace(mesh2, "CG", 1)
F2 = interpolate(F, Q2)

Comments (3)