I was erroneously using a scalar for Q. I made it into a matrix. However, I am pretty sure it is still wrong. Q shouldn't be symmetric.
The name was horrible, we were using gaussian.gaussian() to compute gaussians. stats.gaussian() is much better. Updated all notebooks and scripts to reflect the change.