Release 0.5.5 (#436)

* Update changelog and version number

* Improve speed of multivariate_normality (#437)

docs/changelog.rst

@@ -3,10 +3,18 @@
 What's new
 ##########
 
+v0.5.5 (September 2024)
+-----------------------
+
+This is a minor release with several bugfixes, and major updates to the internal structure and sphinx documentation.
+
+See the full `changelog for 0.5.5 <https://github.com/raphaelvallat/pingouin/releases/tag/v0.5.5>`_.
+
+
 v0.5.4 (January 2024)
 ---------------------
 
-This is a minor release with several bugfixes and no new features. The new version is tested for Python 3.8-3.11 (but should also work with Python 3.12). See `GitHub <https://github.com/raphaelvallat/pingouin/releases>`_ for the full changelog.
+This is a minor release with several bugfixes and no new features. The new version is tested for Python 3.8-3.11 (but should also work with Python 3.12). See the full `changelog for 0.5.4 <https://github.com/raphaelvallat/pingouin/releases/tag/v0.5.4>`_.
 
 This release requires pandas≥1.5. We recommend scipy≥1.11.0.
 

src/pingouin/__init__.py

@@ -20,7 +20,7 @@
 from .config import *
 
 # Current version
-__version__ = "0.5.4"
+__version__ = "0.5.5"
 
 # load default options
 set_default_options()

src/pingouin/multivariate.py

@@ -78,7 +78,8 @@ def multivariate_normality(X, alpha=0.05):
     # Squared-Mahalanobis distances
     Dj = np.diag(np.linalg.multi_dot([difT, S_inv, difT.T]))
     Y = np.linalg.multi_dot([X, S_inv, X.T])
-    Djk = -2 * Y.T + np.repeat(np.diag(Y.T), n).reshape(n, -1) + np.tile(np.diag(Y.T), (n, 1))
+    Y_diag = np.diag(Y)
+    Djk = -2 * Y.T + Y_diag + Y_diag[..., None]
 
     # Smoothing parameter
     b = 1 / (np.sqrt(2)) * ((2 * p + 1) / 4) ** (1 / (p + 4)) * (n ** (1 / (p + 4)))