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