t/Roots.t

   1 package Math::GSL::Roots::Test;
   2 use strict;
   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/;
   9 use Data::Dumper;
  10
  11 BEGIN { gsl_set_error_handler_off(); }
  12
  13 sub make_fixture : Test(setup) {
  14     my $self = shift;
  15     $self->{solver}    = gsl_root_fsolver_alloc($gsl_root_fsolver_bisection);
  16     $self->{fdfsolver} = gsl_root_fdfsolver_alloc($gsl_root_fdfsolver_newton);
  17 }
  18
  19 sub teardown : Test(teardown) {
  20     my $self = shift;
  21 }
  22
  23 sub GSL_FDFSOLVER_BASIC : Tests {
  24     my $self = shift;
  25     local $TODO = 'this blows up';
  26     #ok_status(gsl_root_fdfsolver_set($self->{fdfsolver},
  27     #    sub { my $x=shift; ($x-3.2)**3 },
  28     #    5
  29     #));
  30
  31 }
  32 sub GSL_ROOTS_ALLOC_FREE : Tests {
  33     my $self = shift;
  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');
  38 }
  39
  40 sub GSL_ROOTS_SET : Tests {
  41     my $self = shift;
  42     my $solver = $self->{solver};
  43     ok_status(gsl_root_fsolver_set($solver,
  44         sub { my $x=shift; ($x-3.2)**3 },
  45         0, 5
  46     ));
  47     ok_similar( [$solver->{root} ], [2.5], 'bisection starts of with midpoint as initial guess' );
  48 }
  49
  50 sub GSL_ROOT_ITERATE : Tests {
  51     my $self = shift;
  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 },
  55         0, 5
  56     ));
  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');
  62 }
  63
  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' );
  71 }
  72
  73 sub GSL_ROOTFSOLVER_NAME : Tests {
  74     my $self = shift;
  75     cmp_ok( gsl_root_fsolver_name($self->{solver})   ,'eq','bisection'  );
  76 }
  77
  78 sub GSL_FDFSOLVER_NAME : Tests {
  79     my $self = shift;
  80     cmp_ok( gsl_root_fdfsolver_name($self->{fdfsolver})   ,'eq','newton'  );
  81 }
  82
  83 Test::Class->runtests;