Interatively factorize a matrix into three matrices with linear constraint(Matrix is large)