critical bugfix wish shapes in gaussian_likelihood.c
stats['obs*obs.T'] is being computed correctly now