The help for ldltfact says: Base.ldltfact(A) -> LDLtFactorization
Compute a factorization of a positive definite matrix "A" such
that "A=L*Diagonal(d)*L'" where "L" is a unit lower triangular
matrix and "d" is a vector with non-negative elements.
But for a general positive definite matrix it does not work.
In fact the header of the function is
function ldltfact{T}(M::SymTridiagonal{T})
Or I miss something?
Giovanni
