dig_u             290 random.c           uint64_t x = ((uint64_t)a << dig_u) +
dig_u             291 random.c       	(((uint64_t)b + (a >> dig_u)) >> dig_r64);
dig_u             295 random.c           b = (b >> dig_r64) + (((a >> dig_u) + (b & bmask)) >> dig_r64);
dig_u             296 random.c           r = (double)a * (1 << dig_u) + b;