+
+uint64 uint64_subtract(uint64 x, uint64 y)
+{
+ x.lo = (x.lo - y.lo) & 0xFFFFFFFFU;
+ x.hi = (x.hi - y.hi - (x.lo > (y.lo ^ 0xFFFFFFFFU) ? 1 : 0)) & 0xFFFFFFFFU;
+ return x;
+}
+
+double uint64_to_double(uint64 x)
+{
+ return (4294967296.0 * x.hi) + (double)x.lo;
+}
+
+uint64 uint64_shift_right(uint64 x, int shift)
+{
+ if (shift < 32) {
+ x.lo >>= shift;
+ x.lo |= (x.hi << (32-shift)) & 0xFFFFFFFFU;
+ x.hi >>= shift;
+ } else {
+ x.lo = x.hi >> (shift-32);
+ x.hi = 0;
+ }
+ return x;
+}
+
+uint64 uint64_shift_left(uint64 x, int shift)
+{
+ if (shift < 32) {
+ x.hi = (x.hi << shift) & 0xFFFFFFFFU;
+ x.hi |= (x.lo >> (32-shift));
+ x.lo = (x.lo << shift) & 0xFFFFFFFFU;
+ } else {
+ x.hi = (x.lo << (shift-32)) & 0xFFFFFFFFU;
+ x.lo = 0;
+ }
+ return x;
+}
+
+uint64 uint64_from_decimal(char *str)
+{
+ uint64 ret;
+ ret.hi = ret.lo = 0;
+ while (*str >= '0' && *str <= '9') {
+ ret = uint64_add(uint64_shift_left(ret, 3),
+ uint64_shift_left(ret, 1));
+ ret = uint64_add32(ret, *str - '0');
+ str++;
+ }
+ return ret;
+}
+
+#ifdef TESTMODE
+
+#include <stdio.h>
+
+int main(void)
+{
+ uint64 x, y, z;
+ char buf[80];
+
+ x = uint64_make(0x3456789AUL, 0xDEF01234UL);
+ printf("%08lx.%08lx\n", x.hi, x.lo);
+ uint64_decimal(x, buf);
+ printf("%s\n", buf);
+
+ y = uint64_add32(x, 0xFFFFFFFFU);
+ printf("%08lx.%08lx\n", y.hi, y.lo);
+ uint64_decimal(y, buf);
+ printf("%s\n", buf);
+
+ z = uint64_subtract(y, x);
+ printf("%08lx.%08lx\n", z.hi, z.lo);
+ uint64_decimal(z, buf);
+ printf("%s\n", buf);
+
+ z = uint64_subtract(x, y);
+ printf("%08lx.%08lx\n", z.hi, z.lo);
+ uint64_decimal(z, buf);
+ printf("%s\n", buf);
+
+ y = uint64_shift_right(x, 4);
+ printf("%08lx.%08lx\n", y.hi, y.lo);
+
+ y = uint64_shift_right(x, 36);
+ printf("%08lx.%08lx\n", y.hi, y.lo);
+
+ y = uint64_shift_left(x, 4);
+ printf("%08lx.%08lx\n", x.hi, x.lo);
+
+ y = uint64_shift_left(x, 36);
+ printf("%08lx.%08lx\n", x.hi, x.lo);
+
+ return 0;
+}
+#endif