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