Add note to scipy.optimize.brute_force about the datatype of the evalution.
ENH: scipy.special information theory functions
TST: remove negative inputs
TST: adding test cases for COO sum_duplicates
ENH: adding sum_duplicates method to COO matrix
adding has_canonical_format to COO matrix
fixing build errors
incorporating @jaimefrio's vectorized version
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
DEP: put back veccdf/vecfunc/vec_generic_momemt in stats.distributions.

Removed erroneously in gh-3243 for 0.14.0.  These methods were always
private ones but not marked as such, and as the comments on gh-3243 show
they were used. Hence now deprecated.
MAINT: attempt to improve svds of singular matrices, passes existing tests
MAINT: sparse svd passes new tests
BUG: fix romb in scipy.integrate.quadrature

Routine romb for Romberg integration using 2**k +1 samples of a function
was only using 2**(k-1) +1 samples due to wrong ranges.
Renumbered the values R to the common textbook numbering scheme starting
with 0s. Note that romb with 2**k +1 should give the same values as
romberg with divmax=k. Or with 3 values romb is the same as simps.
ENH: cluster: reimplement the update-step of K-means in Cython
MAINT: cluster: remove the conditional import for _vq
BUG: stub out spurious methods of Akima1DIntgerpolator

from_spline and from_bernstein_basis are inherited
from PPoly, and neither produces an Akima interpolator.
DOC: add See Also sections to Pchip and Akima interpolators
ENH: Added discrete Lyapunov transformation solve
Reverse Cuthill-McKee reordering of sparse graphs

Added the symrcm function that takes a symmetric graph in CSR or CSC
sparse format and returns an array of row and column permutations that
reduces the overall bandwidth of the underlying matrix.