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