Using chol2inv(chol()) for inverting symmetric matrices (faster than solve())