Thanks !

This is indeed the case:
The older system uses libgsl0ldbl version 1.13+dfsg-1.
The newer system has 1.15+dfsg+1build1.

Thus the packaged Ubuntu 12.04 LTS combination has this problem.
Fortunately, I could install the older libgsl0ldbl version from Ubuntu
10.04 LTS, and now it works again.

You are a great help !

Chris Marshall wrote:
> 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
> <[email protected]> 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
>> [email protected]
>> http://mailman.jach.hawaii.edu/mailman/listinfo/perldl
> 


_______________________________________________
Perldl mailing list
[email protected]
http://mailman.jach.hawaii.edu/mailman/listinfo/perldl

Reply via email to