[petsc-users] GAMG issue
John Mousel
john.mousel at gmail.com
Fri Mar 30 10:58:40 CDT 2012
Mark,
I've just run on one core which allows me to get ML to produce a one row
coarse grid. The problem converges. I'm a bit confused.
John
On Fri, Mar 30, 2012 at 10:24 AM, Jed Brown <jedbrown at mcs.anl.gov> wrote:
> -mg_coarse_pc_type svd?
>
> (Use redundant for parallel.)
> On Mar 30, 2012 9:21 AM, "Mark F. Adams" <mark.adams at columbia.edu> wrote:
>
>>
>> On Mar 30, 2012, at 10:52 AM, John Mousel wrote:
>>
>> Mark,
>>
>> I've run GAMG twice with different coarse grid sizes of 2 and 8 with 1
>> sweep of SOR on the coarse grid. For a size of 8 it converges nicely, but
>> for a size of 2, I think the null space is causing too many problems.
>>
>>
>> YEs, the iterative method is seeing the null space because of floating
>> point error.
>>
>> If GAMG were to coarsen to a size of 1, then there would be no hope
>> because only the null space would remain, right? This doesn't ever seem to
>> occur with ML because there are at least as many rows as processors.
>>
>>
>> Yes that seems like a good assumption. The right thing to do here would
>> probably be to do and SVD and filter out the very low modes explicitly.
>> For now I guess tweaking -pc_gamg_coarse_eq_limit n is al that can be
>> done. Not very satisfying. We will think about this ... any thoughts
>> anyone?
>>
>> Mark
>>
>>
>> John
>>
>> On Fri, Mar 30, 2012 at 8:42 AM, Mark F. Adams <mark.adams at columbia.edu>wrote:
>>
>>>
>>> On Mar 29, 2012, at 2:40 PM, John Mousel wrote:
>>>
>>> I'm attempting to solve a non-symmetric discretization of a 3D Poisson
>>> problem. The problem is singular. I've attached the results of KSPView from
>>> runs with ML and GAMG. When I run ML, I get convergence in 30 iterations.
>>> When I attempt to use the same settings with GAMG, I'm not getting
>>> convergence at all. The two things I notice are:
>>>
>>> 1. GAMG is using KSPType preonly, even though I've set it to be
>>> Richardson in my command line options.
>>>
>>>
>>> PETSc seems to switch the coarse grid solver to GMRES in Setup. This
>>> seems to be a bug and I unwisely decide to override this manually. I will
>>> undo this in the next checkin. This should not be the problem however.
>>>
>>> 2. ML only coarsens down to 4 rows while GAMG coarsens to 2. My problem
>>> is singular, and whenever I try to use LU, I get zero pivot problems. To
>>> mitigate this, I've been using Richardson with SOR on the coarse matrix.
>>> Could the smaller coarse grid size of GAMG be causing problems with SOR. If
>>> so, is there a way to put a lower limit on the coarse grid size?
>>>
>>>
>>> I'm thinking that with a 2x2 coarse grid 8 iterations of SOR is picking
>>> up the null space. Maybe try just one SOR iteration on the coarse grid.
>>>
>>> Also, can you run with options_left so that I can see your arguments.
>>> One known bug is the mat_diagaonal_scale breaks GAMG, but it should also
>>> break ML.
>>>
>>> Mark
>>>
>>> John
>>>
>>> On Thu, Mar 29, 2012 at 11:03 AM, Jed Brown <jedbrown at mcs.anl.gov>wrote:
>>>
>>>> On Thu, Mar 29, 2012 at 09:18, John Mousel <john.mousel at gmail.com>wrote:
>>>>
>>>>> [0]PETSC ERROR: Error in external library!
>>>>> [0]PETSC ERROR: Cannot disable floating point exceptions!
>>>>>
>>>>
>>>> Looks like something is strange with your environment because
>>>> fesetenv() is returning an error. I have disabled the call if the trap mode
>>>> is not changing.
>>>>
>>>> http://petsc.cs.iit.edu/petsc/petsc-dev/rev/352b4c19e451
>>>>
>>>
>>> <KSPView_GAMG.txt><KSPView_ML.txt>
>>>
>>>
>>>
>> <KSPView_ML.txt><KSPView_GAMG.txt>
>>
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120330/4d608135/attachment-0001.htm>
-------------- next part --------------
-pres_ksp_type bcgsl -pres_pc_type ml -pres_mg_coarse_ksp_type richardson -pres_mg_coarse_pc_type sor -pres_mg_coarse_pc_sor_its 8 -pres_ksp_monitor_true_residual -pres_ksp_view -pres_pc_ml_Threshold .01 -pres_options_left
Residual norms for pres_ solve.
0 KSP preconditioned resid norm 2.484438808878e+03 true resid norm 1.537048121260e+01 ||r(i)||/||b|| 1.000000000000e+00
2 KSP preconditioned resid norm 4.803121617032e+03 true resid norm 2.730837669491e+02 ||r(i)||/||b|| 1.776676755736e+01
4 KSP preconditioned resid norm 6.585025011351e+02 true resid norm 3.800284284730e+01 ||r(i)||/||b|| 2.472456283031e+00
6 KSP preconditioned resid norm 3.225377737134e+02 true resid norm 1.185212296990e+01 ||r(i)||/||b|| 7.710964156533e-01
8 KSP preconditioned resid norm 4.960746023374e+01 true resid norm 4.118225945036e+00 ||r(i)||/||b|| 2.679308401652e-01
10 KSP preconditioned resid norm 9.030392220992e+00 true resid norm 1.475499217091e+00 ||r(i)||/||b|| 9.599564234079e-02
12 KSP preconditioned resid norm 3.084127960093e+00 true resid norm 1.384958803077e-01 ||r(i)||/||b|| 9.010510366727e-03
14 KSP preconditioned resid norm 4.079306062909e+00 true resid norm 4.367252084412e-01 ||r(i)||/||b|| 2.841324239629e-02
16 KSP preconditioned resid norm 5.156398532433e-01 true resid norm 5.165106549016e-02 ||r(i)||/||b|| 3.360406533519e-03
18 KSP preconditioned resid norm 9.880823156494e-01 true resid norm 4.287333629430e-02 ||r(i)||/||b|| 2.789329475199e-03
20 KSP preconditioned resid norm 2.420500561140e-02 true resid norm 5.178810721229e-03 ||r(i)||/||b|| 3.369322436687e-04
22 KSP preconditioned resid norm 2.874890434196e-02 true resid norm 3.193794305528e-03 ||r(i)||/||b|| 2.077875286631e-04
24 KSP preconditioned resid norm 3.412881480874e-03 true resid norm 1.504951460520e-04 ||r(i)||/||b|| 9.791179857704e-06
26 KSP preconditioned resid norm 4.854602551459e-04 true resid norm 7.508313761033e-05 ||r(i)||/||b|| 4.884891798233e-06
28 KSP preconditioned resid norm 5.821234231049e-06 true resid norm 4.848903629291e-05 ||r(i)||/||b|| 3.154685635552e-06
KSP Object:(pres_) 1 MPI processes
type: bcgsl
BCGSL: Ell = 2
BCGSL: Delta = 0
maximum iterations=1000
tolerances: relative=1e-08, absolute=1e-50, divergence=10000
left preconditioning
has attached null space
using nonzero initial guess
using PRECONDITIONED norm type for convergence test
PC Object:(pres_) 1 MPI processes
type: ml
MG: type is MULTIPLICATIVE, levels=7 cycles=v
Cycles per PCApply=1
Using Galerkin computed coarse grid matrices
Coarse grid solver -- level -------------------------------
KSP Object: (pres_mg_coarse_) 1 MPI processes
type: richardson
Richardson: damping factor=1
maximum iterations=1, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using PRECONDITIONED norm type for convergence test
PC Object: (pres_mg_coarse_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 8, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Matrix Object: 1 MPI processes
type: seqaij
rows=1, cols=1
total: nonzeros=1, allocated nonzeros=1
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Down solver (pre-smoother) on level 1 -------------------------------
KSP Object: (pres_mg_levels_1_) 1 MPI processes
type: richardson
Richardson: damping factor=1
maximum iterations=1
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using PRECONDITIONED norm type for convergence test
PC Object: (pres_mg_levels_1_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Matrix Object: 1 MPI processes
type: seqaij
rows=2, cols=2
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 2 -------------------------------
KSP Object: (pres_mg_levels_2_) 1 MPI processes
type: richardson
Richardson: damping factor=1
maximum iterations=1
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using PRECONDITIONED norm type for convergence test
PC Object: (pres_mg_levels_2_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Matrix Object: 1 MPI processes
type: seqaij
rows=5, cols=5
total: nonzeros=21, allocated nonzeros=21
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 3 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 3 -------------------------------
KSP Object: (pres_mg_levels_3_) 1 MPI processes
type: richardson
Richardson: damping factor=1
maximum iterations=1
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using PRECONDITIONED norm type for convergence test
PC Object: (pres_mg_levels_3_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Matrix Object: 1 MPI processes
type: seqaij
rows=59, cols=59
total: nonzeros=1661, allocated nonzeros=1661
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 4 -------------------------------
KSP Object: (pres_mg_levels_4_) 1 MPI processes
type: richardson
Richardson: damping factor=1
maximum iterations=1
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using PRECONDITIONED norm type for convergence test
PC Object: (pres_mg_levels_4_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Matrix Object: 1 MPI processes
type: seqaij
rows=833, cols=833
total: nonzeros=68657, allocated nonzeros=68657
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 5 -------------------------------
KSP Object: (pres_mg_levels_5_) 1 MPI processes
type: richardson
Richardson: damping factor=1
maximum iterations=1
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using PRECONDITIONED norm type for convergence test
PC Object: (pres_mg_levels_5_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Matrix Object: 1 MPI processes
type: seqaij
rows=16364, cols=16364
total: nonzeros=1501443, allocated nonzeros=1501443
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 6 -------------------------------
KSP Object: (pres_mg_levels_6_) 1 MPI processes
type: richardson
Richardson: damping factor=1
maximum iterations=1
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
has attached null space
using nonzero initial guess
using PRECONDITIONED norm type for convergence test
PC Object: (pres_mg_levels_6_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Matrix Object: 1 MPI processes
type: seqaij
rows=178085, cols=178085
total: nonzeros=4676387, allocated nonzeros=8385165
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:
Matrix Object: 1 MPI processes
type: seqaij
rows=178085, cols=178085
total: nonzeros=4676387, allocated nonzeros=8385165
total number of mallocs used during MatSetValues calls =0
not using I-node routines
total run time = 284.970678000000
WARNING! There are options you set that were not used!
WARNING! could be spelling mistake, etc!
-------------- next part --------------
Options:
-pres_ksp_type bcgsl
-pres_pc_type gamg -pres_pc_gamg_threshold 0.01 -pres_pc_gamg_sym_graph -pres_pc_gamg_coarse_eq_limit 8 -pres_pc_gamg_agg_nsmooths 1
-pres_mg_levels_ksp_type richardson -pres_mg_levels_pc_type sor
-pres_mg_coarse_ksp_type richardson -pres_mg_coarse_pc_type sor -pres_mg_coarse_pc_sor_its 8
-pres_pc_gamg_verbose 2 -pres_ksp_converged_reason -pres_options_left
-pres_ksp_view -pres_ksp_max_it 120 -pres_options_left -pres_ksp_monitor_true_residual
[0]PCSetFromOptions_GAMG threshold set 1.000000e-02
[0]PCSetUp_GAMG level 0 N=178085, n data rows=1, n data cols=1, nnz/row (ave)=26, np=1
[0]PCGAMGFilterGraph 64.0278% nnz after filtering, with threshold 0.01, 26.2593 nnz ave.
[0]maxIndSetAgg removed 0 of 178085 vertices. (0 local) 6889 selected.
[0]PCGAMGProlongator_AGG New grid 6889 nodes
PCGAMGOptprol_AGG smooth P0: max eigen=1.926427e+00 min=4.469907e-02 PC=jacobi
[0]PCSetUp_GAMG 1) N=6889, n data cols=1, nnz/row (ave)=62, 1 active pes
[0]PCGAMGFilterGraph 32.4218% nnz after filtering, with threshold 0.01, 62.979 nnz ave.
[0]maxIndSetAgg removed 0 of 6889 vertices. (0 local) 162 selected.
[0]PCGAMGProlongator_AGG New grid 162 nodes
PCGAMGOptprol_AGG smooth P0: max eigen=1.619990e+00 min=3.528137e-03 PC=jacobi
[0]PCSetUp_GAMG 2) N=162, n data cols=1, nnz/row (ave)=55, 1 active pes
[0]PCGAMGFilterGraph 28.2109% nnz after filtering, with threshold 0.01, 55.2716 nnz ave.
[0]maxIndSetAgg removed 0 of 162 vertices. (0 local) 5 selected.
[0]PCGAMGProlongator_AGG New grid 5 nodes
PCGAMGOptprol_AGG smooth P0: max eigen=1.569842e+00 min=2.465992e-03 PC=jacobi
[0]PCSetUp_GAMG 3) N=5, n data cols=1, nnz/row (ave)=5, 1 active pes
[0]PCGAMGFilterGraph 72% nnz after filtering, with threshold 0.01, 5 nnz ave.
[0]maxIndSetAgg removed 0 of 5 vertices. (0 local) 1 selected.
[0]PCGAMGProlongator_AGG New grid 1 nodes
PCGAMGOptprol_AGG smooth P0: max eigen=3.686042e+00 min=1.557703e-13 PC=jacobi
[0]PCSetUp_GAMG 4) N=1, n data cols=1, nnz/row (ave)=1, 1 active pes
[0]PCSetUp_GAMG 5 levels, grid complexity = 1.0947
Residual norms for pres_ solve.
0 KSP preconditioned resid norm 7.588646719755e+04 true resid norm 1.537048121260e+01 ||r(i)||/||b|| 1.000000000000e+00
2 KSP preconditioned resid norm 1.326222725745e+03 true resid norm 1.669373471357e+01 ||r(i)||/||b|| 1.086090570794e+00
4 KSP preconditioned resid norm 8.694481142514e+02 true resid norm 1.002676416276e+01 ||r(i)||/||b|| 6.523389882247e-01
6 KSP preconditioned resid norm 9.189263949949e+01 true resid norm 6.222035822438e+00 ||r(i)||/||b|| 4.048042306793e-01
8 KSP preconditioned resid norm 3.703995815828e+01 true resid norm 5.983853333761e+00 ||r(i)||/||b|| 3.893081323215e-01
10 KSP preconditioned resid norm 2.183186527735e+01 true resid norm 5.952245622356e+00 ||r(i)||/||b|| 3.872517418307e-01
12 KSP preconditioned resid norm 3.493925445567e+02 true resid norm 7.329828465587e+00 ||r(i)||/||b|| 4.768769672336e-01
14 KSP preconditioned resid norm 1.356658375642e+03 true resid norm 1.739768043232e+01 ||r(i)||/||b|| 1.131889118609e+00
16 KSP preconditioned resid norm 5.412294172351e+02 true resid norm 6.711159067106e+00 ||r(i)||/||b|| 4.366264773548e-01
18 KSP preconditioned resid norm 4.616224894440e+02 true resid norm 7.521105420582e+00 ||r(i)||/||b|| 4.893214022744e-01
20 KSP preconditioned resid norm 2.577285776652e+02 true resid norm 6.858357846531e+00 ||r(i)||/||b|| 4.462031963520e-01
22 KSP preconditioned resid norm 1.231090948563e+02 true resid norm 6.138029315352e+00 ||r(i)||/||b|| 3.993387864994e-01
24 KSP preconditioned resid norm 2.852537299632e+03 true resid norm 1.740135209172e+01 ||r(i)||/||b|| 1.132127995931e+00
26 KSP preconditioned resid norm 2.236778478113e+04 true resid norm 2.346805733410e+02 ||r(i)||/||b|| 1.526826454520e+01
28 KSP preconditioned resid norm 1.869320270064e+04 true resid norm 1.805684545024e+02 ||r(i)||/||b|| 1.174774244247e+01
30 KSP preconditioned resid norm 2.308247715341e+04 true resid norm 2.257355295321e+02 ||r(i)||/||b|| 1.468630203633e+01
32 KSP preconditioned resid norm 2.261854214989e+04 true resid norm 2.331192380366e+02 ||r(i)||/||b|| 1.516668442661e+01
34 KSP preconditioned resid norm 2.497512693081e+04 true resid norm 2.734577063670e+02 ||r(i)||/||b|| 1.779109597055e+01
36 KSP preconditioned resid norm 2.596144404138e+04 true resid norm 3.014136843412e+02 ||r(i)||/||b|| 1.960990551773e+01
38 KSP preconditioned resid norm 2.350126041193e+04 true resid norm 2.252989797587e+02 ||r(i)||/||b|| 1.465790020771e+01
40 KSP preconditioned resid norm 2.477262261324e+04 true resid norm 2.564872285424e+02 ||r(i)||/||b|| 1.668700055611e+01
42 KSP preconditioned resid norm 6.931302940514e+04 true resid norm 8.605328793761e+02 ||r(i)||/||b|| 5.598607275032e+01
44 KSP preconditioned resid norm 8.688325842110e+04 true resid norm 1.094833786804e+03 ||r(i)||/||b|| 7.122963631788e+01
46 KSP preconditioned resid norm 4.479201456138e+04 true resid norm 2.694563200163e+02 ||r(i)||/||b|| 1.753076668773e+01
48 KSP preconditioned resid norm 2.165097859297e+04 true resid norm 3.074224494605e+02 ||r(i)||/||b|| 2.000083440514e+01
50 KSP preconditioned resid norm 2.055902332706e+04 true resid norm 2.890284413681e+02 ||r(i)||/||b|| 1.880412443633e+01
52 KSP preconditioned resid norm 2.226003271579e+04 true resid norm 3.190481467683e+02 ||r(i)||/||b|| 2.075719961888e+01
54 KSP preconditioned resid norm 2.649136916398e+04 true resid norm 3.633403292566e+02 ||r(i)||/||b|| 2.363883890367e+01
56 KSP preconditioned resid norm 3.636407292301e+03 true resid norm 1.770809887959e+02 ||r(i)||/||b|| 1.152084871947e+01
58 KSP preconditioned resid norm 1.911384813961e+03 true resid norm 1.764445605915e+02 ||r(i)||/||b|| 1.147944284574e+01
60 KSP preconditioned resid norm 3.026919479329e+02 true resid norm 1.794750740713e+02 ||r(i)||/||b|| 1.167660736114e+01
62 KSP preconditioned resid norm 3.308470165040e+02 true resid norm 1.771235549660e+02 ||r(i)||/||b|| 1.152361806479e+01
64 KSP preconditioned resid norm 1.051625068085e+02 true resid norm 1.775402024438e+02 ||r(i)||/||b|| 1.155072505461e+01
66 KSP preconditioned resid norm 2.235446953063e+01 true resid norm 1.778283859610e+02 ||r(i)||/||b|| 1.156947420847e+01
68 KSP preconditioned resid norm 2.112406298937e+01 true resid norm 1.777642739041e+02 ||r(i)||/||b|| 1.156530309268e+01
70 KSP preconditioned resid norm 2.326066510932e+01 true resid norm 1.777623800257e+02 ||r(i)||/||b|| 1.156517987738e+01
72 KSP preconditioned resid norm 3.398049334914e+00 true resid norm 1.777679162542e+02 ||r(i)||/||b|| 1.156554006315e+01
74 KSP preconditioned resid norm 5.201069185149e+00 true resid norm 1.777678703221e+02 ||r(i)||/||b|| 1.156553707482e+01
76 KSP preconditioned resid norm 5.533064413464e+00 true resid norm 1.777679308795e+02 ||r(i)||/||b|| 1.156554101466e+01
78 KSP preconditioned resid norm 3.913801853323e+00 true resid norm 1.777680240223e+02 ||r(i)||/||b|| 1.156554707452e+01
80 KSP preconditioned resid norm 2.537859018317e+00 true resid norm 1.777681223249e+02 ||r(i)||/||b|| 1.156555347006e+01
82 KSP preconditioned resid norm 7.741748315959e+00 true resid norm 1.777654007295e+02 ||r(i)||/||b|| 1.156537640369e+01
84 KSP preconditioned resid norm 4.440928932749e+00 true resid norm 1.777677225227e+02 ||r(i)||/||b|| 1.156552745902e+01
86 KSP preconditioned resid norm 3.704557177890e+00 true resid norm 1.777669720986e+02 ||r(i)||/||b|| 1.156547863660e+01
88 KSP preconditioned resid norm 2.046666429633e+00 true resid norm 1.777675304674e+02 ||r(i)||/||b|| 1.156551496395e+01
90 KSP preconditioned resid norm 2.774055639762e+00 true resid norm 1.777673010008e+02 ||r(i)||/||b|| 1.156550003490e+01
92 KSP preconditioned resid norm 2.777963966142e+00 true resid norm 1.777673492779e+02 ||r(i)||/||b|| 1.156550317580e+01
94 KSP preconditioned resid norm 3.047451777681e-01 true resid norm 1.777673215532e+02 ||r(i)||/||b|| 1.156550137203e+01
96 KSP preconditioned resid norm 1.391897161136e-01 true resid norm 1.777673352268e+02 ||r(i)||/||b|| 1.156550226164e+01
98 KSP preconditioned resid norm 4.877455333649e-02 true resid norm 1.777673302262e+02 ||r(i)||/||b|| 1.156550193630e+01
100 KSP preconditioned resid norm 2.319094057406e-02 true resid norm 1.777673333910e+02 ||r(i)||/||b|| 1.156550214220e+01
102 KSP preconditioned resid norm 1.902403543163e-01 true resid norm 1.777673373724e+02 ||r(i)||/||b|| 1.156550240123e+01
104 KSP preconditioned resid norm 2.088352126122e-02 true resid norm 1.777673337586e+02 ||r(i)||/||b|| 1.156550216612e+01
106 KSP preconditioned resid norm 1.557333341580e-03 true resid norm 1.777673349772e+02 ||r(i)||/||b|| 1.156550224540e+01
108 KSP preconditioned resid norm 1.505375301694e-03 true resid norm 1.777673349419e+02 ||r(i)||/||b|| 1.156550224310e+01
110 KSP preconditioned resid norm 1.671356769285e-03 true resid norm 1.777673349043e+02 ||r(i)||/||b|| 1.156550224066e+01
112 KSP preconditioned resid norm 4.450900294941e-04 true resid norm 1.777673343536e+02 ||r(i)||/||b|| 1.156550220483e+01
Linear solve converged due to CONVERGED_RTOL iterations 112
KSP Object:(pres_) 1 MPI processes
type: bcgsl
BCGSL: Ell = 2
BCGSL: Delta = 0
maximum iterations=120
tolerances: relative=1e-08, absolute=1e-50, divergence=10000
left preconditioning
has attached null space
using nonzero initial guess
using PRECONDITIONED norm type for convergence test
PC Object:(pres_) 1 MPI processes
type: gamg
MG: type is MULTIPLICATIVE, levels=5 cycles=v
Cycles per PCApply=1
Using Galerkin computed coarse grid matrices
Coarse grid solver -- level -------------------------------
KSP Object: (pres_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: (pres_mg_coarse_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 8, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Matrix Object: 1 MPI processes
type: seqaij
rows=1, cols=1
total: nonzeros=1, allocated nonzeros=1
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Down solver (pre-smoother) on level 1 -------------------------------
KSP Object: (pres_mg_levels_1_) 1 MPI processes
type: richardson
Richardson: damping factor=1
maximum iterations=1
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (pres_mg_levels_1_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Matrix Object: 1 MPI processes
type: seqaij
rows=5, cols=5
total: nonzeros=25, allocated nonzeros=25
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
Up solver (post-smoother) same as down solver (pre-smoother)
Down solver (pre-smoother) on level 2 -------------------------------
KSP Object: (pres_mg_levels_2_) 1 MPI processes
type: richardson
Richardson: damping factor=1
maximum iterations=1
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (pres_mg_levels_2_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Matrix Object: 1 MPI processes
type: seqaij
rows=162, cols=162
total: nonzeros=8954, allocated nonzeros=8954
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 3 -------------------------------
KSP Object: (pres_mg_levels_3_) 1 MPI processes
type: richardson
Richardson: damping factor=1
maximum iterations=1
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (pres_mg_levels_3_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Matrix Object: 1 MPI processes
type: seqaij
rows=6889, cols=6889
total: nonzeros=433862, allocated nonzeros=433862
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 4 -------------------------------
KSP Object: (pres_mg_levels_4_) 1 MPI processes
type: richardson
Richardson: damping factor=1
maximum iterations=1
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
has attached null space
using nonzero initial guess
using NONE norm type for convergence test
PC Object: (pres_mg_levels_4_) 1 MPI processes
type: sor
SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
linear system matrix = precond matrix:
Matrix Object: 1 MPI processes
type: seqaij
rows=178085, cols=178085
total: nonzeros=4676387, allocated nonzeros=8385165
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:
Matrix Object: 1 MPI processes
type: seqaij
rows=178085, cols=178085
total: nonzeros=4676387, allocated nonzeros=8385165
total number of mallocs used during MatSetValues calls =0
not using I-node routines
total run time = 326.601349000000
WARNING! There are options you set that were not used!
WARNING! could be spelling mistake, etc!
Option left: name:-lspaint_options_left no value
Option left: name:-nullity_options_left no value
Option left: name:-pres_options_left no value
More information about the petsc-users
mailing list