Non-negative Matrix Factorization

I found some Python code for Non-negative Matrix Factorization. I've altered it to run more like the NMF code in Programming Collective Intelligence. Since the factorization depends on two random matrices, it might be best to factorize several times, and then cluster a set of H matrices using Affinity Propagation, or K-Means.