Wrappers for incomplete factorizations from the CUSPARSE library together with convenience functions to use them as valid preconditioners in most Julia iterative solvers.

Based on the tutorial from Krlov.jl https://juliasmoothoptimizers.github.io/Krylov.jl/stable/gpu/#Example-with-a-symmetric-positive-definite-system

ic0(A:<{CuSparseMatrixCSC, CuSparseMatrixCSR})

Compute the incomplete Cholesky factorization of A using the CUSPARSE library, returning a struct for which ldiv! can be used as a preconditioner.

ilu0(A:<{CuSparseMatrixCSC, CuSparseMatrixCSR})

Compute the incomplete LU factorization of A using the CUSPARSE library, returning a struct for which ldiv! can be used as a preconditioner.