Interpolation in parallel


I would like to interpolate a FE function computed in parallel with PETSc from a first mesh onto a second mesh. Is that possible? What is the best way to do that? Interpolate() or something else? Is there an example available?
Also, is it possible to interpolate when the geometry (domain boundary) differs?

Thank you.

I think your question was answered in this this post.

Thanks! This looks promising indeed. I will give it a try.