@@ -36,14 +36,18 @@ using Test
36
36
@test $ A1 * $ x3 ≈ $ b1
37
37
38
38
y = cache ($ x3, $ A1, $ b2) # reuse factorization
39
- @test $ A1 * y ≈ $ b2
39
+ @test $ A1 * y ≈ $ b2 # with different RHS
40
40
@test $ A1 * $ x3 ≈ $ b2
41
41
42
42
y = cache ($ x3, $ A2, $ b3) # new factorization
43
43
@test $ A2 * y ≈ $ b3 # same old cache
44
44
@test $ A2 * $ x3 ≈ $ b3
45
45
end
46
46
47
+ x1 .= 0.0
48
+ x2 .= 0.0
49
+ x3 .= 0.0
50
+
47
51
return
48
52
end
49
53
@@ -52,33 +56,47 @@ using Test
52
56
:LUFactorization ,
53
57
:QRFactorization ,
54
58
:SVDFactorization ,
55
- :DefaultFactorization ,
56
- :DefaultLinSolve
59
+ # :DefaultLinSolve
57
60
)
58
61
test_interface (alg, kwargs, prob1, prob2, prob3)
59
62
end
60
63
64
+ # alg = :DefaultFactorization
65
+ # for fact_alg in (
66
+ # :lu, :lu!,
67
+ # :qr, :qr!,
68
+ # :cholesky, :cholesky!,
69
+ # # :ldlt, :ldlt!,
70
+ # :bunchkaufman, :bunchkaufman!,
71
+ # :lq, :lq!,
72
+ # :svd, :svd!,
73
+ # :(LinearAlgebra.factorize),
74
+ # )
75
+ # kwargs = :(fact_alg=$fact_alg,)
76
+ # test_interface(alg, kwargs, prob1, prob2, prob3)
77
+ # end
78
+
61
79
# KrylovJL
62
- kwargs = :(ifverbose= true , abstol= 1e-8 , reltol= 1e-8 , maxiter= 30 ,
80
+ kwargs = :(ifverbose= false , abstol= 1e-8 , reltol= 1e-8 , maxiter= 30 ,
63
81
gmres_restart= 5 )
64
82
for alg in (
65
83
:KrylovJL ,
66
84
:KrylovJL_CG ,
67
85
:KrylovJL_GMRES ,
68
- # :KrylovJL_BICGSTAB, # fails
86
+ # :KrylovJL_BICGSTAB,
69
87
:KrylovJL_MINRES ,
70
88
)
71
89
test_interface (alg, kwargs, prob1, prob2, prob3)
72
90
end
73
91
74
92
# IterativeSolversJL
75
- kwargs = :(ifverbose= true , abstol= 1e-8 , reltol= 1e-8 , maxiter= 30 ,
93
+ kwargs = :(ifverbose= false , abstol= 1e-8 , reltol= 1e-8 , maxiter= 30 ,
76
94
gmres_restart= 5 )
77
95
for alg in (
78
96
:IterativeSolversJL ,
79
97
:IterativeSolversJL_CG ,
80
98
:IterativeSolversJL_GMRES ,
81
- # :IterativeSolversJL_BICGSTAB, # fails
99
+ # :IterativeSolversJL_BICGSTAB,
82
100
:IterativeSolversJL_MINRES ,
83
101
)
84
102
test_interface (alg, kwargs, prob1, prob2, prob3)
0 commit comments