What if the covariance is a singular matrix?

[zxdawn]

Sometimes, we may get a singular matrix as the covariance. The inv_cov won't work. Is it correct to use np.linalg.pinv instead of np.linalg.inv for this kind of case?

[zxdawn]

Any idea on this? It's easy to fix it in one line ;)