docs: private double pseudoRandom(final long seed, int input) {
      // Default constants from "man drand48"
      final long mult = 0x5DEECE66DL;
      final long add = 0xBL;
      final long mask = (1L << 48) - 1; // 48 bit
      // Produce an initial seed each
      final long i1 = (input ^ seed ^ mult) & mask;
      final long i2 = (input ^ (seed >>> 16) ^ mult) & mask;
      // Compute the first random each
      final long l1 = (i1 * mult + add) & mask;
      final long l2 = (i2 * mult + add) & mask;
      // Use 53 bit total:
      final int r1 = (int) (l1 >>> 22); // 48 - 22 = 26
      final int r2 = (int) (l2 >>> 21); // 48 - 21 = 27
      double random = ((((long) r1) << 27) + r2) / (double) (1L << 53);
      return random;
    }