t/Poly.t

   1 use Test::More 'no_plan';
   2 use Math::GSL;
   3 use Math::GSL::Errno;
   4 use Math::GSL::Poly;
   5 use Data::Dumper;
   6 use strict;
   7 use warnings;
   8
   9
  10 {
  11     my $z = Math::GSL::Poly::gsl_complex->new;
  12
  13     ok( defined $z && $z->isa('Math::GSL::Poly'), 'gsl_complex->new returns correct object');
  14 }
  15
  16 {
  17     my $gsl = Math::GSL->new;
  18
  19     my $vals = Math::GSL::Poly::new_doubleArray(3);
  20
  21     #fugly
  22     Math::GSL::Poly::doubleArray_setitem($vals, 0, 3.14);
  23     Math::GSL::Poly::doubleArray_setitem($vals, 1, 2.72);
  24     Math::GSL::Poly::doubleArray_setitem($vals, 2, 5.55);
  25
  26     my $y = Math::GSL::Poly::gsl_poly_eval( $vals , 3, 1.0);
  27
  28     ok( abs($y - (3.14+2.72+5.55) ) < 1e-8, 'gsl_poly_eval' );
  29 }
  30