jcs.org
Serenity Operating System
1/*
2 * Copyright (c) 2018-2020, Andreas Kling <kling@serenityos.org>
3 * All rights reserved.
4 *
5 * Redistribution and use in source and binary forms, with or without
6 * modification, are permitted provided that the following conditions are met:
7 *
8 * 1. Redistributions of source code must retain the above copyright notice, this
9 * list of conditions and the following disclaimer.
10 *
11 * 2. Redistributions in binary form must reproduce the above copyright notice,
12 * this list of conditions and the following disclaimer in the documentation
13 * and/or other materials provided with the distribution.
14 *
15 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
16 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
17 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
18 * DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE
19 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
20 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
21 * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
22 * CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
23 * OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
24 * OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
25 */
26
27#include <AK/TestSuite.h>
28
29#include <math.h>
30
31#define EXPECT_CLOSE(a, b) { EXPECT(fabs(a - b) < 0.000001); }
32
33TEST_CASE(trig)
34{
35 EXPECT_CLOSE(sin(1234), 0.653316);
36 EXPECT_CLOSE(cos(1234), -0.830914);
37 EXPECT_CLOSE(tan(1234), -0.786262);
38 EXPECT_CLOSE(sqrt(1234), 35.128336)
39 EXPECT_CLOSE(sin(-1), -0.867955);
40 EXPECT_CLOSE(cos(-1), 0.594715);
41 EXPECT_CLOSE(tan(-1), -1.459446);
42 EXPECT(isnan(sqrt(-1)));
43 EXPECT(isnan(asin(1.1)));
44 EXPECT(isnan(asin(-1.1)));
45 EXPECT_CLOSE(asin(0), 0.0);
46 EXPECT_CLOSE(asin(0.01), 0.01);
47 EXPECT_CLOSE(asin(0.1), 0.100167);
48 EXPECT_CLOSE(asin(0.3), 0.304693);
49 EXPECT_CLOSE(asin(0.499), 0.522444);
50 EXPECT_CLOSE(asin(0.5), 0.523599);
51 EXPECT_CLOSE(asin(0.501), 0.524754);
52 EXPECT_CLOSE(asin(0.9), 1.119770);
53 EXPECT_CLOSE(asin(0.99), 1.429246);
54 EXPECT_CLOSE(asin(1.0), 1.570750);
55 EXPECT_CLOSE(atan(0), 0.0)
56 EXPECT_CLOSE(atan(0.5), 0.463648)
57 EXPECT_CLOSE(atan(-0.5), -0.463648)
58 EXPECT_CLOSE(atan(5.5), 1.390943)
59 EXPECT_CLOSE(atan(-5.5), -1.390943)
60 EXPECT_CLOSE(atan(555.5), 1.568996)
61}
62
63TEST_CASE(other)
64{
65 EXPECT_EQ(trunc(9999999999999.5), 9999999999999.0);
66 EXPECT_EQ(trunc(-9999999999999.5), -9999999999999.0);
67}
68
69TEST_CASE(exponents)
70{
71 struct values {
72 double x;
73 double exp;
74 double sinh;
75 double cosh;
76 double tanh;
77 };
78
79 values values[8] {
80 { 1.500000, 4.481626, 2.129246, 2.352379, 0.905148},
81 { 20.990000, 1304956710.432035, 652478355.216017, 652478355.216017, 1.000000},
82 { 20.010000, 490041186.687082, 245020593.343541, 245020593.343541, 1.000000},
83 { 0.000000, 1.000000, 0.000000, 1.000000, 0.000000},
84 { 0.010000, 1.010050, 0.010000, 1.000050, 0.010000},
85 { -0.010000, 0.990050, -0.010000, 1.000050, -0.010000},
86 { -1.000000, 0.367879, -1.175201, 1.543081, -0.761594},
87 { -17.000000, 0.000000, -12077476.376788, 12077476.376788, -1.000000},
88 };
89 for (auto& v : values) {
90 EXPECT_CLOSE(exp(v.x), v.exp);
91 EXPECT_CLOSE(sinh(v.x), v.sinh);
92 EXPECT_CLOSE(cosh(v.x), v.cosh);
93 EXPECT_CLOSE(tanh(v.x), v.tanh);
94 }
95 EXPECT_EQ(exp(1000), std::numeric_limits<double>::infinity());
96}
97
98TEST_CASE(logarithms)
99{
100 EXPECT(isnan(log(-1)));
101 EXPECT(log(0) < -1000000);
102 EXPECT_CLOSE(log(0.5), -0.693233)
103 EXPECT_CLOSE(log(1.1), 0.095310)
104 EXPECT_CLOSE(log(5), 1.609480)
105 EXPECT_CLOSE(log(5.5), 1.704842)
106 EXPECT_CLOSE(log(500), 6.214104)
107 EXPECT_CLOSE(log2(5), 2.321989)
108 EXPECT_CLOSE(log10(5), 0.698988)
109}
110
111TEST_MAIN(Math)