[petsc-users] Multigrid preconditioning of entire linear systems for discretized coupled multiphysics problems
Fabian Gabel
gabel.fabian at gmail.com
Mon Mar 2 21:19:28 CST 2015
On Mo, 2015-03-02 at 19:43 -0600, Barry Smith wrote:
> Do you really want tolerances: relative=1e-90, absolute=1.10423, divergence=10000? That is an absolute tolerance of 1.1? Normally that would be huge.
I started using atol as convergence criterion with -ksp_norm_type
unpreconditioned. The value of atol gets updated every outer iteration.
> You can provide your matrix with a block size that GAMG will use with MatSetBlockSize().
I think something went wrong. Setting the block size to 4 and solving
for (u,v,w,p) the convergence degraded significantly. I attached the
results for a smaller test case that shows the increase of the number of
needed inner iterations when setting the block size via
MatSetBlockSize().
>
> I would use coupledsolve_mg_coarse_sub_pc_type lu it is weird that it is using SOR for 27 points.
>
> So you must have provided a null space since it printed "has attached null space"
The system has indeed a one dimensional null space (from the pressure
equation with Neumann boundary conditions). But now that you mentioned
it: It seems that the outer GMRES doesn't notice that the matrix has an
attached nullspace. Replacing
CALL MatSetNullSpace(CMAT,NULLSP,IERR)
with
CALL KSPSetNullSpace(KRYLOV,NULLSP,IERR)
solves this. What is wrong with using MatSetNullSpace?
Fabian
>
> Barry
>
>
>
> > On Mar 2, 2015, at 6:39 PM, Fabian Gabel <gabel.fabian at gmail.com> wrote:
> >
> > On Mo, 2015-03-02 at 16:29 -0700, Jed Brown wrote:
> >> Fabian Gabel <gabel.fabian at gmail.com> writes:
> >>
> >>> Dear PETSc Team,
> >>>
> >>> I came across the following paragraph in your publication "Composable
> >>> Linear Solvers for Multiphysics" (2012):
> >>>
> >>> "Rather than splitting the matrix into large blocks and
> >>> forming a preconditioner from solvers (for example, multi-
> >>> grid) on each block, one can perform multigrid on the entire
> >>> system, basing the smoother on solves coming from the tiny
> >>> blocks coupling the degrees of freedom at a single point (or
> >>> small number of points). This approach is also handled in
> >>> PETSc, but we will not elaborate on it here."
> >>>
> >>> How would I use a multigrid preconditioner (GAMG)
> >>
> >> The heuristics in GAMG are not appropriate for indefinite/saddle-point
> >> systems such as arise from Navier-Stokes. You can use geometric
> >> multigrid and use the fieldsplit techniques described in the paper as a
> >> smoother, for example.
> >
> > I sadly don't have a solid background on multigrid methods, but as
> > mentioned in a previous thread
> >
> > http://lists.mcs.anl.gov/pipermail/petsc-users/2015-February/024219.html
> >
> > AMG has apparently been used (successfully?) for fully-coupled
> > finite-volume discretizations of Navier-Stokes:
> >
> > http://dx.doi.org/10.1080/10407790.2014.894448
> > http://dx.doi.org/10.1016/j.jcp.2008.08.027
> >
> > I was hoping to achieve something similar with the right configuration
> > of the PETSc preconditioners. So far I have only been using GAMG in a
> > straightforward manner, without providing any details on the structure
> > of the linear system. I attached the output of a test run with GAMG.
> >
> >>
> >>> from PETSc on linear systems of the form (after reordering the
> >>> variables):
> >>>
> >>> [A_uu 0 0 A_up A_uT]
> >>> [0 A_vv 0 A_vp A_vT]
> >>> [0 0 A_ww A_up A_wT]
> >>> [A_pu A_pv A_pw A_pp 0 ]
> >>> [A_Tu A_Tv A_Tw A_Tp A_TT]
> >>>
> >>> where each of the block matrices A_ij, with i,j in {u,v,w,p,T}, results
> >>> directly from a FVM discretization of the incompressible Navier-Stokes
> >>> equations and the temperature equation. The fifth row and column are
> >>> optional, depending on the method I choose to couple the temperature.
> >>> The Matrix is stored as one AIJ Matrix.
> >>>
> >>> Regards,
> >>> Fabian Gabel
> >
> > <cpld_0128.out.578677>
>
-------------- next part --------------
ENTER PROBLEM NAME (SIX CHARACTERS):
****************************************************
NAME OF PROBLEM SOLVED control
****************************************************
***************************************************
CONTROL SETTINGS
***************************************************
LREAD,LWRITE,LPOST,LTEST,LOUTS,LOUTE,LTIME,LGRAD
F F F F F F F F
IMON, JMON, KMON, MMON, RMON, IPR, JPR, KPR, MPR,NPCOR,NIGRAD
2 2 2 1 0 2 2 3 1 1 1
SORMAX, SLARGE, ALFA
0.1000E-07 0.1000E+31 0.9200E+00
(URF(I),I=1,6)
0.9000E+00 0.9000E+00 0.9000E+00 0.1000E+00 0.1000E+01 0.1000E+01
(SOR(I),I=1,6)
0.1000E+00 0.1000E+00 0.1000E+00 0.1000E-01 0.1000E+00 0.1000E+00
(GDS(I),I=1,6) - BLENDING (CDS-UDS)
0.1000E+01 0.1000E+01 0.1000E+01 0.1600E+01 0.1000E+01 0.0000E+00
LSG
10000
***************************************************
START COUPLED ALGORITHM
***************************************************
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 7.680650991441e+01
1 KSP Residual norm 4.185899808477e+01
2 KSP Residual norm 2.202068303753e+01
3 KSP Residual norm 1.741802539838e+01
4 KSP Residual norm 1.364396531738e+01
5 KSP Residual norm 1.172273211505e+01
6 KSP Residual norm 1.162524589201e+01
7 KSP Residual norm 1.011907279997e+01
8 KSP Residual norm 8.685677152271e+00
9 KSP Residual norm 7.396125459146e+00
10 KSP Residual norm 7.296977130121e+00
11 KSP Residual norm 6.959200481951e+00
12 KSP Residual norm 6.905446814136e+00
13 KSP Residual norm 6.880720024552e+00
14 KSP Residual norm 6.824350801593e+00
15 KSP Residual norm 6.824016557799e+00
16 KSP Residual norm 4.724619270852e+00
17 KSP Residual norm 2.668642871524e+00
18 KSP Residual norm 2.191273302044e+00
19 KSP Residual norm 1.944296054954e+00
20 KSP Residual norm 1.900543170284e+00
21 KSP Residual norm 1.564163873693e+00
22 KSP Residual norm 1.506139846488e+00
23 KSP Residual norm 1.505810088310e+00
24 KSP Residual norm 1.503963939770e+00
25 KSP Residual norm 1.423291747222e+00
26 KSP Residual norm 1.396589263932e+00
27 KSP Residual norm 1.382496701568e+00
28 KSP Residual norm 1.111467646867e+00
KSP Object:(coupledsolve_) 1 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=10000, initial guess is zero
tolerances: relative=1e-90, absolute=1.20082, divergence=10000
right preconditioning
has attached null space
using UNPRECONDITIONED norm type for convergence test
PC Object:(coupledsolve_) 1 MPI processes
type: gamg
MG: type is MULTIPLICATIVE, levels=3 cycles=v
Cycles per PCApply=1
Using Galerkin computed coarse grid matrices
Coarse grid solver -- level -------------------------------
KSP Object: (coupledsolve_mg_coarse_) 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: (coupledsolve_mg_coarse_) 1 MPI processes
type: bjacobi
block Jacobi: number of blocks = 1
Local solve is same for all blocks, in the following KSP and PC objects:
KSP Object: (coupledsolve_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: (coupledsolve_mg_coarse_sub_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=224, cols=224, bs=4
total: nonzeros=24768, allocated nonzeros=24768
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 56 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=224, cols=224, bs=4
total: nonzeros=24768, allocated nonzeros=24768
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 56 nodes, limit used is 5
Down solver (pre-smoother) on level 1 -------------------------------
KSP Object: (coupledsolve_mg_levels_1_) 1 MPI processes
type: richardson
Richardson: damping factor=1
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: (coupledsolve_mg_levels_1_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=12260, cols=12260, bs=4
total: nonzeros=1.42581e+06, allocated nonzeros=1.42581e+06
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 3065 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 2 -------------------------------
KSP Object: (coupledsolve_mg_levels_2_) 1 MPI processes
type: richardson
Richardson: damping factor=1
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: (coupledsolve_mg_levels_2_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=256000, cols=256000, bs=4
total: nonzeros=2.35725e+06, allocated nonzeros=2.35725e+06
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Up solver (post-smoother) same as down solver (pre-smoother)
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=256000, cols=256000, bs=4
total: nonzeros=2.35725e+06, allocated nonzeros=2.35725e+06
total number of mallocs used during MatSetValues calls =0
not using I-node routines
0000001 0.1000E+01 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 2.627622580920e+00
1 KSP Residual norm 2.618623510400e+00
2 KSP Residual norm 2.595681404299e+00
3 KSP Residual norm 2.573101727243e+00
4 KSP Residual norm 2.471252938050e+00
5 KSP Residual norm 2.428028929045e+00
6 KSP Residual norm 2.229168066194e+00
7 KSP Residual norm 2.085615414410e+00
8 KSP Residual norm 2.076282297913e+00
9 KSP Residual norm 2.069139125163e+00
10 KSP Residual norm 2.058700104333e+00
11 KSP Residual norm 2.058567796498e+00
12 KSP Residual norm 1.920870596572e+00
13 KSP Residual norm 1.777523400957e+00
14 KSP Residual norm 1.674384000336e+00
15 KSP Residual norm 1.642534139399e+00
16 KSP Residual norm 1.449777567903e+00
17 KSP Residual norm 1.440931159108e+00
18 KSP Residual norm 1.439213065379e+00
19 KSP Residual norm 1.382404265336e+00
20 KSP Residual norm 1.364215277267e+00
21 KSP Residual norm 1.362114300120e+00
22 KSP Residual norm 1.339978879706e+00
23 KSP Residual norm 1.322639467952e+00
24 KSP Residual norm 1.321384304678e+00
25 KSP Residual norm 1.321380759642e+00
26 KSP Residual norm 1.183186737875e+00
27 KSP Residual norm 8.659328686803e-01
28 KSP Residual norm 8.388085654581e-01
29 KSP Residual norm 8.346618418315e-01
30 KSP Residual norm 7.932132874906e-01
31 KSP Residual norm 7.931916049632e-01
32 KSP Residual norm 7.930070151834e-01
33 KSP Residual norm 7.927622711664e-01
34 KSP Residual norm 7.924352125863e-01
35 KSP Residual norm 7.923331320223e-01
36 KSP Residual norm 7.920591938247e-01
37 KSP Residual norm 7.914444071308e-01
38 KSP Residual norm 7.910824991274e-01
39 KSP Residual norm 7.910587420876e-01
40 KSP Residual norm 7.908788515020e-01
41 KSP Residual norm 7.895603769794e-01
42 KSP Residual norm 7.895589633648e-01
43 KSP Residual norm 7.891864300420e-01
44 KSP Residual norm 7.891326664165e-01
45 KSP Residual norm 7.888105721332e-01
46 KSP Residual norm 7.886174398811e-01
47 KSP Residual norm 7.869613760972e-01
48 KSP Residual norm 7.869597845056e-01
49 KSP Residual norm 7.861600492870e-01
50 KSP Residual norm 7.855646154910e-01
51 KSP Residual norm 7.854192275274e-01
52 KSP Residual norm 7.845302897412e-01
53 KSP Residual norm 7.845131211745e-01
54 KSP Residual norm 7.800888294949e-01
55 KSP Residual norm 7.782372196955e-01
56 KSP Residual norm 7.739281372664e-01
57 KSP Residual norm 7.703291682892e-01
58 KSP Residual norm 7.594343557276e-01
59 KSP Residual norm 7.428628814220e-01
60 KSP Residual norm 6.796108855423e-01
61 KSP Residual norm 6.792621836949e-01
62 KSP Residual norm 6.792599058243e-01
63 KSP Residual norm 6.790825225897e-01
64 KSP Residual norm 6.788952262666e-01
65 KSP Residual norm 6.775543713534e-01
66 KSP Residual norm 6.762003547942e-01
67 KSP Residual norm 6.750660083508e-01
68 KSP Residual norm 6.749056179696e-01
69 KSP Residual norm 6.745788595495e-01
70 KSP Residual norm 6.744831909936e-01
71 KSP Residual norm 6.744095168640e-01
72 KSP Residual norm 6.739619993802e-01
73 KSP Residual norm 6.733835789845e-01
74 KSP Residual norm 6.721515742494e-01
75 KSP Residual norm 6.717972280897e-01
76 KSP Residual norm 6.717926885501e-01
77 KSP Residual norm 6.711395282193e-01
78 KSP Residual norm 6.710189393745e-01
79 KSP Residual norm 6.708920353331e-01
80 KSP Residual norm 6.705065744993e-01
81 KSP Residual norm 6.702492254576e-01
82 KSP Residual norm 6.654066904612e-01
83 KSP Residual norm 6.633284243713e-01
84 KSP Residual norm 6.588735035312e-01
85 KSP Residual norm 6.588203269869e-01
86 KSP Residual norm 6.588021175950e-01
87 KSP Residual norm 6.582004015484e-01
88 KSP Residual norm 6.194837793643e-01
89 KSP Residual norm 5.961955972874e-01
90 KSP Residual norm 5.705567260811e-01
91 KSP Residual norm 5.705515967545e-01
92 KSP Residual norm 5.704003341423e-01
93 KSP Residual norm 5.703233684974e-01
94 KSP Residual norm 5.698144671527e-01
95 KSP Residual norm 5.698140906127e-01
96 KSP Residual norm 5.697148973128e-01
97 KSP Residual norm 5.685830320447e-01
98 KSP Residual norm 5.685420277423e-01
99 KSP Residual norm 5.683463895401e-01
100 KSP Residual norm 5.683244725734e-01
101 KSP Residual norm 5.681581335505e-01
102 KSP Residual norm 5.676435098591e-01
103 KSP Residual norm 5.664466011033e-01
104 KSP Residual norm 5.654075348607e-01
105 KSP Residual norm 5.653642207258e-01
106 KSP Residual norm 5.639595704613e-01
107 KSP Residual norm 5.619972750669e-01
108 KSP Residual norm 5.605485908660e-01
109 KSP Residual norm 5.604783703868e-01
110 KSP Residual norm 5.576368973879e-01
111 KSP Residual norm 5.576368717003e-01
112 KSP Residual norm 5.571870008692e-01
113 KSP Residual norm 5.546719813159e-01
114 KSP Residual norm 5.545657244427e-01
115 KSP Residual norm 5.537810648153e-01
116 KSP Residual norm 5.533802839194e-01
117 KSP Residual norm 5.520193889595e-01
118 KSP Residual norm 5.093820108693e-01
119 KSP Residual norm 4.531399469657e-01
120 KSP Residual norm 2.957173921159e-01
121 KSP Residual norm 2.956784310673e-01
122 KSP Residual norm 2.950937537726e-01
123 KSP Residual norm 2.867092781331e-01
124 KSP Residual norm 2.863531544699e-01
125 KSP Residual norm 2.847392979451e-01
126 KSP Residual norm 2.833809327431e-01
127 KSP Residual norm 2.829158304896e-01
128 KSP Residual norm 2.822111361138e-01
129 KSP Residual norm 2.785559895759e-01
130 KSP Residual norm 2.595237994249e-01
0000002 0.2188E+00 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 5.804789475348e-01
1 KSP Residual norm 5.799778769779e-01
2 KSP Residual norm 5.641496919751e-01
3 KSP Residual norm 5.617959304693e-01
4 KSP Residual norm 5.614011665870e-01
5 KSP Residual norm 5.569451426114e-01
6 KSP Residual norm 5.556608778653e-01
7 KSP Residual norm 5.467421792576e-01
8 KSP Residual norm 5.467265438841e-01
9 KSP Residual norm 5.353090266132e-01
10 KSP Residual norm 4.746622832951e-01
11 KSP Residual norm 4.702818488973e-01
12 KSP Residual norm 4.406372118449e-01
13 KSP Residual norm 4.330662669124e-01
14 KSP Residual norm 4.288057798951e-01
15 KSP Residual norm 4.046274693057e-01
16 KSP Residual norm 4.046230837327e-01
17 KSP Residual norm 4.042511516735e-01
18 KSP Residual norm 4.037663498063e-01
19 KSP Residual norm 3.931216729466e-01
20 KSP Residual norm 3.921031432605e-01
21 KSP Residual norm 3.429800353658e-01
22 KSP Residual norm 2.979221673663e-01
23 KSP Residual norm 2.681645180400e-01
24 KSP Residual norm 2.495132152627e-01
25 KSP Residual norm 2.279615743396e-01
26 KSP Residual norm 2.186269030034e-01
27 KSP Residual norm 2.186147575889e-01
28 KSP Residual norm 2.184196067179e-01
29 KSP Residual norm 2.126159051784e-01
30 KSP Residual norm 2.069760657158e-01
31 KSP Residual norm 2.069746584589e-01
32 KSP Residual norm 2.069731656236e-01
33 KSP Residual norm 2.069707108122e-01
34 KSP Residual norm 2.069685614175e-01
35 KSP Residual norm 2.069685610728e-01
36 KSP Residual norm 2.069679499377e-01
37 KSP Residual norm 2.069667560718e-01
38 KSP Residual norm 2.069658387974e-01
39 KSP Residual norm 2.069656414990e-01
40 KSP Residual norm 2.069630365184e-01
41 KSP Residual norm 2.069630254407e-01
42 KSP Residual norm 2.069600524211e-01
43 KSP Residual norm 2.069528577985e-01
44 KSP Residual norm 2.069375150385e-01
45 KSP Residual norm 2.069214325540e-01
46 KSP Residual norm 2.069193513908e-01
47 KSP Residual norm 2.069027188948e-01
48 KSP Residual norm 2.069026823825e-01
49 KSP Residual norm 2.068953523635e-01
50 KSP Residual norm 2.068777919687e-01
51 KSP Residual norm 2.067685327874e-01
52 KSP Residual norm 2.067523582887e-01
53 KSP Residual norm 2.067099161030e-01
54 KSP Residual norm 2.065528777162e-01
55 KSP Residual norm 2.064349053912e-01
56 KSP Residual norm 2.062439881285e-01
57 KSP Residual norm 2.048486269207e-01
58 KSP Residual norm 2.032421705358e-01
59 KSP Residual norm 1.975396841907e-01
60 KSP Residual norm 1.795975675884e-01
61 KSP Residual norm 1.795934380605e-01
62 KSP Residual norm 1.795171646958e-01
63 KSP Residual norm 1.793585200727e-01
64 KSP Residual norm 1.791359937468e-01
65 KSP Residual norm 1.787155890634e-01
66 KSP Residual norm 1.787083647740e-01
67 KSP Residual norm 1.785741702253e-01
68 KSP Residual norm 1.782753198644e-01
69 KSP Residual norm 1.777962163575e-01
70 KSP Residual norm 1.775966385004e-01
71 KSP Residual norm 1.775381320053e-01
72 KSP Residual norm 1.775289200005e-01
73 KSP Residual norm 1.771227389228e-01
74 KSP Residual norm 1.769384491898e-01
75 KSP Residual norm 1.769231925523e-01
76 KSP Residual norm 1.767348805780e-01
77 KSP Residual norm 1.767316532455e-01
78 KSP Residual norm 1.766472311583e-01
79 KSP Residual norm 1.766439185473e-01
80 KSP Residual norm 1.760015558603e-01
81 KSP Residual norm 1.759613105921e-01
82 KSP Residual norm 1.759281035337e-01
83 KSP Residual norm 1.758447163326e-01
84 KSP Residual norm 1.753046343083e-01
85 KSP Residual norm 1.748441698476e-01
86 KSP Residual norm 1.748374509920e-01
87 KSP Residual norm 1.745938060587e-01
88 KSP Residual norm 1.735307552311e-01
89 KSP Residual norm 1.730258903974e-01
90 KSP Residual norm 1.697175326546e-01
91 KSP Residual norm 1.697173566959e-01
92 KSP Residual norm 1.697093289476e-01
93 KSP Residual norm 1.696982044568e-01
94 KSP Residual norm 1.696879274337e-01
95 KSP Residual norm 1.696855069974e-01
96 KSP Residual norm 1.696835982074e-01
97 KSP Residual norm 1.696710777240e-01
98 KSP Residual norm 1.696622260091e-01
99 KSP Residual norm 1.695969278837e-01
100 KSP Residual norm 1.695836616942e-01
101 KSP Residual norm 1.695834190045e-01
102 KSP Residual norm 1.695812895531e-01
103 KSP Residual norm 1.694567934722e-01
104 KSP Residual norm 1.694195781071e-01
105 KSP Residual norm 1.693701486949e-01
106 KSP Residual norm 1.693482350903e-01
107 KSP Residual norm 1.693142066308e-01
108 KSP Residual norm 1.693139366048e-01
109 KSP Residual norm 1.692893155871e-01
110 KSP Residual norm 1.688223587231e-01
111 KSP Residual norm 1.688197233297e-01
112 KSP Residual norm 1.688163962887e-01
113 KSP Residual norm 1.688049236283e-01
114 KSP Residual norm 1.684781906985e-01
115 KSP Residual norm 1.683539285586e-01
116 KSP Residual norm 1.682149256972e-01
117 KSP Residual norm 1.681685812286e-01
118 KSP Residual norm 1.659367994851e-01
119 KSP Residual norm 1.634765294983e-01
120 KSP Residual norm 1.425866729058e-01
121 KSP Residual norm 1.425841891153e-01
122 KSP Residual norm 1.425662768676e-01
123 KSP Residual norm 1.424811258758e-01
124 KSP Residual norm 1.424663772468e-01
125 KSP Residual norm 1.424392146345e-01
126 KSP Residual norm 1.423428247695e-01
127 KSP Residual norm 1.423379715786e-01
128 KSP Residual norm 1.423161428637e-01
129 KSP Residual norm 1.422826009846e-01
130 KSP Residual norm 1.421362715567e-01
131 KSP Residual norm 1.421297151769e-01
132 KSP Residual norm 1.418970488464e-01
133 KSP Residual norm 1.405045151205e-01
134 KSP Residual norm 1.404234190366e-01
135 KSP Residual norm 1.403729297668e-01
136 KSP Residual norm 1.397587138297e-01
137 KSP Residual norm 1.393043951726e-01
138 KSP Residual norm 1.392643702290e-01
139 KSP Residual norm 1.390947717975e-01
140 KSP Residual norm 1.389847569971e-01
141 KSP Residual norm 1.386182745937e-01
142 KSP Residual norm 1.375593051658e-01
143 KSP Residual norm 1.288200479354e-01
144 KSP Residual norm 1.245208235838e-01
145 KSP Residual norm 1.228085303126e-01
146 KSP Residual norm 1.224380746235e-01
147 KSP Residual norm 1.203742807462e-01
148 KSP Residual norm 1.149199456902e-01
149 KSP Residual norm 8.452247950761e-02
150 KSP Residual norm 5.233620210445e-02
0000003 0.4834E-01 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 2.195661259895e-01
1 KSP Residual norm 2.058279904671e-01
2 KSP Residual norm 2.046045866307e-01
3 KSP Residual norm 2.033797944751e-01
4 KSP Residual norm 2.012090188960e-01
5 KSP Residual norm 2.010479779691e-01
6 KSP Residual norm 1.953230724132e-01
7 KSP Residual norm 1.555726024661e-01
8 KSP Residual norm 1.502318724920e-01
9 KSP Residual norm 1.489692702783e-01
10 KSP Residual norm 1.460165310535e-01
11 KSP Residual norm 1.446268056783e-01
12 KSP Residual norm 1.407909493702e-01
13 KSP Residual norm 1.158026701992e-01
14 KSP Residual norm 1.157234255633e-01
15 KSP Residual norm 1.056728059699e-01
16 KSP Residual norm 9.476360181621e-02
17 KSP Residual norm 9.353972455501e-02
18 KSP Residual norm 9.353942709547e-02
19 KSP Residual norm 8.224859685736e-02
20 KSP Residual norm 8.200143406325e-02
21 KSP Residual norm 7.772887684678e-02
22 KSP Residual norm 7.700008764377e-02
23 KSP Residual norm 7.450258745835e-02
24 KSP Residual norm 6.528939332618e-02
25 KSP Residual norm 6.128315999136e-02
26 KSP Residual norm 6.015626596536e-02
27 KSP Residual norm 5.779857199939e-02
28 KSP Residual norm 5.688458134959e-02
29 KSP Residual norm 4.930918193192e-02
30 KSP Residual norm 4.413319071243e-02
31 KSP Residual norm 4.412509475939e-02
32 KSP Residual norm 4.393555086585e-02
33 KSP Residual norm 4.391898171710e-02
34 KSP Residual norm 4.391354959236e-02
35 KSP Residual norm 4.384162021232e-02
36 KSP Residual norm 4.325311444780e-02
37 KSP Residual norm 4.288365140046e-02
38 KSP Residual norm 4.287209988977e-02
39 KSP Residual norm 4.287206142933e-02
40 KSP Residual norm 4.286686909881e-02
41 KSP Residual norm 4.260594081690e-02
42 KSP Residual norm 4.183743120137e-02
43 KSP Residual norm 4.125657765887e-02
44 KSP Residual norm 4.083515451508e-02
45 KSP Residual norm 4.025039016469e-02
46 KSP Residual norm 4.014739619114e-02
47 KSP Residual norm 3.982092444641e-02
48 KSP Residual norm 3.958395232305e-02
49 KSP Residual norm 3.952931641207e-02
50 KSP Residual norm 3.912847460622e-02
51 KSP Residual norm 3.871947928625e-02
52 KSP Residual norm 3.867946020933e-02
53 KSP Residual norm 3.798425553173e-02
54 KSP Residual norm 3.748560407702e-02
55 KSP Residual norm 3.724127504646e-02
56 KSP Residual norm 3.440431176428e-02
57 KSP Residual norm 3.232527531585e-02
58 KSP Residual norm 2.559831144651e-02
59 KSP Residual norm 2.228223272302e-02
60 KSP Residual norm 1.914862438455e-02
0000004 0.1828E-01 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 5.525635412447e-02
1 KSP Residual norm 5.523779213543e-02
2 KSP Residual norm 5.489254753372e-02
3 KSP Residual norm 5.449195079288e-02
4 KSP Residual norm 5.351045474716e-02
5 KSP Residual norm 5.350454984099e-02
6 KSP Residual norm 5.348879006758e-02
7 KSP Residual norm 5.131804406749e-02
8 KSP Residual norm 5.130634818964e-02
9 KSP Residual norm 5.087978009574e-02
10 KSP Residual norm 5.011790353306e-02
11 KSP Residual norm 4.533803591536e-02
12 KSP Residual norm 3.911373694053e-02
13 KSP Residual norm 3.846373237137e-02
14 KSP Residual norm 3.600651859095e-02
15 KSP Residual norm 3.367842601739e-02
16 KSP Residual norm 3.331417578347e-02
17 KSP Residual norm 3.232106678629e-02
18 KSP Residual norm 2.987776906332e-02
19 KSP Residual norm 2.987650978692e-02
20 KSP Residual norm 2.986461389971e-02
21 KSP Residual norm 2.986427312470e-02
22 KSP Residual norm 2.927086107956e-02
23 KSP Residual norm 2.915793209415e-02
24 KSP Residual norm 2.895908986430e-02
25 KSP Residual norm 2.482129709912e-02
26 KSP Residual norm 2.215601003618e-02
27 KSP Residual norm 1.852509765167e-02
28 KSP Residual norm 1.463909533226e-02
29 KSP Residual norm 1.182311494447e-02
30 KSP Residual norm 1.104520437591e-02
31 KSP Residual norm 1.104387877212e-02
32 KSP Residual norm 1.104147164556e-02
33 KSP Residual norm 1.103409135441e-02
34 KSP Residual norm 1.102286965812e-02
35 KSP Residual norm 1.102285643829e-02
36 KSP Residual norm 1.101641432221e-02
37 KSP Residual norm 1.101051494244e-02
38 KSP Residual norm 1.100036449340e-02
39 KSP Residual norm 1.099997781089e-02
40 KSP Residual norm 1.099778407905e-02
41 KSP Residual norm 1.099328811279e-02
42 KSP Residual norm 1.099279011358e-02
43 KSP Residual norm 1.094729805901e-02
44 KSP Residual norm 1.092560978245e-02
45 KSP Residual norm 1.085465993129e-02
46 KSP Residual norm 1.070530862736e-02
47 KSP Residual norm 1.069191452546e-02
48 KSP Residual norm 1.021260176837e-02
49 KSP Residual norm 1.016462435052e-02
50 KSP Residual norm 1.014117683302e-02
51 KSP Residual norm 1.009109769783e-02
52 KSP Residual norm 1.009014269115e-02
53 KSP Residual norm 1.008243622191e-02
54 KSP Residual norm 9.828263449640e-03
55 KSP Residual norm 9.644468912478e-03
56 KSP Residual norm 9.471903670996e-03
57 KSP Residual norm 9.424275132518e-03
58 KSP Residual norm 9.171514955945e-03
59 KSP Residual norm 9.129070363302e-03
60 KSP Residual norm 8.920084490772e-03
61 KSP Residual norm 8.920052881620e-03
62 KSP Residual norm 8.919017185449e-03
63 KSP Residual norm 8.917675663928e-03
64 KSP Residual norm 8.916816143903e-03
65 KSP Residual norm 8.916306367790e-03
66 KSP Residual norm 8.915411548718e-03
67 KSP Residual norm 8.914357204518e-03
68 KSP Residual norm 8.914325973013e-03
69 KSP Residual norm 8.913827907595e-03
70 KSP Residual norm 8.913823013719e-03
71 KSP Residual norm 8.912842938749e-03
72 KSP Residual norm 8.912759077344e-03
73 KSP Residual norm 8.908363036649e-03
74 KSP Residual norm 8.907594393040e-03
75 KSP Residual norm 8.903524069895e-03
76 KSP Residual norm 8.903262538003e-03
77 KSP Residual norm 8.901039140851e-03
78 KSP Residual norm 8.899996065837e-03
79 KSP Residual norm 8.899969808238e-03
80 KSP Residual norm 8.897816058103e-03
81 KSP Residual norm 8.891773515959e-03
82 KSP Residual norm 8.889867086082e-03
83 KSP Residual norm 8.889531100993e-03
84 KSP Residual norm 8.845142763665e-03
85 KSP Residual norm 8.812084759770e-03
86 KSP Residual norm 8.783032596521e-03
87 KSP Residual norm 8.707250471595e-03
88 KSP Residual norm 8.702180431461e-03
89 KSP Residual norm 8.635571960351e-03
90 KSP Residual norm 7.518785765931e-03
91 KSP Residual norm 7.517335145455e-03
92 KSP Residual norm 7.514172655556e-03
93 KSP Residual norm 7.506088977377e-03
94 KSP Residual norm 7.482746241802e-03
95 KSP Residual norm 7.475422608205e-03
96 KSP Residual norm 7.468838236535e-03
97 KSP Residual norm 7.455967440625e-03
98 KSP Residual norm 7.451468344334e-03
99 KSP Residual norm 7.444672451339e-03
100 KSP Residual norm 7.427661463675e-03
101 KSP Residual norm 7.427320924670e-03
102 KSP Residual norm 7.373404855544e-03
103 KSP Residual norm 7.367870254877e-03
104 KSP Residual norm 7.365859749509e-03
105 KSP Residual norm 7.244783986628e-03
106 KSP Residual norm 7.242768509285e-03
107 KSP Residual norm 7.242168100896e-03
108 KSP Residual norm 7.209856849080e-03
109 KSP Residual norm 7.193612184307e-03
110 KSP Residual norm 7.192680079161e-03
111 KSP Residual norm 7.191745723199e-03
112 KSP Residual norm 7.132789915733e-03
113 KSP Residual norm 7.102962868675e-03
114 KSP Residual norm 7.093393180326e-03
115 KSP Residual norm 7.050086753347e-03
116 KSP Residual norm 7.005778938685e-03
117 KSP Residual norm 6.923125763422e-03
118 KSP Residual norm 6.845042283113e-03
119 KSP Residual norm 6.782618558184e-03
120 KSP Residual norm 6.725143356837e-03
121 KSP Residual norm 6.725132622133e-03
122 KSP Residual norm 6.724847388119e-03
123 KSP Residual norm 6.724591741844e-03
124 KSP Residual norm 6.724521076337e-03
125 KSP Residual norm 6.724112824504e-03
126 KSP Residual norm 6.724017992818e-03
127 KSP Residual norm 6.723682776762e-03
128 KSP Residual norm 6.722114858496e-03
129 KSP Residual norm 6.720184071033e-03
130 KSP Residual norm 6.719812160500e-03
131 KSP Residual norm 6.719677445143e-03
132 KSP Residual norm 6.719664963186e-03
133 KSP Residual norm 6.715718760541e-03
134 KSP Residual norm 6.712603576437e-03
135 KSP Residual norm 6.712442649119e-03
136 KSP Residual norm 6.709544343175e-03
137 KSP Residual norm 6.708458187204e-03
138 KSP Residual norm 6.708430822385e-03
139 KSP Residual norm 6.708391457716e-03
140 KSP Residual norm 6.708390792235e-03
141 KSP Residual norm 6.707190601581e-03
142 KSP Residual norm 6.703224755348e-03
143 KSP Residual norm 6.703109804513e-03
144 KSP Residual norm 6.702996838906e-03
145 KSP Residual norm 6.699760767242e-03
146 KSP Residual norm 6.692176493433e-03
147 KSP Residual norm 6.662860649181e-03
148 KSP Residual norm 6.628546635694e-03
149 KSP Residual norm 6.530604292402e-03
150 KSP Residual norm 6.161085281341e-03
151 KSP Residual norm 6.160895992501e-03
152 KSP Residual norm 6.160070254099e-03
153 KSP Residual norm 6.159595725727e-03
154 KSP Residual norm 6.157084924466e-03
155 KSP Residual norm 6.156718079604e-03
156 KSP Residual norm 6.155458765950e-03
157 KSP Residual norm 6.155005765622e-03
158 KSP Residual norm 6.154685720229e-03
159 KSP Residual norm 6.149195573032e-03
160 KSP Residual norm 6.147134547485e-03
161 KSP Residual norm 6.147120074653e-03
162 KSP Residual norm 6.146674768574e-03
163 KSP Residual norm 6.122234963131e-03
164 KSP Residual norm 6.110786892227e-03
165 KSP Residual norm 6.109736500537e-03
166 KSP Residual norm 6.109732765258e-03
167 KSP Residual norm 6.108919009230e-03
168 KSP Residual norm 6.099997144503e-03
169 KSP Residual norm 6.089822903290e-03
170 KSP Residual norm 6.089040042350e-03
171 KSP Residual norm 6.086209099035e-03
172 KSP Residual norm 6.079717411475e-03
173 KSP Residual norm 6.068765961475e-03
174 KSP Residual norm 6.062916731513e-03
175 KSP Residual norm 6.054259223203e-03
176 KSP Residual norm 6.053401430239e-03
177 KSP Residual norm 6.053138364097e-03
178 KSP Residual norm 6.003231859964e-03
179 KSP Residual norm 5.867241012253e-03
180 KSP Residual norm 5.222815812010e-03
0000005 0.4602E-02 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 2.181258502939e-02
1 KSP Residual norm 2.139160683993e-02
2 KSP Residual norm 2.111132578997e-02
3 KSP Residual norm 2.005041182942e-02
4 KSP Residual norm 1.752877255518e-02
5 KSP Residual norm 1.732245712837e-02
6 KSP Residual norm 1.579101376259e-02
7 KSP Residual norm 1.244042616193e-02
8 KSP Residual norm 1.239980884602e-02
9 KSP Residual norm 1.210096906811e-02
10 KSP Residual norm 1.184653943124e-02
11 KSP Residual norm 1.176320516373e-02
12 KSP Residual norm 1.170782473700e-02
13 KSP Residual norm 1.045298473249e-02
14 KSP Residual norm 1.039361338460e-02
15 KSP Residual norm 9.032884384939e-03
16 KSP Residual norm 8.447636755442e-03
17 KSP Residual norm 8.391101877814e-03
18 KSP Residual norm 8.283099686897e-03
-------------- next part --------------
***************************************************
CONTROL SETTINGS
***************************************************
LREAD,LWRITE,LPOST,LTEST,LOUTS,LOUTE,LTIME,LGRAD
F F F F F F F F
IMON, JMON, KMON, MMON, RMON, IPR, JPR, KPR, MPR,NPCOR,NIGRAD
2 2 2 1 0 2 2 3 1 1 1
SORMAX, SLARGE, ALFA
0.1000E-07 0.1000E+31 0.9200E+00
(URF(I),I=1,6)
0.9000E+00 0.9000E+00 0.9000E+00 0.1000E+00 0.1000E+01 0.1000E+01
(SOR(I),I=1,6)
0.1000E+00 0.1000E+00 0.1000E+00 0.1000E-01 0.1000E+00 0.1000E+00
(GDS(I),I=1,6) - BLENDING (CDS-UDS)
0.1000E+01 0.1000E+01 0.1000E+01 0.1600E+01 0.1000E+01 0.0000E+00
LSG
10000
***************************************************
START COUPLED ALGORITHM
***************************************************
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 7.680650991441e+01
1 KSP Residual norm 2.315112892713e+00
2 KSP Residual norm 7.080140933135e-01
KSP Object:(coupledsolve_) 1 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=10000, initial guess is zero
tolerances: relative=1e-90, absolute=1.20082, divergence=10000
right preconditioning
has attached null space
using UNPRECONDITIONED norm type for convergence test
PC Object:(coupledsolve_) 1 MPI processes
type: gamg
MG: type is MULTIPLICATIVE, levels=3 cycles=v
Cycles per PCApply=1
Using Galerkin computed coarse grid matrices
Coarse grid solver -- level -------------------------------
KSP Object: (coupledsolve_mg_coarse_) 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: (coupledsolve_mg_coarse_) 1 MPI processes
type: bjacobi
block Jacobi: number of blocks = 1
Local solve is same for all blocks, in the following KSP and PC objects:
KSP Object: (coupledsolve_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: (coupledsolve_mg_coarse_sub_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=56, cols=56
total: nonzeros=1868, allocated nonzeros=1868
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=56, cols=56
total: nonzeros=1868, allocated nonzeros=1868
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Down solver (pre-smoother) on level 1 -------------------------------
KSP Object: (coupledsolve_mg_levels_1_) 1 MPI processes
type: richardson
Richardson: damping factor=1
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: (coupledsolve_mg_levels_1_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=7783, cols=7783
total: nonzeros=627755, allocated nonzeros=627755
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 2 -------------------------------
KSP Object: (coupledsolve_mg_levels_2_) 1 MPI processes
type: richardson
Richardson: damping factor=1
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: (coupledsolve_mg_levels_2_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=256000, cols=256000
total: nonzeros=2.35725e+06, allocated nonzeros=2.35725e+06
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Up solver (post-smoother) same as down solver (pre-smoother)
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=256000, cols=256000
total: nonzeros=2.35725e+06, allocated nonzeros=2.35725e+06
total number of mallocs used during MatSetValues calls =0
not using I-node routines
0000001 0.1000E+01 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 2.537524950128e+00
1 KSP Residual norm 4.800084858727e-01
2 KSP Residual norm 1.620821084959e-01
0000002 0.2113E+00 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 4.325900644323e-01
1 KSP Residual norm 1.317083025844e-01
2 KSP Residual norm 6.484567311390e-02
3 KSP Residual norm 3.810677528327e-02
0000003 0.3602E-01 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 1.744152280000e-01
1 KSP Residual norm 3.901664369762e-02
2 KSP Residual norm 2.308489607380e-02
3 KSP Residual norm 1.282790886819e-02
0000004 0.1452E-01 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 3.792432575935e-02
1 KSP Residual norm 1.446893562868e-02
2 KSP Residual norm 1.170414412325e-02
3 KSP Residual norm 7.627017584906e-03
4 KSP Residual norm 4.004268384306e-03
5 KSP Residual norm 2.951612023714e-03
0000005 0.3158E-02 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 1.750225257756e-02
1 KSP Residual norm 3.287521188931e-03
2 KSP Residual norm 2.090054722070e-03
3 KSP Residual norm 1.313025221380e-03
0000006 0.1458E-02 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 4.304131293348e-03
1 KSP Residual norm 1.460779398297e-03
2 KSP Residual norm 1.094883111455e-03
3 KSP Residual norm 6.711726958470e-04
4 KSP Residual norm 4.235393183993e-04
0000007 0.3584E-03 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 1.934132941245e-03
1 KSP Residual norm 4.920885681151e-04
2 KSP Residual norm 3.608655529099e-04
3 KSP Residual norm 2.041208433378e-04
4 KSP Residual norm 1.516104297901e-04
0000008 0.1611E-03 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 5.204045108286e-04
1 KSP Residual norm 1.799866485727e-04
2 KSP Residual norm 1.337137496960e-04
3 KSP Residual norm 7.406048243690e-05
4 KSP Residual norm 5.741608689253e-05
5 KSP Residual norm 4.886315540918e-05
0000009 0.4334E-04 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 2.341009109387e-04
1 KSP Residual norm 4.097862927298e-05
2 KSP Residual norm 3.101636775409e-05
3 KSP Residual norm 2.303123486841e-05
0000010 0.1950E-04 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 6.980043269989e-05
1 KSP Residual norm 2.154354863441e-05
2 KSP Residual norm 1.515768659879e-05
3 KSP Residual norm 1.080676666607e-05
4 KSP Residual norm 9.930405731118e-06
5 KSP Residual norm 7.757848714880e-06
6 KSP Residual norm 4.165862290699e-06
0000011 0.5813E-05 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 2.922247023288e-05
1 KSP Residual norm 5.326163951339e-06
2 KSP Residual norm 2.915029700326e-06
0000012 0.2434E-05 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 9.202780292901e-06
1 KSP Residual norm 2.566838156317e-06
2 KSP Residual norm 1.748670689745e-06
3 KSP Residual norm 1.252703127678e-06
4 KSP Residual norm 1.129468798804e-06
5 KSP Residual norm 9.130304621369e-07
0000013 0.7664E-06 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 3.907612884416e-06
1 KSP Residual norm 8.982421498294e-07
2 KSP Residual norm 6.645851329385e-07
3 KSP Residual norm 5.497724351328e-07
4 KSP Residual norm 3.683792060416e-07
0000014 0.3254E-06 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 1.289190862676e-06
1 KSP Residual norm 4.434982001168e-07
2 KSP Residual norm 3.485444131419e-07
3 KSP Residual norm 2.346302677179e-07
4 KSP Residual norm 1.639403498152e-07
5 KSP Residual norm 1.359191571158e-07
6 KSP Residual norm 9.621188205348e-08
0000015 0.1074E-06 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 5.408086146100e-07
1 KSP Residual norm 9.223981776498e-08
2 KSP Residual norm 6.104650170618e-08
3 KSP Residual norm 4.981085059159e-08
0000016 0.4504E-07 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 1.935234341417e-07
1 KSP Residual norm 4.005921857358e-08
2 KSP Residual norm 3.259196479070e-08
3 KSP Residual norm 2.649660480243e-08
4 KSP Residual norm 2.335356743530e-08
5 KSP Residual norm 1.711730147849e-08
0000017 0.1612E-07 0.0000E+00
Residual norms for coupledsolve_ solve.
0 KSP Residual norm 8.100595588300e-08
1 KSP Residual norm 1.486341216211e-08
2 KSP Residual norm 1.322167068956e-08
3 KSP Residual norm 1.133936064325e-08
4 KSP Residual norm 6.626715597309e-09
0000018 0.6746E-08 0.0000E+00
TIME FOR CALCULATION: 0.7428E+02
L2-NORM ERROR U VELOCITY 4.6062628669300923E-004
L2-NORM ERROR V VELOCITY 4.4771789082304948E-004
L2-NORM ERROR W VELOCITY 4.7430150765984163E-004
L2-NORM ERROR ABS. VELOCITY 5.2255204191562464E-004
L2-NORM ERROR PRESSURE 1.0967291493601079E-002
*** CALCULATION FINISHED - SEE RESULTS ***
More information about the petsc-users
mailing list