Fallback for non-square A

SciML · Nov 8, 2023 · 13927ce · 13927ce
1 parent e3c63c9
commit 13927ce
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/factorization.jl b/src/factorization.jl
@@ -12,6 +12,8 @@ _ldiv!(x, A, b) = ldiv!(x, A, b)
 
 function _ldiv!(x::Vector, A::Factorization, b::Vector)
     # workaround https://github.com/JuliaLang/julia/issues/43507
+    # Fallback if working with non-square matrices
+    length(x) != length(b) && return ldiv!(x, A, b)
     copyto!(x, b)
     ldiv!(A, x)
 end