ldiv

Author: vpuri3

src/default.jl

@@ -58,11 +58,11 @@ function SciMLBase.solve(cache::LinearCache, alg::Nothing,
         A = A.A
     end
 
-    if applicable(ldiv!, A, u)
-        alg = FunctionCall(ldiv!, (:A, :u))
+    if applicable(ldiv!, u, A, b)
+        alg = LdivBang3Args()
         SciMLBase.solve(cache, alg, args...; kwargs...)
-    elseif applicable(ldiv!, u, A, b)
-        alg = FunctionCall(ldiv!, (:u, :A, :b))
+    elseif applicable(ldiv!, A, u)
+        alg = LdivBang2Args()
         SciMLBase.solve(cache, alg, args...; kwargs...)
     end
 

src/function_call.jl

@@ -7,6 +7,9 @@ struct FunctionCall{F,A} <: SciMLLinearSolveAlgorithm
     end
 end
 
+LdivBang2Args() = FunctionCall(ldiv!, (:A, :u))
+LdivBang3Args() = FunctionCall(ldiv!, (:u, :A, :b))
+
 function (f::FunctionCall)(cache::LinearCache)
     @unpack func!, argsyms = f
     args = [getproperty(cache,argsym) for argsym in argsyms]