Adding mapping to VectorTools::interpolate. By default it uses MappingQ1.
Adding Rayleigh-Taylor test case. This still needs some work.