Computation of all eigenvalues of matrices used in restricted maximum likelihood estimation of variance components using sparse matrix techniques C. Robert and V. DucrocqGenet. Sel. Evol., 28 1 (1996) 51-65DOI: https://doi.org/10.1051/gse:19960103