Commit | Line | Data |
---|---|---|
8d063cd8 LW |
1 | #!./perl |
2 | ||
b51533f3 JH |
3 | # Simple tests for the basic math functions. |
4 | ||
c057f5bc RGS |
5 | BEGIN { |
6 | chdir 't' if -d 't'; | |
7 | @INC = '../lib'; | |
8 | require './test.pl'; | |
9 | } | |
8d063cd8 | 10 | |
f9abf36c JH |
11 | use Config; |
12 | ||
13 | plan tests => 31; | |
8d063cd8 LW |
14 | |
15 | # compile time evaluation | |
16 | ||
b51533f3 JH |
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 | ||
8d063cd8 | 26 | $s = sqrt(2); |
e759bc31 | 27 | is(substr($s,0,5), '1.414', 'compile time sqrt(2) == 1.414'); |
8d063cd8 | 28 | |
b51533f3 JH |
29 | $s = exp(0); |
30 | is($s, 1, 'compile time exp(0) == 1'); | |
31 | ||
8d063cd8 | 32 | $s = exp(1); |
e759bc31 | 33 | is(substr($s,0,7), '2.71828', 'compile time exp(1) == e'); |
8d063cd8 | 34 | |
b51533f3 JH |
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 | ||
e759bc31 | 44 | cmp_ok(exp(log(1)), '==', 1, 'compile time exp(log(1)) == 1'); |
8d063cd8 | 45 | |
5c7d6202 JH |
46 | cmp_ok(round(atan2(1, 2)), '==', '0.463647609', "atan2(1, 2)"); |
47 | ||
8d063cd8 LW |
48 | # run time evaluation |
49 | ||
b51533f3 | 50 | $x0 = 0; |
8d063cd8 LW |
51 | $x1 = 1; |
52 | $x2 = 2; | |
b51533f3 JH |
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 | ||
8d063cd8 | 63 | $s = sqrt($x2); |
e759bc31 | 64 | is(substr($s,0,5), '1.414', 'run time sqrt(2) == 1.414'); |
8d063cd8 | 65 | |
b51533f3 JH |
66 | $s = exp($x0); |
67 | is($s, 1, 'run time exp(0) = 1'); | |
68 | ||
8d063cd8 | 69 | $s = exp($x1); |
e759bc31 | 70 | is(substr($s,0,7), '2.71828', 'run time exp(1) = e'); |
8d063cd8 | 71 | |
b51533f3 JH |
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 | ||
e759bc31 | 81 | cmp_ok(exp(log($x1)), '==', 1, 'run time exp(log(1)) == 1'); |
0630166f | 82 | |
402bce5f JH |
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. | |
0630166f | 87 | |
402bce5f | 88 | my $pi_2 = 1.5707963267949; |
0630166f SP |
89 | |
90 | sub round { | |
91 | my $result = shift; | |
92 | return sprintf("%.9f", $result); | |
93 | } | |
94 | ||
95 | # sin() tests | |
e759bc31 | 96 | cmp_ok(sin(0), '==', 0.0, 'sin(0) == 0'); |
402bce5f JH |
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'); | |
0630166f | 99 | |
b51533f3 JH |
100 | cmp_ok(round(sin($x1)), '==', '0.841470985', "sin(1)"); |
101 | ||
0630166f | 102 | # cos() tests |
e759bc31 | 103 | cmp_ok(cos(0), '==', 1.0, 'cos(0) == 1'); |
402bce5f JH |
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'); | |
0630166f | 106 | |
b51533f3 JH |
107 | cmp_ok(round(cos($x1)), '==', '0.540302306', "cos(1)"); |
108 | ||
5c7d6202 JH |
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. | |
f9abf36c JH |
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 | } |