A Symmetric Rank-One Quasi-Newton Method for Nonnegative Matrix Factorization

Figure 7

Reconstructed images after 10 iterations.