1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
|
# Test Script for the PDL interface to the GSL library
# This tests mainly that the interface is working, i.e. that the
# functions can be called.
# The GSL library already has a extensive test suite, and we
# do not want to duplicate that effort here.
use strict;
use warnings;
use Test::More;
use PDL::LiteF;
use PDL::GSL::DIFF;
my @res = gsldiff(\&myf,1.5);
ok(abs($res[0]- 28.4632075095177) < 1e-6 );
@res = gsldiff(\&myf,1.5,{Method => 'central'});
ok(abs($res[0]- 28.4632075095177) < 1e-6 );
@res = gsldiff(\&myf,1.5,{Method => 'forward'});
ok(abs($res[0]- 28.4632852673531) < 1e-6 );
@res = gsldiff(\&myf,1.5,{Method => 'backward'});
ok(abs($res[0]-28.4631297516823 ) < 1e-6 );
done_testing;
sub myf{
my ($x) = @_;
return exp($x**2);
}
|
