]>
Commit | Line | Data |
---|---|---|
e4d82761 UD |
1 | /* |
2 | * IBM Accurate Mathematical Library | |
c6c6dd48 | 3 | * Written by International Business Machines Corp. |
b168057a | 4 | * Copyright (C) 2001-2015 Free Software Foundation, Inc. |
e4d82761 UD |
5 | * |
6 | * This program is free software; you can redistribute it and/or modify | |
c6c6dd48 | 7 | * it under the terms of the GNU Lesser General Public License as published by |
cc7375ce | 8 | * the Free Software Foundation; either version 2.1 of the License, or |
e4d82761 | 9 | * (at your option) any later version. |
50944bca | 10 | * |
e4d82761 UD |
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 | |
c6c6dd48 | 14 | * GNU Lesser General Public License for more details. |
e4d82761 UD |
15 | * |
16 | * You should have received a copy of the GNU Lesser General Public License | |
59ba27a6 | 17 | * along with this program; if not, see <http://www.gnu.org/licenses/>. |
e4d82761 | 18 | */ |
c6c6dd48 | 19 | |
e4d82761 UD |
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 | |
a2d61ef8 | 34 | /**/ nZERO = {{0x80000000, 0}}, /* -0.0 */ |
50944bca | 35 | /**/ INF = {{0x7ff00000, 0x00000000}}, /* INF */ |
a2d61ef8 | 36 | /**/ nINF = {{0xfff00000, 0x00000000}}, /* -INF */ |
50944bca UD |
37 | /**/ sqrt_2 = {{0x3ff6a09e, 0x667f3bcc}}, /* sqrt(2) */ |
38 | /**/ ln2a = {{0x3fe62e42, 0xfefa3800}}, /* ln(2) 43 bits */ | |
39 | /**/ ln2b = {{0x3d2ef357, 0x93c76730}}, /* ln(2)-ln2a */ | |
40 | /**/ bigu = {{0x4297ffff, 0xfffffd2c}}, /* 1.5*2**42 -724*2**-10 */ | |
41 | /**/ bigv = {{0x4207ffff, 0xfff8016a}}, /* 1.5*2**33-1+362*2**-19 */ | |
42 | /**/ t52 = {{0x43300000, 0x00000000}}, /* 2**52 */ | |
43 | /**/ two52e = {{0x43300000, 0x000003ff}}; /* 2**52' */ | |
e4d82761 | 44 | |
e4d82761 UD |
45 | #else |
46 | #ifdef LITTLE_ENDI | |
47 | const static mynumber | |
a2d61ef8 | 48 | /**/ nZERO = {{0, 0x80000000}}, /* -0.0 */ |
50944bca | 49 | /**/ INF = {{0x00000000, 0x7ff00000}}, /* INF */ |
a2d61ef8 | 50 | /**/ nINF = {{0x00000000, 0xfff00000}}, /* -INF */ |
50944bca UD |
51 | /**/ sqrt_2 = {{0x667f3bcc, 0x3ff6a09e}}, /* sqrt(2) */ |
52 | /**/ ln2a = {{0xfefa3800, 0x3fe62e42}}, /* ln(2) 43 bits */ | |
53 | /**/ ln2b = {{0x93c76730, 0x3d2ef357}}, /* ln(2)-ln2a */ | |
54 | /**/ bigu = {{0xfffffd2c, 0x4297ffff}}, /* 1.5*2**42 -724*2**-10 */ | |
55 | /**/ bigv = {{0xfff8016a, 0x4207ffff}}, /* 1.5*2**33-1+362*2**-19 */ | |
56 | /**/ t52 = {{0x00000000, 0x43300000}}, /* 2**52 */ | |
57 | /**/ two52e = {{0x000003ff, 0x43300000}}; /* 2**52' */ | |
e4d82761 | 58 | |
e4d82761 UD |
59 | #endif |
60 | #endif | |
61 | ||
62 | const static double p2=-0.5, p3 = 3.3333333333333333333e-1, p4 = -0.25, | |
63 | q2 = -0.5, q3 = 3.3333333333331404e-01, q4 = -2.4999999999996436e-01, | |
64 | q5 = 2.0000010500004459e-01, q6 = -1.6666678916688004e-01, | |
65 | r3 = 3.33333333333333333372884096563030E-01, | |
66 | r4 = -2.50000000000000000213574153875908E-01, | |
67 | r5 = 1.99999999999683593814072199830603E-01, | |
68 | r6 = -1.66666666666065494878165510225378E-01, | |
69 | r7 = 1.42857517857114380606360005067609E-01, | |
70 | r8 = -1.25000449999974370683775964001702E-01, | |
71 | s3 = 0.333251953125000000e0, | |
72 | ss3 = 8.138020833333333333e-05, | |
73 | s4 = -2.500000000000000000e-01, | |
74 | s5 = 1.999999999999960937e-01, | |
75 | s6 = -1.666666666666592447e-01, | |
76 | s7 = 1.428571845238194705e-01, | |
77 | s8 = -1.250000500000149097e-01; | |
78 | #endif |