@ziolai I just tried the example and it just works for me:
julia> using BSplineKit
julia> xs = -1:0.1:1;
julia> ys = cospi.(xs);
julia> S = interpolate(xs, ys, BSplineOrder(4))
SplineInterpolation containing the 21-element Spline{Float64}:
basis: 21-element BSplineBasis of order 4, domain [-1.0, 1.0]
order: 4
knots: [-1.0, -1.0, -1.0, -1.0, -0.8, -0.7, -0.6, -0.5, -0.4, -0.3 ā¦ 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 1.0, 1.0, 1.0, 1.0]
coefficients: [-1.0, -1.00111, -0.8975, -0.597515, -0.314147, 1.3265e-6, 0.314142, 0.597534, 0.822435, 0.96683 ā¦ 0.96683, 0.822435, 0.597534, 0.314142, 1.3265e-6, -0.314147, -0.597515, -0.8975, -1.00111, -1.0]
interpolation points: -1.0:0.1:1.0
julia> (Derivative(1) * S)(-1)
-0.01663433622896893
This is using the latest version of BSplineKit.jl (v0.14) and Julia 1.8.3. Can you tell me which versions you are using?
@rafael.guerra Iām sorry if your scripts no longer work using the latest BSplineKit. Can you provide a minimal example of a non-working code, to try to identify the issue?