A procedure is presented for rapidly computing the diagnonal elements of a large numerator relationship matrix, say $\mathbf {A}$. It also generates a lower triangular matrix, say $\mathbf {L}$, ...