Data class for LLT (Cholesky) decomposition, with rank-1 update/downdate.
Lower-triangular matrix L such that mat = L*L^T
In-place rank-1 downdate: mat := mat - x*x^T
In-place rank-1 update: mat := mat + x*x^T