added fast diagonalgaussian stats gathering, likelihood calculation