The Perl Toolchain Summit needs more sponsors. If your company depends on Perl, please support this very important event.


# Test Script for the PDL interface to the GSL library
#  This tests only that the interface is working, i.e. that the
#   functions can be called. The actual return values are not
#   checked. 
#  The GSL library already has a extensive test suite, and we
#  do not want to duplicate that effort here.

use PDL::LiteF;
use Test::More;

BEGIN
{
   use PDL::Config;
   if ( $PDL::Config{WITH_GSL} ) {
      eval " use PDL::GSLSF::BESSEL; ";
      unless ($@) {
         plan tests => 1;
      } else {
         plan skip_all => "PDL::GSLSF::BESSEL not installed.";
      }
   } else {
      plan skip_all => "PDL::GSLSF::BESSEL not compiled.";
   }
}

$arg = 5.0;
$expected = -0.17759677131433830434739701;

($y,$err) = gsl_sf_bessel_Jn $arg, 0;

print "got $y +- $err\n";

ok abs($y-$expected) < 1e-6;