xref: /petsc/src/snes/tutorials/output/ex12_p4est_solve_bddc.out (revision eaa7cf455101e7391ead982166c6e99faf273403)
1c4762a1bSJed Brown  0 SNES Function norm 17621.4
2f9635d15SStefano Zampini    0 KSP Residual norm 6.602334536848e+00
3f9635d15SStefano Zampini    1 KSP Residual norm 1.629553767711e-01
4f9635d15SStefano Zampini    2 KSP Residual norm 9.722599947211e-03
5ac7f1a8bSStefano Zampini    3 KSP Residual norm 7.839777888694e-04
6ac7f1a8bSStefano Zampini    4 KSP Residual norm 4.516536402495e-05
7ac7f1a8bSStefano Zampini    5 KSP Residual norm 1.551480410721e-06
8c4762a1bSJed Brown  1 SNES Function norm 5222.44
9f9635d15SStefano Zampini    0 KSP Residual norm 3.864965825919e+00
10f9635d15SStefano Zampini    1 KSP Residual norm 1.038768916208e+00
11f9635d15SStefano Zampini    2 KSP Residual norm 2.533940522317e-01
12ac7f1a8bSStefano Zampini    3 KSP Residual norm 1.653538770638e-02
13f9635d15SStefano Zampini    4 KSP Residual norm 8.631538715215e-04
14ac7f1a8bSStefano Zampini    5 KSP Residual norm 1.065288757146e-04
15ac7f1a8bSStefano Zampini    6 KSP Residual norm 7.258558309150e-06
16ac7f1a8bSStefano Zampini    7 KSP Residual norm 3.349242326554e-07
17c4762a1bSJed Brown  2 SNES Function norm 1549.16
18f9635d15SStefano Zampini    0 KSP Residual norm 3.882775661733e+00
19f9635d15SStefano Zampini    1 KSP Residual norm 2.427808286767e-01
20f9635d15SStefano Zampini    2 KSP Residual norm 4.613930408749e-02
21ac7f1a8bSStefano Zampini    3 KSP Residual norm 2.754707602635e-03
22ac7f1a8bSStefano Zampini    4 KSP Residual norm 3.993273630331e-04
23ac7f1a8bSStefano Zampini    5 KSP Residual norm 5.937195329812e-05
24ac7f1a8bSStefano Zampini    6 KSP Residual norm 2.096865040206e-06
25c4762a1bSJed Brown  3 SNES Function norm 459.505
26f9635d15SStefano Zampini    0 KSP Residual norm 2.249583536806e+00
27f9635d15SStefano Zampini    1 KSP Residual norm 1.800531165705e-01
28ac7f1a8bSStefano Zampini    2 KSP Residual norm 2.510122315296e-02
29f9635d15SStefano Zampini    3 KSP Residual norm 1.772932243567e-03
30ac7f1a8bSStefano Zampini    4 KSP Residual norm 1.609940560635e-04
31ac7f1a8bSStefano Zampini    5 KSP Residual norm 1.477658710616e-05
32ac7f1a8bSStefano Zampini    6 KSP Residual norm 6.537408174152e-07
33c4762a1bSJed Brown  4 SNES Function norm 136.479
34f9635d15SStefano Zampini    0 KSP Residual norm 9.548307303131e-01
35ac7f1a8bSStefano Zampini    1 KSP Residual norm 7.841005787762e-02
36f9635d15SStefano Zampini    2 KSP Residual norm 1.051451448518e-02
37ac7f1a8bSStefano Zampini    3 KSP Residual norm 1.022681088477e-03
38ac7f1a8bSStefano Zampini    4 KSP Residual norm 9.681331572457e-05
39ac7f1a8bSStefano Zampini    5 KSP Residual norm 7.875134849433e-06
40ac7f1a8bSStefano Zampini    6 KSP Residual norm 3.526178148335e-07
41c4762a1bSJed Brown  5 SNES Function norm 40.624
42f9635d15SStefano Zampini    0 KSP Residual norm 4.356840816356e-01
43f9635d15SStefano Zampini    1 KSP Residual norm 2.768016354304e-02
44f9635d15SStefano Zampini    2 KSP Residual norm 2.019234551823e-03
45ac7f1a8bSStefano Zampini    3 KSP Residual norm 1.079739269636e-04
46ac7f1a8bSStefano Zampini    4 KSP Residual norm 3.721836486707e-06
47ac7f1a8bSStefano Zampini    5 KSP Residual norm 2.878839455327e-07
48c4762a1bSJed Brown  6 SNES Function norm 12.1062
49f9635d15SStefano Zampini    0 KSP Residual norm 2.413828207078e-01
50f9635d15SStefano Zampini    1 KSP Residual norm 1.271470947688e-02
51ac7f1a8bSStefano Zampini    2 KSP Residual norm 9.046581777608e-04
52ac7f1a8bSStefano Zampini    3 KSP Residual norm 3.689804266215e-05
53ac7f1a8bSStefano Zampini    4 KSP Residual norm 2.906311127719e-06
54ac7f1a8bSStefano Zampini    5 KSP Residual norm 8.933152054376e-08
55c4762a1bSJed Brown  7 SNES Function norm 3.56967
56f9635d15SStefano Zampini    0 KSP Residual norm 1.051513481961e-01
57ac7f1a8bSStefano Zampini    1 KSP Residual norm 5.511357128251e-03
58f9635d15SStefano Zampini    2 KSP Residual norm 4.397748384999e-04
59f9635d15SStefano Zampini    3 KSP Residual norm 1.803541897729e-05
60ac7f1a8bSStefano Zampini    4 KSP Residual norm 1.281155384368e-06
61ac7f1a8bSStefano Zampini    5 KSP Residual norm 2.698307114036e-08
62c4762a1bSJed Brown  8 SNES Function norm 0.948696
63f9635d15SStefano Zampini    0 KSP Residual norm 4.448587457094e-02
64f9635d15SStefano Zampini    1 KSP Residual norm 1.666506843624e-03
65ac7f1a8bSStefano Zampini    2 KSP Residual norm 1.896722722892e-04
66ac7f1a8bSStefano Zampini    3 KSP Residual norm 7.781898834757e-06
67ac7f1a8bSStefano Zampini    4 KSP Residual norm 5.034969807078e-07
68ac7f1a8bSStefano Zampini    5 KSP Residual norm 1.624399106790e-08
69c4762a1bSJed Brown  9 SNES Function norm 0.224801
70f9635d15SStefano Zampini    0 KSP Residual norm 7.502762042127e-03
71f9635d15SStefano Zampini    1 KSP Residual norm 1.656093896390e-04
72f9635d15SStefano Zampini    2 KSP Residual norm 2.498021277745e-05
73ac7f1a8bSStefano Zampini    3 KSP Residual norm 8.800750827469e-07
74ac7f1a8bSStefano Zampini    4 KSP Residual norm 4.269071150712e-08
75ac7f1a8bSStefano Zampini    5 KSP Residual norm 2.587977253298e-09
76ac7f1a8bSStefano Zampini    6 KSP Residual norm 1.033048936067e-10
77c4762a1bSJed Brown 10 SNES Function norm 0.0400413
78ac7f1a8bSStefano Zampini    0 KSP Residual norm 5.771610774235e-04
79ac7f1a8bSStefano Zampini    1 KSP Residual norm 9.630086764651e-06
80ac7f1a8bSStefano Zampini    2 KSP Residual norm 7.184492044802e-07
81ac7f1a8bSStefano Zampini    3 KSP Residual norm 4.700789158564e-08
82ac7f1a8bSStefano Zampini    4 KSP Residual norm 1.719796413714e-09
83ac7f1a8bSStefano Zampini    5 KSP Residual norm 1.204413570031e-10
84ac7f1a8bSStefano Zampini    6 KSP Residual norm 5.983933181073e-12
85c4762a1bSJed Brown 11 SNES Function norm 0.00166356
86ac7f1a8bSStefano Zampini    0 KSP Residual norm 1.521708678812e-05
87ac7f1a8bSStefano Zampini    1 KSP Residual norm 2.816775207578e-07
88ac7f1a8bSStefano Zampini    2 KSP Residual norm 1.415288155243e-08
89ac7f1a8bSStefano Zampini    3 KSP Residual norm 1.055754682927e-09
90ac7f1a8bSStefano Zampini    4 KSP Residual norm 3.964183154450e-11
91ac7f1a8bSStefano Zampini    5 KSP Residual norm 3.034624910197e-12
92ac7f1a8bSStefano Zampini    6 KSP Residual norm 2.147890910700e-13
93ac7f1a8bSStefano Zampini    7 KSP Residual norm 6.319038359390e-15
94c4762a1bSJed Brown 12 SNES Function norm 2.78359e-06
95c4762a1bSJed Brown  Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 12
96c4762a1bSJed BrownSNES Object: 4 MPI processes
97c4762a1bSJed Brown  type: newtonls
98c4762a1bSJed Brown  maximum iterations=20, maximum function evaluations=10000
99c4762a1bSJed Brown  tolerances: relative=1e-08, absolute=1e-50, solution=1e-08
100f9635d15SStefano Zampini  total number of linear solver iterations=69
101c4762a1bSJed Brown  total number of function evaluations=13
102c4762a1bSJed Brown  norm schedule ALWAYS
103c4762a1bSJed Brown  SNESLineSearch Object: 4 MPI processes
104c4762a1bSJed Brown    type: bt
105c4762a1bSJed Brown      interpolation: cubic
106c4762a1bSJed Brown      alpha=1.000000e-04
107a99ef635SJonas Heinzmann    maxlambda=1.000000e+00, minlambda=1.000000e-12
108c4762a1bSJed Brown    tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
109c4762a1bSJed Brown    maximum iterations=40
110c4762a1bSJed Brown  KSP Object: 4 MPI processes
111c4762a1bSJed Brown    type: cg
112c4762a1bSJed Brown    maximum iterations=10000, initial guess is zero
113c4762a1bSJed Brown    tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
114c4762a1bSJed Brown    left preconditioning
115c4762a1bSJed Brown    using PRECONDITIONED norm type for convergence test
116c4762a1bSJed Brown  PC Object: 4 MPI processes
117c4762a1bSJed Brown    type: bddc
118c4762a1bSJed Brown      Use verbose output: 0
119c4762a1bSJed Brown      Use user-defined CSR: 0
120c4762a1bSJed Brown      Use local mat graph: 1
121c4762a1bSJed Brown      Connectivity graph topological dimension: 2
1220dfc91b7SStefano Zampini      Corner selection: 0 (selected 0)
123c4762a1bSJed Brown      Use vertices: 1 (vertex size 1)
124c4762a1bSJed Brown      Use edges: 1
125c4762a1bSJed Brown      Use faces: 0
126c4762a1bSJed Brown      Use true near null space: 0
127c4762a1bSJed Brown      Use QR for single constraints on cc: 0
128c4762a1bSJed Brown      Use change of basis on local edge nodes: 0
129c4762a1bSJed Brown      Use change of basis on local face nodes: 0
130c4762a1bSJed Brown      User defined change of basis matrix: 0
131c4762a1bSJed Brown      Has change of basis matrix: 0
132c4762a1bSJed Brown      Eliminate dirichlet boundary dofs: 0
133c4762a1bSJed Brown      Switch on static condensation ops around the interface preconditioner: 0
134c4762a1bSJed Brown      Use exact dirichlet trick: 1
135c4762a1bSJed Brown      Interface extension: DIRICHLET
136c4762a1bSJed Brown      Multilevel max levels: 0
137c4762a1bSJed Brown      Multilevel coarsening ratio: 8
138c4762a1bSJed Brown      Use estimated eigs for coarse problem: 0
139c4762a1bSJed Brown      Use deluxe scaling: 0
140c4762a1bSJed Brown      Use deluxe zerorows: 0
141c4762a1bSJed Brown      Use deluxe singlemat: 0
142c4762a1bSJed Brown      Rebuild interface graph for Schur principal minors: 0
143c4762a1bSJed Brown      Number of dofs' layers for the computation of principal minors: -1
144c4762a1bSJed Brown      Use user CSR graph to compute successive layers: 0
145c4762a1bSJed Brown      Adaptive constraint selection threshold (active 0, userdefined 0): 0.
146c4762a1bSJed Brown      Min constraints / connected component: 0
147c4762a1bSJed Brown      Max constraints / connected component: 0
148c4762a1bSJed Brown      Invert exact Schur complement for adaptive selection: 0
149c4762a1bSJed Brown      Symmetric computation of primal basis functions: 1
150c4762a1bSJed Brown      Num. Procs. to map coarse adjacency list: 0
151c4762a1bSJed Brown      Coarse eqs per proc (significant at the coarsest level): 1
152c4762a1bSJed Brown      Detect disconnected: 1 (filter 0)
153c4762a1bSJed Brown      Benign subspace trick: 0 (change explicit 0)
154c4762a1bSJed Brown      Benign subspace trick is active: 0
155c4762a1bSJed Brown      Algebraic computation of no-net-flux: 0
156c4762a1bSJed Brown    ********************************** STATISTICS AT LEVEL 0 **********************************
157f9635d15SStefano Zampini      Global dofs sizes: all 1485 interface 92 coarse 20
158f9635d15SStefano Zampini      Coarsening ratios: all/coarse 74 interface/coarse 4
159c4762a1bSJed Brown      Active processes : 4
160c4762a1bSJed Brown      Total subdomains : 5
161c4762a1bSJed Brown      Dofs type        :	MIN	MAX	MEAN
162c4762a1bSJed Brown      Interior  dofs   :	343	361	348
163c4762a1bSJed Brown      Interface dofs   :	33	61	46
164f9635d15SStefano Zampini      Primal    dofs   :	5	16	10
165c4762a1bSJed Brown      Local     dofs   :	376	422	395
166c4762a1bSJed Brown      Local     subs   :	1	2
167c4762a1bSJed Brown    --- Interior solver (rank 0)
1688cc725e6SPierre Jolivet      KSP Object: (pc_bddc_dirichlet_) 1 MPI process
169c4762a1bSJed Brown        type: preonly
170c4762a1bSJed Brown        maximum iterations=10000, initial guess is zero
171c4762a1bSJed Brown        tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
172c4762a1bSJed Brown        left preconditioning
1738405740aSBarry Smith        not checking for convergence
1748cc725e6SPierre Jolivet      PC Object: (pc_bddc_dirichlet_) 1 MPI process
175c4762a1bSJed Brown        type: lu
176c4762a1bSJed Brown          out-of-place factorization
177c4762a1bSJed Brown          tolerance for zero pivot 2.22045e-14
178c4762a1bSJed Brown          matrix ordering: nd
179c4762a1bSJed Brown          factor fill ratio given 5., needed 2.23015
180*ecf3d421SBarry Smith            Factored matrix:
18126cc229bSBarry Smith              Mat Object: (pc_bddc_dirichlet_) 1 MPI process
182c4762a1bSJed Brown                type: seqaij
183c4762a1bSJed Brown                rows=345, cols=345
184c4762a1bSJed Brown                package used to perform factorization: petsc
185c4762a1bSJed Brown                total: nonzeros=10533, allocated nonzeros=10533
186c4762a1bSJed Brown                  not using I-node routines
187*ecf3d421SBarry Smith        linear system matrix, which is also used to construct the preconditioner:
1888cc725e6SPierre Jolivet        Mat Object: (pc_bddc_dirichlet_) 1 MPI process
189c4762a1bSJed Brown          type: seqaij
190c4762a1bSJed Brown          rows=345, cols=345
191c4762a1bSJed Brown          total: nonzeros=4723, allocated nonzeros=4723
192c4762a1bSJed Brown          total number of mallocs used during MatSetValues calls=0
193c4762a1bSJed Brown            not using I-node routines
194c4762a1bSJed Brown    --- Correction solver (rank 0)
1958cc725e6SPierre Jolivet      KSP Object: (pc_bddc_neumann_) 1 MPI process
196c4762a1bSJed Brown        type: preonly
197c4762a1bSJed Brown        maximum iterations=10000, initial guess is zero
198c4762a1bSJed Brown        tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
199c4762a1bSJed Brown        left preconditioning
2008405740aSBarry Smith        not checking for convergence
2018cc725e6SPierre Jolivet      PC Object: (pc_bddc_neumann_) 1 MPI process
202c4762a1bSJed Brown        type: lu
203c4762a1bSJed Brown          out-of-place factorization
204c4762a1bSJed Brown          tolerance for zero pivot 2.22045e-14
205c4762a1bSJed Brown          matrix ordering: nd
206f9635d15SStefano Zampini          factor fill ratio given 5., needed 2.20116
207*ecf3d421SBarry Smith            Factored matrix:
20826cc229bSBarry Smith              Mat Object: (pc_bddc_neumann_) 1 MPI process
209c4762a1bSJed Brown                type: seqaij
210f9635d15SStefano Zampini                rows=382, cols=382
211c4762a1bSJed Brown                package used to perform factorization: petsc
212f9635d15SStefano Zampini                total: nonzeros=12146, allocated nonzeros=12146
213c4762a1bSJed Brown                  not using I-node routines
214*ecf3d421SBarry Smith        linear system matrix, which is also used to construct the preconditioner:
2158cc725e6SPierre Jolivet        Mat Object: (pc_bddc_neumann_) 1 MPI process
216c4762a1bSJed Brown          type: seqaij
217f9635d15SStefano Zampini          rows=382, cols=382
218f9635d15SStefano Zampini          total: nonzeros=5518, allocated nonzeros=5518
219c4762a1bSJed Brown          total number of mallocs used during MatSetValues calls=0
220c4762a1bSJed Brown            not using I-node routines
221c4762a1bSJed Brown    --- Coarse solver
222c4762a1bSJed Brown      KSP Object: (pc_bddc_coarse_) 4 MPI processes
223c4762a1bSJed Brown        type: preonly
224c4762a1bSJed Brown        maximum iterations=1, initial guess is zero
225c4762a1bSJed Brown        tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
226c4762a1bSJed Brown        left preconditioning
2278405740aSBarry Smith        not checking for convergence
228c4762a1bSJed Brown      PC Object: (pc_bddc_coarse_) 4 MPI processes
229c4762a1bSJed Brown        type: redundant
230c4762a1bSJed Brown          First (color=0) of 4 PCs follows
2318cc725e6SPierre Jolivet          KSP Object: (pc_bddc_coarse_redundant_) 1 MPI process
232c4762a1bSJed Brown            type: preonly
233c4762a1bSJed Brown            maximum iterations=10000, initial guess is zero
234c4762a1bSJed Brown            tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
235c4762a1bSJed Brown            left preconditioning
2368405740aSBarry Smith            not checking for convergence
2378cc725e6SPierre Jolivet          PC Object: (pc_bddc_coarse_redundant_) 1 MPI process
238c4762a1bSJed Brown            type: lu
239c4762a1bSJed Brown              out-of-place factorization
240c4762a1bSJed Brown              tolerance for zero pivot 2.22045e-14
241c4762a1bSJed Brown              matrix ordering: nd
242f9635d15SStefano Zampini              factor fill ratio given 5., needed 1.00535
243*ecf3d421SBarry Smith                Factored matrix:
24426cc229bSBarry Smith                  Mat Object: (pc_bddc_coarse_redundant_) 1 MPI process
245c4762a1bSJed Brown                    type: seqaij
246f9635d15SStefano Zampini                    rows=20, cols=20
247c4762a1bSJed Brown                    package used to perform factorization: petsc
248f9635d15SStefano Zampini                    total: nonzeros=376, allocated nonzeros=376
249f9635d15SStefano Zampini                      using I-node routines: found 5 nodes, limit used is 5
250*ecf3d421SBarry Smith            linear system matrix, which is also used to construct the preconditioner:
25126cc229bSBarry Smith            Mat Object: 1 MPI process
252c4762a1bSJed Brown              type: seqaij
253f9635d15SStefano Zampini              rows=20, cols=20
254f9635d15SStefano Zampini              total: nonzeros=374, allocated nonzeros=374
255c4762a1bSJed Brown              total number of mallocs used during MatSetValues calls=0
256f9635d15SStefano Zampini                using I-node routines: found 8 nodes, limit used is 5
257*ecf3d421SBarry Smith        linear system matrix, which is also used to construct the preconditioner:
2585a67bb2bSBarry Smith        Mat Object: (pc_bddc_coarse_) 4 MPI processes
2595a67bb2bSBarry Smith          type: mpiaij
260f9635d15SStefano Zampini          rows=20, cols=20
261ac7f1a8bSStefano Zampini          total: nonzeros=374, allocated nonzeros=374
2625a67bb2bSBarry Smith          total number of mallocs used during MatSetValues calls=0
2635a67bb2bSBarry Smith            using I-node (on process 0) routines: found 1 nodes, limit used is 5
264*ecf3d421SBarry Smith    linear system matrix, which is also used to construct the preconditioner:
265c4762a1bSJed Brown    Mat Object: 4 MPI processes
266c4762a1bSJed Brown      type: is
267c4762a1bSJed Brown      rows=1485, cols=1485
268c4762a1bSJed Brown      total: nonzeros=23046, allocated nonzeros=24772
269c4762a1bSJed Brown      total number of mallocs used during MatSetValues calls=0
270