like algorithm to enhance accuracy. This algorithm begins by a suitable Cholesky-like
factorization of A, A= GT JG. In some applications, A is given implicitly in that form and its
natural Cholesky-like factor G is immediately available, but “tall”, ie, not of full row rank. This
factor G is unsuitable for the Jacobi-like process. To avoid explicit computation of A, and
possible loss of accuracy, the factor has to be preprocessed by a QR-like factorization. In this …