Given a matrix , the problem seeks nonnegative matrices and such that or .

The model is

The algorithm is described in Algorithm 2 of this report.

Matlab function

Synthetic example