1
2/*
3 * IBM Accurate Mathematical Library
4 * Copyright (C) 2001-2022 Free Software Foundation, Inc.
5 *
6 * This program is free software; you can redistribute it and/or modify
7 * it under the terms of the GNU Lesser General Public License as published by
8 * the Free Software Foundation; either version 2.1 of the License, or
9 * (at your option) any later version.
10 *
11 * This program is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14 * GNU Lesser General Public License for more details.
15 *
16 * You should have received a copy of the GNU Lesser General Public License
17 * along with this program; if not, see <https://www.gnu.org/licenses/>.
18 */
19
20/************************************************************************/
21/* MODULE_NAME: atnat2.h */
22/* */
23/* */
24/* common data and variables definition for BIG or LITTLE ENDIAN */
25/************************************************************************/
26
27
28
29#ifndef ATNAT2_H
30#define ATNAT2_H
31
32
33#define MM 5
34#ifdef BIG_ENDI
35
36 static const mynumber
37 /* polynomial I */
38/**/ d3 = {{0xbfd55555, 0x55555555} }, /* -0.333... */
39/**/ d5 = {{0x3fc99999, 0x999997fd} }, /* 0.199... */
40/**/ d7 = {{0xbfc24924, 0x923f7603} }, /* -0.142... */
41/**/ d9 = {{0x3fbc71c6, 0xe5129a3b} }, /* 0.111... */
42/**/ d11 = {{0xbfb74580, 0x22b13c25} }, /* -0.090... */
43/**/ d13 = {{0x3fb375f0, 0x8b31cbce} }, /* 0.076... */
44 /* polynomial II */
45/**/ f3 = {{0xbfd55555, 0x55555555} }, /* -1/3 */
46/**/ ff3 = {{0xbc755555, 0x55555555} }, /* -1/3-f3 */
47/**/ f5 = {{0x3fc99999, 0x9999999a} }, /* 1/5 */
48/**/ ff5 = {{0xbc699999, 0x9999999a} }, /* 1/5-f5 */
49/**/ f7 = {{0xbfc24924, 0x92492492} }, /* -1/7 */
50/**/ ff7 = {{0xbc624924, 0x92492492} }, /* -1/7-f7 */
51/**/ f9 = {{0x3fbc71c7, 0x1c71c71c} }, /* 1/9 */
52/**/ ff9 = {{0x3c5c71c7, 0x1c71c71c} }, /* 1/9-f9 */
53/**/ f11 = {{0xbfb745d1, 0x745d1746} }, /* -1/11 */
54/**/ f13 = {{0x3fb3b13b, 0x13b13b14} }, /* 1/13 */
55/**/ f15 = {{0xbfb11111, 0x11111111} }, /* -1/15 */
56/**/ f17 = {{0x3fae1e1e, 0x1e1e1e1e} }, /* 1/17 */
57/**/ f19 = {{0xbfaaf286, 0xbca1af28} }, /* -1/19 */
58 /* constants */
59/**/ inv16 = {{0x3fb00000, 0x00000000} }, /* 1/16 */
60/**/ opi = {{0x400921fb, 0x54442d18} }, /* pi */
61/**/ opi1 = {{0x3ca1a626, 0x33145c07} }, /* pi-opi */
62/**/ mopi = {{0xc00921fb, 0x54442d18} }, /* -pi */
63/**/ hpi = {{0x3ff921fb, 0x54442d18} }, /* pi/2 */
64/**/ hpi1 = {{0x3c91a626, 0x33145c07} }, /* pi/2-hpi */
65/**/ mhpi = {{0xbff921fb, 0x54442d18} }, /* -pi/2 */
66/**/ qpi = {{0x3fe921fb, 0x54442d18} }, /* pi/4 */
67/**/ mqpi = {{0xbfe921fb, 0x54442d18} }, /* -pi/4 */
68/**/ tqpi = {{0x4002d97c, 0x7f3321d2} }, /* 3pi/4 */
69/**/ mtqpi = {{0xc002d97c, 0x7f3321d2} }, /* -3pi/4 */
70/**/ u1 = {{0x3c314c2a, 0x00000000} }, /* 9.377e-19 */
71/**/ u2 = {{0x3bf955e4, 0x00000000} }, /* 8.584e-20 */
72/**/ u3 = {{0x3bf955e4, 0x00000000} }, /* 8.584e-20 */
73/**/ u4 = {{0x3bf955e4, 0x00000000} }, /* 8.584e-20 */
74/**/ u5 = {{0x3aaef2d1, 0x00000000} }, /* 5e-26 */
75/**/ u6 = {{0x3a6eeb36, 0x00000000} }, /* 3.122e-27 */
76/**/ u7 = {{0x3a6eeb36, 0x00000000} }, /* 3.122e-27 */
77/**/ u8 = {{0x3a6eeb36, 0x00000000} }, /* 3.122e-27 */
78/**/ u91 = {{0x3c6dffc0, 0x00000000} }, /* 1.301e-17 */
79/**/ u92 = {{0x3c527bd0, 0x00000000} }, /* 4.008e-18 */
80/**/ u93 = {{0x3c3cd057, 0x00000000} }, /* 1.562e-18 */
81/**/ u94 = {{0x3c329cdf, 0x00000000} }, /* 1.009e-18 */
82/**/ ua1 = {{0x3c3a1edf, 0x00000000} }, /* 1.416e-18 */
83/**/ ua2 = {{0x3c33f0e1, 0x00000000} }, /* 1.081e-18 */
84/**/ ub = {{0x3a98c56d, 0x00000000} }, /* 2.001e-26 */
85/**/ uc = {{0x3a9375de, 0x00000000} }, /* 1.572e-26 */
86/**/ ud[MM] ={{{0x38c6eddf, 0x00000000} }, /* 3.450e-35 */
87/**/ {{0x35c6ef60, 0x00000000} }, /* 1.226e-49 */
88/**/ {{0x32c6ed2f, 0x00000000} }, /* 4.354e-64 */
89/**/ {{0x23c6eee8, 0x00000000} }, /* 2.465e-136 */
90/**/ {{0x11c6ed16, 0x00000000} }},/* 4.955e-223 */
91/**/ ue = {{0x38900e9d, 0x00000000} }, /* 3.02e-36 */
92/**/ two500 = {{0x5f300000, 0x00000000} }, /* 2**500 */
93/**/ twom500 = {{0x20b00000, 0x00000000} }; /* 2**(-500) */
94
95#else
96#ifdef LITTLE_ENDI
97
98 static const mynumber
99 /* polynomial I */
100/**/ d3 = {{0x55555555, 0xbfd55555} }, /* -0.333... */
101/**/ d5 = {{0x999997fd, 0x3fc99999} }, /* 0.199... */
102/**/ d7 = {{0x923f7603, 0xbfc24924} }, /* -0.142... */
103/**/ d9 = {{0xe5129a3b, 0x3fbc71c6} }, /* 0.111... */
104/**/ d11 = {{0x22b13c25, 0xbfb74580} }, /* -0.090... */
105/**/ d13 = {{0x8b31cbce, 0x3fb375f0} }, /* 0.076... */
106 /* polynomial II */
107/**/ f3 = {{0x55555555, 0xbfd55555} }, /* -1/3 */
108/**/ ff3 = {{0x55555555, 0xbc755555} }, /* -1/3-f3 */
109/**/ f5 = {{0x9999999a, 0x3fc99999} }, /* 1/5 */
110/**/ ff5 = {{0x9999999a, 0xbc699999} }, /* 1/5-f5 */
111/**/ f7 = {{0x92492492, 0xbfc24924} }, /* -1/7 */
112/**/ ff7 = {{0x92492492, 0xbc624924} }, /* -1/7-f7 */
113/**/ f9 = {{0x1c71c71c, 0x3fbc71c7} }, /* 1/9 */
114/**/ ff9 = {{0x1c71c71c, 0x3c5c71c7} }, /* 1/9-f9 */
115/**/ f11 = {{0x745d1746, 0xbfb745d1} }, /* -1/11 */
116/**/ f13 = {{0x13b13b14, 0x3fb3b13b} }, /* 1/13 */
117/**/ f15 = {{0x11111111, 0xbfb11111} }, /* -1/15 */
118/**/ f17 = {{0x1e1e1e1e, 0x3fae1e1e} }, /* 1/17 */
119/**/ f19 = {{0xbca1af28, 0xbfaaf286} }, /* -1/19 */
120 /* constants */
121/**/ inv16 = {{0x00000000, 0x3fb00000} }, /* 1/16 */
122/**/ opi = {{0x54442d18, 0x400921fb} }, /* pi */
123/**/ opi1 = {{0x33145c07, 0x3ca1a626} }, /* pi-opi */
124/**/ mopi = {{0x54442d18, 0xc00921fb} }, /* -pi */
125/**/ hpi = {{0x54442d18, 0x3ff921fb} }, /* pi/2 */
126/**/ hpi1 = {{0x33145c07, 0x3c91a626} }, /* pi/2-hpi */
127/**/ mhpi = {{0x54442d18, 0xbff921fb} }, /* -pi/2 */
128/**/ qpi = {{0x54442d18, 0x3fe921fb} }, /* pi/4 */
129/**/ mqpi = {{0x54442d18, 0xbfe921fb} }, /* -pi/4 */
130/**/ tqpi = {{0x7f3321d2, 0x4002d97c} }, /* 3pi/4 */
131/**/ mtqpi = {{0x7f3321d2, 0xc002d97c} }, /* -3pi/4 */
132/**/ u1 = {{0x00000000, 0x3c314c2a} }, /* 9.377e-19 */
133/**/ u2 = {{0x00000000, 0x3bf955e4} }, /* 8.584e-20 */
134/**/ u3 = {{0x00000000, 0x3bf955e4} }, /* 8.584e-20 */
135/**/ u4 = {{0x00000000, 0x3bf955e4} }, /* 8.584e-20 */
136/**/ u5 = {{0x00000000, 0x3aaef2d1} }, /* 5e-26 */
137/**/ u6 = {{0x00000000, 0x3a6eeb36} }, /* 3.122e-27 */
138/**/ u7 = {{0x00000000, 0x3a6eeb36} }, /* 3.122e-27 */
139/**/ u8 = {{0x00000000, 0x3a6eeb36} }, /* 3.122e-27 */
140/**/ u91 = {{0x00000000, 0x3c6dffc0} }, /* 1.301e-17 */
141/**/ u92 = {{0x00000000, 0x3c527bd0} }, /* 4.008e-18 */
142/**/ u93 = {{0x00000000, 0x3c3cd057} }, /* 1.562e-18 */
143/**/ u94 = {{0x00000000, 0x3c329cdf} }, /* 1.009e-18 */
144/**/ ua1 = {{0x00000000, 0x3c3a1edf} }, /* 1.416e-18 */
145/**/ ua2 = {{0x00000000, 0x3c33f0e1} }, /* 1.081e-18 */
146/**/ ub = {{0x00000000, 0x3a98c56d} }, /* 2.001e-26 */
147/**/ uc = {{0x00000000, 0x3a9375de} }, /* 1.572e-26 */
148/**/ ud[MM] ={{{0x00000000, 0x38c6eddf} }, /* 3.450e-35 */
149/**/ {{0x00000000, 0x35c6ef60} }, /* 1.226e-49 */
150/**/ {{0x00000000, 0x32c6ed2f} }, /* 4.354e-64 */
151/**/ {{0x00000000, 0x23c6eee8} }, /* 2.465e-136 */
152/**/ {{0x00000000, 0x11c6ed16} }},/* 4.955e-223 */
153/**/ ue = {{0x00000000, 0x38900e9d} }, /* 3.02e-36 */
154/**/ two500 = {{0x00000000, 0x5f300000} }, /* 2**500 */
155/**/ twom500 = {{0x00000000, 0x20b00000} }; /* 2**(-500) */
156
157#endif
158#endif
159
160#endif
161

source code of glibc/sysdeps/ieee754/dbl-64/atnat2.h