High order interpolation operators

Hi All,

I want to do interpolation from solution to new set of mesh (not only different polynomial order, but also mesh nodal points). I have searched poly.py and shape.py and also write.vtk.py, but I cannot find interpolation operator via constructing the high order polynomial based on GL points (or other bases) to have the solution at any point in that element. Just to make sure there’s no such tool to do it, otherwise I will write my own one.

Best wishes,
Zhenyang

This is a highly non-trivial problem in the context of curved grids. Specifically, you need to invert the non-linear mapping function inside each element in order to ascertain where in reference space you need to evaluate the nodal basis to obtain a desired value in physical space.

Regards, Freddie.