This is an automated email from the ASF dual-hosted git repository. aherbert pushed a commit to branch master in repository https://gitbox.apache.org/repos/asf/commons-rng.git
commit 4d211d27144e046903f48ed260ee1ca1144c7a85 Author: aherbert <aherb...@apache.org> AuthorDate: Fri Jul 9 10:44:42 2021 +0100 Added missing branch frequencies. --- .../rng/sampling/distribution/ZigguratNormalizedGaussianSampler.java | 2 ++ 1 file changed, 2 insertions(+) diff --git a/commons-rng-sampling/src/main/java/org/apache/commons/rng/sampling/distribution/ZigguratNormalizedGaussianSampler.java b/commons-rng-sampling/src/main/java/org/apache/commons/rng/sampling/distribution/ZigguratNormalizedGaussianSampler.java index 1d1acb2..db3c8c6 100644 --- a/commons-rng-sampling/src/main/java/org/apache/commons/rng/sampling/distribution/ZigguratNormalizedGaussianSampler.java +++ b/commons-rng-sampling/src/main/java/org/apache/commons/rng/sampling/distribution/ZigguratNormalizedGaussianSampler.java @@ -107,6 +107,7 @@ public class ZigguratNormalizedGaussianSampler final long j = rng.nextLong(); final int i = ((int) j) & LAST; if (Math.abs(j) < K[i]) { + // This branch is called about 0.972101 times per sample. return j * W[i]; } return fix(j, i); @@ -152,6 +153,7 @@ public class ZigguratNormalizedGaussianSampler // This branch is called about 0.027323 times per sample. final double x = hz * W[iz]; if (F[iz] + rng.nextDouble() * (F[iz - 1] - F[iz]) < pdf(x)) { + // This branch is called about 0.014961 times per sample. return x; } // Try again.