On 01/mar/08, at 14:23, Carlo de Falco wrote:
Hi,I found a bug in csape from the package splines v1.0.4 that can be reproduced by doing:>> x = linspace(0,1,3); >> pp = csape(x, x, 'variational'); >> xs = linspace(0,1,100); >> ys = ppval(pp,xs); >> plot(xs,ys)Unfortunately I don't have time to debug this myself at the moment, but I attach a function implementing cubic splines with conditions on the 2nd derivatives that I had implemented myself some time ago, maybe it can serve as a reference to someone willing to solve the problem in csape.c.
BTW, I also have the attached implementation of Catmull-Rom spline interpolation which I wish to use in the msh package.
As "splines" is already listed as a dependency of msh I think this would really belong there , is it ok if I check it in? c.
catmullrom.m
Description: Binary data
------------------------------------------------------------------------- This SF.net email is sponsored by: Microsoft Defy all challenges. Microsoft(R) Visual Studio 2008. http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/
_______________________________________________ Octave-dev mailing list Octave-dev@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/octave-dev
