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 #ifdef __mcffpu__
39 #define SKIP_SNAN_CHECKS
40 #endif
41
42 /* pow(x,±0) = 1 for any x (even a zero or quiet NaN) */
43 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x7ff80000, 0x00000000, 0x00000000, 0x00000000 }, /* 1=f(qnan, +0) */
44 #ifndef SKIP_SNAN_CHECKS
45 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff40000, 0x00000000, 0x00000000, 0x00000000 }, /* qnan=f(snan, +0) */
46 #endif
47 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x7ff00000, 0x00000000, 0x00000000, 0x00000000 }, /* 1=f(+inf, +0) */
48 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0xfff00000, 0x00000000, 0x00000000, 0x00000000 }, /* 1=f(-inf, +0) */
49 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x00000000, 0x00000000, 0x00000000, 0x00000000 }, /* 1=f(+0, +0) */
50 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x80000000, 0x00000000, 0x00000000, 0x00000000 }, /* 1=f(-0, +0) */
51 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x7ff80000, 0x00000000, 0x80000000, 0x00000000 }, /* 1=f(qnan, -0) */
52 #ifndef SKIP_SNAN_CHECKS
53 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff40000, 0x00000000, 0x80000000, 0x00000000 }, /* qnan=f(snan, -0) */
54 #endif
55 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x7ff00000, 0x00000000, 0x80000000, 0x00000000 }, /* 1=f(+inf, -0) */
56 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0xfff00000, 0x00000000, 0x80000000, 0x00000000 }, /* 1=f(-inf, -0) */
57 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x00000000, 0x00000000, 0x80000000, 0x00000000 }, /* 1=f(+0, -0) */
58 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x80000000, 0x00000000, 0x80000000, 0x00000000 }, /* 1=f(-0, -0) */
59
60 /* pow(±0, y) = ±∞ and signals divideByZero for y an odd integer < 0 */
61 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0x00000000, 0x00000000, 0xc0080000, 0x00000000 }, /* +inf=f(+0, -3) */
62 {64, 0,123,__LINE__, 0xfff00000, 0x00000000, 0x80000000, 0x00000000, 0xc0080000, 0x00000000 }, /* -inf=f(-0, -3) */
63
64 /* pow(±0, y) = +∞ and signals divideByZero for y < 0 and not an odd integer */
65 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0x00000000, 0x00000000, 0xc0100000, 0x00000000 }, /* +inf=f(+0, -4) */
66 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0x80000000, 0x00000000, 0xc0100000, 0x00000000 }, /* -inf=f(-0, -4) */
67 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0x00000000, 0x00000000, 0xbfe00000, 0x00000000 }, /* +inf=f(+0, -0.5) */
68 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0x80000000, 0x00000000, 0xbfe00000, 0x00000000 }, /* +inf=f(-0, -0.5) */
69
70 /* pow(±0, y)= +0 for y >0 and not an odd integer */
71 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x00000000, 0x00000000, 0x40100000, 0x00000000 }, /* +inf=f(+0, 4) */
72 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x80000000, 0x00000000, 0x40100000, 0x00000000 }, /* -inf=f(-0, 4) */
73 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x00000000, 0x00000000, 0x3fe00000, 0x00000000 }, /* +inf=f(+0, 0.5) */
74 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x80000000, 0x00000000, 0x3fe00000, 0x00000000 }, /* +inf=f(-0, 0.5) */
75
76 /* pow(+1,y)= +1 for any y (even a quiet NaN) */
77 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x3ff00000, 0x00000000, 0x7ff80000, 0x00000000 }, /* 1=f(1, qnan) */
78 #ifndef SKIP_SNAN_CHECKS
79 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x3ff00000, 0x00000000, 0x7ff40000, 0x00000000 }, /* qnan=f(1, snan) */
80 #endif
81 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x3ff00000, 0x00000000, 0x7ff00000, 0x00000000 }, /* 1=f(1, +inf) */
82 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x3ff00000, 0x00000000, 0xfff00000, 0x00000000 }, /* 1=f(1, -inf) */
83 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x3ff00000, 0x00000000, 0x40100000, 0x00000000 }, /* 1=f(1, +4) */
84 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x3ff00000, 0x00000000, 0xc0100000, 0x00000000 }, /* 1=f(1, -4) */
85 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x3ff00000, 0x00000000, 0x3fe00000, 0x00000000 }, /* 1=f(1, +0.5) */
86 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0x3ff00000, 0x00000000, 0xbfe00000, 0x00000000 }, /* 1=f(1, -0.5) */
87
88 /* pow(x,y) returns a quiet NaN and signals invalid for finite x < 0 and finite non-integer y */
89 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0xbff00000, 0x00000000, 0xbfe00000, 0x00000000 }, /* 1=f(-1, -0.5) */
90 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0xc0100000, 0x00000000, 0xbfe00000, 0x00000000 }, /* 1=f(-4, -0.5) */
91
92 /* except for the above cases, nan operands return qnan */
93 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff80000, 0x00000000, 0x7ff80000, 0x00000000 }, /* qnan=f(qnan, qnan) */
94 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff40000, 0x00000000, 0x7ff80000, 0x00000000 }, /* qnan=f(snan, qnan) */
95 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff80000, 0x00000000, 0x7ff40000, 0x00000000 }, /* qnan=f(qnan, snan) */
96 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff40000, 0x00000000, 0x7ff40000, 0x00000000 }, /* qnan=f(snan, snan) */
97 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff80000, 0x00000000, 0x7ff80000, 0x00000000 }, /* qnan=f(qnan, qnan) */
98 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff80000, 0x00000000, 0x3ff00000, 0x00000000 }, /* qnan=f(qnan, 1) */
99 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff40000, 0x00000000, 0x3ff00000, 0x00000000 }, /* qnan=f(snan, 1) */
100 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff80000, 0x00000000, 0xbff00000, 0x00000000 }, /* qnan=f(qnan, -1) */
101 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0x7ff40000, 0x00000000, 0xbff00000, 0x00000000 }, /* qnan=f(snan, -1) */
102 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0xbff00000, 0x00000000, 0x7ff80000, 0x00000000 }, /* qnan=f(-1, qnan) */
103 {64, 0,123,__LINE__, 0x7ff80000, 0x00000000, 0xbff00000, 0x00000000, 0x7ff40000, 0x00000000 }, /* qnan=f(-1, snan) */
104
105 /* If x is +0 (-0), and y is an odd integer greater than 0, the result is +0 (-0). */
106 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x00000000, 0x00000000, 0x40080000, 0x00000000 }, /* +0=f(+0, 3) */
107 {64, 0,123,__LINE__, 0x80000000, 0x00000000, 0x80000000, 0x00000000, 0x40080000, 0x00000000 }, /* -0=f(-0, 3) */
108
109 /* If x is 0, and y greater than 0 and not an odd integer, the result is +0. */
110 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x00000000, 0x00000000, 0x40100000, 0x00000000 }, /* +0=f(+0, 4) */
111 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x80000000, 0x00000000, 0x40100000, 0x00000000 }, /* +0=f(-0, 4) */
112
113 /* If x is -1, and y is positive infinity or negative infinity, the result is 1.0. */
114 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0xbff00000, 0x00000000, 0x7ff00000, 0x00000000 }, /* 1=f(-1, +inf) */
115 {64, 0,123,__LINE__, 0x3ff00000, 0x00000000, 0xbff00000, 0x00000000, 0xfff00000, 0x00000000 }, /* 1=f(-1, -inf) */
116
117 /* If the absolute value of x is less than 1, and y is negative infinity, the result is positive infinity. */
118 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0x3fe00000, 0x00000000, 0xfff00000, 0x00000000 }, /* +inf=f(0.5, -inf) */
119 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0xbfe00000, 0x00000000, 0xfff00000, 0x00000000 }, /* +inf=f(-0.5, -inf) */
120
121 /* If the absolute value of x is greater than 1, and y is negative infinity, the result is +0. */
122 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x40080000, 0x00000000, 0xfff00000, 0x00000000 }, /* +0=f(3.0, -inf) */
123 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0xc0080000, 0x00000000, 0xfff00000, 0x00000000 }, /* +0=f(-3.0, -inf) */
124
125 /* If the absolute value of x is less than 1, and y is positive infinity, the result is +0. */
126 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x3fe00000, 0x00000000, 0x7ff00000, 0x00000000 }, /* +0=f(0.5, +inf) */
127 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0xbfe00000, 0x00000000, 0x7ff00000, 0x00000000 }, /* +0=f(-0.5, +inf) */
128
129 /* If the absolute value of x is greater than 1, and y is positive infinity, the result is positive infinity. */
130 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0x40080000, 0x00000000, 0x7ff00000, 0x00000000 }, /* +inf=f(3.0, inf) */
131 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0xc0080000, 0x00000000, 0x7ff00000, 0x00000000 }, /* +inf=f(-3.0, inf) */
132
133 /* If x is negative infinity, and y is an odd integer less than 0, the result is -0. */
134 {64, 0,123,__LINE__, 0x80000000, 0x00000000, 0xfff00000, 0x00000000, 0xc0080000, 0x00000000 }, /* -0=f(-inf, -3) */
135
136 /* If x is negative infinity, and y less than 0 and not an odd integer, the result is +0. */
137 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0xfff00000, 0x00000000, 0xc0100000, 0x00000000 }, /* +0=f(-inf, -4) */
138 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0xfff00000, 0x00000000, 0xbfe00000, 0x00000000 }, /* +0=f(-inf, -0.5) */
139
140 /* If x is negative infinity, and y is an odd integer greater than 0, the result is negative infinity. */
141 {64, 0,123,__LINE__, 0xfff00000, 0x00000000, 0xfff00000, 0x00000000, 0x40080000, 0x00000000 }, /* -inf=f(-inf, 3) */
142
143 /* If x is negative infinity, and y greater than 0 and not an odd integer, the result is positive infinity. */
144 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0xfff00000, 0x00000000, 0x40100000, 0x00000000 }, /* +inf=f(-inf, 4) */
145 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0xfff00000, 0x00000000, 0x3fe00000, 0x00000000 }, /* +inf=f(-inf, 4) */
146
147 /* If x is positive infinity, and y less than 0, the result is +0. */
148 {64, 0,123,__LINE__, 0x00000000, 0x00000000, 0x7ff00000, 0x00000000, 0xc0080000, 0x00000000 }, /* +0=f(+inf, -3) */
149
150 /* If x is positive infinity, and y greater than 0, the result is positive infinity. */
151 {64, 0,123,__LINE__, 0x7ff00000, 0x00000000, 0x7ff00000, 0x00000000, 0x40100000, 0x00000000 }, /* +inf=f(+inf, 4) */
152
153 /* x is close to one and y is large */
154 {64, 0,123,__LINE__, 0x792ffffe, 0x0bc9e399, 0x3ff00000, 0x2c5e2e99, 0x41ec9eee, 0x35374af6},
155 {64, 0,123,__LINE__, 0x18bfffff, 0xec16bafd, 0x3fefffff, 0xd2e3e669, 0x41f344c9, 0x823eb66c},
156
157 {0},};
test_pow_vec(int m)158 void test_pow_vec(int m) {run_vector_1(m,pow_vec,(char *)(pow),"pow","ddd"); }
159
