Most efficient way of writing specific quadratic form

[spinkney]

Better formulation of X' * diag_matrix(v) * X

I do crossprod(diag_pre_multiply(sqrt(v), X)); assuming v is positive.