1 /*
2 * Copyright (c) 2008-2020 Stefan Krah. All rights reserved.
3 *
4 * Redistribution and use in source and binary forms, with or without
5 * modification, are permitted provided that the following conditions
6 * are met:
7 *
8 * 1. Redistributions of source code must retain the above copyright
9 * notice, this list of conditions and the following disclaimer.
10 *
11 * 2. Redistributions in binary form must reproduce the above copyright
12 * notice, this list of conditions and the following disclaimer in the
13 * documentation and/or other materials provided with the distribution.
14 *
15 * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS "AS IS" AND
16 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
17 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
18 * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
19 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
20 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
21 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
22 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
23 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
24 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
25 * SUCH DAMAGE.
26 */
27
28
29 #include "mpdecimal.h"
30 #include "basearith.h"
31 #include "constants.h"
32
33
34 #if defined(CONFIG_64)
35
36 /* number-theory.c */
37 const mpd_uint_t mpd_moduli[3] = {
38 18446744069414584321ULL, 18446744056529682433ULL, 18446742974197923841ULL
39 };
40 const mpd_uint_t mpd_roots[3] = {7ULL, 10ULL, 19ULL};
41
42 /* crt.c */
43 const mpd_uint_t INV_P1_MOD_P2 = 18446744055098026669ULL;
44 const mpd_uint_t INV_P1P2_MOD_P3 = 287064143708160ULL;
45 const mpd_uint_t LH_P1P2 = 18446744052234715137ULL; /* (P1*P2) % 2^64 */
46 const mpd_uint_t UH_P1P2 = 18446744052234715141ULL; /* (P1*P2) / 2^64 */
47
48 /* transpose.c */
49 const mpd_size_t mpd_bits[64] = {
50 1, 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, 2048, 4096, 8192, 16384,
51 32768, 65536, 131072, 262144, 524288, 1048576, 2097152, 4194304, 8388608,
52 16777216, 33554432, 67108864, 134217728, 268435456, 536870912, 1073741824,
53 2147483648ULL, 4294967296ULL, 8589934592ULL, 17179869184ULL, 34359738368ULL,
54 68719476736ULL, 137438953472ULL, 274877906944ULL, 549755813888ULL,
55 1099511627776ULL, 2199023255552ULL, 4398046511104, 8796093022208ULL,
56 17592186044416ULL, 35184372088832ULL, 70368744177664ULL, 140737488355328ULL,
57 281474976710656ULL, 562949953421312ULL, 1125899906842624ULL,
58 2251799813685248ULL, 4503599627370496ULL, 9007199254740992ULL,
59 18014398509481984ULL, 36028797018963968ULL, 72057594037927936ULL,
60 144115188075855872ULL, 288230376151711744ULL, 576460752303423488ULL,
61 1152921504606846976ULL, 2305843009213693952ULL, 4611686018427387904ULL,
62 9223372036854775808ULL
63 };
64
65 /* mpdecimal.c */
66 const mpd_uint_t mpd_pow10[MPD_RDIGITS+1] = {
67 1,10,100,1000,10000,100000,1000000,10000000,100000000,1000000000,
68 10000000000ULL,100000000000ULL,1000000000000ULL,10000000000000ULL,
69 100000000000000ULL,1000000000000000ULL,10000000000000000ULL,
70 100000000000000000ULL,1000000000000000000ULL,10000000000000000000ULL
71 };
72
73 /* magic number for constant division by MPD_RADIX */
74 const mpd_uint_t mprime_rdx = 15581492618384294730ULL;
75
76 #elif defined(CONFIG_32)
77
78 /* number-theory.c */
79 const mpd_uint_t mpd_moduli[3] = {2113929217UL, 2013265921UL, 1811939329UL};
80 const mpd_uint_t mpd_roots[3] = {5UL, 31UL, 13UL};
81
82 /* PentiumPro modular multiplication: These constants have to be loaded as
83 * 80 bit long doubles, which are not supported by certain compilers. */
84 const uint32_t mpd_invmoduli[3][3] = {
85 {4293885170U, 2181570688U, 16352U}, /* ((long double) 1 / 2113929217UL) */
86 {1698898177U, 2290649223U, 16352U}, /* ((long double) 1 / 2013265921UL) */
87 {2716021846U, 2545165803U, 16352U} /* ((long double) 1 / 1811939329UL) */
88 };
89
90 const float MPD_TWO63 = 9223372036854775808.0; /* 2^63 */
91
92 /* crt.c */
93 const mpd_uint_t INV_P1_MOD_P2 = 2013265901UL;
94 const mpd_uint_t INV_P1P2_MOD_P3 = 54UL;
95 const mpd_uint_t LH_P1P2 = 4127195137UL; /* (P1*P2) % 2^32 */
96 const mpd_uint_t UH_P1P2 = 990904320UL; /* (P1*P2) / 2^32 */
97
98 /* transpose.c */
99 const mpd_size_t mpd_bits[32] = {
100 1, 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, 2048, 4096, 8192, 16384,
101 32768, 65536, 131072, 262144, 524288, 1048576, 2097152, 4194304, 8388608,
102 16777216, 33554432, 67108864, 134217728, 268435456, 536870912, 1073741824,
103 2147483648UL
104 };
105
106 /* mpdecimal.c */
107 const mpd_uint_t mpd_pow10[MPD_RDIGITS+1] = {
108 1,10,100,1000,10000,100000,1000000,10000000,100000000,1000000000
109 };
110
111 #else
112 #error "CONFIG_64 or CONFIG_32 must be defined."
113 #endif
114
115 const char * const mpd_round_string[MPD_ROUND_GUARD] = {
116 "ROUND_UP", /* round away from 0 */
117 "ROUND_DOWN", /* round toward 0 (truncate) */
118 "ROUND_CEILING", /* round toward +infinity */
119 "ROUND_FLOOR", /* round toward -infinity */
120 "ROUND_HALF_UP", /* 0.5 is rounded up */
121 "ROUND_HALF_DOWN", /* 0.5 is rounded down */
122 "ROUND_HALF_EVEN", /* 0.5 is rounded to even */
123 "ROUND_05UP", /* round zero or five away from 0 */
124 "ROUND_TRUNC", /* truncate, but set infinity */
125 };
126
127 const char * const mpd_clamp_string[MPD_CLAMP_GUARD] = {
128 "CLAMP_DEFAULT",
129 "CLAMP_IEEE_754"
130 };