use new Context (evaluation of element of functionspace at a set of points in
a mesh) to speedup eim computation
implement Context
helper functions