1/*
2 * IBM Accurate Mathematical Library
3 * Copyright (C) 2001-2022 Free Software Foundation, Inc.
4 *
5 * This program is free software; you can redistribute it and/or modify
6 * it under the terms of the GNU Lesser General Public License as published by
7 * the Free Software Foundation; either version 2.1 of the License, or
8 * (at your option) any later version.
9 *
10 * This program is distributed in the hope that it will be useful,
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 * GNU Lesser General Public License for more details.
14 *
15 * You should have received a copy of the GNU Lesser General Public License
16 * along with this program; if not, see <https://www.gnu.org/licenses/>.
17 */
18
19/************************************************************************/
20/* MODULE_NAME: atnat.h */
21/* */
22/* */
23/* common data and variables definition for BIG or LITTLE ENDIAN */
24/************************************************************************/
25#ifndef ATNAT_H
26#define ATNAT_H
27
28#define M 4
29
30#ifdef BIG_ENDI
31 static const mynumber
32 /* polynomial I */
33/**/ d3 = {{0xbfd55555, 0x55555555} }, /* -0.333... */
34/**/ d5 = {{0x3fc99999, 0x999997fd} }, /* 0.199... */
35/**/ d7 = {{0xbfc24924, 0x923f7603} }, /* -0.142... */
36/**/ d9 = {{0x3fbc71c6, 0xe5129a3b} }, /* 0.111... */
37/**/ d11 = {{0xbfb74580, 0x22b13c25} }, /* -0.090... */
38/**/ d13 = {{0x3fb375f0, 0x8b31cbce} }, /* 0.076... */
39 /* polynomial II */
40/**/ f3 = {{0xbfd55555, 0x55555555} }, /* -1/3 */
41/**/ ff3 = {{0xbc755555, 0x55555555} }, /* -1/3-f3 */
42/**/ f5 = {{0x3fc99999, 0x9999999a} }, /* 1/5 */
43/**/ ff5 = {{0xbc699999, 0x9999999a} }, /* 1/5-f5 */
44/**/ f7 = {{0xbfc24924, 0x92492492} }, /* -1/7 */
45/**/ ff7 = {{0xbc624924, 0x92492492} }, /* -1/7-f7 */
46/**/ f9 = {{0x3fbc71c7, 0x1c71c71c} }, /* 1/9 */
47/**/ ff9 = {{0x3c5c71c7, 0x1c71c71c} }, /* 1/9-f9 */
48/**/ f11 = {{0xbfb745d1, 0x745d1746} }, /* -1/11 */
49/**/ f13 = {{0x3fb3b13b, 0x13b13b14} }, /* 1/13 */
50/**/ f15 = {{0xbfb11111, 0x11111111} }, /* -1/15 */
51/**/ f17 = {{0x3fae1e1e, 0x1e1e1e1e} }, /* 1/17 */
52/**/ f19 = {{0xbfaaf286, 0xbca1af28} }, /* -1/19 */
53 /* constants */
54/**/ a = {{0x3e4bb67a, 0x00000000} }, /* 1.290e-8 */
55/**/ b = {{0x3fb00000, 0x00000000} }, /* 1/16 */
56/**/ c = {{0x3ff00000, 0x00000000} }, /* 1 */
57/**/ d = {{0x40300000, 0x00000000} }, /* 16 */
58/**/ e = {{0x43349ff2, 0x00000000} }, /* 5.805e15 */
59/**/ hpi = {{0x3ff921fb, 0x54442d18} }, /* pi/2 */
60/**/ mhpi = {{0xbff921fb, 0x54442d18} }, /* -pi/2 */
61/**/ hpi1 = {{0x3c91a626, 0x33145c07} }, /* pi/2-hpi */
62/**/ u1 = {{0x3c2d3382, 0x00000000} }, /* 7.915e-19 */
63/**/ u21 = {{0x3c6dffc0, 0x00000000} }, /* 1.301e-17 */
64/**/ u22 = {{0x3c527bd0, 0x00000000} }, /* 4.008e-18 */
65/**/ u23 = {{0x3c3cd057, 0x00000000} }, /* 1.562e-18 */
66/**/ u24 = {{0x3c329cdf, 0x00000000} }, /* 1.009e-18 */
67/**/ u31 = {{0x3c3a1edf, 0x00000000} }, /* 1.416e-18 */
68/**/ u32 = {{0x3c33f0e1, 0x00000000} }, /* 1.081e-18 */
69/**/ u4 = {{0x3bf955e4, 0x00000000} }, /* 8.584e-20 */
70/**/ u5 = {{0x3aaef2d1, 0x00000000} }, /* 5e-26 */
71/**/ u6 = {{0x3a98c56d, 0x00000000} }, /* 2.001e-26 */
72/**/ u7 = {{0x3a9375de, 0x00000000} }, /* 1.572e-26 */
73/**/ u8 = {{0x3a6eeb36, 0x00000000} }, /* 3.122e-27 */
74/**/ u9[M] ={{{0x38c1aa5b, 0x00000000} }, /* 2.658e-35 */
75/**/ {{0x35c1aa4d, 0x00000000} }, /* 9.443e-50 */
76/**/ {{0x32c1aa88, 0x00000000} }, /* 3.355e-64 */
77/**/ {{0x11c1aa56, 0x00000000} }};/* 3.818e-223 */
78
79#else
80#ifdef LITTLE_ENDI
81 static const mynumber
82 /* polynomial I */
83/**/ d3 = {{0x55555555, 0xbfd55555} }, /* -0.333... */
84/**/ d5 = {{0x999997fd, 0x3fc99999} }, /* 0.199... */
85/**/ d7 = {{0x923f7603, 0xbfc24924} }, /* -0.142... */
86/**/ d9 = {{0xe5129a3b, 0x3fbc71c6} }, /* 0.111... */
87/**/ d11 = {{0x22b13c25, 0xbfb74580} }, /* -0.090... */
88/**/ d13 = {{0x8b31cbce, 0x3fb375f0} }, /* 0.076... */
89 /* polynomial II */
90/**/ f3 = {{0x55555555, 0xbfd55555} }, /* -1/3 */
91/**/ ff3 = {{0x55555555, 0xbc755555} }, /* -1/3-f3 */
92/**/ f5 = {{0x9999999a, 0x3fc99999} }, /* 1/5 */
93/**/ ff5 = {{0x9999999a, 0xbc699999} }, /* 1/5-f5 */
94/**/ f7 = {{0x92492492, 0xbfc24924} }, /* -1/7 */
95/**/ ff7 = {{0x92492492, 0xbc624924} }, /* -1/7-f7 */
96/**/ f9 = {{0x1c71c71c, 0x3fbc71c7} }, /* 1/9 */
97/**/ ff9 = {{0x1c71c71c, 0x3c5c71c7} }, /* 1/9-f9 */
98/**/ f11 = {{0x745d1746, 0xbfb745d1} }, /* -1/11 */
99/**/ f13 = {{0x13b13b14, 0x3fb3b13b} }, /* 1/13 */
100/**/ f15 = {{0x11111111, 0xbfb11111} }, /* -1/15 */
101/**/ f17 = {{0x1e1e1e1e, 0x3fae1e1e} }, /* 1/17 */
102/**/ f19 = {{0xbca1af28, 0xbfaaf286} }, /* -1/19 */
103 /* constants */
104/**/ a = {{0x00000000, 0x3e4bb67a} }, /* 1.290e-8 */
105/**/ b = {{0x00000000, 0x3fb00000} }, /* 1/16 */
106/**/ c = {{0x00000000, 0x3ff00000} }, /* 1 */
107/**/ d = {{0x00000000, 0x40300000} }, /* 16 */
108/**/ e = {{0x00000000, 0x43349ff2} }, /* 5.805e15 */
109/**/ hpi = {{0x54442d18, 0x3ff921fb} }, /* pi/2 */
110/**/ mhpi = {{0x54442d18, 0xbff921fb} }, /* -pi/2 */
111/**/ hpi1 = {{0x33145c07, 0x3c91a626} }, /* pi/2-hpi */
112/**/ u1 = {{0x00000000, 0x3c2d3382} }, /* 7.915e-19 */
113/**/ u21 = {{0x00000000, 0x3c6dffc0} }, /* 1.301e-17 */
114/**/ u22 = {{0x00000000, 0x3c527bd0} }, /* 4.008e-18 */
115/**/ u23 = {{0x00000000, 0x3c3cd057} }, /* 1.562e-18 */
116/**/ u24 = {{0x00000000, 0x3c329cdf} }, /* 1.009e-18 */
117/**/ u31 = {{0x00000000, 0x3c3a1edf} }, /* 1.416e-18 */
118/**/ u32 = {{0x00000000, 0x3c33f0e1} }, /* 1.081e-18 */
119/**/ u4 = {{0x00000000, 0x3bf955e4} }, /* 8.584e-20 */
120/**/ u5 = {{0x00000000, 0x3aaef2d1} }, /* 5e-26 */
121/**/ u6 = {{0x00000000, 0x3a98c56d} }, /* 2.001e-26 */
122/**/ u7 = {{0x00000000, 0x3a9375de} }, /* 1.572e-26 */
123/**/ u8 = {{0x00000000, 0x3a6eeb36} }, /* 3.122e-27 */
124/**/ u9[M] ={{{0x00000000, 0x38c1aa5b} }, /* 2.658e-35 */
125/**/ {{0x00000000, 0x35c1aa4d} }, /* 9.443e-50 */
126/**/ {{0x00000000, 0x32c1aa88} }, /* 3.355e-64 */
127/**/ {{0x00000000, 0x11c1aa56} }};/* 3.818e-223 */
128
129#endif
130#endif
131
132#define A a.d
133#define B b.d
134#define C c.d
135#define D d.d
136#define E e.d
137#define HPI hpi.d
138#define MHPI mhpi.d
139#define HPI1 hpi1.d
140#define U1 u1.d
141#define U21 u21.d
142#define U22 u22.d
143#define U23 u23.d
144#define U24 u24.d
145#define U31 u31.d
146#define U32 u32.d
147#define U4 u4.d
148#define U5 u5.d
149#define U6 u6.d
150#define U7 u7.d
151#define U8 u8.d
152
153#endif
154

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