Radial Basis interpolation multivariable

You just need an array of tuples, via xyz = tuple.(x,y,z), e.g.

julia> x, y, z = rand(1000), rand(1000), rand(1000);

julia> w = @. cos(x^2 - y * sin(z*x));

julia> interp = RadialBasis(tuple.(x,y,z), w, [0,0,0], [1,1,1]);

julia> interp((0.2,0.3,0.4))
0.9999134605897768

julia> cos(0.2^2 - 0.3 * sin(0.4*0.2))
0.9998715929517374
1 Like