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
