command line option "-sub_pc_factor_levels 1" does not work in fortran

[Matthew Knepley]

Works for me

knepley at khan:/PETSc3/petsc/petsc-dev/src/ksp/ksp/examples/tutorials$ ./ex2
-ksp_view -sub_pc_factor_levels 1 -pc_type bjacobi
KSP Object:
  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=0.000138889, absolute=1e-50, divergence=10000
  left preconditioning
  using PRECONDITIONED norm type for convergence test
PC Object:
  type: bjacobi
    block Jacobi: number of blocks = 1
    Local solve is same for all blocks, in the following KSP and PC objects:
    KSP Object:(sub_)
      type: preonly
      maximum iterations=10000, initial guess is zero
      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000
      left preconditioning
      using PRECONDITIONED norm type for convergence test
    PC Object:(sub_)
      type: ilu
        ILU: 1 level of fill
        ILU: factor fill ratio allocated 1
        ILU: tolerance for zero pivot 1e-12
        ILU: using diagonal shift to prevent zero pivot
        ILU: using diagonal shift on blocks to prevent zero pivot
             out-of-place factorization
             matrix ordering: natural
        ILU: factor fill ratio needed 1.336
             Factored matrix follows
            Matrix Object:
              type=seqaij, rows=56, cols=56
              package used to perform factorization: petsc
              total: nonzeros=334, allocated nonzeros=334
                not using I-node routines
      linear system matrix = precond matrix:
      Matrix Object:
        type=seqaij, rows=56, cols=56
        total: nonzeros=250, allocated nonzeros=280
          not using I-node routines
  linear system matrix = precond matrix:
  Matrix Object:
    type=seqaij, rows=56, cols=56
    total: nonzeros=250, allocated nonzeros=280
      not using I-node routines
Norm of error 0.000280658 iterations 4

  Matt

On Tue, Dec 8, 2009 at 4:29 PM, <hxie at umn.edu> wrote:

> Hi,
>
> 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.
>
> Bests,
> Hui
>



-- 
What most experimenters take for granted before they begin their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20091208/62bcdcfe/attachment.htm>