  This article investigates a fast and stable method to solve Henderson's mixed
model equation. The proposed algorithm is stable in that it avoids inverting a
matrix of a large dimension and hence is free from the curse of dimensionality.
This tactic is enabled through row operations performed on the design matrix.
