Add conservative projection functionality to interpolator

[Fuad-HH]

• runs but results are incorrect
• followed the python example
• the petsc mpi comm changes are temp fix for some error

Todos:

• Get the meshes out of the context struct so that only one copy of them are loaded.
• Merge the interpolators in a single struct??