1 package Math::GSL::Min::Test;
2 use base q{Test::Class};
4 use Test::More tests => 23;
5 use Math::GSL qw/:all/;
6 use Math::GSL::Min qw/:all/;
7 use Math::GSL::Test qw/:all/;
8 use Math::GSL::Errno qw/:all/;
9 use Math::GSL::Const qw/$M_PI/;
10 #use Devel::Trace qw/trace/;
15 BEGIN { gsl_set_error_handler_off(); }
19 sub make_fixture : Test(setup) {
21 $self->{min} = gsl_min_fminimizer_alloc($gsl_min_fminimizer_goldensection);
22 $self->{brent} = gsl_min_fminimizer_alloc($gsl_min_fminimizer_brent);
25 sub teardown : Test(teardown) {
27 #gsl_min_fminimizer_free($self->{min});
28 #gsl_min_fminimizer_free($self->{brent});
31 sub GSL_MIN_TYPES : Tests {
33 my $m = gsl_min_fminimizer_alloc($gsl_min_fminimizer_goldensection);
34 isa_ok($m, 'Math::GSL::Min');
36 my $n = gsl_min_fminimizer_alloc($gsl_min_fminimizer_brent);
37 isa_ok($n, 'Math::GSL::Min');
40 sub GSL_MIN_NAME : Tests {
42 cmp_ok( 'goldensection', 'eq', gsl_min_fminimizer_name($self->{min}) );
43 cmp_ok( 'brent', 'eq', gsl_min_fminimizer_name($self->{brent}) );
46 sub GSL_MIN_SET : Tests {
48 my $mini = $self->{min};
50 gsl_min_fminimizer_set($mini,
51 sub { cos($_[0]) }, 3, 0, 2*$M_PI
54 # These are the first guesses first the initial iteration
55 cmp_ok( $mini->{x_minimum}, '==', 3 );
56 cmp_ok( $mini->{x_lower}, '==', 0 );
57 isa_ok( $mini->{function} , 'Math::GSL::Min::gsl_function_struct');
58 ok_similar([$mini->{x_upper}], [2*$M_PI] );
59 ok_similar([$mini->{f_minimum}], [cos(3)] );
60 ok_similar([$mini->{f_lower}], [cos(0)] );
61 ok_similar([$mini->{f_upper}], [cos(2*$M_PI)] );
63 ok_similar([gsl_min_fminimizer_x_lower($mini)],[0],'x_lower');
64 ok_similar([gsl_min_fminimizer_x_upper($mini)],[2*$M_PI],'x_upper');
65 ok_similar([gsl_min_fminimizer_x_minimum($mini)],[3],'x_minimum');
67 ok_similar([gsl_min_fminimizer_f_lower($mini)],[cos(0)],'f_lower');
68 ok_similar([gsl_min_fminimizer_f_upper($mini)],[cos(2*$M_PI)],'f_upper');
69 ok_similar([gsl_min_fminimizer_f_minimum($mini)],[cos(3)],'f_minimum');
72 sub GSL_MIN_TEST_INTERVAL : Tests {
74 my $mini = $self->{min};
76 my ($x_lower, $x_upper, $epsabs, $epsrel) = (0,1e-7, 1e-3,1e-5);
77 ok_status(gsl_min_test_interval ($x_lower, $x_upper, $epsabs, $epsrel),
78 $GSL_SUCCESS, 'gsl_min_test_interval'
81 ($x_lower, $x_upper, $epsabs, $epsrel) = (0,1e-2, 1e-3,1e-5);
82 ok_status(gsl_min_test_interval ($x_lower, $x_upper, $epsabs, $epsrel),
83 $GSL_CONTINUE, 'gsl_min_test_interval'
87 sub GSL_MIN_ITERATE : Tests {
89 my $mini = $self->{min};
90 #warn Dumper [ $mini, $mini->{state} ];
91 ok_status(gsl_min_fminimizer_set_with_values($mini,
97 #warn Dumper [ $mini, $mini->{state} ];
99 #warn Dumper [map { $mini->{$_} } qw(x_minimum x_lower x_upper f_minimum f_lower f_upper) ];
100 local $TODO = 'iterate does not work';
102 #ok_status( gsl_min_fminimizer_iterate($mini));
106 sub GSL_MIN_NEW_FREE : Tests {
108 my $min = $self->{min};
109 isa_ok($min, 'Math::GSL::Min');
111 gsl_min_fminimizer_free($min);
112 ok(!$@, 'gsl_min_fminimizer_free');
115 Test::Class->runtests;