\begin{Schunk}
\begin{Sinput}
 eS <- eigen(Sigma, symmetric = TRUE)
 ev <- eS$values
 if (!all(ev >= -tol * abs(ev[1L])))
          stop("'Sigma' is not positive definite")
\end{Sinput}
\end{Schunk}
