We consider the problem of factorizing matrices into two lower dimensional nonnegative matrices. This problem has gathered lots of attention and has been applied in many disciplines including pattern recognition and data analytics. We discuss approaches to find a global optimal solution of the problem.

(This is a joint work with Henry Wolkowicz.)