lib/Math/GSL/Complex/Test.pm

   1 package Math::GSL::Complex::Test;
   2 use base q{Test::Class};
   3 use Test::More;
   4 use Math::GSL::Complex qw/:all/;
   5 use Math::GSL qw/:all/;
   6 use Data::Dumper;
   7 use strict;
   8
   9 sub make_fixture : Test(setup) {
  10     my $self = shift;
  11     $self->{gsl_complex} = Math::GSL::Complex::gsl_complex->new;
  12 }
  13
  14 sub teardown : Test(teardown) {
  15 }
  16
  17 sub GSL_COMPLEX_NEW : Test {
  18     my $self = shift;
  19     my $x = $self->{gsl_complex};
  20     isa_ok( $x, 'Math::GSL::Complex' );
  21 }
  22
  23 sub GSL_COMPLEX_RECT : Tests {
  24     my $x = gsl_complex_rect(5,3);
  25     ok_similar( [ gsl_parts($x) ], [ 5, 3 ], 'gsl_complex_rect' );
  26 }
  27
  28 sub GSL_COMPLEX_POLAR : Tests {
  29     my $x           = gsl_complex_polar(1.0,0.1);
  30     my ($r,$theta)  = ( gsl_complex_abs($x), gsl_complex_arg($x) );
  31
  32     print "x=$r * e^(i $theta)*I\n" if $ENV{DEBUG};
  33     ok( is_similar($r,1.0),     'gsl_complex_polar r, res=' .($r-1.) . "\n");
  34     ok( is_similar($theta,0.1), 'gsl_complex_polar theta, res=' .($theta-0.1) . "\n");
  35 }
  36
  37 sub GSL_COMPLEX_ABS : Tests {
  38     my $x = gsl_complex_rect(5,3);
  39     my $abs = gsl_complex_abs($x);
  40     ok ( is_similar($abs, sqrt(34)), 'gsl_complex_abs');
  41 }
  42
  43 sub GSL_COMPLEX_ABS2 : Tests {
  44     my $x = gsl_complex_rect(5,3);
  45     my $abs = gsl_complex_abs2($x);
  46     ok( is_similar($abs, 34), 'gsl_complex_abs2' );
  47 }
  48
  49 sub GSL_COMPLEX_LOGABS : Tests {
  50     my $x = gsl_complex_rect(5,3);
  51
  52     my $abs = gsl_complex_logabs($x);
  53     ok ( is_similar($abs, log(sqrt(34)) ), 'gsl_complex_logabs' );
  54 }
  55
  56 sub GSL_COMPLEX_ADD : Tests {
  57     my $x = gsl_complex_rect(5,3);
  58     my $y = gsl_complex_rect(1,2);
  59
  60     my $z = gsl_complex_add($x, $y);
  61     ok_similar( [ gsl_parts($z) ], [ 6, 5 ], 'gsl_complex_add' );
  62 }
  63
  64 sub GSL_COMPLEX_SUB : Tests {
  65     my $x = gsl_complex_rect(5,3);
  66     my $y = gsl_complex_rect(1,2);
  67
  68     my $z = gsl_complex_sub($x, $y);
  69     ok_similar( [ gsl_parts($z) ], [ 4, 1 ], 'gsl_complex_sub' );
  70 }
  71
  72 sub GSL_COMPLEX_MUL : Tests {
  73     my $x = gsl_complex_rect(5,3);
  74     my $y = gsl_complex_rect(1,2);
  75
  76     my $z = gsl_complex_mul($x, $y);
  77     ok_similar( [ gsl_parts($z) ], [ -1, 13 ], 'gsl_complex_mul' );
  78 }
  79
  80 sub GSL_COMPLEX_DIV : Tests {
  81     my $x = gsl_complex_rect(4,5);
  82     my $y = gsl_complex_rect(2,1);
  83
  84     my $z = gsl_complex_div($x, $y);
  85     ok_similar( [ gsl_parts($z) ], [ 13/5, 6/5 ], 'gsl_complex_div' );
  86 }
  87 sub GSL_COMPLEX_EQ : Tests {
  88     my $x = gsl_complex_rect(4,5);
  89     my $y = gsl_complex_rect(2,1);
  90     ok(! gsl_complex_eq($x,$y), 'gsl_complex_eq' );
  91     $y    = gsl_complex_rect(4,5);
  92     ok( gsl_complex_eq($x,$y), 'gsl_complex_eq' );
  93 }
  94 sub GSL_COMPLEX_ADD_REAL : Tests {
  95     my $x = gsl_complex_rect(5,3);
  96
  97     my $z = gsl_complex_add_real($x, 1);
  98     ok( gsl_real($z) ==  6, 'gsl_complex_add_real');
  99 }
 100
 101 sub GSL_COMPLEX_SUB_REAL : Tests {
 102     my $x = gsl_complex_rect(5,3);
 103
 104     my $z = gsl_complex_sub_real($x, 1);
 105     ok( gsl_real($z) == 4, 'gsl_complex_sub_real');
 106 }
 107
 108 sub GSL_COMPLEX_MUL_REAL : Tests {
 109     my $x = gsl_complex_rect(5,3);
 110
 111     my $z = gsl_complex_mul_real($x, 2);
 112     ok( gsl_real($z) == 10, 'gsl_complex_mul_real');
 113 }
 114
 115 sub GSL_COMPLEX_DIV_REAL : Tests {
 116     my $x = gsl_complex_rect(6,3);
 117
 118     my $z = gsl_complex_div_real($x, 2);
 119     ok( gsl_real($z) == 3, 'gsl_complex_div_real');
 120 }
 121
 122 sub GSL_COMPLEX_ADD_IMAG : Tests  {
 123     my $x = gsl_complex_rect(6,3);
 124
 125     my $z = gsl_complex_add_imag($x, 2);
 126     ok( gsl_imag($z) == 5, 'gsl_complex_add_imag');
 127 }
 128
 129 sub GSL_COMPLEX_SUB_IMAG : Tests {
 130     my $x = gsl_complex_rect(6,3);
 131
 132     my $z = gsl_complex_sub_imag($x, 2);
 133     ok( gsl_imag($z) == 1, 'gsl_complex_sub_imag');
 134 }
 135
 136 sub GSL_COMPLEX_MUL_IMAG : Tests {
 137     my $x = gsl_complex_rect(6,3);
 138
 139     my $z = gsl_complex_mul_imag($x, 2);
 140     ok_similar( [ gsl_parts($z) ], [ -6, 12 ], 'gsl_complex_mul_imag' );
 141 }
 142
 143 sub GSL_COMPLEX_DIV_IMAG : Tests {
 144     my $x = gsl_complex_rect(6,4);
 145
 146     my $z = gsl_complex_div_imag($x, 2);
 147     ok_similar( [ gsl_parts($z) ], [ 2, -3 ], 'gsl_complex_div_imag' );
 148 }
 149
 150 sub GSL_COMPLEX_CONJUGATE : Tests {
 151     my $x = gsl_complex_rect(6,4);
 152
 153     my $z = gsl_complex_conjugate($x);
 154     ok_similar( [ gsl_parts($z) ], [6, -4], 'gsl_complex_conjugate');
 155 }
 156
 157
 158 sub GSL_COMPLEX_NEGATIVE : Tests {
 159     my $x = gsl_complex_rect(6,4);
 160
 161     my $z = gsl_complex_negative($x);
 162     ok_similar( [ gsl_parts($z) ], [-6, -4], 'gsl_complex_negative');
 163 }
 164
 165 sub GSL_COMPLEX_SQRT : Tests {
 166     my $x = gsl_complex_rect(-7,24);
 167
 168     my $z = gsl_complex_sqrt($x);
 169     ok_similar( [ gsl_parts($z)  ], [ 3  , 4 ] ,'gsl_complex_sqrt'  );
 170 }
 171
 172 sub GSL_COMPLEX_SQRT_REAL : Tests {
 173     my $x = gsl_complex_rect(9,4);
 174
 175     my $z = gsl_complex_sqrt_real(-4);
 176     ok( gsl_imag($z) == 2, 'gsl_complex_sqrt_real');
 177 }
 178
 179 sub GSL_COMPLEX_POW : Tests {
 180     my $x = gsl_complex_rect(0,1);
 181     my $y = gsl_complex_rect(2,0);
 182     my $z = gsl_complex_pow($x, $y);
 183
 184     ok_similar( [ gsl_parts($z) ], [-1, 0 ], 'gsl_complex_pow' );
 185
 186     $x = gsl_complex_rect(3,4);
 187     $z = gsl_complex_pow($x, $y);
 188
 189     ok_similar( [gsl_parts($z) ], [ 3**2 - 4**2, 2*3*4 ], 'gsl_complex_pow' );
 190 }
 191
 192 sub GSL_COMPLEX_SET_REAL : Tests {
 193     my $x = gsl_complex_rect(3,4);
 194     gsl_set_real($x, 5);
 195     ok_similar(  [ gsl_parts($x) ], [ 5, 4 ], 'gsl_complex_set_real');
 196 }
 197
 198 sub GSL_COMPLEX_SET_IMAG : Tests {
 199     my $x = gsl_complex_rect(3,4);
 200     gsl_set_imag($x, 5);
 201     ok_similar(  [ gsl_parts($x) ], [ 3, 5 ], 'gsl_complex_set_imag');
 202 }
 203
 204 sub GSL_COMPLEX_SET_COMPLEX : Tests {
 205     my $x = gsl_complex_rect(3,4);
 206     gsl_set_complex($x, 5, 3);
 207     ok_similar(  [ gsl_parts($x) ], [ 5, 3 ], 'gsl_complex_set_complex');
 208 }
 209
 210 sub GSL_COMPLEX_SIN : Tests {
 211     my $x = gsl_complex_rect(3,2);
 212     my $z = gsl_complex_sin($x);
 213
 214     ok_similar( [ gsl_parts($z)                       ],
 215                 [ 0.53092108624852, -3.59056458998578 ],
 216                 'gsl_complex_sin'
 217     );
 218 }
 219
 220 sub GSL_COMPLEX_SINH : Tests {
 221     my $x = gsl_complex_rect(0,0);
 222     my $z = gsl_complex_sinh($x);
 223
 224     ok_similar( [ gsl_parts($z) ],
 225                 [ 0.0,      0.0 ],
 226                 'gsl_complex_sinh'
 227     );
 228 }
 229 sub GSL_COMPLEX_COS : Tests {
 230     my $x = gsl_complex_rect(3,2);
 231     my $z = gsl_complex_cos($x);
 232     ok_similar( [ gsl_parts($z)                       ],
 233                 [ -3.72454550491532, -0.511822569987385],
 234                 'gsl_complex_cos'
 235     );
 236 }
 237
 238 sub GSL_COMPLEX_COSH : Tests {
 239     my $x = gsl_complex_rect(0,0);
 240     my $z = gsl_complex_cosh($x);
 241     ok_similar( [ gsl_parts($z)  ],
 242                 [ 1, 0.0         ],
 243                 'gsl_complex_cosh'
 244     );
 245 }
 246 sub GSL_COMPLEX_TAN : Tests {
 247     my $x = gsl_complex_rect(3,2);
 248     my $z = gsl_complex_tan($x);
 249     ok_similar( [ gsl_parts($z)                       ],
 250                 [ -0.0098843750383225, 0.965385879022133 ],
 251                 'gsl_complex_tan'
 252     );
 253 }
 254
 255 sub GSL_COMPLEX_TANH : Tests {
 256     my $x = gsl_complex_rect(0,0);
 257     my $z = gsl_complex_tan($x);
 258     ok_similar( [ gsl_parts($z) ],
 259                 [ 0.0,    0.0   ],
 260                 'gsl_complex_tanh'
 261     );
 262 }
 263 sub GSL_COMPLEX_SEC : Tests {
 264     my $x = gsl_complex_rect(3,2);
 265     my $z = gsl_complex_sec($x);
 266     ok_similar( [ gsl_parts($z)                          ],
 267                 [ -0.263512975158389, 0.0362116365587685 ],
 268                 'gsl_complex_sec'
 269     );
 270 }
 271
 272 sub GSL_COMPLEX_SECH : Tests {
 273     my $x = gsl_complex_rect(1,0);
 274     my $z = gsl_complex_sech($x);
 275     ok_similar( [ gsl_parts($z)             ],
 276                 [ 2/(exp(1)+exp(-1)), 0.0 ],
 277                 'gsl_complex_sech'
 278     );
 279 }
 280
 281 sub GSL_COMPLEX_CSC : Tests {
 282     my $x = gsl_complex_rect(3,2);
 283     my $z = gsl_complex_csc($x);
 284     ok_similar( [ gsl_parts($z)                       ],
 285                 [ 0.0403005788568915, 0.27254866146294],
 286                 'gsl_complex_csc'
 287     );
 288 }
 289
 290 sub GSL_COMPLEX_CSCH : Tests {
 291     my $x = gsl_complex_rect(0,1);
 292     my $z = gsl_complex_csch($x);
 293     ok_similar( [ gsl_parts($z)  ],
 294                 [ 0.0 , -1/sin(1)],
 295                 'gsl_complex_csch'
 296     );
 297 }
 298
 299 sub GSL_COMPLEX_COT : Tests {
 300     my $x = gsl_complex_rect(3,2);
 301     my $z = gsl_complex_cot($x);
 302     ok_similar( [ gsl_parts($z)                       ],
 303                 [ -0.0106047834703371, -1.035746637765],
 304                 'gsl_complex_cot'
 305     );
 306 }
 307
 308 sub GSL_COMPLEX_COTH : Tests {
 309     my $x = gsl_complex_rect(0,1);
 310     my $z = gsl_complex_coth($x);
 311     ok_similar( [ gsl_parts($z)       ],
 312                 [ 0.0, -cos(1)/sin(1) ],
 313                 'gsl_complex_coth'
 314     );
 315 }
 316
 317 42;