Use Cholesky decomposition, and call gemv to compute L*x
Rethrow InternalException when factorization failed

This is needed by TemporalNormalProcess to detect that matrix
is not symmetric definite positive, and try to regularize it.
Add CovarianceHMatrix-AdmissibilityFactor and CovarianceHMatrix-CompressionMethod keys to ResourceMap