[petsc-users] GAMG and zero pivots follow up
David Knezevic
david.knezevic at akselos.com
Fri Nov 13 15:45:50 CST 2015
To follow up on this, I have run some tests on a smaller test case (same
model, coarser mesh). The options I tried are:
8 MPI processes with:
"-ksp_type gmres -pc_type gamg"
"-ksp_type fgmres -pc_type gamg"
"-ksp_type cg -pc_type gamg"
"-ksp_type cg -pc_type gamg -ksp_norm_type natural"
"-ksp_type cg -pc_type gamg -ksp_norm_type natural -pc_mg_type full"
4 MPI processes with:
"-ksp_type cg -pc_type gamg -ksp_norm_type natural -pc_mg_type full"
"-ksp_type cg -pc_type gamg"
"-ksp_type gmres -pc_type gamg"
"-ksp_type fgmres -pc_type gamg"
Let me know if you'd like me to try any other cases.
The "-ksp_monitor -ksp_view" output is attached. As you can see, all of the
"4 MPI processes" cases worked well, and none of the "8 MPI processes"
cases converged.
Best regards,
David
On Fri, Nov 13, 2015 at 2:48 PM, Matthew Knepley <knepley at gmail.com> wrote:
> On Fri, Nov 13, 2015 at 12:28 PM, Jed Brown <jed at jedbrown.org> wrote:
>
>> Matthew Knepley <knepley at gmail.com> writes:
>> > Something very strange is happening here. CG should converge
>> monotonically,
>> > but above it does not. What could be happening?
>>
>> Are you use -ksp_norm_type natural? CG is not monotone in other norms.
>>
>
> Yikes! I did not check that. Why do we have PRECONDITIONED as the default
> for CG?
>
> Matt
>
> Also, if boundary conditions are enforced using a nonsymmetric
>> formulation (for example), then you can get lack of monotonicity with CG
>> that may not be catastrophic.
>>
>
>
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which their
> experiments lead.
> -- Norbert Wiener
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20151113/ec852507/attachment-0001.html>
-------------- next part --------------
----------------------------------------------------------------------------------------------------------
8 MPI processes with "-ksp_type gmres -pc_type gamg -ksp_max_it 200":
0 KSP Residual norm 1.243143122667e-03
1 KSP Residual norm 6.790564176716e-04
2 KSP Residual norm 2.863642167759e-04
3 KSP Residual norm 2.169460015654e-04
4 KSP Residual norm 1.623823395363e-04
5 KSP Residual norm 1.115831456626e-04
6 KSP Residual norm 7.905550178072e-05
7 KSP Residual norm 6.676506348708e-05
8 KSP Residual norm 5.411318965562e-05
9 KSP Residual norm 4.946393535811e-05
10 KSP Residual norm 4.873373148542e-05
11 KSP Residual norm 4.873372878673e-05
12 KSP Residual norm 4.832393150132e-05
13 KSP Residual norm 4.682644869071e-05
14 KSP Residual norm 4.177965048741e-05
15 KSP Residual norm 3.561738393315e-05
16 KSP Residual norm 3.183178450210e-05
17 KSP Residual norm 2.820849441834e-05
18 KSP Residual norm 2.411305833934e-05
19 KSP Residual norm 2.073531106031e-05
20 KSP Residual norm 1.832253875945e-05
21 KSP Residual norm 1.613725457732e-05
22 KSP Residual norm 1.447115239529e-05
23 KSP Residual norm 1.332661204650e-05
24 KSP Residual norm 1.248919278483e-05
25 KSP Residual norm 1.196188151016e-05
26 KSP Residual norm 1.161737695363e-05
27 KSP Residual norm 1.153632298642e-05
28 KSP Residual norm 1.153553211867e-05
29 KSP Residual norm 1.150729646659e-05
30 KSP Residual norm 1.134640588584e-05
31 KSP Residual norm 1.125651483355e-05
32 KSP Residual norm 1.121985823782e-05
33 KSP Residual norm 1.121682797994e-05
34 KSP Residual norm 1.120526536096e-05
35 KSP Residual norm 1.112698441144e-05
36 KSP Residual norm 1.099161361515e-05
37 KSP Residual norm 1.077160664786e-05
38 KSP Residual norm 1.046319447066e-05
39 KSP Residual norm 1.002732866634e-05
40 KSP Residual norm 9.687406818053e-06
41 KSP Residual norm 9.291736292845e-06
42 KSP Residual norm 8.787280517217e-06
43 KSP Residual norm 8.323595238657e-06
44 KSP Residual norm 7.891080867185e-06
45 KSP Residual norm 7.537064831605e-06
46 KSP Residual norm 7.316511381129e-06
47 KSP Residual norm 7.185951262668e-06
48 KSP Residual norm 7.117216131634e-06
49 KSP Residual norm 7.104770082988e-06
50 KSP Residual norm 7.099139305978e-06
51 KSP Residual norm 7.038487040610e-06
52 KSP Residual norm 6.861458611935e-06
53 KSP Residual norm 6.625293689513e-06
54 KSP Residual norm 6.329429663991e-06
55 KSP Residual norm 5.919056214664e-06
56 KSP Residual norm 5.507182558921e-06
57 KSP Residual norm 5.196803884679e-06
58 KSP Residual norm 5.002188092285e-06
59 KSP Residual norm 4.880759791404e-06
60 KSP Residual norm 4.770150595855e-06
61 KSP Residual norm 4.724469615685e-06
62 KSP Residual norm 4.673829760077e-06
63 KSP Residual norm 4.629705280910e-06
64 KSP Residual norm 4.601474765626e-06
65 KSP Residual norm 4.593132269745e-06
66 KSP Residual norm 4.593013889961e-06
67 KSP Residual norm 4.587628601477e-06
68 KSP Residual norm 4.552820908762e-06
69 KSP Residual norm 4.477855982146e-06
70 KSP Residual norm 4.405304333703e-06
71 KSP Residual norm 4.330447444642e-06
72 KSP Residual norm 4.237398563528e-06
73 KSP Residual norm 4.138174613148e-06
74 KSP Residual norm 4.031940389494e-06
75 KSP Residual norm 3.924707157992e-06
76 KSP Residual norm 3.802185445933e-06
77 KSP Residual norm 3.721305730027e-06
78 KSP Residual norm 3.679963259079e-06
79 KSP Residual norm 3.667845615364e-06
80 KSP Residual norm 3.667179799479e-06
81 KSP Residual norm 3.662313644020e-06
82 KSP Residual norm 3.638833884448e-06
83 KSP Residual norm 3.598532435205e-06
84 KSP Residual norm 3.535852321120e-06
85 KSP Residual norm 3.456701541505e-06
86 KSP Residual norm 3.365433403050e-06
87 KSP Residual norm 3.271989106911e-06
88 KSP Residual norm 3.176356037348e-06
89 KSP Residual norm 3.034288347877e-06
90 KSP Residual norm 2.938417615975e-06
91 KSP Residual norm 2.900711516997e-06
92 KSP Residual norm 2.869381676648e-06
93 KSP Residual norm 2.855261464067e-06
94 KSP Residual norm 2.850036755015e-06
95 KSP Residual norm 2.849803273543e-06
96 KSP Residual norm 2.849494112992e-06
97 KSP Residual norm 2.846681415018e-06
98 KSP Residual norm 2.838146221432e-06
99 KSP Residual norm 2.823475351161e-06
100 KSP Residual norm 2.804320481647e-06
101 KSP Residual norm 2.775770817616e-06
102 KSP Residual norm 2.740056888907e-06
103 KSP Residual norm 2.691643851812e-06
104 KSP Residual norm 2.631661051625e-06
105 KSP Residual norm 2.579149456214e-06
106 KSP Residual norm 2.544055525393e-06
107 KSP Residual norm 2.525575063547e-06
108 KSP Residual norm 2.511478767787e-06
109 KSP Residual norm 2.505297851010e-06
110 KSP Residual norm 2.504876255779e-06
111 KSP Residual norm 2.504461173319e-06
112 KSP Residual norm 2.500656666173e-06
113 KSP Residual norm 2.492291225825e-06
114 KSP Residual norm 2.473685710456e-06
115 KSP Residual norm 2.444819997807e-06
116 KSP Residual norm 2.412359679152e-06
117 KSP Residual norm 2.381936043848e-06
118 KSP Residual norm 2.360962621352e-06
119 KSP Residual norm 2.339143019400e-06
120 KSP Residual norm 2.312025372293e-06
121 KSP Residual norm 2.295275489137e-06
122 KSP Residual norm 2.280756556273e-06
123 KSP Residual norm 2.267434577949e-06
124 KSP Residual norm 2.258143104977e-06
125 KSP Residual norm 2.252859925246e-06
126 KSP Residual norm 2.250650586520e-06
127 KSP Residual norm 2.250303886635e-06
128 KSP Residual norm 2.250181880096e-06
129 KSP Residual norm 2.248447177342e-06
130 KSP Residual norm 2.244592613442e-06
131 KSP Residual norm 2.237578766797e-06
132 KSP Residual norm 2.225002721330e-06
133 KSP Residual norm 2.208977931816e-06
134 KSP Residual norm 2.190972169698e-06
135 KSP Residual norm 2.163454716687e-06
136 KSP Residual norm 2.113637340502e-06
137 KSP Residual norm 2.056404468594e-06
138 KSP Residual norm 2.018194032501e-06
139 KSP Residual norm 2.002756982750e-06
140 KSP Residual norm 1.998776984922e-06
141 KSP Residual norm 1.998751658408e-06
142 KSP Residual norm 1.998070253266e-06
143 KSP Residual norm 1.995347469782e-06
144 KSP Residual norm 1.989518461703e-06
145 KSP Residual norm 1.977152184939e-06
146 KSP Residual norm 1.956268435730e-06
147 KSP Residual norm 1.933741835071e-06
148 KSP Residual norm 1.903282487554e-06
149 KSP Residual norm 1.852992801238e-06
150 KSP Residual norm 1.805399988610e-06
151 KSP Residual norm 1.780994347161e-06
152 KSP Residual norm 1.758724640678e-06
153 KSP Residual norm 1.741785139789e-06
154 KSP Residual norm 1.729960287252e-06
155 KSP Residual norm 1.722639037882e-06
156 KSP Residual norm 1.717886600510e-06
157 KSP Residual norm 1.716690823426e-06
158 KSP Residual norm 1.716687901763e-06
159 KSP Residual norm 1.715728009211e-06
160 KSP Residual norm 1.711931771052e-06
161 KSP Residual norm 1.705427293500e-06
162 KSP Residual norm 1.694134402206e-06
163 KSP Residual norm 1.681925884623e-06
164 KSP Residual norm 1.669192420569e-06
165 KSP Residual norm 1.655373298279e-06
166 KSP Residual norm 1.643029516607e-06
167 KSP Residual norm 1.633824784916e-06
168 KSP Residual norm 1.625200980894e-06
169 KSP Residual norm 1.618026284854e-06
170 KSP Residual norm 1.614455149899e-06
171 KSP Residual norm 1.613267262776e-06
172 KSP Residual norm 1.613263468733e-06
173 KSP Residual norm 1.612824286477e-06
174 KSP Residual norm 1.611353951523e-06
175 KSP Residual norm 1.608719287421e-06
176 KSP Residual norm 1.604804091012e-06
177 KSP Residual norm 1.599631030973e-06
178 KSP Residual norm 1.593875036132e-06
179 KSP Residual norm 1.587155532692e-06
180 KSP Residual norm 1.579858618170e-06
181 KSP Residual norm 1.575395248383e-06
182 KSP Residual norm 1.569657005101e-06
183 KSP Residual norm 1.562953726886e-06
184 KSP Residual norm 1.558094493575e-06
185 KSP Residual norm 1.554659285749e-06
186 KSP Residual norm 1.552539529764e-06
187 KSP Residual norm 1.551536559097e-06
188 KSP Residual norm 1.551451922610e-06
189 KSP Residual norm 1.551156446567e-06
190 KSP Residual norm 1.549585509867e-06
191 KSP Residual norm 1.546047672843e-06
192 KSP Residual norm 1.540124267251e-06
193 KSP Residual norm 1.533684533355e-06
194 KSP Residual norm 1.525913861531e-06
195 KSP Residual norm 1.514318734030e-06
196 KSP Residual norm 1.496433144904e-06
197 KSP Residual norm 1.479608675862e-06
198 KSP Residual norm 1.465580795626e-06
199 KSP Residual norm 1.457272640901e-06
200 KSP Residual norm 1.454275659455e-06
KSP Object: 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=200
tolerances: relative=1e-12, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using PRECONDITIONED norm type for convergence test
PC Object: 8 MPI processes
type: gamg
MG: type is MULTIPLICATIVE, levels=4 cycles=v
Cycles per PCApply=1
Using Galerkin computed coarse grid matrices
GAMG specific options
Threshold for dropping small values from graph 0
AGG specific options
Symmetric graph false
Coarse grid solver -- level -------------------------------
KSP Object: (mg_coarse_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_) 8 MPI processes
type: bjacobi
block Jacobi: number of blocks = 8
Local solve is same for all blocks, in the following KSP and PC objects:
KSP Object: (mg_coarse_sub_) 1 MPI processes
type: preonly
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_sub_) 1 MPI processes
type: lu
LU: out-of-place factorization
tolerance for zero pivot 2.22045e-14
using diagonal shift on blocks to prevent zero pivot [INBLOCKS]
matrix ordering: nd
factor fill ratio given 5, needed 1
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
package used to perform factorization: petsc
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 8 MPI processes
type: mpiaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 8 nodes, limit used is 5
Down solver (pre-smoother) on level 1 -------------------------------
KSP Object: (mg_levels_1_) 8 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.164941, max = 1.81436
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_1_esteig_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_1_) 8 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 8 MPI processes
type: mpiaij
rows=672, cols=672, bs=6
total: nonzeros=255600, allocated nonzeros=255600
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 10 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 2 -------------------------------
KSP Object: (mg_levels_2_) 8 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.138297, max = 1.52126
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_2_esteig_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_2_) 8 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 8 MPI processes
type: mpiaij
rows=8940, cols=8940, bs=6
total: nonzeros=1.99829e+06, allocated nonzeros=1.99829e+06
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 364 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 3 -------------------------------
KSP Object: (mg_levels_3_) 8 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.1, max = 1.1
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_3_esteig_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_3_) 8 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: () 8 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 9163 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
linear system matrix = precond matrix:
Mat Object: () 8 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 9163 nodes, limit used is 5
Error, conv_flag < 0!
----------------------------------------------------------------------------------------------------------
8 MPI processes with "-ksp_type fgmres -pc_type gamg -ksp_max_it 200":
0 KSP Residual norm 1.000000000000e+00
1 KSP Residual norm 4.670997735323e-03
2 KSP Residual norm 4.107753808617e-03
3 KSP Residual norm 2.889337330348e-03
4 KSP Residual norm 2.116306152456e-03
5 KSP Residual norm 1.627169532601e-03
6 KSP Residual norm 1.218768613562e-03
7 KSP Residual norm 8.120536960099e-04
8 KSP Residual norm 6.255941040493e-04
9 KSP Residual norm 4.871101195846e-04
10 KSP Residual norm 3.715553641296e-04
11 KSP Residual norm 3.186530628626e-04
12 KSP Residual norm 2.916461487976e-04
13 KSP Residual norm 2.617076737616e-04
14 KSP Residual norm 2.434026947079e-04
15 KSP Residual norm 2.279904817094e-04
16 KSP Residual norm 2.094017928661e-04
17 KSP Residual norm 1.920994891443e-04
18 KSP Residual norm 1.824925891815e-04
19 KSP Residual norm 1.734956192133e-04
20 KSP Residual norm 1.641803577415e-04
21 KSP Residual norm 1.588154081656e-04
22 KSP Residual norm 1.551465552751e-04
23 KSP Residual norm 1.527361686091e-04
24 KSP Residual norm 1.507178096592e-04
25 KSP Residual norm 1.497023666544e-04
26 KSP Residual norm 1.493164041697e-04
27 KSP Residual norm 1.492075267955e-04
28 KSP Residual norm 1.492068301125e-04
29 KSP Residual norm 1.489442512775e-04
30 KSP Residual norm 1.481048645778e-04
31 KSP Residual norm 1.474912906430e-04
32 KSP Residual norm 1.472011893427e-04
33 KSP Residual norm 1.471630138354e-04
34 KSP Residual norm 1.471550027371e-04
35 KSP Residual norm 1.470274749944e-04
36 KSP Residual norm 1.466892693934e-04
37 KSP Residual norm 1.460358831019e-04
38 KSP Residual norm 1.446742442148e-04
39 KSP Residual norm 1.431445836508e-04
40 KSP Residual norm 1.412836143884e-04
41 KSP Residual norm 1.388559210259e-04
42 KSP Residual norm 1.362542327411e-04
43 KSP Residual norm 1.340755156773e-04
44 KSP Residual norm 1.320627198519e-04
45 KSP Residual norm 1.304003591661e-04
46 KSP Residual norm 1.296389081606e-04
47 KSP Residual norm 1.293033497828e-04
48 KSP Residual norm 1.291691278421e-04
49 KSP Residual norm 1.291623043886e-04
50 KSP Residual norm 1.291167471680e-04
51 KSP Residual norm 1.289206750819e-04
52 KSP Residual norm 1.285100433935e-04
53 KSP Residual norm 1.279213654135e-04
54 KSP Residual norm 1.267775022913e-04
55 KSP Residual norm 1.254444474248e-04
56 KSP Residual norm 1.232791465833e-04
57 KSP Residual norm 1.208718136472e-04
58 KSP Residual norm 1.188202291365e-04
59 KSP Residual norm 1.172567876791e-04
60 KSP Residual norm 1.161519332344e-04
61 KSP Residual norm 1.157179339808e-04
62 KSP Residual norm 1.154472502196e-04
63 KSP Residual norm 1.153127504460e-04
64 KSP Residual norm 1.152777949860e-04
65 KSP Residual norm 1.152770118788e-04
66 KSP Residual norm 1.152598598481e-04
67 KSP Residual norm 1.151629129170e-04
68 KSP Residual norm 1.149661894772e-04
69 KSP Residual norm 1.146467643611e-04
70 KSP Residual norm 1.141866143599e-04
71 KSP Residual norm 1.137025615373e-04
72 KSP Residual norm 1.131870901828e-04
73 KSP Residual norm 1.125975338408e-04
74 KSP Residual norm 1.121327693566e-04
75 KSP Residual norm 1.117787078570e-04
76 KSP Residual norm 1.115324968645e-04
77 KSP Residual norm 1.114305998751e-04
78 KSP Residual norm 1.113970100465e-04
79 KSP Residual norm 1.113957597622e-04
80 KSP Residual norm 1.113847867416e-04
81 KSP Residual norm 1.113271032325e-04
82 KSP Residual norm 1.112054965331e-04
83 KSP Residual norm 1.110712297565e-04
84 KSP Residual norm 1.108472210129e-04
85 KSP Residual norm 1.104980395243e-04
86 KSP Residual norm 1.100251900688e-04
87 KSP Residual norm 1.093770551681e-04
88 KSP Residual norm 1.087264856232e-04
89 KSP Residual norm 1.079864655032e-04
90 KSP Residual norm 1.069761018146e-04
91 KSP Residual norm 1.063544715487e-04
92 KSP Residual norm 1.059754177887e-04
93 KSP Residual norm 1.057041553935e-04
94 KSP Residual norm 1.055091860555e-04
95 KSP Residual norm 1.054083360521e-04
96 KSP Residual norm 1.053609855313e-04
97 KSP Residual norm 1.053537755585e-04
98 KSP Residual norm 1.053427425561e-04
99 KSP Residual norm 1.052632351147e-04
100 KSP Residual norm 1.050672404418e-04
101 KSP Residual norm 1.047076764607e-04
102 KSP Residual norm 1.042881690190e-04
103 KSP Residual norm 1.038070303216e-04
104 KSP Residual norm 1.032844486594e-04
105 KSP Residual norm 1.027143454602e-04
106 KSP Residual norm 1.020711306007e-04
107 KSP Residual norm 1.016161730578e-04
108 KSP Residual norm 1.012557310931e-04
109 KSP Residual norm 1.009893064643e-04
110 KSP Residual norm 1.008471358431e-04
111 KSP Residual norm 1.007996592947e-04
112 KSP Residual norm 1.007934604418e-04
113 KSP Residual norm 1.007885339167e-04
114 KSP Residual norm 1.007508141131e-04
115 KSP Residual norm 1.006505753610e-04
116 KSP Residual norm 1.004519339301e-04
117 KSP Residual norm 1.000280738585e-04
118 KSP Residual norm 9.944350404969e-05
119 KSP Residual norm 9.877759176521e-05
120 KSP Residual norm 9.818389395388e-05
121 KSP Residual norm 9.773577436930e-05
122 KSP Residual norm 9.728467175021e-05
123 KSP Residual norm 9.694926350524e-05
124 KSP Residual norm 9.672180646719e-05
125 KSP Residual norm 9.660114484763e-05
126 KSP Residual norm 9.654419729617e-05
127 KSP Residual norm 9.652671659533e-05
128 KSP Residual norm 9.652657055221e-05
129 KSP Residual norm 9.651524413669e-05
130 KSP Residual norm 9.647086564822e-05
131 KSP Residual norm 9.636753857657e-05
132 KSP Residual norm 9.618894386546e-05
133 KSP Residual norm 9.593619806440e-05
134 KSP Residual norm 9.569428638515e-05
135 KSP Residual norm 9.546096106785e-05
136 KSP Residual norm 9.525162010664e-05
137 KSP Residual norm 9.513447840595e-05
138 KSP Residual norm 9.505399070083e-05
139 KSP Residual norm 9.500332960990e-05
140 KSP Residual norm 9.497709591387e-05
141 KSP Residual norm 9.496352888228e-05
142 KSP Residual norm 9.495949083864e-05
143 KSP Residual norm 9.495948564734e-05
144 KSP Residual norm 9.495712057209e-05
145 KSP Residual norm 9.494127387935e-05
146 KSP Residual norm 9.489993555506e-05
147 KSP Residual norm 9.480961347934e-05
148 KSP Residual norm 9.469034802471e-05
149 KSP Residual norm 9.450844275660e-05
150 KSP Residual norm 9.417598101754e-05
151 KSP Residual norm 9.390512254846e-05
152 KSP Residual norm 9.376400510459e-05
153 KSP Residual norm 9.367492079299e-05
154 KSP Residual norm 9.361909990247e-05
155 KSP Residual norm 9.358335650237e-05
156 KSP Residual norm 9.356353195699e-05
157 KSP Residual norm 9.355891777491e-05
158 KSP Residual norm 9.355716324420e-05
159 KSP Residual norm 9.353449945065e-05
160 KSP Residual norm 9.346172941840e-05
161 KSP Residual norm 9.332283377696e-05
162 KSP Residual norm 9.315361593823e-05
163 KSP Residual norm 9.295337730885e-05
164 KSP Residual norm 9.270022939874e-05
165 KSP Residual norm 9.239090400786e-05
166 KSP Residual norm 9.198176650477e-05
167 KSP Residual norm 9.156845906463e-05
168 KSP Residual norm 9.110592161428e-05
169 KSP Residual norm 9.066318157752e-05
170 KSP Residual norm 9.045691253321e-05
171 KSP Residual norm 9.040420523474e-05
172 KSP Residual norm 9.040265102252e-05
173 KSP Residual norm 9.037594896598e-05
174 KSP Residual norm 9.027023066846e-05
175 KSP Residual norm 9.005602932590e-05
176 KSP Residual norm 8.974010911762e-05
177 KSP Residual norm 8.909184657206e-05
178 KSP Residual norm 8.820223874376e-05
179 KSP Residual norm 8.713852748360e-05
180 KSP Residual norm 8.640949295202e-05
181 KSP Residual norm 8.583932475736e-05
182 KSP Residual norm 8.510819872706e-05
183 KSP Residual norm 8.457549405131e-05
184 KSP Residual norm 8.420203270956e-05
185 KSP Residual norm 8.400296218178e-05
186 KSP Residual norm 8.389188464879e-05
187 KSP Residual norm 8.383544110933e-05
188 KSP Residual norm 8.382332577117e-05
189 KSP Residual norm 8.382244057144e-05
190 KSP Residual norm 8.379390563248e-05
191 KSP Residual norm 8.369470603151e-05
192 KSP Residual norm 8.349180033405e-05
193 KSP Residual norm 8.318078193378e-05
194 KSP Residual norm 8.288342685543e-05
195 KSP Residual norm 8.260634837066e-05
196 KSP Residual norm 8.240897737972e-05
197 KSP Residual norm 8.232131544913e-05
198 KSP Residual norm 8.226983935342e-05
199 KSP Residual norm 8.223729915616e-05
200 KSP Residual norm 8.222055281446e-05
KSP Object: 8 MPI processes
type: fgmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=200
tolerances: relative=1e-12, absolute=1e-50, divergence=10000
right preconditioning
using nonzero initial guess
using UNPRECONDITIONED norm type for convergence test
PC Object: 8 MPI processes
type: gamg
MG: type is MULTIPLICATIVE, levels=4 cycles=v
Cycles per PCApply=1
Using Galerkin computed coarse grid matrices
GAMG specific options
Threshold for dropping small values from graph 0
AGG specific options
Symmetric graph false
Coarse grid solver -- level -------------------------------
KSP Object: (mg_coarse_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_) 8 MPI processes
type: bjacobi
block Jacobi: number of blocks = 8
Local solve is same for all blocks, in the following KSP and PC objects:
KSP Object: (mg_coarse_sub_) 1 MPI processes
type: preonly
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_sub_) 1 MPI processes
type: lu
LU: out-of-place factorization
tolerance for zero pivot 2.22045e-14
using diagonal shift on blocks to prevent zero pivot [INBLOCKS]
matrix ordering: nd
factor fill ratio given 5, needed 1
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
package used to perform factorization: petsc
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 8 MPI processes
type: mpiaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 8 nodes, limit used is 5
Down solver (pre-smoother) on level 1 -------------------------------
KSP Object: (mg_levels_1_) 8 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.164941, max = 1.81436
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_1_esteig_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_1_) 8 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 8 MPI processes
type: mpiaij
rows=672, cols=672, bs=6
total: nonzeros=255600, allocated nonzeros=255600
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 10 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 2 -------------------------------
KSP Object: (mg_levels_2_) 8 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.138297, max = 1.52126
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_2_esteig_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_2_) 8 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 8 MPI processes
type: mpiaij
rows=8940, cols=8940, bs=6
total: nonzeros=1.99829e+06, allocated nonzeros=1.99829e+06
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 364 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 3 -------------------------------
KSP Object: (mg_levels_3_) 8 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.1, max = 1.1
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_3_esteig_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_3_) 8 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: () 8 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 9163 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
linear system matrix = precond matrix:
Mat Object: () 8 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 9163 nodes, limit used is 5
Error, conv_flag < 0!
----------------------------------------------------------------------------------------------------------
8 MPI processes with "-ksp_type cg -pc_type gamg":
0 KSP Residual norm 1.243143122667e-03
1 KSP Residual norm 7.792397481334e-04
KSP Object: 8 MPI processes
type: cg
maximum iterations=200
tolerances: relative=1e-12, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using PRECONDITIONED norm type for convergence test
PC Object: 8 MPI processes
type: gamg
MG: type is MULTIPLICATIVE, levels=4 cycles=v
Cycles per PCApply=1
Using Galerkin computed coarse grid matrices
GAMG specific options
Threshold for dropping small values from graph 0
AGG specific options
Symmetric graph false
Coarse grid solver -- level -------------------------------
KSP Object: (mg_coarse_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_) 8 MPI processes
type: bjacobi
block Jacobi: number of blocks = 8
Local solve is same for all blocks, in the following KSP and PC objects:
KSP Object: (mg_coarse_sub_) 1 MPI processes
type: preonly
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_sub_) 1 MPI processes
type: lu
LU: out-of-place factorization
tolerance for zero pivot 2.22045e-14
using diagonal shift on blocks to prevent zero pivot [INBLOCKS]
matrix ordering: nd
factor fill ratio given 5, needed 1
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
package used to perform factorization: petsc
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 8 MPI processes
type: mpiaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 8 nodes, limit used is 5
Down solver (pre-smoother) on level 1 -------------------------------
KSP Object: (mg_levels_1_) 8 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.164941, max = 1.81436
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_1_esteig_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_1_) 8 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 8 MPI processes
type: mpiaij
rows=672, cols=672, bs=6
total: nonzeros=255600, allocated nonzeros=255600
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 10 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 2 -------------------------------
KSP Object: (mg_levels_2_) 8 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.138297, max = 1.52126
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_2_esteig_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_2_) 8 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 8 MPI processes
type: mpiaij
rows=8940, cols=8940, bs=6
total: nonzeros=1.99829e+06, allocated nonzeros=1.99829e+06
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 364 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 3 -------------------------------
KSP Object: (mg_levels_3_) 8 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.1, max = 1.1
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_3_esteig_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_3_) 8 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: () 8 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 9163 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
linear system matrix = precond matrix:
Mat Object: () 8 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 9163 nodes, limit used is 5
Error, conv_flag < 0!
----------------------------------------------------------------------------------------------------------
8 MPI processes with "-ksp_type cg -pc_type gamg -ksp_norm_type natural":
0 KSP Residual norm 3.519959598679e-02
1 KSP Residual norm 6.406057376702e-04
KSP Object: 8 MPI processes
type: cg
maximum iterations=200
tolerances: relative=1e-12, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NATURAL norm type for convergence test
PC Object: 8 MPI processes
type: gamg
MG: type is MULTIPLICATIVE, levels=4 cycles=v
Cycles per PCApply=1
Using Galerkin computed coarse grid matrices
GAMG specific options
Threshold for dropping small values from graph 0
AGG specific options
Symmetric graph false
Coarse grid solver -- level -------------------------------
KSP Object: (mg_coarse_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_) 8 MPI processes
type: bjacobi
block Jacobi: number of blocks = 8
Local solve is same for all blocks, in the following KSP and PC objects:
KSP Object: (mg_coarse_sub_) 1 MPI processes
type: preonly
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_sub_) 1 MPI processes
type: lu
LU: out-of-place factorization
tolerance for zero pivot 2.22045e-14
using diagonal shift on blocks to prevent zero pivot [INBLOCKS]
matrix ordering: nd
factor fill ratio given 5, needed 1
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
package used to perform factorization: petsc
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 8 MPI processes
type: mpiaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 8 nodes, limit used is 5
Down solver (pre-smoother) on level 1 -------------------------------
KSP Object: (mg_levels_1_) 8 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.164941, max = 1.81436
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_1_esteig_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_1_) 8 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 8 MPI processes
type: mpiaij
rows=672, cols=672, bs=6
total: nonzeros=255600, allocated nonzeros=255600
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 10 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 2 -------------------------------
KSP Object: (mg_levels_2_) 8 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.138297, max = 1.52126
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_2_esteig_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_2_) 8 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 8 MPI processes
type: mpiaij
rows=8940, cols=8940, bs=6
total: nonzeros=1.99829e+06, allocated nonzeros=1.99829e+06
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 364 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 3 -------------------------------
KSP Object: (mg_levels_3_) 8 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.1, max = 1.1
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_3_esteig_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_3_) 8 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: () 8 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 9163 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
linear system matrix = precond matrix:
Mat Object: () 8 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 9163 nodes, limit used is 5
Error, conv_flag < 0!
----------------------------------------------------------------------------------------------------------
8 MPI processes with "-ksp_type cg -pc_type gamg -ksp_norm_type natural -pc_mg_type full":
0 KSP Residual norm 3.520284638135e-02
1 KSP Residual norm 1.694661905493e-02
KSP Object: 8 MPI processes
type: cg
maximum iterations=200
tolerances: relative=1e-12, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NATURAL norm type for convergence test
PC Object: 8 MPI processes
type: gamg
MG: type is FULL, levels=4 cycles=v
Using Galerkin computed coarse grid matrices
GAMG specific options
Threshold for dropping small values from graph 0
AGG specific options
Symmetric graph false
Coarse grid solver -- level -------------------------------
KSP Object: (mg_coarse_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_) 8 MPI processes
type: bjacobi
block Jacobi: number of blocks = 8
Local solve is same for all blocks, in the following KSP and PC objects:
KSP Object: (mg_coarse_sub_) 1 MPI processes
type: preonly
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_sub_) 1 MPI processes
type: lu
LU: out-of-place factorization
tolerance for zero pivot 2.22045e-14
using diagonal shift on blocks to prevent zero pivot [INBLOCKS]
matrix ordering: nd
factor fill ratio given 5, needed 1
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
package used to perform factorization: petsc
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 8 MPI processes
type: mpiaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 8 nodes, limit used is 5
Down solver (pre-smoother) on level 1 -------------------------------
KSP Object: (mg_levels_1_) 8 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.162676, max = 1.78943
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_1_esteig_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_1_) 8 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 8 MPI processes
type: mpiaij
rows=672, cols=672, bs=6
total: nonzeros=255600, allocated nonzeros=255600
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 10 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 2 -------------------------------
KSP Object: (mg_levels_2_) 8 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.138297, max = 1.52126
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_2_esteig_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_2_) 8 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 8 MPI processes
type: mpiaij
rows=8940, cols=8940, bs=6
total: nonzeros=1.99829e+06, allocated nonzeros=1.99829e+06
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 364 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 3 -------------------------------
KSP Object: (mg_levels_3_) 8 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.1, max = 1.1
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_3_esteig_) 8 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_3_) 8 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: () 8 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 9163 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
linear system matrix = precond matrix:
Mat Object: () 8 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 9163 nodes, limit used is 5
Error, conv_flag < 0!
----------------------------------------------------------------------------------------------------------
4 MPI processes with "-ksp_type cg -pc_type gamg -ksp_norm_type natural -pc_mg_type full":
0 KSP Residual norm 3.848292390006e-01
1 KSP Residual norm 3.394562853380e-02
2 KSP Residual norm 1.709369896514e-02
3 KSP Residual norm 6.553998141868e-03
4 KSP Residual norm 1.571187862685e-03
5 KSP Residual norm 5.542725674501e-04
6 KSP Residual norm 2.190719066108e-04
7 KSP Residual norm 8.746263629667e-05
8 KSP Residual norm 2.888658373789e-05
9 KSP Residual norm 1.070029893769e-05
10 KSP Residual norm 4.344983313819e-06
11 KSP Residual norm 1.580078249323e-06
12 KSP Residual norm 6.433222755357e-07
13 KSP Residual norm 3.060196103782e-07
14 KSP Residual norm 1.370666779262e-07
15 KSP Residual norm 5.152436428051e-08
16 KSP Residual norm 2.069305024182e-08
17 KSP Residual norm 8.324752474842e-09
18 KSP Residual norm 3.288645476001e-09
19 KSP Residual norm 1.304811823492e-09
20 KSP Residual norm 4.871732366599e-10
21 KSP Residual norm 2.006752778499e-10
22 KSP Residual norm 8.721524661501e-11
23 KSP Residual norm 3.839304119068e-11
24 KSP Residual norm 1.592613431818e-11
25 KSP Residual norm 5.759010833611e-12
26 KSP Residual norm 2.300499231327e-12
27 KSP Residual norm 9.338357091753e-13
28 KSP Residual norm 3.271574791340e-13
KSP Object: 4 MPI processes
type: cg
maximum iterations=200
tolerances: relative=1e-12, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NATURAL norm type for convergence test
PC Object: 4 MPI processes
type: gamg
MG: type is FULL, levels=4 cycles=v
Using Galerkin computed coarse grid matrices
GAMG specific options
Threshold for dropping small values from graph 0
AGG specific options
Symmetric graph false
Coarse grid solver -- level -------------------------------
KSP Object: (mg_coarse_) 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_) 4 MPI processes
type: bjacobi
block Jacobi: number of blocks = 4
Local solve is same for all blocks, in the following KSP and PC objects:
KSP Object: (mg_coarse_sub_) 1 MPI processes
type: preonly
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_sub_) 1 MPI processes
type: lu
LU: out-of-place factorization
tolerance for zero pivot 2.22045e-14
using diagonal shift on blocks to prevent zero pivot [INBLOCKS]
matrix ordering: nd
factor fill ratio given 5, needed 1
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
package used to perform factorization: petsc
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 4 MPI processes
type: mpiaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 8 nodes, limit used is 5
Down solver (pre-smoother) on level 1 -------------------------------
KSP Object: (mg_levels_1_) 4 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.15311, max = 1.6842
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_1_esteig_) 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_1_) 4 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 4 MPI processes
type: mpiaij
rows=690, cols=690, bs=6
total: nonzeros=271044, allocated nonzeros=271044
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 47 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 2 -------------------------------
KSP Object: (mg_levels_2_) 4 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.140891, max = 1.5498
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_2_esteig_) 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_2_) 4 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 4 MPI processes
type: mpiaij
rows=8892, cols=8892, bs=6
total: nonzeros=1.9759e+06, allocated nonzeros=1.9759e+06
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 742 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 3 -------------------------------
KSP Object: (mg_levels_3_) 4 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.152642, max = 1.67907
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_3_esteig_) 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_3_) 4 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: () 4 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 17386 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
linear system matrix = precond matrix:
Mat Object: () 4 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 17386 nodes, limit used is 5
----------------------------------------------------------------------------------------------------------
4 MPI processes with "-ksp_type cg -pc_type gamg":
0 KSP Residual norm 2.475992358942e-01
1 KSP Residual norm 1.010115669246e-01
2 KSP Residual norm 5.483162946279e-02
3 KSP Residual norm 3.603244951862e-02
4 KSP Residual norm 1.685453948015e-02
5 KSP Residual norm 6.509593131622e-03
6 KSP Residual norm 2.632449026853e-03
7 KSP Residual norm 1.241917834713e-03
8 KSP Residual norm 5.224606030158e-04
9 KSP Residual norm 2.415967606517e-04
10 KSP Residual norm 1.260567680797e-04
11 KSP Residual norm 6.128703145411e-05
12 KSP Residual norm 2.679437838882e-05
13 KSP Residual norm 1.365494948139e-05
14 KSP Residual norm 7.038019548649e-06
15 KSP Residual norm 3.080810878376e-06
16 KSP Residual norm 9.636609709684e-07
17 KSP Residual norm 4.213041474760e-07
18 KSP Residual norm 2.289578124495e-07
19 KSP Residual norm 1.263375775690e-07
20 KSP Residual norm 5.703870904870e-08
21 KSP Residual norm 2.485007763011e-08
22 KSP Residual norm 1.237773305875e-08
23 KSP Residual norm 5.385311109191e-09
24 KSP Residual norm 2.329545414721e-09
25 KSP Residual norm 8.675622114080e-10
26 KSP Residual norm 3.570856507066e-10
27 KSP Residual norm 1.734555265152e-10
28 KSP Residual norm 8.579357543855e-11
29 KSP Residual norm 3.633034228722e-11
30 KSP Residual norm 1.400577345736e-11
31 KSP Residual norm 6.244917840390e-12
32 KSP Residual norm 3.389807223488e-12
33 KSP Residual norm 1.737000616285e-12
34 KSP Residual norm 8.070277344018e-13
35 KSP Residual norm 3.849555117241e-13
36 KSP Residual norm 1.738497314772e-13
KSP Object: 4 MPI processes
type: cg
maximum iterations=200
tolerances: relative=1e-12, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using PRECONDITIONED norm type for convergence test
PC Object: 4 MPI processes
type: gamg
MG: type is MULTIPLICATIVE, levels=4 cycles=v
Cycles per PCApply=1
Using Galerkin computed coarse grid matrices
GAMG specific options
Threshold for dropping small values from graph 0
AGG specific options
Symmetric graph false
Coarse grid solver -- level -------------------------------
KSP Object: (mg_coarse_) 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_) 4 MPI processes
type: bjacobi
block Jacobi: number of blocks = 4
Local solve is same for all blocks, in the following KSP and PC objects:
KSP Object: (mg_coarse_sub_) 1 MPI processes
type: preonly
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_sub_) 1 MPI processes
type: lu
LU: out-of-place factorization
tolerance for zero pivot 2.22045e-14
using diagonal shift on blocks to prevent zero pivot [INBLOCKS]
matrix ordering: nd
factor fill ratio given 5, needed 1
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
package used to perform factorization: petsc
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 4 MPI processes
type: mpiaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 8 nodes, limit used is 5
Down solver (pre-smoother) on level 1 -------------------------------
KSP Object: (mg_levels_1_) 4 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.155999, max = 1.71599
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_1_esteig_) 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_1_) 4 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 4 MPI processes
type: mpiaij
rows=690, cols=690, bs=6
total: nonzeros=271044, allocated nonzeros=271044
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 47 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 2 -------------------------------
KSP Object: (mg_levels_2_) 4 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.140448, max = 1.54493
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_2_esteig_) 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_2_) 4 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 4 MPI processes
type: mpiaij
rows=8892, cols=8892, bs=6
total: nonzeros=1.9759e+06, allocated nonzeros=1.9759e+06
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 742 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 3 -------------------------------
KSP Object: (mg_levels_3_) 4 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.152642, max = 1.67907
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_3_esteig_) 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_3_) 4 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: () 4 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 17386 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
linear system matrix = precond matrix:
Mat Object: () 4 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 17386 nodes, limit used is 5
----------------------------------------------------------------------------------------------------------
4 MPI processes with "-ksp_type gmres -pc_type gamg":
0 KSP Residual norm 2.475992358942e-01
1 KSP Residual norm 9.788251866410e-02
2 KSP Residual norm 5.360475672741e-02
3 KSP Residual norm 3.205714918186e-02
4 KSP Residual norm 1.499738421848e-02
5 KSP Residual norm 5.656223448764e-03
6 KSP Residual norm 2.237512723165e-03
7 KSP Residual norm 1.072453752902e-03
8 KSP Residual norm 4.596044917514e-04
9 KSP Residual norm 2.133364000371e-04
10 KSP Residual norm 1.112111466773e-04
11 KSP Residual norm 5.012424812027e-05
12 KSP Residual norm 2.154120596339e-05
13 KSP Residual norm 1.113658201430e-05
14 KSP Residual norm 5.834290931458e-06
15 KSP Residual norm 2.584285286138e-06
16 KSP Residual norm 8.085531298248e-07
17 KSP Residual norm 3.730338029295e-07
18 KSP Residual norm 1.913619404427e-07
19 KSP Residual norm 1.034223785467e-07
20 KSP Residual norm 4.752281229701e-08
21 KSP Residual norm 2.194488512677e-08
22 KSP Residual norm 1.127906237520e-08
23 KSP Residual norm 4.848127907381e-09
24 KSP Residual norm 2.073720505286e-09
25 KSP Residual norm 7.549640791079e-10
26 KSP Residual norm 3.159362670485e-10
27 KSP Residual norm 1.480081790670e-10
28 KSP Residual norm 7.315280619947e-11
29 KSP Residual norm 3.096498631652e-11
30 KSP Residual norm 1.189973158428e-11
31 KSP Residual norm 6.479291407172e-12
32 KSP Residual norm 3.492865137025e-12
33 KSP Residual norm 1.805801366839e-12
34 KSP Residual norm 8.167217552481e-13
35 KSP Residual norm 3.806289288690e-13
36 KSP Residual norm 1.782376313089e-13
KSP Object: 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=200
tolerances: relative=1e-12, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using PRECONDITIONED norm type for convergence test
PC Object: 4 MPI processes
type: gamg
MG: type is MULTIPLICATIVE, levels=4 cycles=v
Cycles per PCApply=1
Using Galerkin computed coarse grid matrices
GAMG specific options
Threshold for dropping small values from graph 0
AGG specific options
Symmetric graph false
Coarse grid solver -- level -------------------------------
KSP Object: (mg_coarse_) 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_) 4 MPI processes
type: bjacobi
block Jacobi: number of blocks = 4
Local solve is same for all blocks, in the following KSP and PC objects:
KSP Object: (mg_coarse_sub_) 1 MPI processes
type: preonly
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_sub_) 1 MPI processes
type: lu
LU: out-of-place factorization
tolerance for zero pivot 2.22045e-14
using diagonal shift on blocks to prevent zero pivot [INBLOCKS]
matrix ordering: nd
factor fill ratio given 5, needed 1
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
package used to perform factorization: petsc
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 4 MPI processes
type: mpiaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 8 nodes, limit used is 5
Down solver (pre-smoother) on level 1 -------------------------------
KSP Object: (mg_levels_1_) 4 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.155999, max = 1.71599
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_1_esteig_) 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_1_) 4 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 4 MPI processes
type: mpiaij
rows=690, cols=690, bs=6
total: nonzeros=271044, allocated nonzeros=271044
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 47 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 2 -------------------------------
KSP Object: (mg_levels_2_) 4 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.140448, max = 1.54493
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_2_esteig_) 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_2_) 4 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 4 MPI processes
type: mpiaij
rows=8892, cols=8892, bs=6
total: nonzeros=1.9759e+06, allocated nonzeros=1.9759e+06
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 742 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 3 -------------------------------
KSP Object: (mg_levels_3_) 4 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.152642, max = 1.67907
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_3_esteig_) 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_3_) 4 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: () 4 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 17386 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
linear system matrix = precond matrix:
Mat Object: () 4 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 17386 nodes, limit used is 5
----------------------------------------------------------------------------------------------------------
4 MPI processes with "-ksp_type fgmres -pc_type gamg":
0 KSP Residual norm 1.000000000000e+00
1 KSP Residual norm 4.455458968821e-02
2 KSP Residual norm 1.488836455945e-02
3 KSP Residual norm 9.951640321741e-03
4 KSP Residual norm 4.935202239340e-03
5 KSP Residual norm 2.196694493952e-03
6 KSP Residual norm 8.473972720318e-04
7 KSP Residual norm 3.682116362426e-04
8 KSP Residual norm 1.810781420824e-04
9 KSP Residual norm 8.751313625489e-05
10 KSP Residual norm 4.184578916542e-05
11 KSP Residual norm 2.187519733560e-05
12 KSP Residual norm 9.952092586537e-06
13 KSP Residual norm 4.670524263053e-06
14 KSP Residual norm 2.273179771211e-06
15 KSP Residual norm 1.077537786223e-06
16 KSP Residual norm 5.567436295073e-07
17 KSP Residual norm 2.877739699547e-07
18 KSP Residual norm 1.308592111530e-07
19 KSP Residual norm 5.712099664096e-08
20 KSP Residual norm 2.699430333773e-08
21 KSP Residual norm 1.305214840427e-08
22 KSP Residual norm 5.720307231306e-09
23 KSP Residual norm 2.481505127019e-09
24 KSP Residual norm 1.111010201625e-09
25 KSP Residual norm 4.839255906413e-10
26 KSP Residual norm 2.109362228309e-10
27 KSP Residual norm 9.986868872181e-11
28 KSP Residual norm 4.448162227505e-11
29 KSP Residual norm 2.209195701027e-11
30 KSP Residual norm 9.636290761114e-12
31 KSP Residual norm 5.491604445745e-12
32 KSP Residual norm 2.816561935838e-12
33 KSP Residual norm 1.430060877293e-12
34 KSP Residual norm 7.093552021344e-13
KSP Object: 4 MPI processes
type: fgmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=200
tolerances: relative=1e-12, absolute=1e-50, divergence=10000
right preconditioning
using nonzero initial guess
using UNPRECONDITIONED norm type for convergence test
PC Object: 4 MPI processes
type: gamg
MG: type is MULTIPLICATIVE, levels=4 cycles=v
Cycles per PCApply=1
Using Galerkin computed coarse grid matrices
GAMG specific options
Threshold for dropping small values from graph 0
AGG specific options
Symmetric graph false
Coarse grid solver -- level -------------------------------
KSP Object: (mg_coarse_) 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_) 4 MPI processes
type: bjacobi
block Jacobi: number of blocks = 4
Local solve is same for all blocks, in the following KSP and PC objects:
KSP Object: (mg_coarse_sub_) 1 MPI processes
type: preonly
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (mg_coarse_sub_) 1 MPI processes
type: lu
LU: out-of-place factorization
tolerance for zero pivot 2.22045e-14
using diagonal shift on blocks to prevent zero pivot [INBLOCKS]
matrix ordering: nd
factor fill ratio given 5, needed 1
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
package used to perform factorization: petsc
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 8 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 4 MPI processes
type: mpiaij
rows=36, cols=36, bs=6
total: nonzeros=1296, allocated nonzeros=1296
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 8 nodes, limit used is 5
Down solver (pre-smoother) on level 1 -------------------------------
KSP Object: (mg_levels_1_) 4 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.155999, max = 1.71599
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_1_esteig_) 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_1_) 4 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 4 MPI processes
type: mpiaij
rows=690, cols=690, bs=6
total: nonzeros=271044, allocated nonzeros=271044
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 47 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 2 -------------------------------
KSP Object: (mg_levels_2_) 4 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.140448, max = 1.54493
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_2_esteig_) 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_2_) 4 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 4 MPI processes
type: mpiaij
rows=8892, cols=8892, bs=6
total: nonzeros=1.9759e+06, allocated nonzeros=1.9759e+06
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 742 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 3 -------------------------------
KSP Object: (mg_levels_3_) 4 MPI processes
type: chebyshev
Chebyshev: eigenvalue estimates: min = 0.152642, max = 1.67907
Chebyshev: eigenvalues estimated using gmres with translations [0 0.1; 0 1.1]
KSP Object: (mg_levels_3_esteig_) 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
maximum iterations=2
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (mg_levels_3_) 4 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: () 4 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 17386 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
linear system matrix = precond matrix:
Mat Object: () 4 MPI processes
type: mpiaij
rows=199752, cols=199752, bs=3
total: nonzeros=1.54794e+07, allocated nonzeros=1.54794e+07
total number of mallocs used during MatSetValues calls =0
has attached near null space
using I-node (on process 0) routines: found 17386 nodes, limit used is 5
----------------------------------------------------------------------------------------------------------
More information about the petsc-users
mailing list