Removed computation of xoyo at each step and computed it before entering for loop
Note Armadillo Chol(E)=R such that R.t()*R=X, so Chol(xaxa)=xa such that xa.t()*xa=xaxa BUT NOTE GENERATION OF MVNORMAL, Y=m+LZ REQUIRES LOWER TRIANGULAR