ScatteredInterpolation - basic example from docs not working

Specific Domains Numerics
1

I have tried to use the ScatteredInterpolation package.
When i duplicate the first simple example from the documentation Home · ScatteredInterpolation.jl

using ScatteredInterpolation
samples = [0.0; 0.5; 0.5; 0.5; 1.0];
points = [0.0 0.0; 0.0 1.0; 1.0 0.0; 0.5 0.5; 1.0 1.0];
itp = interpolate(Multiquadratic(), points, samples);

I get

DimensionMismatch: arguments must have the same number of rows

Stacktrace:
 [1] \(F::LinearAlgebra.LU{Float64, Matrix{Float64}, Vector{Int64}}, B::Vector{Float64})
   @ LinearAlgebra ~/bin/julia-1.8.4/share/julia/stdlib/v1.8/LinearAlgebra/src/LinearAlgebra.jl:472
 [2] \(A::Matrix{Float64}, B::Vector{Float64})
   @ LinearAlgebra ~/bin/julia-1.8.4/share/julia/stdlib/v1.8/LinearAlgebra/src/generic.jl:1110
 [3] solveForWeights
   @ ~/.julia/packages/ScatteredInterpolation/tlSER/src/rbf.jl:249 [inlined]
 [4] interpolate(rbf::Multiquadratic{Int64}, points::Matrix{Float64}, samples::Vector{Float64}; metric::Distances.Euclidean, returnRBFmatrix::Bool, smooth::Bool)
   @ ScatteredInterpolation ~/.julia/packages/ScatteredInterpolation/tlSER/src/rbf.jl:201
 [5] interpolate(rbf::Multiquadratic{Int64}, points::Matrix{Float64}, samples::Vector{Float64})
   @ ScatteredInterpolation ~/.julia/packages/ScatteredInterpolation/tlSER/src/rbf.jl:185
 [6] top-level scope
   @ In[11]:3

I wonder what i am missing here. Am i doing something wrong? Is the documentation showing some deprecated syntax here? Is the whole package broken?

2

Works fine for me:

julia> using ScatteredInterpolation

julia> samples = [0.0; 0.5; 0.5; 0.5; 1.0];

julia> points = [0.0 0.0; 0.0 1.0; 1.0 0.0; 0.5 0.5; 1.0 1.0]'
2×5 adjoint(::Matrix{Float64}) with eltype Float64:
 0.0  0.0  1.0  0.5  1.0
 0.0  1.0  0.0  0.5  1.0

julia> itp = interpolate(Multiquadratic(), points, samples)
ScatteredInterpolation.RBFInterpolant{Vector{Float64}, LinearAlgebra.Adjoint{Float64, Matrix{Float64}}, Multiquadratic{Int64}, Distances.Euclidean}([0.9386823120387296, 0.2556696101465102, 0.2556696101465102, -0.7525201751905015, -0.42734309174570895], [0.0 0.0 … 0.5 1.0; 0.0 1.0 … 0.5 1.0], Multiquadratic{Int64}(1), Distances.Euclidean(0.0))

julia> interpolated = evaluate(itp, [0.6; 0.6])
1-element Vector{Float64}:
 0.6105036860019829

It looks to me like you’ve just made an error copy pasting from the docs - you are missing a transpose (') in the definition of points.

3

This is the working solution, thanks and sorry for missing this simple error.