Hello!
My question is relatively straightforward.
I’d like to use the LU factorization of the adjoint of a sparse matrix of type CuSparseMatrixCSR with the package CUDSS.jl, e.g. lu(A'), by only having access to lu(A).
I’m looking for a method such as
ldiv!(X::StridedVecOrMat{T}, adjlu::AdjointFactorization{T,<:UmfpackLU{T}}, B::StridedVecOrMat{T}) where {T<:UMFVTypes},
which exists for “conventional” sparse matrices. In this case, everything is abstracted away so that lu(A')=lu(A)'.
I understand this high-level feature may not be implemented just yet, but I wonder if there is a simple workaround to achieve the same.