emmreml function

Solver for Gaussian mixed model with known covariance structure.