Calculation of the inverse of the error variance-covariance matrix is required for both feasible generalized least squares and maximum-likelihood estimation of the regression parameters in the two-way error-components model. Since in many applications this matrix can be quite large, efficient...