docs: private static final int getR(final double theta, final double alpha, final double p) {
      final double split1 = (p * (alpha + 1.0)) / 2.0;
      if (theta > split1) { return 0; }
      if (theta > (alpha * split1)) { return 1; }
      return 2;
    }