The content of factorization.jl is already part of the base julia, but we haven't had the cholpfact method for a while. Instead you can use a keyword to get the pivoted cholesky. E.g.
cholfact(A, pivot = true, tol = 1e-8)
Where tol is the tolerance for the rank determination.