Helper function for updating the precision matrices in the EMFA algorithm.
updatePrec(m, nc, omega, w, p, wNew, pNew, cNew, het, maxDiag)An integer, the order of the model.
An integer, the number of traits or genotypes.
A computed matrix for the current step in the algoritm.
A model matrix for the current step in the algorithm.
A model matrix for the current step in the algorithm.
A pointer to the updated model matrix for w.
A pointer to the updated model matrix for p.
A pointer to the updated matrix c.
Should an extra diagonal part be added in the model for the precision matrix.
A numerical value for the maximum value of sigma2.