xref: /petsc/src/snes/tutorials/output/ex19_tut_2.out (revision 70646cd191a02c3aba559ba717dac5da7a8a1e20)
1c4762a1bSJed Brownlid velocity = 0.000106281, prandtl # = 1., grashof # = 1.
2c4762a1bSJed Brown  0 SNES Function norm 1.036007954337e-02
3c4762a1bSJed Brown    0 KSP Residual norm 2.388589583549e+00
4c4762a1bSJed Brown    1 KSP Residual norm 5.715829806981e-01
5c4762a1bSJed Brown    2 KSP Residual norm 4.623679005936e-02
6c4762a1bSJed Brown    3 KSP Residual norm 1.143381177646e-02
7c4762a1bSJed Brown    4 KSP Residual norm 2.015139840224e-03
8c4762a1bSJed Brown    5 KSP Residual norm 4.356196119798e-04
9c4762a1bSJed Brown    6 KSP Residual norm 4.240953066710e-05
10c4762a1bSJed Brown    7 KSP Residual norm 8.848315297175e-06
11c4762a1bSJed Brown  1 SNES Function norm 9.854304971115e-06
12c4762a1bSJed Brown    0 KSP Residual norm 3.868049496775e-05
13c4762a1bSJed Brown    1 KSP Residual norm 7.693574326868e-06
14c4762a1bSJed Brown    2 KSP Residual norm 1.059429116239e-06
15c4762a1bSJed Brown    3 KSP Residual norm 4.004524784804e-07
16c4762a1bSJed Brown    4 KSP Residual norm 1.050186948327e-07
17c4762a1bSJed Brown    5 KSP Residual norm 5.073180513583e-08
18c4762a1bSJed Brown    6 KSP Residual norm 2.510513776297e-08
19c4762a1bSJed Brown    7 KSP Residual norm 1.211886495400e-08
20c4762a1bSJed Brown    8 KSP Residual norm 1.911963112131e-09
21c4762a1bSJed Brown    9 KSP Residual norm 3.005260864225e-10
22c4762a1bSJed Brown  2 SNES Function norm 3.117674497824e-10
23c4762a1bSJed Brown    0 KSP Residual norm 3.005042584730e-10
24c4762a1bSJed Brown    1 KSP Residual norm 1.120673922713e-10
25c4762a1bSJed Brown    2 KSP Residual norm 3.288439453292e-11
26c4762a1bSJed Brown    3 KSP Residual norm 5.822504321413e-12
27c4762a1bSJed Brown    4 KSP Residual norm 2.486684466178e-12
28c4762a1bSJed Brown    5 KSP Residual norm 1.198858055503e-12
29c4762a1bSJed Brown    6 KSP Residual norm 6.255669709502e-13
30c4762a1bSJed Brown    7 KSP Residual norm 1.544647758005e-13
31c4762a1bSJed Brown    8 KSP Residual norm 4.592122224907e-14
32c4762a1bSJed Brown    9 KSP Residual norm 4.984149547392e-15
33c4762a1bSJed Brown   10 KSP Residual norm 8.905129652955e-16
34c4762a1bSJed Brown  3 SNES Function norm 1.045594761851e-14
35c4762a1bSJed BrownSNES Object: 4 MPI processes
36c4762a1bSJed Brown  type: newtonls
37c4762a1bSJed Brown  maximum iterations=50, maximum function evaluations=10000
38c4762a1bSJed Brown  tolerances: relative=1e-08, absolute=1e-50, solution=1e-08
39c4762a1bSJed Brown  total number of linear solver iterations=26
40c4762a1bSJed Brown  total number of function evaluations=4
41c4762a1bSJed Brown  norm schedule ALWAYS
42c4762a1bSJed Brown  Jacobian is built using colored finite differences on a DMDA
43c4762a1bSJed Brown  SNESLineSearch Object: 4 MPI processes
44c4762a1bSJed Brown    type: bt
45c4762a1bSJed Brown      interpolation: cubic
46c4762a1bSJed Brown      alpha=1.000000e-04
47a99ef635SJonas Heinzmann    maxlambda=1.000000e+00, minlambda=1.000000e-12
48c4762a1bSJed Brown    tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
49c4762a1bSJed Brown    maximum iterations=40
50c4762a1bSJed Brown  KSP Object: 4 MPI processes
51c4762a1bSJed Brown    type: gmres
52f971d498SPierre Jolivet      restart=30, using classical (unmodified) Gram-Schmidt orthogonalization with no iterative refinement
53*143f2514SPierre Jolivet      happy breakdown tolerance=1e-30
54c4762a1bSJed Brown    maximum iterations=10000, initial guess is zero
55c4762a1bSJed Brown    tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
56c4762a1bSJed Brown    left preconditioning
57c4762a1bSJed Brown    using PRECONDITIONED norm type for convergence test
58c4762a1bSJed Brown  PC Object: 4 MPI processes
59c4762a1bSJed Brown    type: mg
60c4762a1bSJed Brown      type is MULTIPLICATIVE, levels=6 cycles=v
61c4762a1bSJed Brown        Cycles per PCApply=1
62c4762a1bSJed Brown        Not using Galerkin computed coarse grid matrices
6363a3b9bcSJacob Faibussowitsch    Coarse grid solver -- level 0 -------------------------------
64c4762a1bSJed Brown      KSP Object: (mg_coarse_) 4 MPI processes
65c4762a1bSJed Brown        type: preonly
66c4762a1bSJed Brown        maximum iterations=10000, initial guess is zero
67c4762a1bSJed Brown        tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
68c4762a1bSJed Brown        left preconditioning
698405740aSBarry Smith        not checking for convergence
70c4762a1bSJed Brown      PC Object: (mg_coarse_) 4 MPI processes
71c4762a1bSJed Brown        type: redundant
72c4762a1bSJed Brown          First (color=0) of 4 PCs follows
738cc725e6SPierre Jolivet          KSP Object: (mg_coarse_redundant_) 1 MPI process
74c4762a1bSJed Brown            type: preonly
75c4762a1bSJed Brown            maximum iterations=10000, initial guess is zero
76c4762a1bSJed Brown            tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
77c4762a1bSJed Brown            left preconditioning
788405740aSBarry Smith            not checking for convergence
798cc725e6SPierre Jolivet          PC Object: (mg_coarse_redundant_) 1 MPI process
80c4762a1bSJed Brown            type: lu
81c4762a1bSJed Brown              out-of-place factorization
82c4762a1bSJed Brown              tolerance for zero pivot 2.22045e-14
83c4762a1bSJed Brown              using diagonal shift on blocks to prevent zero pivot [INBLOCKS]
84c4762a1bSJed Brown              matrix ordering: nd
85c4762a1bSJed Brown              factor fill ratio given 5., needed 1.875
86ecf3d421SBarry Smith                Factored matrix:
8726cc229bSBarry Smith                  Mat Object: (mg_coarse_redundant_) 1 MPI process
88c4762a1bSJed Brown                    type: seqaij
89c4762a1bSJed Brown                    rows=64, cols=64, bs=4
90c4762a1bSJed Brown                    package used to perform factorization: petsc
91c4762a1bSJed Brown                    total: nonzeros=1920, allocated nonzeros=1920
92c4762a1bSJed Brown                      using I-node routines: found 16 nodes, limit used is 5
93ecf3d421SBarry Smith            linear system matrix, which is also used to construct the preconditioner:
9426cc229bSBarry Smith            Mat Object: 1 MPI process
95c4762a1bSJed Brown              type: seqaij
96c4762a1bSJed Brown              rows=64, cols=64, bs=4
97c4762a1bSJed Brown              total: nonzeros=1024, allocated nonzeros=1024
98c4762a1bSJed Brown              total number of mallocs used during MatSetValues calls=0
99c4762a1bSJed Brown                using I-node routines: found 16 nodes, limit used is 5
100ecf3d421SBarry Smith        linear system matrix, which is also used to construct the preconditioner:
101c4762a1bSJed Brown        Mat Object: 4 MPI processes
102c4762a1bSJed Brown          type: mpiaij
103c4762a1bSJed Brown          rows=64, cols=64, bs=4
104c4762a1bSJed Brown          total: nonzeros=1024, allocated nonzeros=1024
105c4762a1bSJed Brown          total number of mallocs used during MatSetValues calls=0
106c4762a1bSJed Brown    Down solver (pre-smoother) on level 1 -------------------------------
107c4762a1bSJed Brown      KSP Object: (mg_levels_1_) 4 MPI processes
108c4762a1bSJed Brown        type: chebyshev
109f2edd1f0SMalachi Phillips          Chebyshev polynomial of first kind
11073f7197eSJed Brown          eigenvalue targets used: min 0.148269, max 1.63095
11173f7197eSJed Brown          eigenvalues estimated via gmres: min 0.144902, max 1.48269
11273f7197eSJed Brown          eigenvalues estimated using gmres with transform: [0. 0.1; 0. 1.1]
113c4762a1bSJed Brown          KSP Object: (mg_levels_1_esteig_) 4 MPI processes
114c4762a1bSJed Brown            type: gmres
115f971d498SPierre Jolivet              restart=30, using classical (unmodified) Gram-Schmidt orthogonalization with no iterative refinement
116*143f2514SPierre Jolivet              happy breakdown tolerance=1e-30
117c4762a1bSJed Brown            maximum iterations=10, initial guess is zero
118c4762a1bSJed Brown            tolerances: relative=1e-12, absolute=1e-50, divergence=10000.
119c4762a1bSJed Brown            left preconditioning
120c4762a1bSJed Brown            using PRECONDITIONED norm type for convergence test
121dd8e379bSPierre Jolivet          estimating eigenvalues using a noisy random number generated right-hand side
122c4762a1bSJed Brown        maximum iterations=2, nonzero initial guess
123c4762a1bSJed Brown        tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
124c4762a1bSJed Brown        left preconditioning
1258405740aSBarry Smith        not checking for convergence
126c4762a1bSJed Brown      PC Object: (mg_levels_1_) 4 MPI processes
127c4762a1bSJed Brown        type: sor
128c4762a1bSJed Brown          type = local_symmetric, iterations = 1, local iterations = 1, omega = 1.
129ecf3d421SBarry Smith        linear system matrix, which is also used to construct the preconditioner:
130c4762a1bSJed Brown        Mat Object: 4 MPI processes
131c4762a1bSJed Brown          type: mpiaij
132c4762a1bSJed Brown          rows=196, cols=196, bs=4
133c4762a1bSJed Brown          total: nonzeros=3472, allocated nonzeros=3472
134c4762a1bSJed Brown          total number of mallocs used during MatSetValues calls=0
135c4762a1bSJed Brown    Up solver (post-smoother) same as down solver (pre-smoother)
136c4762a1bSJed Brown    Down solver (pre-smoother) on level 2 -------------------------------
137c4762a1bSJed Brown      KSP Object: (mg_levels_2_) 4 MPI processes
138c4762a1bSJed Brown        type: chebyshev
139f2edd1f0SMalachi Phillips          Chebyshev polynomial of first kind
14073f7197eSJed Brown          eigenvalue targets used: min 0.149178, max 1.64096
14173f7197eSJed Brown          eigenvalues estimated via gmres: min 0.0843938, max 1.49178
14273f7197eSJed Brown          eigenvalues estimated using gmres with transform: [0. 0.1; 0. 1.1]
143c4762a1bSJed Brown          KSP Object: (mg_levels_2_esteig_) 4 MPI processes
144c4762a1bSJed Brown            type: gmres
145f971d498SPierre Jolivet              restart=30, using classical (unmodified) Gram-Schmidt orthogonalization with no iterative refinement
146*143f2514SPierre Jolivet              happy breakdown tolerance=1e-30
147c4762a1bSJed Brown            maximum iterations=10, initial guess is zero
148c4762a1bSJed Brown            tolerances: relative=1e-12, absolute=1e-50, divergence=10000.
149c4762a1bSJed Brown            left preconditioning
150c4762a1bSJed Brown            using PRECONDITIONED norm type for convergence test
151dd8e379bSPierre Jolivet          estimating eigenvalues using a noisy random number generated right-hand side
152c4762a1bSJed Brown        maximum iterations=2, nonzero initial guess
153c4762a1bSJed Brown        tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
154c4762a1bSJed Brown        left preconditioning
1558405740aSBarry Smith        not checking for convergence
156c4762a1bSJed Brown      PC Object: (mg_levels_2_) 4 MPI processes
157c4762a1bSJed Brown        type: sor
158c4762a1bSJed Brown          type = local_symmetric, iterations = 1, local iterations = 1, omega = 1.
159ecf3d421SBarry Smith        linear system matrix, which is also used to construct the preconditioner:
160c4762a1bSJed Brown        Mat Object: 4 MPI processes
161c4762a1bSJed Brown          type: mpiaij
162c4762a1bSJed Brown          rows=676, cols=676, bs=4
163c4762a1bSJed Brown          total: nonzeros=12688, allocated nonzeros=12688
164c4762a1bSJed Brown          total number of mallocs used during MatSetValues calls=0
165c4762a1bSJed Brown    Up solver (post-smoother) same as down solver (pre-smoother)
166c4762a1bSJed Brown    Down solver (pre-smoother) on level 3 -------------------------------
167c4762a1bSJed Brown      KSP Object: (mg_levels_3_) 4 MPI processes
168c4762a1bSJed Brown        type: chebyshev
169f2edd1f0SMalachi Phillips          Chebyshev polynomial of first kind
17073f7197eSJed Brown          eigenvalue targets used: min 0.146454, max 1.61099
17173f7197eSJed Brown          eigenvalues estimated via gmres: min 0.0659153, max 1.46454
17273f7197eSJed Brown          eigenvalues estimated using gmres with transform: [0. 0.1; 0. 1.1]
173c4762a1bSJed Brown          KSP Object: (mg_levels_3_esteig_) 4 MPI processes
174c4762a1bSJed Brown            type: gmres
175f971d498SPierre Jolivet              restart=30, using classical (unmodified) Gram-Schmidt orthogonalization with no iterative refinement
176*143f2514SPierre Jolivet              happy breakdown tolerance=1e-30
177c4762a1bSJed Brown            maximum iterations=10, initial guess is zero
178c4762a1bSJed Brown            tolerances: relative=1e-12, absolute=1e-50, divergence=10000.
179c4762a1bSJed Brown            left preconditioning
180c4762a1bSJed Brown            using PRECONDITIONED norm type for convergence test
181dd8e379bSPierre Jolivet          estimating eigenvalues using a noisy random number generated right-hand side
182c4762a1bSJed Brown        maximum iterations=2, nonzero initial guess
183c4762a1bSJed Brown        tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
184c4762a1bSJed Brown        left preconditioning
1858405740aSBarry Smith        not checking for convergence
186c4762a1bSJed Brown      PC Object: (mg_levels_3_) 4 MPI processes
187c4762a1bSJed Brown        type: sor
188c4762a1bSJed Brown          type = local_symmetric, iterations = 1, local iterations = 1, omega = 1.
189ecf3d421SBarry Smith        linear system matrix, which is also used to construct the preconditioner:
190c4762a1bSJed Brown        Mat Object: 4 MPI processes
191c4762a1bSJed Brown          type: mpiaij
192c4762a1bSJed Brown          rows=2500, cols=2500, bs=4
193c4762a1bSJed Brown          total: nonzeros=48400, allocated nonzeros=48400
194c4762a1bSJed Brown          total number of mallocs used during MatSetValues calls=0
195c4762a1bSJed Brown    Up solver (post-smoother) same as down solver (pre-smoother)
196c4762a1bSJed Brown    Down solver (pre-smoother) on level 4 -------------------------------
197c4762a1bSJed Brown      KSP Object: (mg_levels_4_) 4 MPI processes
198c4762a1bSJed Brown        type: chebyshev
199f2edd1f0SMalachi Phillips          Chebyshev polynomial of first kind
20073f7197eSJed Brown          eigenvalue targets used: min 0.141089, max 1.55197
20173f7197eSJed Brown          eigenvalues estimated via gmres: min 0.044097, max 1.41089
20273f7197eSJed Brown          eigenvalues estimated using gmres with transform: [0. 0.1; 0. 1.1]
203c4762a1bSJed Brown          KSP Object: (mg_levels_4_esteig_) 4 MPI processes
204c4762a1bSJed Brown            type: gmres
205f971d498SPierre Jolivet              restart=30, using classical (unmodified) Gram-Schmidt orthogonalization with no iterative refinement
206*143f2514SPierre Jolivet              happy breakdown tolerance=1e-30
207c4762a1bSJed Brown            maximum iterations=10, initial guess is zero
208c4762a1bSJed Brown            tolerances: relative=1e-12, absolute=1e-50, divergence=10000.
209c4762a1bSJed Brown            left preconditioning
210c4762a1bSJed Brown            using PRECONDITIONED norm type for convergence test
211dd8e379bSPierre Jolivet          estimating eigenvalues using a noisy random number generated right-hand side
212c4762a1bSJed Brown        maximum iterations=2, nonzero initial guess
213c4762a1bSJed Brown        tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
214c4762a1bSJed Brown        left preconditioning
2158405740aSBarry Smith        not checking for convergence
216c4762a1bSJed Brown      PC Object: (mg_levels_4_) 4 MPI processes
217c4762a1bSJed Brown        type: sor
218c4762a1bSJed Brown          type = local_symmetric, iterations = 1, local iterations = 1, omega = 1.
219ecf3d421SBarry Smith        linear system matrix, which is also used to construct the preconditioner:
220c4762a1bSJed Brown        Mat Object: 4 MPI processes
221c4762a1bSJed Brown          type: mpiaij
222c4762a1bSJed Brown          rows=9604, cols=9604, bs=4
223c4762a1bSJed Brown          total: nonzeros=188944, allocated nonzeros=188944
224c4762a1bSJed Brown          total number of mallocs used during MatSetValues calls=0
225c4762a1bSJed Brown    Up solver (post-smoother) same as down solver (pre-smoother)
226c4762a1bSJed Brown    Down solver (pre-smoother) on level 5 -------------------------------
227c4762a1bSJed Brown      KSP Object: (mg_levels_5_) 4 MPI processes
228c4762a1bSJed Brown        type: chebyshev
229f2edd1f0SMalachi Phillips          Chebyshev polynomial of first kind
23073f7197eSJed Brown          eigenvalue targets used: min 0.127956, max 1.40751
23173f7197eSJed Brown          eigenvalues estimated via gmres: min 0.0380398, max 1.27956
23273f7197eSJed Brown          eigenvalues estimated using gmres with transform: [0. 0.1; 0. 1.1]
233c4762a1bSJed Brown          KSP Object: (mg_levels_5_esteig_) 4 MPI processes
234c4762a1bSJed Brown            type: gmres
235f971d498SPierre Jolivet              restart=30, using classical (unmodified) Gram-Schmidt orthogonalization with no iterative refinement
236*143f2514SPierre Jolivet              happy breakdown tolerance=1e-30
237c4762a1bSJed Brown            maximum iterations=10, initial guess is zero
238c4762a1bSJed Brown            tolerances: relative=1e-12, absolute=1e-50, divergence=10000.
239c4762a1bSJed Brown            left preconditioning
240c4762a1bSJed Brown            using PRECONDITIONED norm type for convergence test
241dd8e379bSPierre Jolivet          estimating eigenvalues using a noisy random number generated right-hand side
242c4762a1bSJed Brown        maximum iterations=2, nonzero initial guess
243c4762a1bSJed Brown        tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
244c4762a1bSJed Brown        left preconditioning
2458405740aSBarry Smith        not checking for convergence
246c4762a1bSJed Brown      PC Object: (mg_levels_5_) 4 MPI processes
247c4762a1bSJed Brown        type: sor
248c4762a1bSJed Brown          type = local_symmetric, iterations = 1, local iterations = 1, omega = 1.
249ecf3d421SBarry Smith        linear system matrix, which is also used to construct the preconditioner:
250c4762a1bSJed Brown        Mat Object: 4 MPI processes
251c4762a1bSJed Brown          type: mpiaij
252c4762a1bSJed Brown          rows=37636, cols=37636, bs=4
253c4762a1bSJed Brown          total: nonzeros=746512, allocated nonzeros=746512
254c4762a1bSJed Brown          total number of mallocs used during MatSetValues calls=0
255c4762a1bSJed Brown    Up solver (post-smoother) same as down solver (pre-smoother)
256ecf3d421SBarry Smith    linear system matrix, which is also used to construct the preconditioner:
257c4762a1bSJed Brown    Mat Object: 4 MPI processes
258c4762a1bSJed Brown      type: mpiaij
259c4762a1bSJed Brown      rows=37636, cols=37636, bs=4
260c4762a1bSJed Brown      total: nonzeros=746512, allocated nonzeros=746512
261c4762a1bSJed Brown      total number of mallocs used during MatSetValues calls=0
262c4762a1bSJed BrownNumber of SNES iterations = 3
263