Test performance of interpolations

interpolations_speed.jl

# based on https://github.com/tlycken/Interpolations.jl/pull/37

Pkg.checkout("Interpolations", "master")
Pkg.checkout("Ratios", "checked")

using Interpolations, Grid

A = rand(300,300,30);
itp = interpolate!(copy(A), BSpline(Quadratic(Flat)), OnCell);
yi = InterpGrid(copy(A), BCreflect, InterpQuadratic);

function interp_all!(dest, itp)
    for k = 1:size(itp,3)
        for j = 1:size(itp,2)
            for i = 1:size(itp,1)
                dest[i,j,k] = itp[i,j,k]
            end
        end
    end
    dest
end

B = similar(A);

# warmup
@time interp_all!(B, itp);
@time interp_all!(B, yi);

@time interp_all!(B, itp);
@time interp_all!(B, yi);
@time interp_all!(B, itp);
@time interp_all!(B, yi);