handle a few more cases

Author: ChrisRackauckas

src/default.jl

@@ -46,6 +46,15 @@ function defaultalg(A, b::GPUArraysCore.AbstractGPUArray, ::OperatorAssumptions{
     end
 end
 
+function defaultalg(A::SciMLBase.AbstractDiffEqOperator, b, assumptions::OperatorAssumptions)
+    KrylovJL_GMRES()
+end
+
+# Ambiguity handling
+function defaultalg(A::SciMLBase.AbstractDiffEqOperator, b, assumptions::OperatorAssumptions{nothing})
+    KrylovJL_GMRES()
+end
+
 # Handle ambiguity
 function defaultalg(A::GPUArraysCore.AbstractGPUArray, b::GPUArraysCore.AbstractGPUArray,
                     ::OperatorAssumptions{true})
@@ -70,6 +79,8 @@ function defaultalg(A::GPUArraysCore.AbstractGPUArray, b::GPUArraysCore.Abstract
     QRFactorization()
 end
 
+
+
 # Allows A === nothing as a stand-in for dense matrix
 function defaultalg(A, b, ::OperatorAssumptions{true})
     # Special case on Arrays: avoid BLAS for RecursiveFactorization.jl when
@@ -119,6 +130,6 @@ end
 
 function init_cacheval(alg::Nothing, A, b, u, Pl, Pr, maxiters::Int, abstol, reltol,
                        verbose::Bool, assumptions::OperatorAssumptions)
-    init_cacheval(default_alg(A, b), A, b, u, Pl, Pr, maxiters, abstol, reltol, verbose,
+    init_cacheval(default_alg(A, b, assumptions), A, b, u, Pl, Pr, maxiters, abstol, reltol, verbose,
                   assumptions)
 end