This is trying to implement some algorithms for non-negative matrix factorization.