Works for me<br><br>knepley@khan:/PETSc3/petsc/petsc-dev/src/ksp/ksp/examples/tutorials$ ./ex2 -ksp_view -sub_pc_factor_levels 1 -pc_type bjacobi<br>KSP Object:<br> type: gmres<br> GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement<br>
GMRES: happy breakdown tolerance 1e-30<br> maximum iterations=10000, initial guess is zero<br> tolerances: relative=0.000138889, absolute=1e-50, divergence=10000<br> left preconditioning<br> using PRECONDITIONED norm type for convergence test<br>
PC Object:<br> type: bjacobi<br> block Jacobi: number of blocks = 1<br> Local solve is same for all blocks, in the following KSP and PC objects:<br> KSP Object:(sub_)<br> type: preonly<br> maximum iterations=10000, initial guess is zero<br>
tolerances: relative=1e-05, absolute=1e-50, divergence=10000<br> left preconditioning<br> using PRECONDITIONED norm type for convergence test<br> PC Object:(sub_)<br> type: ilu<br> ILU: 1 level of fill<br>
ILU: factor fill ratio allocated 1<br> ILU: tolerance for zero pivot 1e-12<br> ILU: using diagonal shift to prevent zero pivot<br> ILU: using diagonal shift on blocks to prevent zero pivot<br>
out-of-place factorization<br> matrix ordering: natural<br> ILU: factor fill ratio needed 1.336<br> Factored matrix follows<br> Matrix Object:<br> type=seqaij, rows=56, cols=56<br>
package used to perform factorization: petsc<br> total: nonzeros=334, allocated nonzeros=334<br> not using I-node routines<br> linear system matrix = precond matrix:<br> Matrix Object:<br>
type=seqaij, rows=56, cols=56<br> total: nonzeros=250, allocated nonzeros=280<br> not using I-node routines<br> linear system matrix = precond matrix:<br> Matrix Object:<br> type=seqaij, rows=56, cols=56<br>
total: nonzeros=250, allocated nonzeros=280<br> not using I-node routines<br>Norm of error 0.000280658 iterations 4<br><br> Matt<br><br><div class="gmail_quote">On Tue, Dec 8, 2009 at 4:29 PM, <span dir="ltr"><<a href="mailto:hxie@umn.edu">hxie@umn.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi,<br>
<br>
I want to change the ilu levels for block jacobi method. What is the command line option for it? Or should I use "PCFactorSetLevels" in the code? The "-sub_pc_factor_levels 1" does not work in fortran. Thanks.<br>
<br>
Bests,<br>
Hui<br>
</blockquote></div><br><br clear="all"><br>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener<br>