IEEE math for the masses
[p5sagit/p5-mst-13.2.git] / t / op / exp.t
1 #!./perl
2
3 BEGIN {
4     chdir 't' if -d 't';
5     @INC = '../lib';
6     require './test.pl';
7 }
8
9 plan tests => 22;
10
11 # compile time evaluation
12
13 $s = sqrt(2);
14 is(substr($s,0,5), '1.414');
15
16 $s = exp(1);
17 is(substr($s,0,7), '2.71828');
18
19 ok(exp(log(1)) == 1);
20
21 # run time evaluation
22
23 $x1 = 1;
24 $x2 = 2;
25 $s = sqrt($x2);
26 is(substr($s,0,5), '1.414');
27
28 $s = exp($x1);
29 is(substr($s,0,7), '2.71828');
30
31 ok(exp(log($x1)) == 1);
32
33 # tests for transcendental functions
34
35 my $pi = 3.1415926535897931160;
36 my $pi_2 = 1.5707963267948965580;
37
38 sub round {
39    my $result = shift;
40    return sprintf("%.9f", $result);
41 }
42
43 # sin() tests
44 ok(sin(0) == 0.0);
45 ok(round(sin($pi)) == 0.0);
46 ok(round(sin(-1 * $pi)) == 0.0);
47 ok(round(sin($pi_2)) == 1.0);
48 ok(round(sin(-1 * $pi_2)) == -1.0);
49
50 # cos() tests
51 ok(cos(0) == 1.0);
52 ok(round(cos($pi)) == -1.0);
53 ok(round(cos(-1 * $pi)) == -1.0);
54 ok(round(cos($pi_2)) == 0.0);
55 ok(round(cos(-1 * $pi_2)) == 0.0);
56
57 # atan2() tests
58 ok(round(atan2(-0.0, 0.0)) == 0);
59 ok(round(atan2(0.0, 0.0)) == 0);
60 ok(round(atan2(-0.0, -0.0)) == round(-1 * $pi));
61 ok(round(atan2(0.0, -0.0)) == round($pi));
62 ok(round(atan2(-1.0, 0.0)) == round(-1 * $pi_2));
63 ok(round(atan2(1.0, 0.0)) == round($pi_2));