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
44
45
|
# 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::More;
BEGIN{
# would be nice to use the 'use_ok' routine here
# instead of the following logic, but that needs
# an easy way to find out of the GSL lib is
# available
#
eval " use PDL::GSL::MROOT; ";
if ($@) {
plan skip_all => "PDL::GSL::MROOT not installed";
} else {
plan tests => 2;
}
}
my $init = pdl (-10.00, -5.0);
my $epsabs = 1e-7;
$res = gslmroot_fsolver($init, \&rosenbrock,{Method => 0, EpsAbs => $epsabs});
my @res = list ($res);
ok(abs($res[0]- 1) < 1e-6 );
ok(abs($res[1]- 1) < 1e-6 );
sub rosenbrock{
my ($x) = @_;
my $a = 1;
my $b = 10;
my $y = zeroes($x);
my $tmp; # work around perl -d "feature"
($tmp = $y->slice(0)) .= $a * (1 - $x->slice(0));
($tmp = $y->slice(1)) .= $b * ($x->slice(1) - $x->slice(0)**2);
return $y;
}
|
