Efficient algorithms for regularized Poisson Non-negative Matrix Factorization