Convex non-negative matrix factorization for massive datasets