separate gauss utility function
maintain axis order for eig, change values to diag matrix