Block coordinate descent for Sparse NMF