BUG: sparse: fix operator overloading with custom types

See issue #3580.
ENH: sparse: inplace ops should fall back to x = x op y if not implemented
BUG: sparse: itruediv should be allowed on sparse matrix
Adds default params to pdfs of multivariate_norm

Adds default params to pdf and logpdf to match
`norm` and the documentation.
BUG: signal.welch: ensure floating point dtype regardless of input dtype

Fixes gh-3655.
Merge pull request #3667 from tshauck/multivariate-norm-docs

ENH: Adds default params to pdfs of multivariate_norm
Merge pull request #3657 from ewmoore/specint

BUG: signal.welch: ensure floating point dtype regardless of input dtype
TST: speed up TravisCI tests of scipy.special, by installing gmpy2.

Gain is ~230 seconds on my machine (will be more on TravisCI).
DOC: remove `special.jn` from refguide and mention it's an alias of jv.

Closes gh-2564.
Merge pull request #3681 from rgommers/speedup-tests

MAINT: install gmpy2 in travis-ci to speed up tests
Merge pull request #3589 from larsmans/sparse-binop

BUG: sparse: allow operator overloading

Make binop methods should return NotImplemented rather than raising errors.
DOC: add some docstrings for special functions in the nctdtr family.

These were not yet in the built docs, see gh-2746.
DOC: add some docstrings for special functions in the ncfdtr family.

These were not yet in the built docs, see gh-2746.
DOC: add docstrings for special.nrdtrisd/nrdtrimn.
GEN: regenerate special/_ufuncs.pyx
DEP: deprecate special.bessel_diff_formula, it's a private function.
DOC: add docstring to special.assoc_laguerre and add that to built docs.

Closes gh-2746.
Merge pull request #3684 from rgommers/special-docs

Update special docs
DOC: spacing in special docstrings
GEN: regenerate special/_ufuncs.pyx