Re: [Help-gsl] The spline polynomial coefficients are unique (????)

[Athanasios Anastasiou]

Dear Brian

Thank you for replying to my message.

Well, i have done the same thing that you are suggesting but using 
MATLAB, testing the output of my implementation at different stages.

In the case of long and very long time series interpolations (more than 
15 and more than 500 x,y pairs respectively) there are errors only in 
the begining and at the end of the interpolation.

The only source of errors i can think of would be round off errors at 
the stage of the solution of the tridiagonal system and the calculations 
that follow this to derive the rest of the coefficients. But then again 
the rest of the coefficients should have some error :-?

I have selected the implementation where a tridiagonal system of 
equations is set up and solved and from this solution all the rest of 
the coefficients are derived.. However, i ran on a site today and there 
was an algorithm in pseudocode which did not set up a system of 
equations at all and it did not look like a recursive based method too. 
Any more information on that? (It was odd because every implementation i 
have seen so far uses a linear system of some sort).

All the best.
thanOS

Brian Gough wrote:

> Athanasios Anastasiou writes:
> > The problem is that i get three different implementations for three
> > different spline interpolation "providers" and honestly i dont know
> > which one to assume true!!! (so much for the uniqueness of the polyonym
> > :-) [only joking])
>
> Hell,
>
> I would suggest using GNU octave to generate the matrix from the
> recurrence relations and solve it, then you can see at what stage the
> discrepancy arises (by looking at the intermediate results in GSL with
> gdb, for example).
>
>