if (k - x <= s || u >= hIntegral(k + 0.5) - h(k)) {
if (k - x <= s || u >= hIntegral(k + 0.5) - h(k)) {
/** Simple constructor. * @param numberOfElements number of elements * @param exponent exponent parameter of the distribution */ ZipfRejectionInversionSampler(final int numberOfElements, final double exponent) { this.exponent = exponent; this.numberOfElements = numberOfElements; this.hIntegralX1 = hIntegral(1.5) - 1d; this.hIntegralNumberOfElements = hIntegral(numberOfElements + 0.5); this.s = 2d - hIntegralInverse(hIntegral(2.5) - h(2)); }
/** Simple constructor. * @param numberOfElements number of elements * @param exponent exponent parameter of the distribution */ ZipfRejectionInversionSampler(final int numberOfElements, final double exponent) { this.exponent = exponent; this.numberOfElements = numberOfElements; this.hIntegralX1 = hIntegral(1.5) - 1d; this.hIntegralNumberOfElements = hIntegral(numberOfElements + 0.5); this.s = 2d - hIntegralInverse(hIntegral(2.5) - h(2)); }
if (k - x <= s || u >= hIntegral(k + 0.5) - h(k)) {
/** Simple constructor. * @param numberOfElements number of elements * @param exponent exponent parameter of the distribution */ ZipfRejectionInversionSampler(final int numberOfElements, final double exponent) { this.exponent = exponent; this.numberOfElements = numberOfElements; this.hIntegralX1 = hIntegral(1.5) - 1d; this.hIntegralNumberOfElements = hIntegral(numberOfElements + 0.5); this.s = 2d - hIntegralInverse(hIntegral(2.5) - h(2)); }