We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 7ccbae7 commit 8dd1489Copy full SHA for 8dd1489
test/lalqmr.jl
@@ -57,6 +57,7 @@ end
57
# Guaranteed to create blocks during Lanczos process
58
# This satisfies the condition that in the V-W sequence, the first
59
# iterates are orthogonal: <Av - v<A, v>, Atv - v<At, v>> under transpose inner product
60
+ # These do _not_ work for regular qmr
61
dl = fill(one(T), n-1)
62
du = fill(one(T), n-1)
63
d = fill(one(T), n)
@@ -72,9 +73,9 @@ end
72
73
end
74
75
@testset "Linear operator defined as a function" begin
- A = LinearMap(cumsum!, 100; ismutating=true)
76
+ A = LinearMap(identity, identity, 100; ismutating=false)
77
b = rand(100)
- reltol = 1e-5
78
+ reltol = 1e-6
79
80
x = lalqmr(A, b; reltol=reltol, maxiter=2000)
81
@test norm(A * x - b) / norm(b) ≤ reltol
0 commit comments