Set size of zero vector to be size(A, 2)

Wimmerer · Wimmerer · commit 6c9cbc8f7f04 · 2022-01-29T08:43:53.000-05:00
diff --git a/src/common.jl b/src/common.jl
@@ -87,7 +87,7 @@ function SciMLBase.init(prob::LinearProblem, alg::Union{SciMLLinearSolveAlgorith
                        )
     @unpack A, b, u0, p = prob
 
-    u0 = (u0 !== nothing) ? u0 : zero(b)
+    u0 = (u0 !== nothing) ? u0 : zeros(size(A, 2))
 
     cacheval = init_cacheval(alg, A, b, u0, Pl, Pr, maxiters, abstol, reltol, verbose)
     isfresh = true

Original file line number	Diff line number	Diff line change
`@@ -87,7 +87,7 @@ function SciMLBase.init(prob::LinearProblem, alg::Union{SciMLLinearSolveAlgorith`
`87`	`87`	`)`
`88`	`88`	`@unpack A, b, u0, p = prob`
`89`	`89`
`90`		`- u0 = (u0 !== nothing) ? u0 : zero(b)`
	`90`	`+ u0 = (u0 !== nothing) ? u0 : zeros(size(A, 2))`
`91`	`91`
`92`	`92`	`cacheval = init_cacheval(alg, A, b, u0, Pl, Pr, maxiters, abstol, reltol, verbose)`
`93`	`93`	`isfresh = true`