Branch data     Line data    Source code

       1                 :            : // Copyright (c) 2009-2010 Satoshi Nakamoto
       2                 :            : // Copyright (c) 2009-2022 The Bitcoin Core developers
       3                 :            : // Distributed under the MIT software license, see the accompanying
       4                 :            : // file COPYING or http://www.opensource.org/licenses/mit-license.php.
       5                 :            : 
       6                 :            : #include <arith_uint256.h>
       7                 :            : 
       8                 :            : #include <uint256.h>
       9                 :            : #include <crypto/common.h>
      10                 :            : 
      11                 :            : 
      12                 :            : template <unsigned int BITS>
      13                 :          0 : base_uint<BITS>::base_uint(const std::string& str)
      14                 :            : {
      15                 :            :     static_assert(BITS/32 > 0 && BITS%32 == 0, "Template parameter BITS must be a positive multiple of 32.");
      16                 :            : 
      17                 :          0 :     SetHex(str);
      18                 :          0 : }
      19                 :            : 
      20                 :            : template <unsigned int BITS>
      21                 :       1042 : base_uint<BITS>& base_uint<BITS>::operator<<=(unsigned int shift)
      22                 :            : {
      23                 :       1042 :     base_uint<BITS> a(*this);
      24         [ +  + ]:       9378 :     for (int i = 0; i < WIDTH; i++)
      25                 :       8336 :         pn[i] = 0;
      26                 :       1042 :     int k = shift / 32;
      27                 :       1042 :     shift = shift % 32;
      28         [ +  + ]:       9378 :     for (int i = 0; i < WIDTH; i++) {
      29 [ +  + ][ -  + ]:       8336 :         if (i + k + 1 < WIDTH && shift != 0)
      30                 :       1057 :             pn[i + k + 1] |= (a.pn[i] >> (32 - shift));
      31         [ +  + ]:       8336 :         if (i + k < WIDTH)
      32                 :       2099 :             pn[i + k] |= (a.pn[i] << shift);
      33                 :       8336 :     }
      34                 :       1042 :     return *this;
      35                 :            : }
      36                 :            : 
      37                 :            : template <unsigned int BITS>
      38                 :        902 : base_uint<BITS>& base_uint<BITS>::operator>>=(unsigned int shift)
      39                 :            : {
      40                 :        902 :     base_uint<BITS> a(*this);
      41         [ +  + ]:       8118 :     for (int i = 0; i < WIDTH; i++)
      42                 :       7216 :         pn[i] = 0;
      43                 :        902 :     int k = shift / 32;
      44                 :        902 :     shift = shift % 32;
      45         [ +  + ]:       8118 :     for (int i = 0; i < WIDTH; i++) {
      46 [ +  + ][ -  + ]:       7216 :         if (i - k - 1 >= 0 && shift != 0)
      47                 :       2114 :             pn[i - k - 1] |= (a.pn[i] << (32 - shift));
      48         [ +  + ]:       7216 :         if (i - k >= 0)
      49                 :       3016 :             pn[i - k] |= (a.pn[i] >> shift);
      50                 :       7216 :     }
      51                 :        902 :     return *this;
      52                 :            : }
      53                 :            : 
      54                 :            : template <unsigned int BITS>
      55                 :          0 : base_uint<BITS>& base_uint<BITS>::operator*=(uint32_t b32)
      56                 :            : {
      57                 :          0 :     uint64_t carry = 0;
      58         [ #  # ]:          0 :     for (int i = 0; i < WIDTH; i++) {
      59                 :          0 :         uint64_t n = carry + (uint64_t)b32 * pn[i];
      60                 :          0 :         pn[i] = n & 0xffffffff;
      61                 :          0 :         carry = n >> 32;
      62                 :          0 :     }
      63                 :          0 :     return *this;
      64                 :            : }
      65                 :            : 
      66                 :            : template <unsigned int BITS>
      67                 :          0 : base_uint<BITS>& base_uint<BITS>::operator*=(const base_uint& b)
      68                 :            : {
      69                 :          0 :     base_uint<BITS> a;
      70         [ #  # ]:          0 :     for (int j = 0; j < WIDTH; j++) {
      71                 :          0 :         uint64_t carry = 0;
      72         [ #  # ]:          0 :         for (int i = 0; i + j < WIDTH; i++) {
      73                 :          0 :             uint64_t n = carry + a.pn[i + j] + (uint64_t)pn[j] * b.pn[i];
      74                 :          0 :             a.pn[i + j] = n & 0xffffffff;
      75                 :          0 :             carry = n >> 32;
      76                 :          0 :         }
      77                 :          0 :     }
      78                 :          0 :     *this = a;
      79                 :          0 :     return *this;
      80                 :            : }
      81                 :            : 
      82                 :            : template <unsigned int BITS>
      83                 :        151 : base_uint<BITS>& base_uint<BITS>::operator/=(const base_uint& b)
      84                 :            : {
      85                 :        151 :     base_uint<BITS> div = b;     // make a copy, so we can shift.
      86                 :        151 :     base_uint<BITS> num = *this; // make a copy, so we can subtract.
      87                 :        151 :     *this = 0;                   // the quotient.
      88                 :        151 :     int num_bits = num.bits();
      89                 :        151 :     int div_bits = div.bits();
      90         [ +  - ]:        151 :     if (div_bits == 0)
      91 [ #  # ][ #  # ]:          0 :         throw uint_error("Division by zero");
         [ #  # ][ #  # ]
      92         [ -  + ]:        151 :     if (div_bits > num_bits) // the result is certainly 0.
      93                 :          0 :         return *this;
      94                 :        151 :     int shift = num_bits - div_bits;
      95                 :        151 :     div <<= shift; // shift so that div and num align.
      96         [ +  + ]:        453 :     while (shift >= 0) {
      97         [ +  + ]:        302 :         if (num >= div) {
      98                 :        151 :             num -= div;
      99                 :        151 :             pn[shift / 32] |= (1U << (shift & 31)); // set a bit of the result.
     100                 :        151 :         }
     101                 :        302 :         div >>= 1; // shift back.
     102                 :        302 :         shift--;
     103                 :            :     }
     104                 :            :     // num now contains the remainder of the division.
     105                 :        151 :     return *this;
     106                 :        151 : }
     107                 :            : 
     108                 :            : template <unsigned int BITS>
     109                 :     210299 : int base_uint<BITS>::CompareTo(const base_uint<BITS>& b) const
     110                 :            : {
     111         [ +  + ]:    1673527 :     for (int i = WIDTH - 1; i >= 0; i--) {
     112         [ +  + ]:    1669918 :         if (pn[i] < b.pn[i])
     113                 :     137544 :             return -1;
     114         [ +  + ]:    1532374 :         if (pn[i] > b.pn[i])
     115                 :      69146 :             return 1;
     116                 :    1463228 :     }
     117                 :       3609 :     return 0;
     118                 :     210299 : }
     119                 :            : 
     120                 :            : template <unsigned int BITS>
     121                 :        891 : bool base_uint<BITS>::EqualTo(uint64_t b) const
     122                 :            : {
     123         [ +  - ]:        891 :     for (int i = WIDTH - 1; i >= 2; i--) {
     124         [ -  + ]:        891 :         if (pn[i])
     125                 :        891 :             return false;
     126                 :          0 :     }
     127         [ #  # ]:          0 :     if (pn[1] != (b >> 32))
     128                 :          0 :         return false;
     129         [ #  # ]:          0 :     if (pn[0] != (b & 0xfffffffful))
     130                 :          0 :         return false;
     131                 :          0 :     return true;
     132                 :        891 : }
     133                 :            : 
     134                 :            : template <unsigned int BITS>
     135                 :        151 : double base_uint<BITS>::getdouble() const
     136                 :            : {
     137                 :        151 :     double ret = 0.0;
     138                 :        151 :     double fact = 1.0;
     139         [ +  + ]:       1359 :     for (int i = 0; i < WIDTH; i++) {
     140                 :       1208 :         ret += fact * pn[i];
     141                 :       1208 :         fact *= 4294967296.0;
     142                 :       1208 :     }
     143                 :        151 :     return ret;
     144                 :            : }
     145                 :            : 
     146                 :            : template <unsigned int BITS>
     147                 :          1 : std::string base_uint<BITS>::GetHex() const
     148                 :            : {
     149                 :          1 :     base_blob<BITS> b;
     150         [ +  + ]:          9 :     for (int x = 0; x < this->WIDTH; ++x) {
     151                 :          8 :         WriteLE32(b.begin() + x*4, this->pn[x]);
     152                 :          8 :     }
     153                 :          1 :     return b.GetHex();
     154                 :            : }
     155                 :            : 
     156                 :            : template <unsigned int BITS>
     157                 :          0 : void base_uint<BITS>::SetHex(const char* psz)
     158                 :            : {
     159                 :          0 :     base_blob<BITS> b;
     160                 :          0 :     b.SetHex(psz);
     161         [ #  # ]:          0 :     for (int x = 0; x < this->WIDTH; ++x) {
     162                 :          0 :         this->pn[x] = ReadLE32(b.begin() + x*4);
     163                 :          0 :     }
     164                 :          0 : }
     165                 :            : 
     166                 :            : template <unsigned int BITS>
     167                 :          0 : void base_uint<BITS>::SetHex(const std::string& str)
     168                 :            : {
     169                 :          0 :     SetHex(str.c_str());
     170                 :          0 : }
     171                 :            : 
     172                 :            : template <unsigned int BITS>
     173                 :          0 : std::string base_uint<BITS>::ToString() const
     174                 :            : {
     175                 :          0 :     return GetHex();
     176                 :            : }
     177                 :            : 
     178                 :            : template <unsigned int BITS>
     179                 :        902 : unsigned int base_uint<BITS>::bits() const
     180                 :            : {
     181         [ +  - ]:        902 :     for (int pos = WIDTH - 1; pos >= 0; pos--) {
     182         [ -  + ]:        902 :         if (pn[pos]) {
     183         [ +  - ]:       1653 :             for (int nbits = 31; nbits > 0; nbits--) {
     184         [ +  + ]:       1653 :                 if (pn[pos] & 1U << nbits)
     185                 :        902 :                     return 32 * pos + nbits + 1;
     186                 :        751 :             }
     187                 :          0 :             return 32 * pos + 1;
     188                 :            :         }
     189                 :          0 :     }
     190                 :          0 :     return 0;
     191                 :        902 : }
     192                 :            : 
     193                 :            : // Explicit instantiations for base_uint<256>
     194                 :            : template class base_uint<256>;
     195                 :            : 
     196                 :            : // This implementation directly uses shifts instead of going
     197                 :            : // through an intermediate MPI representation.
     198                 :        891 : arith_uint256& arith_uint256::SetCompact(uint32_t nCompact, bool* pfNegative, bool* pfOverflow)
     199                 :            : {
     200                 :        891 :     int nSize = nCompact >> 24;
     201                 :        891 :     uint32_t nWord = nCompact & 0x007fffff;
     202         [ -  + ]:        891 :     if (nSize <= 3) {
     203                 :          0 :         nWord >>= 8 * (3 - nSize);
     204                 :          0 :         *this = nWord;
     205                 :          0 :     } else {
     206                 :        891 :         *this = nWord;
     207                 :        891 :         *this <<= 8 * (nSize - 3);
     208                 :            :     }
     209         [ -  + ]:        891 :     if (pfNegative)
     210         [ -  + ]:        891 :         *pfNegative = nWord != 0 && (nCompact & 0x00800000) != 0;
     211         [ -  + ]:        891 :     if (pfOverflow)
     212 [ -  + ][ +  - ]:       1782 :         *pfOverflow = nWord != 0 && ((nSize > 34) ||
     213 [ +  - ][ -  + ]:       1782 :                                      (nWord > 0xff && nSize > 33) ||
     214         [ -  + ]:        891 :                                      (nWord > 0xffff && nSize > 32));
     215                 :        891 :     return *this;
     216                 :            : }
     217                 :            : 
     218                 :        600 : uint32_t arith_uint256::GetCompact(bool fNegative) const
     219                 :            : {
     220                 :        600 :     int nSize = (bits() + 7) / 8;
     221                 :        600 :     uint32_t nCompact = 0;
     222         [ -  + ]:        600 :     if (nSize <= 3) {
     223                 :          0 :         nCompact = GetLow64() << 8 * (3 - nSize);
     224                 :          0 :     } else {
     225                 :        600 :         arith_uint256 bn = *this >> 8 * (nSize - 3);
     226                 :        600 :         nCompact = bn.GetLow64();
     227                 :            :     }
     228                 :            :     // The 0x00800000 bit denotes the sign.
     229                 :            :     // Thus, if it is already set, divide the mantissa by 256 and increase the exponent.
     230         [ +  - ]:        600 :     if (nCompact & 0x00800000) {
     231                 :          0 :         nCompact >>= 8;
     232                 :          0 :         nSize++;
     233                 :          0 :     }
     234         [ +  - ]:        600 :     assert((nCompact & ~0x007fffffU) == 0);
     235         [ -  + ]:        600 :     assert(nSize < 256);
     236                 :        600 :     nCompact |= nSize << 24;
     237         [ +  - ]:        600 :     nCompact |= (fNegative && (nCompact & 0x007fffff) ? 0x00800000 : 0);
     238                 :        600 :     return nCompact;
     239                 :            : }
     240                 :            : 
     241                 :          0 : uint256 ArithToUint256(const arith_uint256 &a)
     242                 :            : {
     243                 :          0 :     uint256 b;
     244         [ #  # ]:          0 :     for(int x=0; x<a.WIDTH; ++x)
     245                 :          0 :         WriteLE32(b.begin() + x*4, a.pn[x]);
     246                 :          0 :     return b;
     247                 :            : }
     248                 :       2082 : arith_uint256 UintToArith256(const uint256 &a)
     249                 :            : {
     250                 :       2082 :     arith_uint256 b;
     251         [ +  + ]:      18738 :     for(int x=0; x<b.WIDTH; ++x)
     252                 :      16656 :         b.pn[x] = ReadLE32(a.begin() + x*4);
     253                 :       2082 :     return b;
     254                 :            : }