Dear friends, I have been doing some stuff with a DG1 x CG2 field, as we use it for temperature in the compatible FEM framework. I'm getting some funny output when trying to interpolate. In the following code: https://gist.github.com/colinjcotter/96e23745d40a73303b05 I interpolate a function to a DG1 x CG2 field. The function represented should be continuous. Then, I remap it using a Firedrake kernel to a DG1 x DG1 field. The function should still be continuous. When I visualise the field in Paraview, I see that it has jumps in the vertical direction. What went wrong? cheers --cjc -- http://www.imperial.ac.uk/people/colin.cotter www.cambridge.org/9781107663916