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
MAINT: avoid redundant copying, now that the trmm wrapper has been fixed
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
TST: compare values of pdf with precomputed values
MAINT: create new multivariate_normal_gen instance in frozen class