-
Notifications
You must be signed in to change notification settings - Fork 5
/
1.5.0-DEV-7206b56e94.log
201 lines (196 loc) · 6.86 KB
/
1.5.0-DEV-7206b56e94.log
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
Julia Version 1.5.0-DEV.485
Commit 7206b56e94 (2020-03-18 17:25 UTC)
Platform Info:
OS: Linux (x86_64-linux-gnu)
CPU: Intel(R) Xeon(R) Silver 4114 CPU @ 2.20GHz
WORD_SIZE: 64
LIBM: libopenlibm
LLVM: libLLVM-9.0.1 (ORCJIT, skylake)
Environment:
JULIA_DEPOT_PATH = ::/usr/local/share/julia
JULIA_NUM_THREADS = 2
Resolving package versions...
Installed KrylovMethods ─ v0.6.0
Updating `~/.julia/environments/v1.5/Project.toml`
9a2cd570 + KrylovMethods v0.6.0
Updating `~/.julia/environments/v1.5/Manifest.toml`
9a2cd570 + KrylovMethods v0.6.0
8f399da3 + Libdl
37e2e46d + LinearAlgebra
de0858da + Printf
9a3f8284 + Random
9e88b42a + Serialization
2f01184e + SparseArrays
4ec0a83e + Unicode
Testing KrylovMethods
Status `/tmp/jl_hCs73n/Project.toml`
9a2cd570 KrylovMethods v0.6.0
5c8ed15e LinearOperators v1.0.1
37e2e46d LinearAlgebra
de0858da Printf
2f01184e SparseArrays
8dfed614 Test
Status `/tmp/jl_hCs73n/Manifest.toml`
9aa1b823 FastClosures v0.3.2
9a2cd570 KrylovMethods v0.6.0
5c8ed15e LinearOperators v1.0.1
a759f4b9 TimerOutputs v0.5.3
2a0f44e3 Base64
8ba89e20 Distributed
b77e0a4c InteractiveUtils
8f399da3 Libdl
37e2e46d LinearAlgebra
56ddb016 Logging
d6f4376e Markdown
de0858da Printf
9a3f8284 Random
9e88b42a Serialization
6462fe0b Sockets
2f01184e SparseArrays
8dfed614 Test
4ec0a83e Unicode
=== bicgstb ===
iter relres
1 1.00e-01
2 1.12e-02
3 8.94e-04
bicgstb iterated maxIter (=3) times but reached only residual norm 1.12e-02 instead of tol=1.00e-06.
bcgstb achieved desired tolerance at iteration 6. Residual norm is 7.46e-07.
=== block bicgstb ===
iter relres
1 9.49e-02
2 1.52e-02
3 4.03e-03
blockBiCGSTB iterated maxIter (=3) times but reached only residual norm 1.52e-02 instead of tol=1.00e-06.
blockBiCGSTB achieved desired tolerance at iteration 6.5. Residual norm is 8.89e-07.
=== cg ===
iter relres
1 5.55e-17
cg achieved desired tolerance at iteration 1. Residual norm is 5.55e-17.
=== cg ===
iter relres
1 2.51e-01
cg iterated maxIter (=1) times but reached only residual norm 2.51e-01 instead of tol=1.00e-15.
Matrix A in cg has to be positive definite.
cg achieved desired tolerance at iteration 24. Residual norm is 9.34e-06.
=== blockcg ===
iter max(relres)
1 3.06e-15
2 2.48e-28
blockCG achieved desired tolerance at iteration 2. Residual norm is 2.48e-28.
=== blockcg ===
iter max(relres)
1 1.56e-16
blockCG iterated maxIter (=1) times but reached only residual norm 1.56e-16 instead of tol=1.00e-30.
=== blockcg ===
iter max(relres)
1 3.97e-01
2 2.29e-01
3 1.68e-01
4 1.22e-01
5 8.93e-02
6 7.05e-02
7 4.36e-02
8 2.34e-02
9 1.19e-02
10 5.95e-03
11 2.92e-03
12 1.28e-03
13 5.23e-04
blockCG achieved desired tolerance at iteration 13. Residual norm is 5.23e-04.
=== blockcg ===
iter max(relres)
1 3.97e-01
2 2.29e-01
3 1.68e-01
4 1.22e-01
5 8.93e-02
6 7.05e-02
7 4.36e-02
8 2.34e-02
9 1.19e-02
10 5.95e-03
11 2.92e-03
12 1.28e-03
13 5.23e-04
blockCG achieved desired tolerance at iteration 13. Residual norm is 5.23e-04.
=== cgls ===
iter |A'r| norm(r) norm(x)
1 1.25e+01 5.48e-01 9.33e+00
2 2.85e+00 6.97e-01 9.27e+00
cgls iterated maxIter (=2) times witout reaching tolerance. Returned result has residual norm 2.85e+00.
cgls relative gradient norm below desired tolerance at iteration 10. Returned result has residual norm 5.62e-11.
=== gmres ===
iter relres
1 3.1e-01 1.0e-01 3.0e-02 3.0e-02
2 9.4e-03 2.8e-03 8.0e-04 8.0e-04
3 2.2e-04 6.6e-05 1.8e-05 1.8e-05
gmres iterated maxIter (=3) times without achieving the desired tolerance.
gmres achieved desired tolerance at iteration 3. Residual norm is 5.68e-07.
*******************************************************
=== fgmres ===
iter relres
1 3.1e-01 1.0e-01 3.0e-02
2 9.4e-03 2.8e-03 8.0e-04
3 2.2e-04 6.6e-05 1.8e-05
fgmres iterated maxIter (=3) outer iterations without achieving the desired tolerance. Acheived: 1.80e-05
fgmres achieved desired tolerance at inner iteration 12. Residual norm is 5.68e-07.
*******************************************************
=== blockFGMRES ===
iter relres
1 2.0e-01 4.2e-02 8.3e-03
2 1.6e-03 3.1e-04 5.8e-05
3 1.2e-05 2.6e-06 4.7e-07
blockFGMRES iterated maxIter (=3) outer iterations without achieving the desired tolerance. Acheived: 4.70e-07
-----------------------------------
blockFGMRES achieved desired tolerance at inner iteration 9. Residual norm is 5.18e-07.
*****************************************************************************
Givens rotation for v=[0.0,0.5973099101664067]
Givens rotation for v=[-0.1466461410642105,0.0]
Givens rotation for v=[-1.542348888485669,0.027894366541275903]
=== Testing Lanczos Bidiagonalization ===
ssor achieved desired tolerance at iteration 34. Residual norm is 9.81e-03.
ssor achieved desired tolerance at iteration 50. Residual norm is 1.00e-02.
cg achieved desired tolerance at iteration 54. Residual norm is 6.99e-13.
=== sor ===
iter relres
1 1.44e-01
2 6.52e-02
3 4.08e-02
ssor iterated maxIter (=3) times but reached only residual norm 4.08e-02 instead of tol=1.00e-02.
ssor achieved desired tolerance at iteration 16. Residual norm is 9.66e-03.
ssor achieved desired tolerance at iteration 24. Residual norm is 9.54e-03.
iter |r_k| |A'*r_k| |A| cond(A) STOP
001 4.629e-01 4.792e-01 1.673e+00 1.000e+00 [0,0,0]
002 2.979e-16 4.965e-16 2.000e+00 2.309e+00 [1,0,0]
lsqr achieved desired tolerance at iteration 2. The system should be compatible. Residual norm is 2.98e-16.
lsqr achieved desired tolerance at iteration 2. The system should be incompatible. Gradient norm is 1.15e+00.
lsqr achieved desired tolerance at iteration 10. The system should be incompatible. Gradient norm is 9.26e+00.
minres iterated maxIter (=1) times but reached only residual norm 3.48e-01 instead of tol=1.00e-12.
=== minres ===
iter |r_k| |A*r_k| |A| cond(A)
1 3.48e-01 1.98e+00 5.99e+00 1.19e+00
2 8.45e-02 5.91e-01 1.08e+01 2.14e+00
3 1.52e-02 6.43e-02 1.08e+01 2.14e+00
4 4.53e-03 3.43e-02 1.08e+01 2.14e+00
5 1.04e-03 4.94e-03 1.08e+01 2.14e+00
6 1.62e-04 4.51e-04 1.08e+01 2.14e+00
7 5.56e-19 3.38e-32 1.08e+01 2.14e+00
minres converged at iteration 7. |A r_k|=3.38e-32 and |r_k|=5.56e-19.
=== minres ===
iter |r_k| |A*r_k| |A| cond(A)
1 1.42e-01 1.78e+00 9.43e+00 1.01e+00
2 3.14e-02 1.17e-01 9.43e+00 1.67e+00
3 6.43e-03 3.35e-02 9.43e+00 1.67e+00
4 1.06e-03 4.72e-03 9.43e+00 1.67e+00
5 2.46e-04 9.77e-04 9.43e+00 1.77e+00
6 1.27e-06 1.15e-07 9.43e+00 1.77e+00
7 6.80e-18 6.42e-28 9.43e+00 1.77e+00
minres converged at iteration 7. |A r_k|=6.42e-28 and |r_k|=6.80e-18.
test stop by condlim
test stop by beta
minres stopped because beta (=4.65e-14) became smaller than the tolerance (btol=1.00e-10) at iteration 5.
test if anything is done for zero right hand side
Test Summary: | Pass Total
KrylovMethods | 155 155
Testing KrylovMethods tests passed