split along highest variance, faster diag llpg
faster implementation of rand()
now use fast version for diag llpg
include/require problem, moved all includes back.
Changed full covariance to list of matrices.
fix avll(), add em! after kmeans init.