diag-formula-method function

Diagonalize the random-effect covariance structure, possibly assisting convergence