Fix defaultalg method ambiguity for SciMLOp + GPUArray

should fix this error
```julia
Candidates:                                                                                                                                                        
  defaultalg(A::SciMLOperators.AbstractSciMLOperator, b, assump::LinearSolve.OperatorAssumptions)
    @ LinearSolve ~/.julia/packages/LinearSolve/Rq2MY/src/default.jl:132                                                                                           
  defaultalg(A, b::GPUArraysCore.AbstractGPUArray, assump::LinearSolve.OperatorAssumptions)                                   
    @ LinearSolve ~/.julia/packages/LinearSolve/Rq2MY/src/default.jl:106

Possible fix, define
  defaultalg(::SciMLOperators.AbstractSciMLOperator, ::GPUArraysCore.AbstractGPUArray, ::LinearSolve.OperatorAssumptions)

```

src/default.jl

@@ -129,7 +129,7 @@ function defaultalg(A::GPUArraysCore.AbstractGPUArray, b::GPUArraysCore.Abstract
     end
 end
 
-function defaultalg(A::SciMLBase.AbstractSciMLOperator, b,
+function _defaultalg_scimlops(A::SciMLBase.AbstractSciMLOperator, b,
     assump::OperatorAssumptions)
     if has_ldiv!(A)
         return DefaultLinearSolver(DefaultAlgorithmChoice.DirectLdiv!)
@@ -145,6 +145,9 @@ function defaultalg(A::SciMLBase.AbstractSciMLOperator, b,
     end
 end
 
+defaultalg(A::SciMLBase.AbstractSciMLOperator, b, assump::OperatorAssumptions) = _defaultalg_scimlops(A, b, assump)
+defaultalg(A::SciMLBase.AbstractSciMLOperator, b::GPUArraysCore.AbstractGPUArray, assump::OperatorAssumptions) = _defaultalg_scimlops(A, b, assump)
+
 # Allows A === nothing as a stand-in for dense matrix
 function defaultalg(A, b, assump::OperatorAssumptions)
     alg = if assump.issq