I would like to use the preconditioner in this package with solvers in Krylov.jl. The operator that I would like to use only supports `mul!(y, A, x)`.