The old approach used a `std::map` and compared coordinates
in the inner most loop to avoid duplicate evaluations
of `LocalBasis`. Instead the new approach determines the
interpolation points first and caches the values in advaced
based on the evaluation order. This improved the performance
significantly.