1 2 /* 3 * IBM Accurate Mathematical Library 4 * Written by International Business Machines Corp. 5 * Copyright (C) 2001 Free Software Foundation, Inc. 6 * 7 * This program is free software; you can redistribute it and/or modify 8 * it under the terms of the GNU Lesser General Public License as published by 9 * the Free Software Foundation; either version 2.1 of the License, or 10 * (at your option) any later version. 11 * 12 * This program is distributed in the hope that it will be useful, 13 * but WITHOUT ANY WARRANTY; without even the implied warranty of 14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 15 * GNU Lesser General Public License for more details. 16 * 17 * You should have received a copy of the GNU Lesser General Public License 18 * along with this program; if not, write to the Free Software 19 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. 20 */ 21 22 /******************************************************************/ 23 /* */ 24 /* MODULE_NAME:mpatan.h */ 25 /* */ 26 /* common data and variables prototype and definition */ 27 /******************************************************************/ 28 29 #ifndef MPATAN_H 30 #define MPATAN_H 31 32 #ifdef BIG_ENDI 33 static const number 34 xm[8] = { /* x[m] */ 35 /**/ {{0x00000000, 0x00000000} }, /* 0.0 */ 36 /**/ {{0x3f8930be, 0x00000000} }, /* 0.0123 */ 37 /**/ {{0x3f991687, 0x00000000} }, /* 0.0245 */ 38 /**/ {{0x3fa923a2, 0x00000000} }, /* 0.0491 */ 39 /**/ {{0x3fb930be, 0x00000000} }, /* 0.0984 */ 40 /**/ {{0x3fc95810, 0x00000000} }, /* 0.198 */ 41 /**/ {{0x3fda7ef9, 0x00000000} }, /* 0.414 */ 42 /**/ {{0x3ff00000, 0x00000000} }, /* 1.0 */ 43 }; 44 static const number 45 twonm1[33] = { /* 2n-1 */ 46 /**/ {{0x00000000, 0x00000000} }, /* 0 */ 47 /**/ {{0x00000000, 0x00000000} }, /* 0 */ 48 /**/ {{0x00000000, 0x00000000} }, /* 0 */ 49 /**/ {{0x00000000, 0x00000000} }, /* 0 */ 50 /**/ {{0x40260000, 0x00000000} }, /* 11 */ 51 /**/ {{0x402e0000, 0x00000000} }, /* 15 */ 52 /**/ {{0x40330000, 0x00000000} }, /* 19 */ 53 /**/ {{0x40350000, 0x00000000} }, /* 21 */ 54 /**/ {{0x40390000, 0x00000000} }, /* 25 */ 55 /**/ {{0x403d0000, 0x00000000} }, /* 29 */ 56 /**/ {{0x40408000, 0x00000000} }, /* 33 */ 57 /**/ {{0x40428000, 0x00000000} }, /* 37 */ 58 /**/ {{0x40448000, 0x00000000} }, /* 41 */ 59 /**/ {{0x40468000, 0x00000000} }, /* 45 */ 60 /**/ {{0x40488000, 0x00000000} }, /* 49 */ 61 /**/ {{0x404a8000, 0x00000000} }, /* 53 */ 62 /**/ {{0x404b8000, 0x00000000} }, /* 55 */ 63 /**/ {{0x404d8000, 0x00000000} }, /* 59 */ 64 /**/ {{0x404f8000, 0x00000000} }, /* 63 */ 65 /**/ {{0x4050c000, 0x00000000} }, /* 67 */ 66 /**/ {{0x4051c000, 0x00000000} }, /* 71 */ 67 /**/ {{0x4052c000, 0x00000000} }, /* 75 */ 68 /**/ {{0x4053c000, 0x00000000} }, /* 79 */ 69 /**/ {{0x4054c000, 0x00000000} }, /* 83 */ 70 /**/ {{0x40554000, 0x00000000} }, /* 85 */ 71 /**/ {{0x40564000, 0x00000000} }, /* 89 */ 72 /**/ {{0x40574000, 0x00000000} }, /* 93 */ 73 /**/ {{0x40584000, 0x00000000} }, /* 97 */ 74 /**/ {{0x40594000, 0x00000000} }, /* 101 */ 75 /**/ {{0x405a4000, 0x00000000} }, /* 105 */ 76 /**/ {{0x405b4000, 0x00000000} }, /* 109 */ 77 /**/ {{0x405c4000, 0x00000000} }, /* 113 */ 78 /**/ {{0x405d4000, 0x00000000} }, /* 117 */ 79 }; 80 81 static const number 82 twom[8] = { /* 2**m */ 83 /**/ {{0x3ff00000, 0x00000000} }, /* 1.0 */ 84 /**/ {{0x40000000, 0x00000000} }, /* 2.0 */ 85 /**/ {{0x40100000, 0x00000000} }, /* 4.0 */ 86 /**/ {{0x40200000, 0x00000000} }, /* 8.0 */ 87 /**/ {{0x40300000, 0x00000000} }, /* 16.0 */ 88 /**/ {{0x40400000, 0x00000000} }, /* 32.0 */ 89 /**/ {{0x40500000, 0x00000000} }, /* 64.0 */ 90 /**/ {{0x40600000, 0x00000000} }, /* 128.0 */ 91 }; 92 93 static const number 94 /**/ one = {{0x3ff00000, 0x00000000} }, /* 1 */ 95 /**/ two = {{0x40000000, 0x00000000} }; /* 2 */ 96 97 #else 98 #ifdef LITTLE_ENDI 99 100 static const number 101 xm[8] = { /* x[m] */ 102 /**/ {{0x00000000, 0x00000000} }, /* 0.0 */ 103 /**/ {{0x00000000, 0x3f8930be} }, /* 0.0123 */ 104 /**/ {{0x00000000, 0x3f991687} }, /* 0.0245 */ 105 /**/ {{0x00000000, 0x3fa923a2} }, /* 0.0491 */ 106 /**/ {{0x00000000, 0x3fb930be} }, /* 0.0984 */ 107 /**/ {{0x00000000, 0x3fc95810} }, /* 0.198 */ 108 /**/ {{0x00000000, 0x3fda7ef9} }, /* 0.414 */ 109 /**/ {{0x00000000, 0x3ff00000} }, /* 1.0 */ 110 }; 111 static const number 112 twonm1[33] = { /* 2n-1 */ 113 /**/ {{0x00000000, 0x00000000} }, /* 0 */ 114 /**/ {{0x00000000, 0x00000000} }, /* 0 */ 115 /**/ {{0x00000000, 0x00000000} }, /* 0 */ 116 /**/ {{0x00000000, 0x00000000} }, /* 0 */ 117 /**/ {{0x00000000, 0x40260000} }, /* 11 */ 118 /**/ {{0x00000000, 0x402e0000} }, /* 15 */ 119 /**/ {{0x00000000, 0x40330000} }, /* 19 */ 120 /**/ {{0x00000000, 0x40350000} }, /* 21 */ 121 /**/ {{0x00000000, 0x40390000} }, /* 25 */ 122 /**/ {{0x00000000, 0x403d0000} }, /* 29 */ 123 /**/ {{0x00000000, 0x40408000} }, /* 33 */ 124 /**/ {{0x00000000, 0x40428000} }, /* 37 */ 125 /**/ {{0x00000000, 0x40448000} }, /* 41 */ 126 /**/ {{0x00000000, 0x40468000} }, /* 45 */ 127 /**/ {{0x00000000, 0x40488000} }, /* 49 */ 128 /**/ {{0x00000000, 0x404a8000} }, /* 53 */ 129 /**/ {{0x00000000, 0x404b8000} }, /* 55 */ 130 /**/ {{0x00000000, 0x404d8000} }, /* 59 */ 131 /**/ {{0x00000000, 0x404f8000} }, /* 63 */ 132 /**/ {{0x00000000, 0x4050c000} }, /* 67 */ 133 /**/ {{0x00000000, 0x4051c000} }, /* 71 */ 134 /**/ {{0x00000000, 0x4052c000} }, /* 75 */ 135 /**/ {{0x00000000, 0x4053c000} }, /* 79 */ 136 /**/ {{0x00000000, 0x4054c000} }, /* 83 */ 137 /**/ {{0x00000000, 0x40554000} }, /* 85 */ 138 /**/ {{0x00000000, 0x40564000} }, /* 89 */ 139 /**/ {{0x00000000, 0x40574000} }, /* 93 */ 140 /**/ {{0x00000000, 0x40584000} }, /* 97 */ 141 /**/ {{0x00000000, 0x40594000} }, /* 101 */ 142 /**/ {{0x00000000, 0x405a4000} }, /* 105 */ 143 /**/ {{0x00000000, 0x405b4000} }, /* 109 */ 144 /**/ {{0x00000000, 0x405c4000} }, /* 113 */ 145 /**/ {{0x00000000, 0x405d4000} }, /* 117 */ 146 }; 147 148 static const number 149 twom[8] = { /* 2**m */ 150 /**/ {{0x00000000, 0x3ff00000} }, /* 1.0 */ 151 /**/ {{0x00000000, 0x40000000} }, /* 2.0 */ 152 /**/ {{0x00000000, 0x40100000} }, /* 4.0 */ 153 /**/ {{0x00000000, 0x40200000} }, /* 8.0 */ 154 /**/ {{0x00000000, 0x40300000} }, /* 16.0 */ 155 /**/ {{0x00000000, 0x40400000} }, /* 32.0 */ 156 /**/ {{0x00000000, 0x40500000} }, /* 64.0 */ 157 /**/ {{0x00000000, 0x40600000} }, /* 128.0 */ 158 }; 159 160 static const number 161 /**/ one = {{0x00000000, 0x3ff00000} }, /* 1 */ 162 /**/ two = {{0x00000000, 0x40000000} }; /* 2 */ 163 164 #endif 165 #endif 166 167 #define ONE one.d 168 #define TWO two.d 169 170 static const int 171 np[33] = { 0, 0, 0, 0, 6, 8,10,11,13,15,17,19,21,23,25,27,28, 172 30,32,34,36,38,40,42,43,45,47,49,51,53,55,57,59}; 173 174 #endif 175