GaussianMixture in scikit-learn triggers deadlocks

[fpavogt]

Describe the bug

The ``GaussianMixture'' routine triggers deadlocks, which is a know feature apparently: scikit-learn/scikit-learn#19272

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

[fpavogt]

Some more info about this underlying OpenMP issue from the scikit-learn crowd.

They tried to submit a patch to OpenMP apparently, but it got rejected. Seems that the real issue is with python's "multiprocessing" package.

Interesting quote from the OpenMP patch thread.

The really nasty case is when the code using omp and the code using fork() are in entirely different pieces of code, that don't know about each other (e.g., different shared libraries).

Which I believe is exactly what is happening to us.

[fpavogt]

This turns out to be a feature of Python. As such, I don't plan to fix it anytime soon.