GaussianBP algorithm does not converge for some specific precision matrices (inverse covariances)

[mehta-a]

site-packages/numpy/linalg/linalg.py", line 88, in _raise_linalgerror_singular
raise LinAlgError("Singular matrix")
numpy.linalg.LinAlgError: Singular matrix