Hi Jan-
>From the git log there appear to be no functional changes
to PDL/GSL/INTERP that would affect the result between
PDL-2.4.5 and PDL-2.4.7_001. I believe the problem is that
the underlying GSL changed at some point and a default
error handler was called which caused the Extrapolate option
to fail. There is a commit on 24-Oct-2010, just after the
PDL-2.4.7_001 release that addresses the problem.
If you use the same GSL library version for both systems,
I believe you will get essentially the same output.
--Chris
On Fri, Jul 25, 2014 at 11:55 AM, Jan Hoogenraad
<jan-per...@hoogenraad.net> wrote:
> I have 2 machines on different versions of PDL and GSL.
> Both are pre-packed Ubuntu versions.
> The older one supports PDL::GSL::INTERP with {Extrapolate => 1}
> The somewhat newer one doesn't. Is this known behavior ?
>
> I modified the example code to ask for an extrapolated version.
>
> ===code
>
> use PDL;
> use PDL::GSL::INTERP;
>
> my $x = sequence(10);
> my $y = exp($x);
>
> my $spl = PDL::GSL::INTERP->init('cspline',$x,$y);
>
> my $res = $spl->eval(14.35,{Extrapolate => 1});
> print $res."\n";
>
> === result1
>
> PDL v2.4.5 (supports bad values)
>
> -81198.810683691
>
>
> === result2
>
>
> PDL v2.4.7_001 (supports bad values)
>
> ./tstipol.pl
> gsl: interp.c:150: ERROR: interpolation error
> Default GSL error handler invoked.
> Aborted (core dumped)
>
> _______________________________________________
> Perldl mailing list
> Perldl@jach.hawaii.edu
> http://mailman.jach.hawaii.edu/mailman/listinfo/perldl
_______________________________________________
Perldl mailing list
Perldl@jach.hawaii.edu
http://mailman.jach.hawaii.edu/mailman/listinfo/perldl
