1 | /* |
2 | * IBM Accurate Mathematical Library |
3 | * Written by International Business Machines Corp. |
4 | * Copyright (C) 2001-2018 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 <http://www.gnu.org/licenses/>. |
18 | */ |
19 | |
20 | /******************************************************************/ |
21 | /* */ |
22 | /* MODULE_NAME:upow.h */ |
23 | /* */ |
24 | /* common data and variables prototype and definition */ |
25 | /******************************************************************/ |
26 | |
27 | #ifndef UPOW_H |
28 | #define UPOW_H |
29 | |
30 | #include "mydefs.h" |
31 | |
32 | #ifdef BIG_ENDI |
33 | const static mynumber |
34 | /**/ nZERO = {{0x80000000, 0}}, /* -0.0 */ |
35 | /**/ INF = {{0x7ff00000, 0x00000000}}, /* INF */ |
36 | /**/ nINF = {{0xfff00000, 0x00000000}}, /* -INF */ |
37 | /**/ ln2a = {{0x3fe62e42, 0xfefa3800}}, /* ln(2) 43 bits */ |
38 | /**/ ln2b = {{0x3d2ef357, 0x93c76730}}, /* ln(2)-ln2a */ |
39 | /**/ bigu = {{0x4297ffff, 0xfffffd2c}}, /* 1.5*2**42 -724*2**-10 */ |
40 | /**/ bigv = {{0x4207ffff, 0xfff8016a}}, /* 1.5*2**33-1+362*2**-19 */ |
41 | /**/ t52 = {{0x43300000, 0x00000000}}, /* 2**52 */ |
42 | /**/ two52e = {{0x43300000, 0x000003ff}}; /* 2**52' */ |
43 | |
44 | #else |
45 | #ifdef LITTLE_ENDI |
46 | const static mynumber |
47 | /**/ nZERO = {{0, 0x80000000}}, /* -0.0 */ |
48 | /**/ INF = {{0x00000000, 0x7ff00000}}, /* INF */ |
49 | /**/ nINF = {{0x00000000, 0xfff00000}}, /* -INF */ |
50 | /**/ ln2a = {{0xfefa3800, 0x3fe62e42}}, /* ln(2) 43 bits */ |
51 | /**/ ln2b = {{0x93c76730, 0x3d2ef357}}, /* ln(2)-ln2a */ |
52 | /**/ bigu = {{0xfffffd2c, 0x4297ffff}}, /* 1.5*2**42 -724*2**-10 */ |
53 | /**/ bigv = {{0xfff8016a, 0x4207ffff}}, /* 1.5*2**33-1+362*2**-19 */ |
54 | /**/ t52 = {{0x00000000, 0x43300000}}, /* 2**52 */ |
55 | /**/ two52e = {{0x000003ff, 0x43300000}}; /* 2**52' */ |
56 | |
57 | #endif |
58 | #endif |
59 | |
60 | const static double p2=-0.5, p3 = 3.3333333333333333333e-1, p4 = -0.25, |
61 | q2 = -0.5, q3 = 3.3333333333331404e-01, q4 = -2.4999999999996436e-01, |
62 | q5 = 2.0000010500004459e-01, q6 = -1.6666678916688004e-01, |
63 | r3 = 3.33333333333333333372884096563030E-01, |
64 | r4 = -2.50000000000000000213574153875908E-01, |
65 | r5 = 1.99999999999683593814072199830603E-01, |
66 | r6 = -1.66666666666065494878165510225378E-01, |
67 | r7 = 1.42857517857114380606360005067609E-01, |
68 | r8 = -1.25000449999974370683775964001702E-01, |
69 | s3 = 0.333251953125000000e0, |
70 | ss3 = 8.138020833333333333e-05, |
71 | s4 = -2.500000000000000000e-01, |
72 | s5 = 1.999999999999960937e-01, |
73 | s6 = -1.666666666666592447e-01, |
74 | s7 = 1.428571845238194705e-01, |
75 | s8 = -1.250000500000149097e-01; |
76 | #endif |
77 | |