Using CHOLMOD.Factor as preconditionner

Dear all, I have recently tried to use Cholesky factors (of simpler problem) as a preconditionner to a Krylov solve, e.g.: ```julia PC_fact = cholesky(PC) gmres!(x, A, b; Pl=PC_fact) ``` This fails as `ldiv!` is not defined for `CHOLMOD.Factors`: ``` ERROR: MethodError: no method matching ldiv!(::SparseArrays.CHOLMOD.Factor{Float64}, ``` The issue is discussed [there] in more details (https://discourse.julialang.org/t/defining-a-preconditionner-for-iterativesolvers/86977). A temporary solution (defining `ldiv!` for `CHOLMOD.Factors` manually) is presented in that post but is there now a better approach to overcomes this? thanks!