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
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
|
package Math::GSL::Min::Test;
use base q{Test::Class};
use strict;
use Test::Most;
use Math::GSL qw/:all/;
use Math::GSL::Min qw/:all/;
use Math::GSL::Test qw/:all/;
use Math::GSL::Errno qw/:all/;
use Math::GSL::Const qw/$M_PI/;
#use Devel::Trace qw/trace/;
use Data::Dumper;
BEGIN { gsl_set_error_handler_off(); }
sub trace { }
sub make_fixture : Test(setup) {
my $self = shift;
$self->{min} = gsl_min_fminimizer_alloc($gsl_min_fminimizer_goldensection);
$self->{brent} = gsl_min_fminimizer_alloc($gsl_min_fminimizer_brent);
}
sub teardown : Test(teardown) {
my $self = shift;
}
sub GSL_MIN_TYPES : Tests {
my $m = gsl_min_fminimizer_alloc($gsl_min_fminimizer_goldensection);
isa_ok($m, 'Math::GSL::Min');
my $n = gsl_min_fminimizer_alloc($gsl_min_fminimizer_brent);
isa_ok($n, 'Math::GSL::Min');
}
sub GSL_MIN_NAME : Tests {
my $self = shift;
cmp_ok( 'goldensection', 'eq', gsl_min_fminimizer_name($self->{min}) );
cmp_ok( 'brent', 'eq', gsl_min_fminimizer_name($self->{brent}) );
}
sub GSL_MIN_SET : Tests {
my $self = shift;
my $mini = $self->{min};
ok_status(
gsl_min_fminimizer_set($mini,
sub { cos($_[0]) }, 3, 0, 2*$M_PI
),
);
# These are the first guesses first the initial iteration
cmp_ok( $mini->{x_minimum}, '==', 3 );
cmp_ok( $mini->{x_lower}, '==', 0 );
isa_ok( $mini->{function} , 'Math::GSL::Min::gsl_function_struct');
ok_similar([$mini->{x_upper}], [2*$M_PI] );
ok_similar([$mini->{f_minimum}], [cos(3)] );
ok_similar([$mini->{f_lower}], [cos(0)] );
ok_similar([$mini->{f_upper}], [cos(2*$M_PI)] );
ok_similar([gsl_min_fminimizer_x_lower($mini)],[0],'x_lower');
ok_similar([gsl_min_fminimizer_x_upper($mini)],[2*$M_PI],'x_upper');
ok_similar([gsl_min_fminimizer_x_minimum($mini)],[3],'x_minimum');
ok_similar([gsl_min_fminimizer_f_lower($mini)],[cos(0)],'f_lower');
ok_similar([gsl_min_fminimizer_f_upper($mini)],[cos(2*$M_PI)],'f_upper');
ok_similar([gsl_min_fminimizer_f_minimum($mini)],[cos(3)],'f_minimum');
}
sub GSL_MIN_TEST_INTERVAL : Tests {
my $self = shift;
my $mini = $self->{min};
my ($x_lower, $x_upper, $epsabs, $epsrel) = (0,1e-7, 1e-3,1e-5);
ok_status(gsl_min_test_interval ($x_lower, $x_upper, $epsabs, $epsrel),
$GSL_SUCCESS, 'gsl_min_test_interval'
);
($x_lower, $x_upper, $epsabs, $epsrel) = (0,1e-2, 1e-3,1e-5);
ok_status(gsl_min_test_interval ($x_lower, $x_upper, $epsabs, $epsrel),
$GSL_CONTINUE, 'gsl_min_test_interval'
);
}
sub GSL_MIN_ITERATE : Tests {
my $self = shift;
my $m = $self->{min};
ok_status(gsl_min_fminimizer_set_with_values($m,
sub { cos($_[0]) },
3, cos(3),
0, cos(0),
2*$M_PI, cos(2*$M_PI)
));
my $status = gsl_min_fminimizer_iterate($m);
ok_status($status);
lives_ok(sub{
my $min = gsl_min_fminimizer_x_minimum($m);
my $lower = gsl_min_fminimizer_x_lower($m);
my $upper = gsl_min_fminimizer_x_upper($m);
}, 'getting the min, lower and upper lives');
}
Test::Class->runtests;
|
