I would call this only preliminary, since I'm not at all sure it's the
most efficient code. but this can be changed in the future...
also, only interpolation of vector fields is implemented.