1 package Math::GSL::Roots::Test;
3 use base q{Test::Class};
4 use Test::More tests => 14;
5 use Math::GSL qw/:all/;
6 use Math::GSL::Roots qw/:all/;
7 use Math::GSL::Test qw/:all/;
8 use Math::GSL::Errno qw/:all/;
11 BEGIN { gsl_set_error_handler_off(); }
13 sub make_fixture : Test(setup) {
15 $self->{solver} = gsl_root_fsolver_alloc($gsl_root_fsolver_bisection);
16 $self->{fdfsolver} = gsl_root_fdfsolver_alloc($gsl_root_fdfsolver_newton);
19 sub teardown : Test(teardown) {
23 sub GSL_FDFSOLVER_BASIC : Tests {
25 local $TODO = 'this blows up';
26 #ok_status(gsl_root_fdfsolver_set($self->{fdfsolver},
27 # sub { my $x=shift; ($x-3.2)**3 },
32 sub GSL_ROOTS_ALLOC_FREE : Tests {
34 my $x = $self->{solver};
35 isa_ok($x, 'Math::GSL::Roots', 'gsl_root_fsolver_alloc' );
36 gsl_root_fsolver_free($x);
37 ok(!$@, 'gsl_root_fsolver_free');
40 sub GSL_ROOTS_SET : Tests {
42 my $solver = $self->{solver};
43 ok_status(gsl_root_fsolver_set($solver,
44 sub { my $x=shift; ($x-3.2)**3 },
47 ok_similar( [$solver->{root} ], [2.5], 'bisection starts of with midpoint as initial guess' );
50 sub GSL_ROOT_ITERATE : Tests {
52 my $solver = gsl_root_fsolver_alloc($gsl_root_fsolver_brent);
53 ok_status(gsl_root_fsolver_set($solver,
54 sub { my $x=shift; ($x-3.2)**3 },
57 # This currently blows up
58 #local $TODO = q{???};
59 #ok_status( gsl_root_fsolver_iterate($solver));
60 my $root = gsl_root_fsolver_root($solver);
61 ok_similar([$root], [2.5], 'gsl_root_fsolver_root');
64 sub SOlVER_TYPES : Tests {
65 cmp_ok( $gsl_root_fsolver_bisection->{name} ,'eq','bisection' );
66 cmp_ok( $gsl_root_fsolver_brent->{name} ,'eq','brent' );
67 cmp_ok( $gsl_root_fsolver_falsepos->{name} ,'eq','falsepos' );
68 cmp_ok( $gsl_root_fdfsolver_newton->{name} ,'eq','newton' );
69 cmp_ok( $gsl_root_fdfsolver_secant->{name} ,'eq','secant' );
70 cmp_ok( $gsl_root_fdfsolver_steffenson->{name},'eq','steffenson' );
73 sub GSL_ROOTFSOLVER_NAME : Tests {
75 cmp_ok( gsl_root_fsolver_name($self->{solver}) ,'eq','bisection' );
78 sub GSL_FDFSOLVER_NAME : Tests {
80 cmp_ok( gsl_root_fdfsolver_name($self->{fdfsolver}) ,'eq','newton' );
83 Test::Class->runtests;