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
35
36
37
38
39
40
41
42
43
|
# 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 PDL;
use Test;
BEGIN{
eval " use PDL::GSL::DIFF; ";
unless ($@){
plan tests => 4;
}
else {
plan tests => 1;
print "ok 1 # Skipped: PDL::GSL::DIFF not installed\n";
exit;
}
}
@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 );
sub myf{
my ($x) = @_;
return exp($x**2);
}
|
