MAINT: spatial.distance: generate boilerplate with macros
BUG: spatial.distance: fix unchecked malloc
MAINT: spatial.distance: further refactoring

* static on all functions in distance_impl.h
* commented out unused compute_mean_vector
* replaced &&, || with &, | in tight loops
* reindented code according to NumPy conventions
* removed some useless variables
ENH: convert the cryptically named scipy.special information theory functions into ufuncs
TST: add negative values to scipy.special information theory function tests
DOC: improve scipy.special information theory ufunc docstrings
MAINT: add .pxd files
MAINT: simplify kl_div
MAINT: combine multiple ufuncs into the same .pxd file
DOC: remove docstring from UnivariateSpline.__init__

same for InterpolatedUnivariateSpline and LSQUnivariateSpline.
These simply duplicated the class docstrings.
DOC: improve the docstrings of *UnivariateSpline classes

Tweak examples to make plots actually render.
DOC: add a nan masking example to the UnivariateSpline docstring

example is by @juliantaylor,
https://github.com/scipy/scipy/issues/3531#issuecomment-40233093
ENH: add huber loss function to the collection of convexity-related ufuncs
PEP8: pacify the pep8 overseer
ENH: add pseudo-huber loss function ufunc
BUG: stats: correct value of powerlaw.logpdf(0, 1) is 0. Closes gh-3986.
BUG: stats: Better handling of x=0 in the 'exponpow' distribution.

exponpow.logpdf(0, 1) now returns 0.0 instead of nan.

When b < 1, exponpow.pdf(0, b) and exponpow.logpdf(0, b) now
return inf instead of nan.

Closes gh-3982.
BUG: stats: gengamma.pdf(0,1,1) now returns 1.0 instead of nan.

Closes gh-3985.
ENH: using python's warnings instead of prints in fitpack.

When fitpack detects a recoverable problem, it will signal it with
a RuntimeWarning instead of printing to stdout. This way the output
can be programatically silenced or made exceptions.
MAINT: optimize: remove a couple unused variables in zeros.c