Add note to scipy.optimize.brute_force about the datatype of the evalution.
ENH: Small speedup of FFT size check

Use a bit-twiddling power of 2 check instead of repeatedly dividing by 2
ENH: Use floor division instead of float division

This uses from __future__ import division, but if n % c == 0, n / c == n // c
ENH: add the limit shape c=0 to genpareto distribution

For c=0, genpareto is equivalent to the exponential
distribution.
BUG: use special.log1p instead of np.log1p

np.log1p(np.inf) is reported to produce nans on some platforms,
see numpy issue gh-4225
BUG: use scipy.special.expm1 instead of np.expm1
ENH: vectorize genpareto._munp

implementation is by @pbrod
BUG: special-case genpareto(c=-1)
ENH/BUG: support negative m in elliptic integrals

Prior, ellipkinc, and ellipeinc silently returned garbage for
negative m.  Both functions are still inaccurate due at extremely
large phi.
BUG: ellipkinc(p, 1) imprecise.

The old expression was imprecise for small phi and for phi very
close to pi/2.
TST: special: enable some mpmath test for incomplete elliptic integrals E and K
TST: enable ellipkinc tests against boost data
TST: Add tables driven tests of ellip(k,e)inc values with neg m
ENH: fix failing argument reduction for large phi. (ellip(e,k)inc)

The floating point multiples of pi/2 won't fit into an int.
Argument reduce into a double instead.
ENH: Improve precision of ellipeinc for small p, m in [0,1].
DOC: Improve docstrings for elliptic integrals.
scipy.special.orthogonal.p_roots with cache
the mu parameter can also be a float
Coerce args into tuple if necessary.

Same logic as in scipy.integrate.quad.

Fixes #3503 and Fixes #3785.
Added regression test.