contrib/putty/INT64.C

   1 /*\r
   2  * Handling of the int64 and uint64 types. Done in 32-bit integers,\r
   3  * for (pre-C99) portability. Hopefully once C99 becomes widespread\r
   4  * we can kiss this lot goodbye...\r
   5  */\r
   6 \r
   7 #include <assert.h>\r
   8 #include <string.h>\r
   9 \r
  10 #include "int64.h"\r
  11 \r
  12 uint64 uint64_div10(uint64 x, int *remainder)\r
  13 {\r
  14     uint64 y;\r
  15     unsigned int rem, r2;\r
  16     y.hi = x.hi / 10;\r
  17     y.lo = x.lo / 10;\r
  18     rem = x.lo % 10;\r
  19     /*\r
  20      * Now we have to add in the remainder left over from x.hi.\r
  21      */\r
  22     r2 = x.hi % 10;\r
  23     y.lo += r2 * 429496729;\r
  24     rem += r2 * 6;\r
  25     y.lo += rem / 10;\r
  26     rem %= 10;\r
  27 \r
  28     if (remainder)\r
  29         *remainder = rem;\r
  30     return y;\r
  31 }\r
  32 \r
  33 void uint64_decimal(uint64 x, char *buffer)\r
  34 {\r
  35     char buf[20];\r
  36     int start = 20;\r
  37     int d;\r
  38 \r
  39     do {\r
  40         x = uint64_div10(x, &d);\r
  41         assert(start > 0);\r
  42         buf[--start] = d + '0';\r
  43     } while (x.hi || x.lo);\r
  44 \r
  45     memcpy(buffer, buf + start, sizeof(buf) - start);\r
  46     buffer[sizeof(buf) - start] = '\0';\r
  47 }\r
  48 \r
  49 uint64 uint64_make(unsigned long hi, unsigned long lo)\r
  50 {\r
  51     uint64 y;\r
  52     y.hi = hi & 0xFFFFFFFFU;\r
  53     y.lo = lo & 0xFFFFFFFFU;\r
  54     return y;\r
  55 }\r
  56 \r
  57 uint64 uint64_add(uint64 x, uint64 y)\r
  58 {\r
  59     x.lo = (x.lo + y.lo) & 0xFFFFFFFFU;\r
  60     x.hi += y.hi + (x.lo < y.lo ? 1 : 0);\r
  61     return x;\r
  62 }\r
  63 \r
  64 uint64 uint64_add32(uint64 x, unsigned long y)\r
  65 {\r
  66     uint64 yy;\r
  67     yy.hi = 0;\r
  68     yy.lo = y;\r
  69     return uint64_add(x, yy);\r
  70 }\r
  71 \r
  72 int uint64_compare(uint64 x, uint64 y)\r
  73 {\r
  74     if (x.hi != y.hi)\r
  75         return x.hi < y.hi ? -1 : +1;\r
  76     if (x.lo != y.lo)\r
  77         return x.lo < y.lo ? -1 : +1;\r
  78     return 0;\r
  79 }\r
  80 \r
  81 uint64 uint64_subtract(uint64 x, uint64 y)\r
  82 {\r
  83     x.lo = (x.lo - y.lo) & 0xFFFFFFFFU;\r
  84     x.hi = (x.hi - y.hi - (x.lo > (y.lo ^ 0xFFFFFFFFU) ? 1 : 0)) & 0xFFFFFFFFU;\r
  85     return x;\r
  86 }\r
  87 \r
  88 double uint64_to_double(uint64 x)\r
  89 {\r
  90     return (4294967296.0 * x.hi) + (double)x.lo;\r
  91 }\r
  92 \r
  93 uint64 uint64_shift_right(uint64 x, int shift)\r
  94 {\r
  95     if (shift < 32) {\r
  96         x.lo >>= shift;\r
  97         x.lo |= (x.hi << (32-shift)) & 0xFFFFFFFFU;\r
  98         x.hi >>= shift;\r
  99     } else {\r
 100         x.lo = x.hi >> (shift-32);\r
 101         x.hi = 0;\r
 102     }\r
 103     return x;\r
 104 }\r
 105 \r
 106 uint64 uint64_shift_left(uint64 x, int shift)\r
 107 {\r
 108     if (shift < 32) {\r
 109         x.hi = (x.hi << shift) & 0xFFFFFFFFU;\r
 110         x.hi |= (x.lo >> (32-shift));\r
 111         x.lo = (x.lo << shift) & 0xFFFFFFFFU;\r
 112     } else {\r
 113         x.hi = (x.lo << (shift-32)) & 0xFFFFFFFFU;\r
 114         x.lo = 0;\r
 115     }\r
 116     return x;\r
 117 }\r
 118 \r
 119 uint64 uint64_from_decimal(char *str)\r
 120 {\r
 121     uint64 ret;\r
 122     ret.hi = ret.lo = 0;\r
 123     while (*str >= '0' && *str <= '9') {\r
 124         ret = uint64_add(uint64_shift_left(ret, 3),\r
 125                          uint64_shift_left(ret, 1));\r
 126         ret = uint64_add32(ret, *str - '0');\r
 127         str++;\r
 128     }\r
 129     return ret;\r
 130 }\r
 131 \r
 132 #ifdef TESTMODE\r
 133 \r
 134 #include <stdio.h>\r
 135 \r
 136 int main(void)\r
 137 {\r
 138     uint64 x, y, z;\r
 139     char buf[80];\r
 140 \r
 141     x = uint64_make(0x3456789AUL, 0xDEF01234UL);\r
 142     printf("%08lx.%08lx\n", x.hi, x.lo);\r
 143     uint64_decimal(x, buf);\r
 144     printf("%s\n", buf);\r
 145 \r
 146     y = uint64_add32(x, 0xFFFFFFFFU);\r
 147     printf("%08lx.%08lx\n", y.hi, y.lo);\r
 148     uint64_decimal(y, buf);\r
 149     printf("%s\n", buf);\r
 150 \r
 151     z = uint64_subtract(y, x);\r
 152     printf("%08lx.%08lx\n", z.hi, z.lo);\r
 153     uint64_decimal(z, buf);\r
 154     printf("%s\n", buf);\r
 155 \r
 156     z = uint64_subtract(x, y);\r
 157     printf("%08lx.%08lx\n", z.hi, z.lo);\r
 158     uint64_decimal(z, buf);\r
 159     printf("%s\n", buf);\r
 160 \r
 161     y = uint64_shift_right(x, 4);\r
 162     printf("%08lx.%08lx\n", y.hi, y.lo);\r
 163 \r
 164     y = uint64_shift_right(x, 36);\r
 165     printf("%08lx.%08lx\n", y.hi, y.lo);\r
 166 \r
 167     y = uint64_shift_left(x, 4);\r
 168     printf("%08lx.%08lx\n", x.hi, x.lo);\r
 169 \r
 170     y = uint64_shift_left(x, 36);\r
 171     printf("%08lx.%08lx\n", x.hi, x.lo);\r
 172 \r
 173     return 0;\r
 174 }\r
 175 #endif\r