element interpolate takes nodal coords for hermite elements,
operator must call element.interpolate with nodal coords

Additionally, added link to arxiv paper