docs: def tau_rand_int(state):
      """A fast (pseudo)-random number generator.
  
      Parameters
      ----------
      state: array of int64, shape (3,)
          The internal state of the rng
  
      Returns
      -------
      A (pseudo)-random int32 value
      """
      state[0] = (((state[0] & 4294967294) << 12) & 0xffffffff) ^ (
          (((state[0] << 13) & 0xffffffff) ^ state[0]) >> 19
      )
      state[1] = (((state[1] & 4294967288) << 4) & 0xffffffff) ^ (
          (((state[1] << 2) & 0xffffffff) ^ state[1]) >> 25
      )
      state[2] = (((state[2] & 4294967280) << 17) & 0xffffffff) ^ (
          (((state[2] << 3) & 0xffffffff) ^ state[2]) >> 11
      )
  
      return state[0] ^ state[1] ^ state[2]