Regrid by applying the regridder object to dst and src arrays

[simone-silvestri]

@navidcy what do you think about this solution?

[simone-silvestri]

This pattern is the same pattern as for Interpolations.jl and Dierckz.jl where you first build an intepolator object and then apply it as a functor (from the Dierckz.jl README):

Fit a 1-d spline to some input data (points can be unevenly spaced):

x = [0., 1., 2., 3., 4.]
y = [-1., 0., 7., 26., 63.]  # x.^3 - 1.
spl = Spline1D(x, y)

Evaluate the spline at some new points:

spl([1.5, 2.5])  # result = [2.375, 14.625]
spl(1.5)  # result = 2.375 

[navidcy]

does LinearAlgebra.mul! return its first arg?

[simone-silvestri]

I think so:

julia> A=[1.0 2.0; 3.0 4.0]; B=[1.0 1.0; 1.0 1.0]; Y = similar(B); mul!(Y, A, B)
2×2 Matrix{Float64}:
 3.0  3.0
 7.0  7.0

[navidcy]

I'll test this with Oceananigans and merge if all seems good