numerical stability in matrix factorization