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