Update optimize.rst
ENH: allow sparse __rmul__ to work with custom types
BUG: special: fix asymptotic expansion in rgamma
BUG: special: fix several issues in poch

Handle better the cases where the gamma functions in the definition hit
poles or over/underflow.

Fixes gh-1437
GEN: special: regenerate files
DOC: enhanced fftpack tutorial
DOC: further changes in fftpack tutorial
enhanced one example
Fixed issue #2720 in tutorial
adressed review comments, used [] throughout the document
DOC: Add Notes section to chebwin docstring
ENH integrate: quadrature accepts non-tuple args argument

Copy quad behavior accepting a single element (integer) instead of the
expected tuple

Closes gh-2044
BUG: fix the 'lower' attribute of dpotri to make dpotrf/dpotri consistent.
Should close gh-2691

Also remove (deprecated) clapack dpotri/dpotrf functions.
changing the doc strings of measurements.sum
correcting the measurements.sum() doc strings
BUG: special: fix overflow at low order in eval_hermite

Hermite polynomials are very large and using the relationship to
the laguerre polynomials overflow earlier than necessary (around
N = 67 or so.)  Using the Clenshaw recursion we can evaluate to
much higher orders.

From now we evaluate the hermite polynomial as a stretched and
scaled normalized hermite (He) polynomial since He is slightly
smaller at any given order.
BUG: special: improve h_roots accuracy at high order

The roots of the hermite polynomials are still computed by using the
Golub-Welsch algorithm.  However a single pass of newton's method is
applied to improve the roots.  The symmetry of the Hermite roots
about zero is used as well.  The weights are calculated and then
rescaled to ensure that the integral of 1 is correct.

These changes were largely copied from the hermgauss implementation
in numpy (by Charles Harris).

Fixes #1991
ENH: expose control truncation of ndimage.gaussian_filter().

This change gives the user control over the size at which the Gaussian filter is truncated,
with default behavior to truncate at 4 standard deviations (the previous fixed value).
ENH: Pass extra keyword arguments to gaussian_laplace() and gaussian_gradient_magnitude() to gaussian_filter()

This change allows control over the truncation radius of the gaussian filter used in
gaussian_laplace() and gaussian_gradient_magnitude().
TST: ndimage: add more tests in test_gaussian_truncate()