A flexible method is introduced to model the structure of a covariance matrix C and study the dependence of the covariances on explanatory variables by observing that for any real symmetric matrix A, ...

If the matrix is symmetric, this splitting is combined with the conjugate-gradient method to provide a fast iterative solution algorithm. Comparisons have been made with other well-known methods.