--- /dev/null
+/*\r
+ * Handling of the int64 and uint64 types. Done in 32-bit integers,\r
+ * for (pre-C99) portability. Hopefully once C99 becomes widespread\r
+ * we can kiss this lot goodbye...\r
+ */\r
+\r
+#include <assert.h>\r
+#include <string.h>\r
+\r
+#include "int64.h"\r
+\r
+uint64 uint64_div10(uint64 x, int *remainder)\r
+{\r
+ uint64 y;\r
+ unsigned int rem, r2;\r
+ y.hi = x.hi / 10;\r
+ y.lo = x.lo / 10;\r
+ rem = x.lo % 10;\r
+ /*\r
+ * Now we have to add in the remainder left over from x.hi.\r
+ */\r
+ r2 = x.hi % 10;\r
+ y.lo += r2 * 429496729;\r
+ rem += r2 * 6;\r
+ y.lo += rem / 10;\r
+ rem %= 10;\r
+\r
+ if (remainder)\r
+ *remainder = rem;\r
+ return y;\r
+}\r
+\r
+void uint64_decimal(uint64 x, char *buffer)\r
+{\r
+ char buf[20];\r
+ int start = 20;\r
+ int d;\r
+\r
+ do {\r
+ x = uint64_div10(x, &d);\r
+ assert(start > 0);\r
+ buf[--start] = d + '0';\r
+ } while (x.hi || x.lo);\r
+\r
+ memcpy(buffer, buf + start, sizeof(buf) - start);\r
+ buffer[sizeof(buf) - start] = '\0';\r
+}\r
+\r
+uint64 uint64_make(unsigned long hi, unsigned long lo)\r
+{\r
+ uint64 y;\r
+ y.hi = hi & 0xFFFFFFFFU;\r
+ y.lo = lo & 0xFFFFFFFFU;\r
+ return y;\r
+}\r
+\r
+uint64 uint64_add(uint64 x, uint64 y)\r
+{\r
+ x.lo = (x.lo + y.lo) & 0xFFFFFFFFU;\r
+ x.hi += y.hi + (x.lo < y.lo ? 1 : 0);\r
+ return x;\r
+}\r
+\r
+uint64 uint64_add32(uint64 x, unsigned long y)\r
+{\r
+ uint64 yy;\r
+ yy.hi = 0;\r
+ yy.lo = y;\r
+ return uint64_add(x, yy);\r
+}\r
+\r
+int uint64_compare(uint64 x, uint64 y)\r
+{\r
+ if (x.hi != y.hi)\r
+ return x.hi < y.hi ? -1 : +1;\r
+ if (x.lo != y.lo)\r
+ return x.lo < y.lo ? -1 : +1;\r
+ return 0;\r
+}\r
+\r
+uint64 uint64_subtract(uint64 x, uint64 y)\r
+{\r
+ x.lo = (x.lo - y.lo) & 0xFFFFFFFFU;\r
+ x.hi = (x.hi - y.hi - (x.lo > (y.lo ^ 0xFFFFFFFFU) ? 1 : 0)) & 0xFFFFFFFFU;\r
+ return x;\r
+}\r
+\r
+double uint64_to_double(uint64 x)\r
+{\r
+ return (4294967296.0 * x.hi) + (double)x.lo;\r
+}\r
+\r
+uint64 uint64_shift_right(uint64 x, int shift)\r
+{\r
+ if (shift < 32) {\r
+ x.lo >>= shift;\r
+ x.lo |= (x.hi << (32-shift)) & 0xFFFFFFFFU;\r
+ x.hi >>= shift;\r
+ } else {\r
+ x.lo = x.hi >> (shift-32);\r
+ x.hi = 0;\r
+ }\r
+ return x;\r
+}\r
+\r
+uint64 uint64_shift_left(uint64 x, int shift)\r
+{\r
+ if (shift < 32) {\r
+ x.hi = (x.hi << shift) & 0xFFFFFFFFU;\r
+ x.hi |= (x.lo >> (32-shift));\r
+ x.lo = (x.lo << shift) & 0xFFFFFFFFU;\r
+ } else {\r
+ x.hi = (x.lo << (shift-32)) & 0xFFFFFFFFU;\r
+ x.lo = 0;\r
+ }\r
+ return x;\r
+}\r
+\r
+uint64 uint64_from_decimal(char *str)\r
+{\r
+ uint64 ret;\r
+ ret.hi = ret.lo = 0;\r
+ while (*str >= '0' && *str <= '9') {\r
+ ret = uint64_add(uint64_shift_left(ret, 3),\r
+ uint64_shift_left(ret, 1));\r
+ ret = uint64_add32(ret, *str - '0');\r
+ str++;\r
+ }\r
+ return ret;\r
+}\r
+\r
+#ifdef TESTMODE\r
+\r
+#include <stdio.h>\r
+\r
+int main(void)\r
+{\r
+ uint64 x, y, z;\r
+ char buf[80];\r
+\r
+ x = uint64_make(0x3456789AUL, 0xDEF01234UL);\r
+ printf("%08lx.%08lx\n", x.hi, x.lo);\r
+ uint64_decimal(x, buf);\r
+ printf("%s\n", buf);\r
+\r
+ y = uint64_add32(x, 0xFFFFFFFFU);\r
+ printf("%08lx.%08lx\n", y.hi, y.lo);\r
+ uint64_decimal(y, buf);\r
+ printf("%s\n", buf);\r
+\r
+ z = uint64_subtract(y, x);\r
+ printf("%08lx.%08lx\n", z.hi, z.lo);\r
+ uint64_decimal(z, buf);\r
+ printf("%s\n", buf);\r
+\r
+ z = uint64_subtract(x, y);\r
+ printf("%08lx.%08lx\n", z.hi, z.lo);\r
+ uint64_decimal(z, buf);\r
+ printf("%s\n", buf);\r
+\r
+ y = uint64_shift_right(x, 4);\r
+ printf("%08lx.%08lx\n", y.hi, y.lo);\r
+\r
+ y = uint64_shift_right(x, 36);\r
+ printf("%08lx.%08lx\n", y.hi, y.lo);\r
+\r
+ y = uint64_shift_left(x, 4);\r
+ printf("%08lx.%08lx\n", x.hi, x.lo);\r
+\r
+ y = uint64_shift_left(x, 36);\r
+ printf("%08lx.%08lx\n", x.hi, x.lo);\r
+\r
+ return 0;\r
+}\r
+#endif\r