Add note to scipy.optimize.brute_force about the datatype of the evalution.
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
BUG: FITPACK conditions were not properly checked in LSQUnivariateSpline.

fixes #2611
STY: PEP8 for LSQUnivariateSpline
ENH: release GIL for filter and interpolation of ndimage

allows using python threads for parallization of indepent operations.
Revert "the mu parameter can also be a float"

This reverts commit 48674c2835861a2ce6a663f94fbc95f5cccdd663.