Add RecursiveFactorization options

Requires JuliaLinearAlgebra/RecursiveFactorization.jl#39

Author: ChrisRackauckas

Project.toml

@@ -27,7 +27,7 @@ IterativeSolvers = "0.9.2"
 KLU = "0.2.3"
 Krylov = "0.7.9"
 KrylovKit = "0.5"
-RecursiveFactorization = "0.2"
+RecursiveFactorization = "0.2.8"
 Reexport = "1"
 Requires = "1"
 SciMLBase = "1.25"

src/factorization.jl

@@ -224,6 +224,12 @@ end
 
 ## RFLUFactorization
 
-RFLUFactorization() = GenericFactorization(;fact_alg=RecursiveFactorization.lu!)
+struct RFWrapper{P,T}
+    RFWrapper(::Val{P},::Val{T}) where {P,T} = new{P,T}()
+end
+(::RFWrapper{P,T})(A) where {P,T} = RecursiveFactorization.lu!(A,val(P),val(T))
+
+RFLUFactorization(;pivot = Val(true), thread = Val(true)) = GenericFactorization(;fact_alg=RFWrapper(pivot,thread))
+
 init_cacheval(alg::GenericFactorization{typeof(RecursiveFactorization.lu!)}, A, b, u, Pl, Pr, maxiters, abstol, reltol, verbose) = ArrayInterface.lu_instance(convert(AbstractMatrix,A))
 init_cacheval(alg::GenericFactorization{typeof(RecursiveFactorization.lu!)}, A::StridedMatrix{<:LinearAlgebra.BlasFloat}, b, u, Pl, Pr, maxiters, abstol, reltol, verbose) = ArrayInterface.lu_instance(convert(AbstractMatrix,A))