GaussianMixture in scikit-learn triggers deadlocks

**Describe the bug**

The ``GaussianMixture'' routine triggers deadlocks, which is a know feature apparently: https://github.com/scikit-learn/scikit-learn/issues/19272

Can we find an alternative method (possibly from another package) that is not subject to this ?

Start by reviewing how GaussianMixture is used and the referenced scikit-learn issue #19272 about deadlocks. Determine whether a viable alternative method or package can replace it, then verify that the replacement avoids the reported deadlock behavior.

Written by the indexing model from the issue text.