Default non-square operators to LSMR

Author: ChrisRackauckas

src/default.jl

@@ -69,6 +69,11 @@ function defaultalg(A::SciMLBase.AbstractDiffEqOperator, b,
     KrylovJL_GMRES()
 end
 
+function defaultalg(A::SciMLBase.AbstractDiffEqOperator, b,
+    assumptions::OperatorAssumptions{false})
+    KrylovJL_LSMR()
+end
+
 # Handle ambiguity
 function defaultalg(A::GPUArraysCore.AbstractGPUArray, b::GPUArraysCore.AbstractGPUArray,
                     ::OperatorAssumptions{true})

src/iterative_wrappers.jl

@@ -24,6 +24,10 @@ function KrylovJL_MINRES(args...; kwargs...)
     KrylovJL(args...; KrylovAlg = Krylov.minres!, kwargs...)
 end
 
+function KrylovJL_LSMR(args...; kwargs...)
+    KrylovJL(args...; KrylovAlg = Krylov.lsmr!, kwargs...)
+end
+
 function get_KrylovJL_solver(KrylovAlg)
     KS = if (KrylovAlg === Krylov.lsmr!)
         Krylov.LsmrSolver