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.
BUG: FITPACK conditions were not properly checked in LSQUnivariateSpline.

fixes #2611
STY: PEP8 for LSQUnivariateSpline
TST: sparse: replace __numpy_ufunc__ version checks by functionality check
DOC: Updated integrate.ode doc (lband, uband)
TST: integrate: Added CoupledDecay to test_integrate.py
MAINT: integrate: Fix documentation of the ode 'with_jacobian' argument.

* Fixed the documentation of the 'with_jacobian' argument of the vode
  and lsoda solvers.
* Changed the default value of 'with_jacobian' from 0 to False.