| 790 | } |
| 791 | |
| 792 | int64_t random_binomial(bitgen_t *bitgen_state, double p, int64_t n, |
| 793 | binomial_t *binomial) { |
| 794 | double q; |
| 795 | |
| 796 | if ((n == 0LL) || (p == 0.0f)) |
| 797 | return 0; |
| 798 | |
| 799 | if (p <= 0.5) { |
| 800 | if (p * n <= 30.0) { |
| 801 | return random_binomial_inversion(bitgen_state, n, p, binomial); |
| 802 | } else { |
| 803 | return random_binomial_btpe(bitgen_state, n, p, binomial); |
| 804 | } |
| 805 | } else { |
| 806 | q = 1.0 - p; |
| 807 | if (q * n <= 30.0) { |
| 808 | return n - random_binomial_inversion(bitgen_state, n, q, binomial); |
| 809 | } else { |
| 810 | return n - random_binomial_btpe(bitgen_state, n, q, binomial); |
| 811 | } |
| 812 | } |
| 813 | } |
| 814 | |
| 815 | double random_noncentral_chisquare(bitgen_t *bitgen_state, double df, |
| 816 | double nonc) { |
no test coverage detected