improved speed of NMF by using crossprod (see R News, 4/1:17-20)