|
152 | 152 | A = A' * A + I |
153 | 153 | tol = √eps(real(T)) |
154 | 154 | b = rand(T, n, 1) |
155 | | - itr = LOBPCGIterator(A, b, largest) |
| 155 | + itr = LOBPCGIterator(A, largest, b) |
156 | 156 |
|
157 | 157 | r = lobpcg!(itr; tol=tol, maxiter=Inf, log=false) |
158 | 158 | λ, X = r.λ, r.X |
|
169 | 169 | B = B' * B + I |
170 | 170 | b = rand(T, n, 1) |
171 | 171 | tol = √eps(real(T)) |
172 | | - itr = LOBPCGIterator(A, B, b, largest) |
| 172 | + itr = LOBPCGIterator(A, B, largest, b) |
173 | 173 |
|
174 | 174 | r = lobpcg!(itr; tol=tol, maxiter=Inf, log=true) |
175 | 175 | λ, X = r.λ, r.X |
|
221 | 221 | r = lobpcg(A, largest, 1; C=copy(r.X), tol=tol, maxiter=Inf, log=false) |
222 | 222 | λ2, X2 = r.λ, r.X |
223 | 223 | @test norm(A*X2 - X2*λ2) ≤ tol |
224 | | - @test isapprox(real(Ac_mul_B(X1, X2)[1,1]), 0, atol=n*tol) |
| 224 | + @test isapprox(real(Ac_mul_B(X1, X2)[1,1]), 0, atol=2*n*tol) |
225 | 225 | end |
226 | 226 | end |
227 | 227 | end |
|
238 | 238 | r = lobpcg(A, B, largest, 1; C=copy(r.X), tol=tol, maxiter=Inf, log=false) |
239 | 239 | λ2, X2 = r.λ, r.X |
240 | 240 | @test norm(A*X2 - B*X2*λ2) ≤ tol |
241 | | - @test isapprox(real(Ac_mul_B(X1, B*X2)[1,1]), 0, atol=n*tol) |
| 241 | + @test isapprox(real(Ac_mul_B(X1, B*X2)[1,1]), 0, atol=2*n*tol) |
242 | 242 | end |
243 | 243 | end |
244 | 244 | end |
|
298 | 298 | r = lobpcg(A, largest, X0, 3, tol=tol, maxiter=Inf, log=true) |
299 | 299 | λ, X = r.λ, r.X |
300 | 300 | @test max_err(A*X - X*diagm(λ)) ≤ tol |
301 | | - @test all(isapprox.(Ac_mul_B(X, X), eye(3), atol=n*tol)) |
| 301 | + @test all(isapprox.(Ac_mul_B(X, X), eye(3), atol=2*n*tol)) |
302 | 302 | end |
303 | 303 | end |
304 | 304 | end |
|
315 | 315 | r = lobpcg(A, B, largest, X0, 3, tol=tol, maxiter=Inf, log=true) |
316 | 316 | λ, X = r.λ, r.X |
317 | 317 | @test max_err(A*X - B*X*diagm(λ)) ≤ tol |
318 | | - @test all(isapprox.(Ac_mul_B(X, B*X), eye(3), atol=n*tol)) |
| 318 | + @test all(isapprox.(Ac_mul_B(X, B*X), eye(3), atol=2*n*tol)) |
319 | 319 | end |
320 | 320 | end |
321 | 321 | end |
|
0 commit comments