powered by
Returns XdiagyXT.
fastVDiag(X, y)
A symmetric, square matrix of dimension n by n
n
A vector of length n
A fast implementation of XdiagyXT. This helps speed up the I-prior EM algorithm.