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