1 /*
2 * SPDX-License-Identifier: BSD-3-Clause
3 *
4 * Copyright © 2019 Keith Packard
5 *
6 * Redistribution and use in source and binary forms, with or without
7 * modification, are permitted provided that the following conditions
8 * are met:
9 *
10 * 1. Redistributions of source code must retain the above copyright
11 * notice, this list of conditions and the following disclaimer.
12 *
13 * 2. Redistributions in binary form must reproduce the above
14 * copyright notice, this list of conditions and the following
15 * disclaimer in the documentation and/or other materials provided
16 * with the distribution.
17 *
18 * 3. Neither the name of the copyright holder nor the names of its
19 * contributors may be used to endorse or promote products derived
20 * from this software without specific prior written permission.
21 *
22 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
23 * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
24 * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
25 * FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
26 * COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT,
27 * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
28 * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
29 * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
30 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,
31 * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
32 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED
33 * OF THE POSSIBILITY OF SUCH DAMAGE.
34 */
35 #include "test.h"
36 one_line_type pow_vec[] = {
37
38 /* x 1.021918654441833 y 4092.023193359375 result is 3.4027371058175e38 */
39 {64, 0,123,__LINE__, 0x47efffca, 0xa6b471f3, 0x3ff059c7, 0x60000000, 0x40aff80b, 0xe0000000},
40
41 #ifdef __mcffpu__
42 #define SKIP_SNAN_CHECKS
43 #endif
44
45 /* pow(x,±0) = 1 for any x (even a zero or quiet NaN) */
46 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x7ff80000, 0x00000000, 0x00000000, 0x00000000 }, /* 1=f(qnan, +0) */
47 #ifndef SKIP_SNAN_CHECKS
48 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff40000, 0x00000000, 0x00000000, 0x00000000 }, /* qnan=f(snan, +0) */
49 #endif
50 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x7ff00000, 0x00000000, 0x00000000, 0x00000000 }, /* 1=f(+inf, +0) */
51 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0xfff00000, 0x00000000, 0x00000000, 0x00000000 }, /* 1=f(-inf, +0) */
52 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x00000000, 0x00000000, 0x00000000, 0x00000000 }, /* 1=f(+0, +0) */
53 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x80000000, 0x00000000, 0x00000000, 0x00000000 }, /* 1=f(-0, +0) */
54 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x7ff80000, 0x00000000, 0x80000000, 0x00000000 }, /* 1=f(qnan, -0) */
55 #ifndef SKIP_SNAN_CHECKS
56 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff40000, 0x00000000, 0x80000000, 0x00000000 }, /* qnan=f(snan, -0) */
57 #endif
58 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x7ff00000, 0x00000000, 0x80000000, 0x00000000 }, /* 1=f(+inf, -0) */
59 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0xfff00000, 0x00000000, 0x80000000, 0x00000000 }, /* 1=f(-inf, -0) */
60 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x00000000, 0x00000000, 0x80000000, 0x00000000 }, /* 1=f(+0, -0) */
61 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x80000000, 0x00000000, 0x80000000, 0x00000000 }, /* 1=f(-0, -0) */
62
63 /* pow(±0, y) = ±∞ and signals divideByZero for y an odd integer < 0 */
64 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0x00000000, 0x00000000, 0xc0080000, 0x00000000 }, /* +inf=f(+0, -3) */
65 {64, 0,123,__LINE__, 0xfff00000, 0x00000000, 0x80000000, 0x00000000, 0xc0080000, 0x00000000 }, /* -inf=f(-0, -3) */
66
67 /* pow(±0, y) = +∞ and signals divideByZero for y < 0 and not an odd integer */
68 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0x00000000, 0x00000000, 0xc0100000, 0x00000000 }, /* +inf=f(+0, -4) */
69 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0x80000000, 0x00000000, 0xc0100000, 0x00000000 }, /* -inf=f(-0, -4) */
70 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0x00000000, 0x00000000, 0xbfe00000, 0x00000000 }, /* +inf=f(+0, -0.5) */
71 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0x80000000, 0x00000000, 0xbfe00000, 0x00000000 }, /* +inf=f(-0, -0.5) */
72
73 /* pow(±0, y)= +0 for y >0 and not an odd integer */
74 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x00000000, 0x00000000, 0x40100000, 0x00000000 }, /* +inf=f(+0, 4) */
75 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x80000000, 0x00000000, 0x40100000, 0x00000000 }, /* -inf=f(-0, 4) */
76 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x00000000, 0x00000000, 0x3fe00000, 0x00000000 }, /* +inf=f(+0, 0.5) */
77 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x80000000, 0x00000000, 0x3fe00000, 0x00000000 }, /* +inf=f(-0, 0.5) */
78
79 /* pow(+1,y)= +1 for any y (even a quiet NaN) */
80 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x3ff00000, 0x00000000, 0x7ff80000, 0x00000000 }, /* 1=f(1, qnan) */
81 #ifndef SKIP_SNAN_CHECKS
82 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x3ff00000, 0x00000000, 0x7ff40000, 0x00000000 }, /* qnan=f(1, snan) */
83 #endif
84 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x3ff00000, 0x00000000, 0x7ff00000, 0x00000000 }, /* 1=f(1, +inf) */
85 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x3ff00000, 0x00000000, 0xfff00000, 0x00000000 }, /* 1=f(1, -inf) */
86 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x3ff00000, 0x00000000, 0x40100000, 0x00000000 }, /* 1=f(1, +4) */
87 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x3ff00000, 0x00000000, 0xc0100000, 0x00000000 }, /* 1=f(1, -4) */
88 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x3ff00000, 0x00000000, 0x3fe00000, 0x00000000 }, /* 1=f(1, +0.5) */
89 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x3ff00000, 0x00000000, 0xbfe00000, 0x00000000 }, /* 1=f(1, -0.5) */
90
91 /* pow(x,y) returns a quiet NaN and signals invalid for finite x < 0 and finite non-integer y */
92 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0xbff00000, 0x00000000, 0xbfe00000, 0x00000000 }, /* 1=f(-1, -0.5) */
93 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0xc0100000, 0x00000000, 0xbfe00000, 0x00000000 }, /* 1=f(-4, -0.5) */
94
95 /* except for the above cases, nan operands return qnan */
96 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff80000, 0x00000000, 0x7ff80000, 0x00000000 }, /* qnan=f(qnan, qnan) */
97 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff40000, 0x00000000, 0x7ff80000, 0x00000000 }, /* qnan=f(snan, qnan) */
98 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff80000, 0x00000000, 0x7ff40000, 0x00000000 }, /* qnan=f(qnan, snan) */
99 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff40000, 0x00000000, 0x7ff40000, 0x00000000 }, /* qnan=f(snan, snan) */
100 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff80000, 0x00000000, 0x7ff80000, 0x00000000 }, /* qnan=f(qnan, qnan) */
101 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff80000, 0x00000000, 0x3ff00000, 0x00000000 }, /* qnan=f(qnan, 1) */
102 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff40000, 0x00000000, 0x3ff00000, 0x00000000 }, /* qnan=f(snan, 1) */
103 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff80000, 0x00000000, 0xbff00000, 0x00000000 }, /* qnan=f(qnan, -1) */
104 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff40000, 0x00000000, 0xbff00000, 0x00000000 }, /* qnan=f(snan, -1) */
105 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0xbff00000, 0x00000000, 0x7ff80000, 0x00000000 }, /* qnan=f(-1, qnan) */
106 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0xbff00000, 0x00000000, 0x7ff40000, 0x00000000 }, /* qnan=f(-1, snan) */
107
108 /* If x is +0 (-0), and y is an odd integer greater than 0, the result is +0 (-0). */
109 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x00000000, 0x00000000, 0x40080000, 0x00000000 }, /* +0=f(+0, 3) */
110 {64, 0,123,__LINE__, 0x80000000, 0x00000000, 0x80000000, 0x00000000, 0x40080000, 0x00000000 }, /* -0=f(-0, 3) */
111
112 /* If x is 0, and y greater than 0 and not an odd integer, the result is +0. */
113 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x00000000, 0x00000000, 0x40100000, 0x00000000 }, /* +0=f(+0, 4) */
114 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x80000000, 0x00000000, 0x40100000, 0x00000000 }, /* +0=f(-0, 4) */
115
116 /* If x is -1, and y is positive infinity or negative infinity, the result is 1.0. */
117 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0xbff00000, 0x00000000, 0x7ff00000, 0x00000000 }, /* 1=f(-1, +inf) */
118 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0xbff00000, 0x00000000, 0xfff00000, 0x00000000 }, /* 1=f(-1, -inf) */
119
120 /* If the absolute value of x is less than 1, and y is negative infinity, the result is positive infinity. */
121 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0x3fe00000, 0x00000000, 0xfff00000, 0x00000000 }, /* +inf=f(0.5, -inf) */
122 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0xbfe00000, 0x00000000, 0xfff00000, 0x00000000 }, /* +inf=f(-0.5, -inf) */
123
124 /* If the absolute value of x is greater than 1, and y is negative infinity, the result is +0. */
125 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x40080000, 0x00000000, 0xfff00000, 0x00000000 }, /* +0=f(3.0, -inf) */
126 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0xc0080000, 0x00000000, 0xfff00000, 0x00000000 }, /* +0=f(-3.0, -inf) */
127
128 /* If the absolute value of x is less than 1, and y is positive infinity, the result is +0. */
129 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x3fe00000, 0x00000000, 0x7ff00000, 0x00000000 }, /* +0=f(0.5, +inf) */
130 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0xbfe00000, 0x00000000, 0x7ff00000, 0x00000000 }, /* +0=f(-0.5, +inf) */
131
132 /* If the absolute value of x is greater than 1, and y is positive infinity, the result is positive infinity. */
133 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0x40080000, 0x00000000, 0x7ff00000, 0x00000000 }, /* +inf=f(3.0, inf) */
134 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0xc0080000, 0x00000000, 0x7ff00000, 0x00000000 }, /* +inf=f(-3.0, inf) */
135
136 /* If x is negative infinity, and y is an odd integer less than 0, the result is -0. */
137 {64, 0,123,__LINE__, 0x80000000, 0x00000000, 0xfff00000, 0x00000000, 0xc0080000, 0x00000000 }, /* -0=f(-inf, -3) */
138
139 /* If x is negative infinity, and y less than 0 and not an odd integer, the result is +0. */
140 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0xfff00000, 0x00000000, 0xc0100000, 0x00000000 }, /* +0=f(-inf, -4) */
141 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0xfff00000, 0x00000000, 0xbfe00000, 0x00000000 }, /* +0=f(-inf, -0.5) */
142
143 /* If x is negative infinity, and y is an odd integer greater than 0, the result is negative infinity. */
144 {64, 0,123,__LINE__, 0xfff00000, 0x00000000, 0xfff00000, 0x00000000, 0x40080000, 0x00000000 }, /* -inf=f(-inf, 3) */
145
146 /* If x is negative infinity, and y greater than 0 and not an odd integer, the result is positive infinity. */
147 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0xfff00000, 0x00000000, 0x40100000, 0x00000000 }, /* +inf=f(-inf, 4) */
148 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0xfff00000, 0x00000000, 0x3fe00000, 0x00000000 }, /* +inf=f(-inf, 4) */
149
150 /* If x is positive infinity, and y less than 0, the result is +0. */
151 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x7ff00000, 0x00000000, 0xc0080000, 0x00000000 }, /* +0=f(+inf, -3) */
152
153 /* If x is positive infinity, and y greater than 0, the result is positive infinity. */
154 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0x7ff00000, 0x00000000, 0x40100000, 0x00000000 }, /* +inf=f(+inf, 4) */
155
156 /* x is close to one and y is large */
157 {64, 0,123,__LINE__, 0x792ffffe, 0x0bc9e399, 0x3ff00000, 0x2c5e2e99, 0x41ec9eee, 0x35374af6},
158 {64, 0,123,__LINE__, 0x18bfffff, 0xec16bafd, 0x3fefffff, 0xd2e3e669, 0x41f344c9, 0x823eb66c},
159
160 /* x -1.1 y 101. float result is -15158.707 */
161 {64, 0,123,__LINE__, 0xc0cd9b5a, 0x770b46a4, 0xbff19999, 0xa0000000, 0x40594000, 0x00000000},
162
163 {0},};
test_pow_vec(int m)164 void test_pow_vec(int m) {run_vector_1(m,pow_vec,(char *)(pow),"pow","ddd"); }
165
