Non-negative Matrix Factorization code C and CUDA implementations of NMF. Refererence Lin, Chih-Jen. "Projected gradient methods for nonnegative matrix factorization." Neural computation 19.10 (2007): 2756-2779.