1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
|
0 SNES Function norm 0.370399
KSP Object:
type: fgmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
right preconditioning
using UNPRECONDITIONED norm type for convergence test
PC Object:
type: fieldsplit
FieldSplit with MULTIPLICATIVE composition: total splits = 2, blocksize = -1
Solver info for each split is in the following KSP objects:
Split number 0 Defined by IS
KSP Object:(fieldsplit_0_)
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using PRECONDITIONED norm type for convergence test
PC Object:(fieldsplit_0_)
type: ilu
ILU: out-of-place factorization
0 levels of fill
tolerance for zero pivot 1e-12
using diagonal shift to prevent zero pivot
matrix ordering: natural
factor fill ratio given 1, needed 1
Factored matrix follows:
Matrix Object:
type=seqaij, rows=48, cols=48
package used to perform factorization: petsc
total: nonzeros=576, allocated nonzeros=576
using I-node routines: found 16 nodes, limit used is 5
linear system matrix = precond matrix:
Matrix Object:
type=seqaij, rows=48, cols=48
total: nonzeros=576, allocated nonzeros=576
using I-node routines: found 16 nodes, limit used is 5
Split number 1 Defined by IS
KSP Object:(fieldsplit_1_)
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using PRECONDITIONED norm type for convergence test
PC Object:(fieldsplit_1_)
type: ilu
ILU: out-of-place factorization
0 levels of fill
tolerance for zero pivot 1e-12
using diagonal shift to prevent zero pivot
matrix ordering: natural
factor fill ratio given 1, needed 1
Factored matrix follows:
Matrix Object:
type=seqaij, rows=16, cols=16
package used to perform factorization: petsc
total: nonzeros=64, allocated nonzeros=64
not using I-node routines
linear system matrix = precond matrix:
Matrix Object:
type=seqaij, rows=16, cols=16
total: nonzeros=64, allocated nonzeros=64
not using I-node routines
linear system matrix followed by preconditioner matrix:
Matrix Object:
type=mffd, rows=64, cols=64
matrix-free approximation:
err=1e-07 (relative error in function evaluation)
Using wp compute h routine
Computes normA
Does not compute normU
Matrix Object:
type=seqaij, rows=64, cols=64
total: nonzeros=696, allocated nonzeros=696
using I-node routines: found 48 nodes, limit used is 5
1 SNES Function norm 6.81478e-05
KSP Object:
type: fgmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
right preconditioning
using UNPRECONDITIONED norm type for convergence test
PC Object:
type: fieldsplit
FieldSplit with MULTIPLICATIVE composition: total splits = 2, blocksize = -1
Solver info for each split is in the following KSP objects:
Split number 0 Defined by IS
KSP Object:(fieldsplit_0_)
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using PRECONDITIONED norm type for convergence test
PC Object:(fieldsplit_0_)
type: ilu
ILU: out-of-place factorization
0 levels of fill
tolerance for zero pivot 1e-12
using diagonal shift to prevent zero pivot
matrix ordering: natural
factor fill ratio given 1, needed 1
Factored matrix follows:
Matrix Object:
type=seqaij, rows=48, cols=48
package used to perform factorization: petsc
total: nonzeros=576, allocated nonzeros=576
using I-node routines: found 16 nodes, limit used is 5
linear system matrix = precond matrix:
Matrix Object:
type=seqaij, rows=48, cols=48
total: nonzeros=576, allocated nonzeros=576
using I-node routines: found 16 nodes, limit used is 5
Split number 1 Defined by IS
KSP Object:(fieldsplit_1_)
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using PRECONDITIONED norm type for convergence test
PC Object:(fieldsplit_1_)
type: ilu
ILU: out-of-place factorization
0 levels of fill
tolerance for zero pivot 1e-12
using diagonal shift to prevent zero pivot
matrix ordering: natural
factor fill ratio given 1, needed 1
Factored matrix follows:
Matrix Object:
type=seqaij, rows=16, cols=16
package used to perform factorization: petsc
total: nonzeros=64, allocated nonzeros=64
not using I-node routines
linear system matrix = precond matrix:
Matrix Object:
type=seqaij, rows=16, cols=16
total: nonzeros=64, allocated nonzeros=64
not using I-node routines
linear system matrix followed by preconditioner matrix:
Matrix Object:
type=mffd, rows=64, cols=64
matrix-free approximation:
err=1e-07 (relative error in function evaluation)
Using wp compute h routine
Computes normA
Does not compute normU
Matrix Object:
type=seqaij, rows=64, cols=64
total: nonzeros=696, allocated nonzeros=696
using I-node routines: found 48 nodes, limit used is 5
2 SNES Function norm 5.968e-11
Composite Physics: Number of Newton iterations = 2
|
