Merge pull request #102 from Wimmerer/nonsquare_sizes

Nonsquare matrix fix
SciML · Jan 29, 2022 · acbc6d8 · acbc6d8
2 parents 3ba3bc3 + e4d13ce
commit acbc6d8
Showing 1 changed file with 6 additions and 1 deletion.
diff --git a/src/common.jl b/src/common.jl
@@ -87,7 +87,12 @@ function SciMLBase.init(prob::LinearProblem, alg::Union{SciMLLinearSolveAlgorith
                        )
     @unpack A, b, u0, p = prob
 
-    u0 = (u0 !== nothing) ? u0 : zero(b)
+    u0 = if u0 !== nothing
+        u0 
+    else
+        u0 = similar(b, size(A, 2))
+        fill!(u0,false)
+    end
 
     cacheval = init_cacheval(alg, A, b, u0, Pl, Pr, maxiters, abstol, reltol, verbose)
     isfresh = true