Jakub Wilk is now a Perl author.
[perl.git] / t / op / exp.t
1 #!./perl
2
3 # Simple tests for the basic math functions.
4
5 BEGIN {
6     chdir 't' if -d 't';
7     require './test.pl';
8     set_up_inc('../lib');
9 }
10
11 use Config;
12
13 plan tests => 31;
14
15 # compile time evaluation
16
17 eval { $s = sqrt(-1) }; # Kind of compile time.
18 like($@, qr/sqrt of -1/, 'compile time sqrt(-1) fails');
19
20 $s = sqrt(0);
21 is($s, 0, 'compile time sqrt(0)');
22
23 $s = sqrt(1);
24 is($s, 1, 'compile time sqrt(1)');
25
26 $s = sqrt(2);
27 is(substr($s,0,5), '1.414', 'compile time sqrt(2) == 1.414');
28
29 $s = exp(0);
30 is($s, 1, 'compile time exp(0) == 1');
31
32 $s = exp(1);
33 is(substr($s,0,7), '2.71828', 'compile time exp(1) == e');
34
35 eval { $s = log(0) };  # Kind of compile time.
36 like($@, qr/log of 0/, 'compile time log(0) fails');
37
38 $s = log(1);
39 is($s, 0, 'compile time log(1) == 0');
40
41 $s = log(2);
42 is(substr($s,0,5), '0.693', 'compile time log(2)');
43
44 cmp_ok(exp(log(1)), '==', 1, 'compile time exp(log(1)) == 1');
45
46 cmp_ok(round(atan2(1, 2)), '==', '0.463647609', "atan2(1, 2)");
47
48 # run time evaluation
49
50 $x0 = 0;
51 $x1 = 1;
52 $x2 = 2;
53
54 eval { $s = sqrt(-$x1) };
55 like($@, qr/sqrt of -1/, 'run time sqrt(-1) fails');
56
57 $s = sqrt($x0);
58 is($s, 0, 'run time sqrt(0)');
59
60 $s = sqrt($x1);
61 is($s, 1, 'run time sqrt(1)');
62
63 $s = sqrt($x2);
64 is(substr($s,0,5), '1.414', 'run time sqrt(2) == 1.414');
65
66 $s = exp($x0);
67 is($s, 1, 'run time exp(0) = 1');
68
69 $s = exp($x1);
70 is(substr($s,0,7), '2.71828', 'run time exp(1) = e');
71
72 eval { $s = log($x0) };
73 like($@, qr/log of 0/, 'run time log(0) fails');
74
75 $s = log($x1);
76 is($s, 0, 'compile time log(1) == 0');
77
78 $s = log($x2);
79 is(substr($s,0,5), '0.693', 'run time log(2)');
80
81 cmp_ok(exp(log($x1)), '==', 1, 'run time exp(log(1)) == 1');
82
83 # NOTE: do NOT test the trigonometric functions at [+-]Pi
84 # and expect to get exact results like 0, 1, -1, because
85 # you may not be able to feed them exactly [+-]Pi given
86 # all the variations of different long doubles.
87
88 my $pi_2 = 1.5707963267949;
89
90 sub round {
91    my $result = shift;
92    return sprintf("%.9f", $result);
93 }
94
95 # sin() tests
96 cmp_ok(sin(0), '==', 0.0, 'sin(0) == 0');
97 cmp_ok(abs(sin($pi_2) - 1), '<', 1e-9, 'sin(pi/2) == 1');
98 cmp_ok(abs(sin(-1 * $pi_2) - -1), '<', 1e-9, 'sin(-pi/2) == -1');
99
100 cmp_ok(round(sin($x1)), '==', '0.841470985', "sin(1)");
101
102 # cos() tests
103 cmp_ok(cos(0), '==', 1.0, 'cos(0) == 1');
104 cmp_ok(abs(cos($pi_2)), '<', 1e-9, 'cos(pi/2) == 0');
105 cmp_ok(abs(cos(-1 * $pi_2)), '<', 1e-9, 'cos(-pi/2) == 0');
106
107 cmp_ok(round(cos($x1)), '==', '0.540302306', "cos(1)");
108
109 cmp_ok(round(atan2($x1, $x2)), '==', '0.463647609', "atan2($x1, $x2)");
110
111 # atan2() tests testing with -0.0, 0.0, -1.0, 1.0 were removed due to
112 # differing results from calls to atan2() on various OS's and
113 # architectures.  See perlport.pod for more information.
114
115 SKIP: {
116     unless ($Config{usequadmath}) {
117         skip "need usequadmath", 1;
118     }
119     # For quadmath we have a known precision.  
120     is(sqrt(2), '1.4142135623730950488016887242097', "quadmath sqrt");
121 }