I've computed a displacement field for my simulation using warp.fem, but I'm unclear about the correct way to apply this displacement to update my undeformed mesh geometry.

Initially, I tried using fem.interpolate with PicQuadrature to map displacement fields back onto mesh vertices. However, I realised this interpolation approach only works for structured grids (fem.Grid2D, fem.Grid3D), not for general meshes.

Given that there isn't necessarily a one-to-one correspondence between mesh vertices and polynomial space nodes, what's the recommended way to interpolate or otherwise apply displacement fields computed in arbitrary polynomial spaces back onto a general (possibly unstructured) mesh?

Any guidance or example workflows would be greatly appreciated.

Thanks,
Chris