Solved by al3x in post #3 
              
              
                You could try working in log space (i.e. log(det(Sigma))), which typically results in much better stability when working with products of small numbers. 
In the case of a diagonal covariance matrix, the log determinant is the sum of the log of each diagonal element.