1 package Math
::GSL
::Complex
::Test
;
2 use Math
::GSL
::Test qw
/:all/;
3 use base
q{Test::Class};
5 use Math
::GSL
::Complex qw
/:all/;
6 use Math
::GSL
::Const qw
/:all/;
7 use Math
::GSL qw
/:all/;
11 sub make_fixture
: Test
(setup
) {
13 $self->{gsl_complex
} = Math
::GSL
::Complex
::gsl_complex
->new;
16 sub teardown
: Test
(teardown
) {
19 sub GSL_COMPLEX_NEW
: Test
{
21 my $x = $self->{gsl_complex
};
22 isa_ok
( $x, 'Math::GSL::Complex' );
25 sub GSL_COMPLEX_RECT
: Tests
{
26 my $x = gsl_complex_rect
(5,3);
27 ok_similar
( [ gsl_parts
($x) ], [ 5, 3 ], 'gsl_complex_rect' );
30 sub GSL_COMPLEX_POLAR
: Tests
{
31 my $x = gsl_complex_polar
(1.0,0.1);
32 my ($r,$theta) = ( gsl_complex_abs
($x), gsl_complex_arg
($x) );
34 print "x=$r * e^(i $theta)*I\n" if $ENV{DEBUG
};
35 ok
( is_similar
($r,1.0), 'gsl_complex_polar r, res=' .($r-1.) . "\n");
36 ok
( is_similar
($theta,0.1), 'gsl_complex_polar theta, res=' .($theta-0.1) . "\n");
39 sub GSL_COMPLEX_ABS
: Tests
{
40 my $x = gsl_complex_rect
(5,3);
41 my $abs = gsl_complex_abs
($x);
42 ok
( is_similar
($abs, sqrt(34)), 'gsl_complex_abs');
45 sub GSL_COMPLEX_ABS2
: Tests
{
46 my $x = gsl_complex_rect
(5,3);
47 my $abs = gsl_complex_abs2
($x);
48 ok
( is_similar
($abs, 34), 'gsl_complex_abs2' );
51 sub GSL_COMPLEX_LOGABS
: Tests
{
52 my $x = gsl_complex_rect
(5,3);
54 my $abs = gsl_complex_logabs
($x);
55 ok
( is_similar
($abs, log(sqrt(34)) ), 'gsl_complex_logabs' );
58 sub GSL_COMPLEX_LOG
: Tests
{
59 my $x = gsl_complex_rect
($M_E,0);
60 $x = gsl_complex_log
($x);
61 ok_similar
([ gsl_parts
($x) ] ,
66 sub GSL_COMPLEX_LOG10
: Tests
{
67 my $x = gsl_complex_rect
(10,0);
68 $x = gsl_complex_log10
($x);
69 ok_similar
([ gsl_parts
($x) ] ,
75 sub GSL_COMPLEX_LOG_B
: Tests
{
76 my $z = gsl_complex_rect
(5.0,0);
77 my $base = gsl_complex_rect
(5.0,0);
78 $z = gsl_complex_log_b
($z, $base);
79 ok_similar
([ gsl_parts
($z) ] ,
85 sub GSL_COMPLEX_EXP
: Tests
{
86 my $z = gsl_complex_rect
(2, 5);
87 $z = gsl_complex_exp
($z);
88 ok_similar
( [ gsl_parts
($z) ],
89 [ $M_E**2*cos(5), $M_E**2*sin(5) ]
93 sub GSL_COMPLEX_EXP_EULERS_IDENTITY
: Tests
{ # e^(i pi) + 1 = 0
94 my $z = gsl_complex_rect
(0, $M_PI);
95 $z = gsl_complex_exp
($z);
96 $z = gsl_complex_add_real
($z, 1);
97 ok_similar
( [ gsl_parts
($z) ],
102 sub GSL_COMPLEX_ADD
: Tests
{
103 my $x = gsl_complex_rect
(5,3);
104 my $y = gsl_complex_rect
(1,2);
106 my $z = gsl_complex_add
($x, $y);
107 ok_similar
( [ gsl_parts
($z) ], [ 6, 5 ], 'gsl_complex_add' );
110 sub GSL_COMPLEX_SUB
: Tests
{
111 my $x = gsl_complex_rect
(5,3);
112 my $y = gsl_complex_rect
(1,2);
114 my $z = gsl_complex_sub
($x, $y);
115 ok_similar
( [ gsl_parts
($z) ], [ 4, 1 ], 'gsl_complex_sub' );
118 sub GSL_COMPLEX_MUL
: Tests
{
119 my $x = gsl_complex_rect
(5,3);
120 my $y = gsl_complex_rect
(1,2);
122 my $z = gsl_complex_mul
($x, $y);
123 ok_similar
( [ gsl_parts
($z) ], [ -1, 13 ], 'gsl_complex_mul' );
126 sub GSL_COMPLEX_DIV
: Tests
{
127 my $x = gsl_complex_rect
(4,5);
128 my $y = gsl_complex_rect
(2,1);
130 my $z = gsl_complex_div
($x, $y);
131 ok_similar
( [ gsl_parts
($z) ], [ 13/5, 6/5 ], 'gsl_complex_div' );
133 sub GSL_COMPLEX_EQ
: Tests
{
134 my $x = gsl_complex_rect
(4,5);
135 my $y = gsl_complex_rect
(2,1);
136 ok
(! gsl_complex_eq
($x,$y), 'gsl_complex_eq' );
137 $y = gsl_complex_rect
(4,5);
138 ok
( gsl_complex_eq
($x,$y), 'gsl_complex_eq' );
140 sub GSL_COMPLEX_ADD_REAL
: Tests
{
141 my $x = gsl_complex_rect
(5,3);
143 my $z = gsl_complex_add_real
($x, 1);
144 ok
( gsl_real
($z) == 6, 'gsl_complex_add_real');
147 sub GSL_COMPLEX_SUB_REAL
: Tests
{
148 my $x = gsl_complex_rect
(5,3);
150 my $z = gsl_complex_sub_real
($x, 1);
151 ok
( gsl_real
($z) == 4, 'gsl_complex_sub_real');
154 sub GSL_COMPLEX_MUL_REAL
: Tests
{
155 my $x = gsl_complex_rect
(5,3);
157 my $z = gsl_complex_mul_real
($x, 2);
158 ok
( gsl_real
($z) == 10, 'gsl_complex_mul_real');
161 sub GSL_COMPLEX_DIV_REAL
: Tests
{
162 my $x = gsl_complex_rect
(6,3);
164 my $z = gsl_complex_div_real
($x, 2);
165 ok
( gsl_real
($z) == 3, 'gsl_complex_div_real');
168 sub GSL_COMPLEX_ADD_IMAG
: Tests
{
169 my $x = gsl_complex_rect
(6,3);
171 my $z = gsl_complex_add_imag
($x, 2);
172 ok
( gsl_imag
($z) == 5, 'gsl_complex_add_imag');
175 sub GSL_COMPLEX_SUB_IMAG
: Tests
{
176 my $x = gsl_complex_rect
(6,3);
178 my $z = gsl_complex_sub_imag
($x, 2);
179 ok
( gsl_imag
($z) == 1, 'gsl_complex_sub_imag');
182 sub GSL_COMPLEX_MUL_IMAG
: Tests
{
183 my $x = gsl_complex_rect
(6,3);
185 my $z = gsl_complex_mul_imag
($x, 2);
186 ok_similar
( [ gsl_parts
($z) ], [ -6, 12 ], 'gsl_complex_mul_imag' );
189 sub GSL_COMPLEX_DIV_IMAG
: Tests
{
190 my $x = gsl_complex_rect
(6,4);
192 my $z = gsl_complex_div_imag
($x, 2);
193 ok_similar
( [ gsl_parts
($z) ], [ 2, -3 ], 'gsl_complex_div_imag' );
196 sub GSL_COMPLEX_CONJUGATE
: Tests
{
197 my $x = gsl_complex_rect
(6,4);
199 my $z = gsl_complex_conjugate
($x);
200 ok_similar
( [ gsl_parts
($z) ], [6, -4], 'gsl_complex_conjugate');
204 sub GSL_COMPLEX_NEGATIVE
: Tests
{
205 my $x = gsl_complex_rect
(6,4);
207 my $z = gsl_complex_negative
($x);
208 ok_similar
( [ gsl_parts
($z) ], [-6, -4], 'gsl_complex_negative');
211 sub GSL_COMPLEX_SQRT
: Tests
{
212 my $x = gsl_complex_rect
(-7,24);
214 my $z = gsl_complex_sqrt
($x);
215 ok_similar
( [ gsl_parts
($z) ], [ 3 , 4 ] ,'gsl_complex_sqrt' );
218 sub GSL_COMPLEX_SQRT_REAL
: Tests
{
219 my $x = gsl_complex_rect
(9,4);
221 my $z = gsl_complex_sqrt_real
(-4);
222 ok
( gsl_imag
($z) == 2, 'gsl_complex_sqrt_real');
225 sub GSL_COMPLEX_POW
: Tests
{
226 my $x = gsl_complex_rect
(0,1);
227 my $y = gsl_complex_rect
(2,0);
228 my $z = gsl_complex_pow
($x, $y);
230 ok_similar
( [ gsl_parts
($z) ], [-1, 0 ], 'gsl_complex_pow' );
232 $x = gsl_complex_rect
(3,4);
233 $z = gsl_complex_pow
($x, $y);
235 ok_similar
( [gsl_parts
($z) ], [ 3**2 - 4**2, 2*3*4 ], 'gsl_complex_pow' );
238 sub GSL_COMPLEX_SET_REAL
: Tests
{
239 my $x = gsl_complex_rect
(3,4);
241 ok_similar
( [ gsl_parts
($x) ], [ 5, 4 ], 'gsl_complex_set_real');
244 sub GSL_COMPLEX_SET_IMAG
: Tests
{
245 my $x = gsl_complex_rect
(3,4);
247 ok_similar
( [ gsl_parts
($x) ], [ 3, 5 ], 'gsl_complex_set_imag');
250 sub GSL_COMPLEX_SET_COMPLEX
: Tests
{
251 my $x = gsl_complex_rect
(3,4);
252 gsl_set_complex
($x, 5, 3);
253 ok_similar
( [ gsl_parts
($x) ], [ 5, 3 ], 'gsl_complex_set_complex');
256 sub GSL_COMPLEX_SIN
: Tests
{
257 my $x = gsl_complex_rect
(3,2);
258 my $z = gsl_complex_sin
($x);
260 ok_similar
( [ gsl_parts
($z) ],
261 [ 0.53092108624852, -3.59056458998578 ],
266 sub GSL_COMPLEX_ARCSIN
: Tests
{
267 my $x = gsl_complex_rect
(1,0);
268 my $z = gsl_complex_arcsin
($x);
270 ok_similar
( [ gsl_parts
($z) ],
275 sub GSL_COMPLEX_SINH
: Tests
{
276 my $x = gsl_complex_rect
(0,0);
277 my $z = gsl_complex_sinh
($x);
279 ok_similar
( [ gsl_parts
($z) ],
285 sub GSL_COMPLEX_ARCSINH
: Tests
{
286 my $x = gsl_complex_rect
(0,0);
287 my $z = gsl_complex_arcsinh
($x);
289 ok_similar
( [ gsl_parts
($z) ],
291 'gsl_complex_arcsinh'
294 sub GSL_COMPLEX_COS
: Tests
{
295 my $x = gsl_complex_rect
(3,2);
296 my $z = gsl_complex_cos
($x);
297 ok_similar
( [ gsl_parts
($z) ],
298 [ -3.72454550491532, -0.511822569987385],
303 sub GSL_COMPLEX_ARCCOS
: Tests
{
304 my $x = gsl_complex_rect
(0,0);
305 my $z = gsl_complex_arccos
($x);
306 ok_similar
( [ gsl_parts
($z)],
311 sub GSL_COMPLEX_COSH
: Tests
{
312 my $x = gsl_complex_rect
(0,0);
313 my $z = gsl_complex_cosh
($x);
314 ok_similar
( [ gsl_parts
($z) ],
320 sub GSL_COMPLEX_ARCCOSH
: Tests
{
321 my $x = gsl_complex_rect
(1,0);
322 my $z = gsl_complex_arccosh
($x);
323 ok_similar
( [ gsl_parts
($z) ],
325 'gsl_complex_arccosh'
329 sub GSL_COMPLEX_TAN
: Tests
{
330 my $x = gsl_complex_rect
(3,2);
331 my $z = gsl_complex_tan
($x);
332 ok_similar
( [ gsl_parts
($z) ],
333 [ -0.0098843750383225, 0.965385879022133 ],
338 sub GSL_COMPLEX_ARCTAN
: Tests
{
339 my $x = gsl_complex_rect
(1,0);
340 my $z = gsl_complex_arctan
($x);
341 ok_similar
( [ gsl_parts
($z) ],
347 sub GSL_COMPLEX_TANH
: Tests
{
348 my $x = gsl_complex_rect
(0,0);
349 my $z = gsl_complex_tan
($x);
350 ok_similar
( [ gsl_parts
($z) ],
356 sub GSL_COMPLEX_ARCTANH
: Tests
{
357 my $x = gsl_complex_rect
(0,0);
358 my $z = gsl_complex_arctan
($x);
359 ok_similar
( [ gsl_parts
($z) ],
361 'gsl_complex_arctanh'
364 sub GSL_COMPLEX_SEC
: Tests
{
365 my $x = gsl_complex_rect
(3,2);
366 my $z = gsl_complex_sec
($x);
367 ok_similar
( [ gsl_parts
($z) ],
368 [ -0.263512975158389, 0.0362116365587685 ],
373 sub GSL_COMPLEX_ARCSEC
: Tests
{
374 my $x = gsl_complex_rect
(0,0);
375 my $z = gsl_complex_arcsec
($x);
376 ok_similar
( [ gsl_parts
($z) ],
381 sub GSL_COMPLEX_SECH
: Tests
{
382 my $x = gsl_complex_rect
(1,0);
383 my $z = gsl_complex_sech
($x);
384 ok_similar
( [ gsl_parts
($z) ],
385 [ 2/(exp(1)+exp(-1)), 0.0 ],
390 sub GSL_COMPLEX_ARCSECH
: Tests
{
391 my $x = gsl_complex_rect
(2/(exp(1)+exp(-1)),0);
392 my $z = gsl_complex_arcsech
($x);
393 ok_similar
( [ gsl_parts
($z) ],
395 'gsl_complex_arcsech'
399 sub GSL_COMPLEX_CSC
: Tests
{
400 my $x = gsl_complex_rect
(3,2);
401 my $z = gsl_complex_csc
($x);
402 ok_similar
( [ gsl_parts
($z) ],
403 [ 0.0403005788568915, 0.27254866146294],
408 sub GSL_COMPLEX_ARCCSC
: Tests
{
409 my $x = gsl_complex_rect
(0,0);
410 my $z = gsl_complex_arccsc
($x);
411 ok_similar
( [ gsl_parts
($z)],
417 sub GSL_COMPLEX_CSCH
: Tests
{
418 my $x = gsl_complex_rect
(0,1);
419 my $z = gsl_complex_csch
($x);
420 ok_similar
( [ gsl_parts
($z) ],
426 sub GSL_COMPLEX_ARCCSCH
: Tests
{
427 my $x = gsl_complex_rect
(0,0);
428 my $z = gsl_complex_arccsch
($x);
429 ok_similar
( [ gsl_parts
($z) ],
431 'gsl_complex_arccsch'
435 sub GSL_COMPLEX_COT
: Tests
{
436 my $x = gsl_complex_rect
(3,2);
437 my $z = gsl_complex_cot
($x);
438 ok_similar
( [ gsl_parts
($z) ],
439 [ -0.0106047834703371, -1.035746637765],
444 sub GSL_COMPLEX_ARCCOT
: Tests
{
445 my $x = gsl_complex_rect
(1,0);
446 my $z = gsl_complex_arccot
($x);
447 ok_similar
( [ gsl_parts
($z) ],
453 sub GSL_COMPLEX_COTH
: Tests
{
454 my $x = gsl_complex_rect
(0,1);
455 my $z = gsl_complex_coth
($x);
456 ok_similar
( [ gsl_parts
($z) ],
457 [ 0.0, -cos(1)/sin(1) ],
462 sub GSL_COMPLEX_ARCCOTH
: Tests
{
463 my $x = gsl_complex_rect
(0,0);
464 my $z = gsl_complex_arccoth
($x);
465 ok_similar
( [ gsl_parts
($z) ],
467 'gsl_complex_arccoth'