Add L-bfgs solver for NMF