MAINT: speed up expm a bit
MAINT: expm now works for numpy matrices as input, including a test
MAINT: move the code that computes the one-norm of a power of a non-negative matrix into its own function, and add a test
MAINT: fblas -> blas
MAINT: blas -> get_blas_funcs
BUG: optimize: mark constraint violations in cobyla as failed solution
MAINT: avoid redundant copying, now that the trmm wrapper has been fixed
DOC: optimize: approx_fprime missing from documentation
ENH: stats/_multivariate: implementation of multivariate normal variable
MAINT: mvn distn -- do not choke on large determinants, and make pinv and pdet agree on which singular values are small
TST: remove a call to a numpy function that is too new to use
MAINT: use a more standard python version compatibility preamble
MAINT: simplified parameter specification of multivariate_normal_frozen class

Suggestion by bburlo.
MAINT: changed wording in docstring to reflect order of return values
MAINT: replace einsum, using a decomposed precision matrix
MAINT: use eigh instead of svd
MAINT: trivial changes to docstring formatting
MAINT: implemented more careful handling of cutoff parameter
DOC: more explicit description for parameters in _logpdf
MAINT: explicit parameters in construction frozen class