1 /* 2 * Copyright (c) 2008-2016 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 <stdio.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 *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 *mpd_clamp_string[MPD_CLAMP_GUARD] = { 128 "CLAMP_DEFAULT", 129 "CLAMP_IEEE_754" 130 }; 131 132 133