Efficiently sample multivariate normal using precision matrix from x N(Q−1a,Q−1)x ~ N(Q^{-1}a, Q^{-1})x N(Q−1a,Q−1), where Q−1Q^{-1}Q−1 is the precision matrix
Useful links