Using solve() with default precision rather than call to
chol2inv(chol()) for consistency (mixing can halt with error)