Skip to content

Commit a219756

Browse files
committed
Fix minor bug in solve cycle
1 parent 08f9724 commit a219756

File tree

1 file changed

+3
-2
lines changed

1 file changed

+3
-2
lines changed

src/multilevel.jl

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -66,11 +66,12 @@ function solve{T}(ml::MultiLevel, b::Vector{T},
6666
tol = 1e-5;
6767
verbose = false,
6868
log = false)
69-
V = promote_type(eltype(ml.levels[1].A), eltype(b))
69+
70+
A = length(ml) == 1 ? ml.final_A : ml.levels[1].A
71+
V = promote_type(eltype(A), eltype(b))
7072
x = zeros(V, size(b))
7173
tol = eltype(b)(tol)
7274
residuals = Vector{V}()
73-
A = length(ml) == 1 ? ml.final_A : ml.levels[1].A
7475
normb = norm(b)
7576
if normb != 0
7677
tol *= normb

0 commit comments

Comments
 (0)