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