Code cleanup
[novacoin.git] / src / uint256.h
1 // Copyright (c) 2009-2010 Satoshi Nakamoto
2 // Copyright (c) 2009-2012 The Bitcoin developers
3 // Distributed under the MIT/X11 software license, see the accompanying
4 // file COPYING or http://www.opensource.org/licenses/mit-license.php.
5 #ifndef BITCOIN_UINT256_H
6 #define BITCOIN_UINT256_H
7
8 #include <limits.h>
9 #include <stdio.h>
10 #include <string.h>
11 #include <string>
12 #include <vector>
13 #include <stdint.h>
14
15 inline int Testuint256AdHoc(std::vector<std::string> vArg);
16
17
18
19 /** Base class without constructors for uint256 and uint160.
20  * This makes the compiler let u use it in a union.
21  */
22 template<unsigned int BITS>
23 class base_uint
24 {
25 protected:
26     enum { WIDTH=BITS/32 };
27     unsigned int pn[WIDTH];
28 public:
29
30     bool operator!() const
31     {
32         for (int i = 0; i < WIDTH; i++)
33             if (pn[i] != 0)
34                 return false;
35         return true;
36     }
37
38     const base_uint operator~() const
39     {
40         base_uint ret;
41         for (int i = 0; i < WIDTH; i++)
42             ret.pn[i] = ~pn[i];
43         return ret;
44     }
45
46     const base_uint operator-() const
47     {
48         base_uint ret;
49         for (int i = 0; i < WIDTH; i++)
50             ret.pn[i] = ~pn[i];
51         ret++;
52         return ret;
53     }
54
55     double getdouble() const
56     {
57         double ret = 0.0;
58         double fact = 1.0;
59         for (int i = 0; i < WIDTH; i++) {
60             ret += fact * pn[i];
61             fact *= 4294967296.0;
62         }
63         return ret;
64     }
65
66     base_uint& operator=(uint64_t b)
67     {
68         pn[0] = (unsigned int)b;
69         pn[1] = (unsigned int)(b >> 32);
70         for (int i = 2; i < WIDTH; i++)
71             pn[i] = 0;
72         return *this;
73     }
74
75     base_uint& operator^=(const base_uint& b)
76     {
77         for (int i = 0; i < WIDTH; i++)
78             pn[i] ^= b.pn[i];
79         return *this;
80     }
81
82     base_uint& operator&=(const base_uint& b)
83     {
84         for (int i = 0; i < WIDTH; i++)
85             pn[i] &= b.pn[i];
86         return *this;
87     }
88
89     base_uint& operator|=(const base_uint& b)
90     {
91         for (int i = 0; i < WIDTH; i++)
92             pn[i] |= b.pn[i];
93         return *this;
94     }
95
96     base_uint& operator^=(uint64_t b)
97     {
98         pn[0] ^= (unsigned int)b;
99         pn[1] ^= (unsigned int)(b >> 32);
100         return *this;
101     }
102
103     base_uint& operator|=(uint64_t b)
104     {
105         pn[0] |= (unsigned int)b;
106         pn[1] |= (unsigned int)(b >> 32);
107         return *this;
108     }
109
110     base_uint& operator<<=(unsigned int shift)
111     {
112         base_uint a(*this);
113         for (int i = 0; i < WIDTH; i++)
114             pn[i] = 0;
115         int k = shift / 32;
116         shift = shift % 32;
117         for (int i = 0; i < WIDTH; i++)
118         {
119             if (i+k+1 < WIDTH && shift != 0)
120                 pn[i+k+1] |= (a.pn[i] >> (32-shift));
121             if (i+k < WIDTH)
122                 pn[i+k] |= (a.pn[i] << shift);
123         }
124         return *this;
125     }
126
127     base_uint& operator>>=(unsigned int shift)
128     {
129         base_uint a(*this);
130         for (int i = 0; i < WIDTH; i++)
131             pn[i] = 0;
132         int k = shift / 32;
133         shift = shift % 32;
134         for (int i = 0; i < WIDTH; i++)
135         {
136             if (i-k-1 >= 0 && shift != 0)
137                 pn[i-k-1] |= (a.pn[i] << (32-shift));
138             if (i-k >= 0)
139                 pn[i-k] |= (a.pn[i] >> shift);
140         }
141         return *this;
142     }
143
144     base_uint& operator+=(const base_uint& b)
145     {
146         uint64_t carry = 0;
147         for (int i = 0; i < WIDTH; i++)
148         {
149             uint64_t n = carry + pn[i] + b.pn[i];
150             pn[i] = n & 0xffffffff;
151             carry = n >> 32;
152         }
153         return *this;
154     }
155
156     base_uint& operator-=(const base_uint& b)
157     {
158         *this += -b;
159         return *this;
160     }
161
162     base_uint& operator+=(uint64_t b64)
163     {
164         base_uint b;
165         b = b64;
166         *this += b;
167         return *this;
168     }
169
170     base_uint& operator-=(uint64_t b64)
171     {
172         base_uint b;
173         b = b64;
174         *this += -b;
175         return *this;
176     }
177
178
179     base_uint& operator++()
180     {
181         // prefix operator
182         int i = 0;
183         while (++pn[i] == 0 && i < WIDTH-1)
184             i++;
185         return *this;
186     }
187
188     const base_uint operator++(int)
189     {
190         // postfix operator
191         const base_uint ret = *this;
192         ++(*this);
193         return ret;
194     }
195
196     base_uint& operator--()
197     {
198         // prefix operator
199         int i = 0;
200         while (--pn[i] == -1 && i < WIDTH-1)
201             i++;
202         return *this;
203     }
204
205     const base_uint operator--(int)
206     {
207         // postfix operator
208         const base_uint ret = *this;
209         --(*this);
210         return ret;
211     }
212
213
214     friend inline bool operator<(const base_uint& a, const base_uint& b)
215     {
216         for (int i = base_uint::WIDTH-1; i >= 0; i--)
217         {
218             if (a.pn[i] < b.pn[i])
219                 return true;
220             else if (a.pn[i] > b.pn[i])
221                 return false;
222         }
223         return false;
224     }
225
226     friend inline bool operator<=(const base_uint& a, const base_uint& b)
227     {
228         for (int i = base_uint::WIDTH-1; i >= 0; i--)
229         {
230             if (a.pn[i] < b.pn[i])
231                 return true;
232             else if (a.pn[i] > b.pn[i])
233                 return false;
234         }
235         return true;
236     }
237
238     friend inline bool operator>(const base_uint& a, const base_uint& b)
239     {
240         for (int i = base_uint::WIDTH-1; i >= 0; i--)
241         {
242             if (a.pn[i] > b.pn[i])
243                 return true;
244             else if (a.pn[i] < b.pn[i])
245                 return false;
246         }
247         return false;
248     }
249
250     friend inline bool operator>=(const base_uint& a, const base_uint& b)
251     {
252         for (int i = base_uint::WIDTH-1; i >= 0; i--)
253         {
254             if (a.pn[i] > b.pn[i])
255                 return true;
256             else if (a.pn[i] < b.pn[i])
257                 return false;
258         }
259         return true;
260     }
261
262     friend inline bool operator==(const base_uint& a, const base_uint& b)
263     {
264         for (int i = 0; i < base_uint::WIDTH; i++)
265             if (a.pn[i] != b.pn[i])
266                 return false;
267         return true;
268     }
269
270     friend inline bool operator==(const base_uint& a, uint64_t b)
271     {
272         if (a.pn[0] != (unsigned int)b)
273             return false;
274         if (a.pn[1] != (unsigned int)(b >> 32))
275             return false;
276         for (int i = 2; i < base_uint::WIDTH; i++)
277             if (a.pn[i] != 0)
278                 return false;
279         return true;
280     }
281
282     friend inline bool operator!=(const base_uint& a, const base_uint& b)
283     {
284         return (!(a == b));
285     }
286
287     friend inline bool operator!=(const base_uint& a, uint64_t b)
288     {
289         return (!(a == b));
290     }
291
292
293
294     std::string GetHex() const
295     {
296         char psz[sizeof(pn)*2 + 1];
297         for (unsigned int i = 0; i < sizeof(pn); i++)
298             sprintf(psz + i*2, "%02x", ((unsigned char*)pn)[sizeof(pn) - i - 1]);
299         return std::string(psz, psz + sizeof(pn)*2);
300     }
301
302     void SetHex(const char* psz)
303     {
304         for (int i = 0; i < WIDTH; i++)
305             pn[i] = 0;
306
307         // skip leading spaces
308         while (isspace(*psz))
309             psz++;
310
311         // skip 0x
312         if (psz[0] == '0' && tolower(psz[1]) == 'x')
313             psz += 2;
314
315         // hex string to uint
316         static const unsigned char phexdigit[256] = { 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,1,2,3,4,5,6,7,8,9,0,0,0,0,0,0, 0,0xa,0xb,0xc,0xd,0xe,0xf,0,0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,0xa,0xb,0xc,0xd,0xe,0xf,0,0,0,0,0,0,0,0,0 };
317         const char* pbegin = psz;
318         while (phexdigit[(unsigned char)*psz] || *psz == '0')
319             psz++;
320         psz--;
321         unsigned char* p1 = (unsigned char*)pn;
322         unsigned char* pend = p1 + WIDTH * 4;
323         while (psz >= pbegin && p1 < pend)
324         {
325             *p1 = phexdigit[(unsigned char)*psz--];
326             if (psz >= pbegin)
327             {
328                 *p1 |= (phexdigit[(unsigned char)*psz--] << 4);
329                 p1++;
330             }
331         }
332     }
333
334     void SetHex(const std::string& str)
335     {
336         SetHex(str.c_str());
337     }
338
339     std::string ToString() const
340     {
341         return (GetHex());
342     }
343
344     unsigned char* begin()
345     {
346         return (unsigned char*)&pn[0];
347     }
348
349     unsigned char* end()
350     {
351         return (unsigned char*)&pn[WIDTH];
352     }
353
354     std::vector<unsigned char> getBytes()
355     {
356         return std::vector<unsigned char>(begin(), end());
357     }
358
359     unsigned int size()
360     {
361         return sizeof(pn);
362     }
363
364     uint64_t Get64(int n=0) const
365     {
366         return pn[2*n] | (uint64_t)pn[2*n+1] << 32;
367     }
368
369     unsigned int GetSerializeSize(int nType, int nVersion) const
370     {
371         return sizeof(pn);
372     }
373
374     template<typename Stream>
375     void Serialize(Stream& s, int nType, int nVersion) const
376     {
377         s.write((char*)pn, sizeof(pn));
378     }
379
380     template<typename Stream>
381     void Unserialize(Stream& s, int nType, int nVersion)
382     {
383         s.read((char*)pn, sizeof(pn));
384     }
385
386     friend class uint160;
387     friend class uint256;
388     friend inline int Testuint256AdHoc(std::vector<std::string> vArg);
389 };
390
391 typedef base_uint<160> base_uint160;
392 typedef base_uint<256> base_uint256;
393
394 //
395 // uint160 and uint256 could be implemented as templates, but to keep
396 // compile errors and debugging cleaner, they're copy and pasted.
397 //
398
399
400
401 //////////////////////////////////////////////////////////////////////////////
402 //
403 // uint160
404 //
405
406 /** 160-bit unsigned integer */
407 class uint160 : public base_uint160
408 {
409 public:
410     typedef base_uint160 basetype;
411
412     uint160()
413     {
414         for (int i = 0; i < WIDTH; i++)
415             pn[i] = 0;
416     }
417
418     uint160(const basetype& b)
419     {
420         for (int i = 0; i < WIDTH; i++)
421             pn[i] = b.pn[i];
422     }
423
424     uint160& operator=(const basetype& b)
425     {
426         for (int i = 0; i < WIDTH; i++)
427             pn[i] = b.pn[i];
428         return *this;
429     }
430
431     uint160(uint64_t b)
432     {
433         pn[0] = (unsigned int)b;
434         pn[1] = (unsigned int)(b >> 32);
435         for (int i = 2; i < WIDTH; i++)
436             pn[i] = 0;
437     }
438
439     uint160& operator=(uint64_t b)
440     {
441         pn[0] = (unsigned int)b;
442         pn[1] = (unsigned int)(b >> 32);
443         for (int i = 2; i < WIDTH; i++)
444             pn[i] = 0;
445         return *this;
446     }
447
448     explicit uint160(const std::string& str)
449     {
450         SetHex(str);
451     }
452
453     explicit uint160(const std::vector<unsigned char>& vch)
454     {
455         if (vch.size() == sizeof(pn))
456             memcpy(pn, &vch[0], sizeof(pn));
457         else
458             *this = 0;
459     }
460 };
461
462 inline bool operator==(const uint160& a, uint64_t b)                           { return (base_uint160)a == b; }
463 inline bool operator!=(const uint160& a, uint64_t b)                           { return (base_uint160)a != b; }
464 inline const uint160 operator<<(const base_uint160& a, unsigned int shift)   { return uint160(a) <<= shift; }
465 inline const uint160 operator>>(const base_uint160& a, unsigned int shift)   { return uint160(a) >>= shift; }
466 inline const uint160 operator<<(const uint160& a, unsigned int shift)        { return uint160(a) <<= shift; }
467 inline const uint160 operator>>(const uint160& a, unsigned int shift)        { return uint160(a) >>= shift; }
468
469 inline const uint160 operator^(const base_uint160& a, const base_uint160& b) { return uint160(a) ^= b; }
470 inline const uint160 operator&(const base_uint160& a, const base_uint160& b) { return uint160(a) &= b; }
471 inline const uint160 operator|(const base_uint160& a, const base_uint160& b) { return uint160(a) |= b; }
472 inline const uint160 operator+(const base_uint160& a, const base_uint160& b) { return uint160(a) += b; }
473 inline const uint160 operator-(const base_uint160& a, const base_uint160& b) { return uint160(a) -= b; }
474
475 inline bool operator<(const base_uint160& a, const uint160& b)          { return (base_uint160)a <  (base_uint160)b; }
476 inline bool operator<=(const base_uint160& a, const uint160& b)         { return (base_uint160)a <= (base_uint160)b; }
477 inline bool operator>(const base_uint160& a, const uint160& b)          { return (base_uint160)a >  (base_uint160)b; }
478 inline bool operator>=(const base_uint160& a, const uint160& b)         { return (base_uint160)a >= (base_uint160)b; }
479 inline bool operator==(const base_uint160& a, const uint160& b)         { return (base_uint160)a == (base_uint160)b; }
480 inline bool operator!=(const base_uint160& a, const uint160& b)         { return (base_uint160)a != (base_uint160)b; }
481 inline const uint160 operator^(const base_uint160& a, const uint160& b) { return (base_uint160)a ^  (base_uint160)b; }
482 inline const uint160 operator&(const base_uint160& a, const uint160& b) { return (base_uint160)a &  (base_uint160)b; }
483 inline const uint160 operator|(const base_uint160& a, const uint160& b) { return (base_uint160)a |  (base_uint160)b; }
484 inline const uint160 operator+(const base_uint160& a, const uint160& b) { return (base_uint160)a +  (base_uint160)b; }
485 inline const uint160 operator-(const base_uint160& a, const uint160& b) { return (base_uint160)a -  (base_uint160)b; }
486
487 inline bool operator<(const uint160& a, const base_uint160& b)          { return (base_uint160)a <  (base_uint160)b; }
488 inline bool operator<=(const uint160& a, const base_uint160& b)         { return (base_uint160)a <= (base_uint160)b; }
489 inline bool operator>(const uint160& a, const base_uint160& b)          { return (base_uint160)a >  (base_uint160)b; }
490 inline bool operator>=(const uint160& a, const base_uint160& b)         { return (base_uint160)a >= (base_uint160)b; }
491 inline bool operator==(const uint160& a, const base_uint160& b)         { return (base_uint160)a == (base_uint160)b; }
492 inline bool operator!=(const uint160& a, const base_uint160& b)         { return (base_uint160)a != (base_uint160)b; }
493 inline const uint160 operator^(const uint160& a, const base_uint160& b) { return (base_uint160)a ^  (base_uint160)b; }
494 inline const uint160 operator&(const uint160& a, const base_uint160& b) { return (base_uint160)a &  (base_uint160)b; }
495 inline const uint160 operator|(const uint160& a, const base_uint160& b) { return (base_uint160)a |  (base_uint160)b; }
496 inline const uint160 operator+(const uint160& a, const base_uint160& b) { return (base_uint160)a +  (base_uint160)b; }
497 inline const uint160 operator-(const uint160& a, const base_uint160& b) { return (base_uint160)a -  (base_uint160)b; }
498
499 inline bool operator<(const uint160& a, const uint160& b)               { return (base_uint160)a <  (base_uint160)b; }
500 inline bool operator<=(const uint160& a, const uint160& b)              { return (base_uint160)a <= (base_uint160)b; }
501 inline bool operator>(const uint160& a, const uint160& b)               { return (base_uint160)a >  (base_uint160)b; }
502 inline bool operator>=(const uint160& a, const uint160& b)              { return (base_uint160)a >= (base_uint160)b; }
503 inline bool operator==(const uint160& a, const uint160& b)              { return (base_uint160)a == (base_uint160)b; }
504 inline bool operator!=(const uint160& a, const uint160& b)              { return (base_uint160)a != (base_uint160)b; }
505 inline const uint160 operator^(const uint160& a, const uint160& b)      { return (base_uint160)a ^  (base_uint160)b; }
506 inline const uint160 operator&(const uint160& a, const uint160& b)      { return (base_uint160)a &  (base_uint160)b; }
507 inline const uint160 operator|(const uint160& a, const uint160& b)      { return (base_uint160)a |  (base_uint160)b; }
508 inline const uint160 operator+(const uint160& a, const uint160& b)      { return (base_uint160)a +  (base_uint160)b; }
509 inline const uint160 operator-(const uint160& a, const uint160& b)      { return (base_uint160)a -  (base_uint160)b; }
510
511
512
513
514
515
516 //////////////////////////////////////////////////////////////////////////////
517 //
518 // uint256
519 //
520
521 /** 256-bit unsigned integer */
522 class uint256 : public base_uint256
523 {
524 public:
525     typedef base_uint256 basetype;
526
527     uint256()
528     {
529         for (int i = 0; i < WIDTH; i++)
530             pn[i] = 0;
531     }
532
533     uint256(const basetype& b)
534     {
535         for (int i = 0; i < WIDTH; i++)
536             pn[i] = b.pn[i];
537     }
538
539     uint256& operator=(const basetype& b)
540     {
541         for (int i = 0; i < WIDTH; i++)
542             pn[i] = b.pn[i];
543         return *this;
544     }
545
546     uint256(uint64_t b)
547     {
548         pn[0] = (unsigned int)b;
549         pn[1] = (unsigned int)(b >> 32);
550         for (int i = 2; i < WIDTH; i++)
551             pn[i] = 0;
552     }
553
554     uint256& operator=(uint64_t b)
555     {
556         pn[0] = (unsigned int)b;
557         pn[1] = (unsigned int)(b >> 32);
558         for (int i = 2; i < WIDTH; i++)
559             pn[i] = 0;
560         return *this;
561     }
562
563     explicit uint256(const std::string& str)
564     {
565         SetHex(str);
566     }
567
568     explicit uint256(const std::vector<unsigned char>& vch)
569     {
570         if (vch.size() == sizeof(pn))
571             memcpy(pn, &vch[0], sizeof(pn));
572         else
573             *this = 0;
574     }
575 };
576
577 inline bool operator==(const uint256& a, uint64_t b)                           { return (base_uint256)a == b; }
578 inline bool operator!=(const uint256& a, uint64_t b)                           { return (base_uint256)a != b; }
579 inline const uint256 operator<<(const base_uint256& a, unsigned int shift)   { return uint256(a) <<= shift; }
580 inline const uint256 operator>>(const base_uint256& a, unsigned int shift)   { return uint256(a) >>= shift; }
581 inline const uint256 operator<<(const uint256& a, unsigned int shift)        { return uint256(a) <<= shift; }
582 inline const uint256 operator>>(const uint256& a, unsigned int shift)        { return uint256(a) >>= shift; }
583
584 inline const uint256 operator^(const base_uint256& a, const base_uint256& b) { return uint256(a) ^= b; }
585 inline const uint256 operator&(const base_uint256& a, const base_uint256& b) { return uint256(a) &= b; }
586 inline const uint256 operator|(const base_uint256& a, const base_uint256& b) { return uint256(a) |= b; }
587 inline const uint256 operator+(const base_uint256& a, const base_uint256& b) { return uint256(a) += b; }
588 inline const uint256 operator-(const base_uint256& a, const base_uint256& b) { return uint256(a) -= b; }
589
590 inline bool operator<(const base_uint256& a, const uint256& b)          { return (base_uint256)a <  (base_uint256)b; }
591 inline bool operator<=(const base_uint256& a, const uint256& b)         { return (base_uint256)a <= (base_uint256)b; }
592 inline bool operator>(const base_uint256& a, const uint256& b)          { return (base_uint256)a >  (base_uint256)b; }
593 inline bool operator>=(const base_uint256& a, const uint256& b)         { return (base_uint256)a >= (base_uint256)b; }
594 inline bool operator==(const base_uint256& a, const uint256& b)         { return (base_uint256)a == (base_uint256)b; }
595 inline bool operator!=(const base_uint256& a, const uint256& b)         { return (base_uint256)a != (base_uint256)b; }
596 inline const uint256 operator^(const base_uint256& a, const uint256& b) { return (base_uint256)a ^  (base_uint256)b; }
597 inline const uint256 operator&(const base_uint256& a, const uint256& b) { return (base_uint256)a &  (base_uint256)b; }
598 inline const uint256 operator|(const base_uint256& a, const uint256& b) { return (base_uint256)a |  (base_uint256)b; }
599 inline const uint256 operator+(const base_uint256& a, const uint256& b) { return (base_uint256)a +  (base_uint256)b; }
600 inline const uint256 operator-(const base_uint256& a, const uint256& b) { return (base_uint256)a -  (base_uint256)b; }
601
602 inline bool operator<(const uint256& a, const base_uint256& b)          { return (base_uint256)a <  (base_uint256)b; }
603 inline bool operator<=(const uint256& a, const base_uint256& b)         { return (base_uint256)a <= (base_uint256)b; }
604 inline bool operator>(const uint256& a, const base_uint256& b)          { return (base_uint256)a >  (base_uint256)b; }
605 inline bool operator>=(const uint256& a, const base_uint256& b)         { return (base_uint256)a >= (base_uint256)b; }
606 inline bool operator==(const uint256& a, const base_uint256& b)         { return (base_uint256)a == (base_uint256)b; }
607 inline bool operator!=(const uint256& a, const base_uint256& b)         { return (base_uint256)a != (base_uint256)b; }
608 inline const uint256 operator^(const uint256& a, const base_uint256& b) { return (base_uint256)a ^  (base_uint256)b; }
609 inline const uint256 operator&(const uint256& a, const base_uint256& b) { return (base_uint256)a &  (base_uint256)b; }
610 inline const uint256 operator|(const uint256& a, const base_uint256& b) { return (base_uint256)a |  (base_uint256)b; }
611 inline const uint256 operator+(const uint256& a, const base_uint256& b) { return (base_uint256)a +  (base_uint256)b; }
612 inline const uint256 operator-(const uint256& a, const base_uint256& b) { return (base_uint256)a -  (base_uint256)b; }
613
614 inline bool operator<(const uint256& a, const uint256& b)               { return (base_uint256)a <  (base_uint256)b; }
615 inline bool operator<=(const uint256& a, const uint256& b)              { return (base_uint256)a <= (base_uint256)b; }
616 inline bool operator>(const uint256& a, const uint256& b)               { return (base_uint256)a >  (base_uint256)b; }
617 inline bool operator>=(const uint256& a, const uint256& b)              { return (base_uint256)a >= (base_uint256)b; }
618 inline bool operator==(const uint256& a, const uint256& b)              { return (base_uint256)a == (base_uint256)b; }
619 inline bool operator!=(const uint256& a, const uint256& b)              { return (base_uint256)a != (base_uint256)b; }
620 inline const uint256 operator^(const uint256& a, const uint256& b)      { return (base_uint256)a ^  (base_uint256)b; }
621 inline const uint256 operator&(const uint256& a, const uint256& b)      { return (base_uint256)a &  (base_uint256)b; }
622 inline const uint256 operator|(const uint256& a, const uint256& b)      { return (base_uint256)a |  (base_uint256)b; }
623 inline const uint256 operator+(const uint256& a, const uint256& b)      { return (base_uint256)a +  (base_uint256)b; }
624 inline const uint256 operator-(const uint256& a, const uint256& b)      { return (base_uint256)a -  (base_uint256)b; }
625
626 #endif