trunk/lib/Math/GSL.pm

   1 package Math::GSL;
   2
   3 use warnings;
   4 use Config;
   5 use Data::Dumper;
   6 use Test::More;
   7 use strict;
   8
   9 our $VERSION = 0.01;
  10
  11 =head1 NAME
  12
  13 Math::GSL - OO Perl interface to the  GNU Scientific Library (GSL)
  14
  15 =head1 VERSION
  16
  17 Version 0.01
  18
  19 =cut
  20
  21 =head1 SYNOPSIS
  22
  23
  24     use Math::GSL::Sf qw/ ... /;
  25     ...
  26
  27 =head1 EXPORT
  28
  29 Nothing is exported by default. This is alpha software, the API could change at any second.
  30
  31
  32 =head1 AUTHOR
  33
  34 Jonathan Leto, C<< <jonathan at leto.net> >>
  35
  36 =head1 BUGS
  37
  38 Please report any bugs or feature requests to C<bug-math-gsl at rt.cpan.org>, or through
  39 the web interface at L<http://rt.cpan.org/NoAuth/ReportBug.html?Queue=Math::GSL>.  I will be notified, and then you'll
  40 automatically be notified of progress on your bug as I make changes.
  41
  42
  43 =head1 SUPPORT
  44
  45 You can find documentation for this module with the perldoc command.
  46
  47     perldoc Math::GSL
  48
  49 or online at L<http://leto.net/code/Math-GSL/>
  50
  51
  52 You can also look for information at:
  53
  54 =over 4
  55
  56 =item * RT: CPAN's request tracker
  57
  58 L<http://rt.cpan.org/NoAuth/Bugs.html?Dist=Math::GSL>
  59
  60 =item * AnnoCPAN: Annotated CPAN documentation
  61
  62 L<http://annocpan.org/dist/Math::GSL>
  63
  64 =item * CPAN Ratings
  65
  66 L<http://cpanratings.perl.org/d/Math::GSL>
  67
  68 =item * Search CPAN
  69
  70 L<http://search.cpan.org/dist/Math::GSL>
  71
  72 =back
  73
  74
  75 =head1 ACKNOWLEDGEMENTS
  76
  77
  78 =head1 COPYRIGHT & LICENSE
  79
  80 Copyright 2008 Jonathan Leto, all rights reserved.
  81
  82 This program is free software; you can redistribute it and/or modify it
  83 under the same terms as Perl itself.
  84
  85 =cut
  86
  87 sub new
  88 {
  89     my ($self,$args) = @_;
  90     my $class = ref $self || $self || 'Math::GSL';
  91     my $this = { };
  92     bless $this, $class;
  93 }
  94
  95 sub subsystems
  96 {
  97     return qw/ CDF Errno Fit Machine PowInt Randist SF Types /;
  98 }
  99
 100 sub verify_results
 101 {
 102     my ($self,$results,$eps) = @_;
 103     $eps ||= 1e-8;
 104     while (my($k,$v)=each %$results){
 105         my $x = eval $k;
 106         print "got $x for $k\n" if defined $ENV{DEBUG};
 107         if($x =~ /nan/i){
 108                 ok( $v eq $x, "'$v'?='$x'" );
 109         } else {
 110             my $res = abs($x-$v);
 111             $@ ? ok(0)
 112                : ok( $res < $eps, "$k ?= $x, +- $res" );
 113         }
 114     }
 115 }
 116
 117 use constant MIN_DOUBLE => 2.2250738585072014e-308;
 118 use constant MAX_DOUBLE => 1.7976931348623157E+308;
 119 use constant MAX_FLOAT  => 3.40282347E+38;
 120 use constant MIN_FLOAT  => 1.175494351e-38;
 121
 122 sub _has_quads          { $Config{use64bitint} eq 'define' || ($Config{longsize} >= 8) }
 123 sub _has_long_doubles                 { $Config{d_longdbl}     eq 'define'             }
 124 sub _has_long_doubles_as_default      { $Config{uselongdouble} eq 'define'             }
 125 sub _has_long_doubles_same_as_doubles { $Config{doublesize}    == $Config{longdblsize} }
 126
 127 sub _assert_dies($;$)
 128 {
 129     my ($code,$msg) = @_;
 130     my $status = eval { &$code };
 131     print "status=||$status||\n\$\?=$?\n\$\!=$!\n" if 0;
 132     $@  ?  ok(1, $msg) : ok (0, join "\n", $@,  $msg );
 133 }
 134
 135 42;