Mercurial > octave-nkf
comparison scripts/statistics/distributions/binoinv.m @ 20644:4e307c55a2b5
Use isempty () rather than any () for faster code in inverse statistical distributions.
betainv.m, binoinv.m, gaminv.m, poissinv.m: Use '! isempty (k)' rather than
'any (k)' for faster code.
| author | Rik <rik@octave.org> |
|---|---|
| date | Sun, 11 Oct 2015 21:09:41 -0700 |
| parents | 9d2023d1a63c |
| children |
comparison
equal
deleted
inserted
replaced
| 20643:d6d04088ac9e | 20644:4e307c55a2b5 |
|---|---|
| 53 !(p >= 0) | !(p <= 1)); | 53 !(p >= 0) | !(p <= 1)); |
| 54 inv(k) = NaN; | 54 inv(k) = NaN; |
| 55 | 55 |
| 56 k = find ((x >= 0) & (x <= 1) & (n >= 0) & (n == fix (n) | 56 k = find ((x >= 0) & (x <= 1) & (n >= 0) & (n == fix (n) |
| 57 & (p >= 0) & (p <= 1))); | 57 & (p >= 0) & (p <= 1))); |
| 58 if (any (k)) | 58 if (! isempty (k)) |
| 59 x = x(k); | 59 x = x(k); |
| 60 if (isscalar (n) && isscalar (p)) | 60 if (isscalar (n) && isscalar (p)) |
| 61 [inv(k), unfinished] = scalar_binoinv (x(:), n, p); | 61 [inv(k), unfinished] = scalar_binoinv (x(:), n, p); |
| 62 k = k(unfinished); | 62 k = k(unfinished); |
| 63 if (! isempty (k)) | 63 if (! isempty (k)) |