Function rk_ulong

numpy/random/src/mt19937/randomkit.c:439–445

long rk_long(rk_state *state) { return rk_ulong(state) >> 1; }

unsigned long rk_ulong(rk_state *state) {
#if ULONG_MAX <= 0xffffffffUL
  return rk_random(state);
#else
  return (rk_random(state) << 32) | (rk_random(state));
#endif
}

unsigned long rk_interval(unsigned long max, rk_state *state) {
  unsigned long mask = max, value;