Generic LDL' factorization

[ViralBShah]

From: JuliaLang/julia#10953

@mschauer wrote:

It would be nice to have a the semidefinite Cholesky or LDL' decomposition, which allows semidefinite matrices, is more stable and avoids computation of square roots. This makes most sense as generic julia algorithm, because for BLAS floats pivoted Cholesky does the trick. I remember that pivoting was difficult for abstract matrices.