t/FFT.t

   1 package Math::GSL::FFT::Test;
   2 use Math::GSL::Test qw/:all/;
   3 use base q{Test::Class};
   4 use Test::More;
   5 use Math::GSL::FFT qw/:all/;
   6 use Math::GSL qw/:all/;
   7 use Data::Dumper;
   8 use Math::GSL::Errno qw/:all/;
   9 use strict;
  10
  11 BEGIN { gsl_set_error_handler_off() }
  12
  13 sub make_fixture : Test(setup) {
  14 }
  15
  16 sub teardown : Test(teardown) {
  17     unlink 'fft' if -f 'fft';
  18 }
  19 sub FFT_REAL_RADIX2_TRANSFORM : Tests
  20 {
  21     my $data = [ (0) x 5, (1) x 22, (0) x 5 ];
  22     my ($status, $pass ) = gsl_fft_real_radix2_transform ($data, 1, 32);
  23     ok_status($status);
  24     ok_similar( $pass, [
  25             22, -8.44205264582682, -4.64465605976076, -0.643126602526688, 1.70710678118655, 1.8349201998544,
  26             0.572726230154202, -0.676964287646119, -1, -0.455944707054924, 0.255700894591988, 0.524240654352147,
  27             0.292893218813453, -0.059180002187481, -0.183771064985432, -0.0818926089645147, 0, -0.831469612302545,
  28             -0.923879532511287, -0.195090322016128, 0.707106781186547, 0.98078528040323, 0.38268343236509,
  29             -0.555570233019602, -1, -0.555570233019602, 0.38268343236509, 0.980785280403231, 0.707106781186547,
  30             -0.195090322016128, -0.923879532511287, -0.831469612302545
  31             ]);
  32 }
  33
  34 sub FFT_COMPLEX_RADIX2_FORWARD : Tests
  35 {
  36     my $data = [ (1) x 10, (0) x 236, (1) x 10 ];
  37     my ($status, $fft) = gsl_fft_complex_radix2_forward ($data, 1, 128);
  38     ok_status($status);
  39     local $TODO = 'fft of complex_packed_array does not work';
  40     # we should propably use the gsl_fft_real_unpack function here to create a suitable complex_packed_array
  41     ok( defined $fft );
  42 }
  43 sub FFT_VARS : Tests {
  44     cmp_ok( $gsl_fft_forward, '==', -1, 'gsl_fft_forward' );
  45     cmp_ok( $gsl_fft_backward, '==', +1, 'gsl_fft_backward' );
  46 }
  47
  48 sub WAVETABLE_ALLOC_FREE: Tests {
  49     my $wavetable = gsl_fft_complex_wavetable_alloc(42);
  50     isa_ok($wavetable, 'Math::GSL::FFT' );
  51     gsl_fft_complex_wavetable_free($wavetable);
  52     ok(!$@, 'gsl_fft_complex_wavetable_free');
  53
  54     $wavetable = gsl_fft_halfcomplex_wavetable_alloc(42);
  55     isa_ok($wavetable, 'Math::GSL::FFT' );
  56     gsl_fft_halfcomplex_wavetable_free($wavetable);
  57     ok(!$@, 'gsl_fft_halfcomplex_wavetable_free');
  58
  59     $wavetable = gsl_fft_real_wavetable_alloc(42);
  60     isa_ok($wavetable, 'Math::GSL::FFT' );
  61     gsl_fft_real_wavetable_free($wavetable);
  62     ok(!$@, 'gsl_fft_real_wavetable_free');
  63
  64 }
  65
  66 sub WORKSPACE_ALLOC_FREE: Tests {
  67     my $workspace = gsl_fft_complex_workspace_alloc(42);
  68     isa_ok($workspace, 'Math::GSL::FFT' );
  69     gsl_fft_complex_workspace_free($workspace);
  70     ok(!$@, 'gsl_fft_complex_workspace_free');
  71
  72     # there are no gsl_fft_halfcomplex_workspace_* functions
  73
  74     $workspace = gsl_fft_real_workspace_alloc(42);
  75     isa_ok($workspace, 'Math::GSL::FFT' );
  76     gsl_fft_real_workspace_free($workspace);
  77     ok(!$@, 'gsl_fft_real_workspace_free');
  78 }
  79 Test::Class->runtests;