Hi, Barry,<br>Thanks a lot.<br>So this is equivalent to using ksp richardson +pclu with 1 iteration on the coarse grid, right?<br><br>Yan<br><br><div class="gmail_quote">On Mon, Dec 14, 2009 at 11:16 PM, Barry Smith <span dir="ltr">&lt;<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>&gt;</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;"><br>
  preonly means apply the preconditioner to the right hand side only, don&#39;t do any iterations. Since LU is the preconditioner in this case preonly means it is using LU as a direct solver.<div><div></div><div class="h5">
<br>
<br>
<br>
On Dec 14, 2009, at 10:09 PM, Ryan Yan wrote:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi All,<br>
I have the following output from -snes_view with the default snes setting., but I don&#39;t know what the KSPPREONLY does on the coarse grid. Can anyone provide some hints on this?<br>
<br>
Assuming:<br>
A_{coarse} * e_{coarse}  =r _{coarse}<br>
<br>
<br>
Thanks a lot,<br>
<br>
Yan<br>
<br>
<br>
------------------------------------------------------------------------------<br>
SNES Object:<br>
  type: ls<br>
    line search variant: SNESLineSearchCubic<br>
    alpha=0.0001, maxstep=1e+08<br>
  maximum iterations=50, maximum function evaluations=10000<br>
  tolerances: relative=1e-08, absolute=1e-50, solution=1e-08<br>
  total number of linear solver iterations=32<br>
  total number of function evaluations=5<br>
  KSP Object:<br>
    type: fgmres<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=1e-05, absolute=1e-50, divergence=10000<br>
    right preconditioning<br>
  PC Object:<br>
    type: mg<br>
      MG: type is FULL, levels=2 cycles=v<br>
    Coarse gride solver -- level 0 presmooths=1 postsmooths=1 -----<br>
      KSP Object:(mg_coarse_)<br>
        type: preonly<br>
        maximum iterations=1, initial guess is zero<br>
        tolerances:  relative=1e-05, absolute=1e-50, divergence=10000<br>
        left preconditioning<br>
      PC Object:(mg_coarse_)<br>
        type: lu<br>
          LU: out-of-place factorization<br>
            matrix ordering: nd<br>
          LU: tolerance for zero pivot 1e-12<br>
          LU: factor fill ratio needed 5.02244<br>
               Factored matrix follows<br>
              Matrix Object:<br>
                type=seqaij, rows=4096, cols=4096<br>
                package used to perform factorization: petsc<br>
                total: nonzeros=401152, allocated nonzeros=401152<br>
                  using I-node routines: found 1024 nodes, limit used is 5<br>
        linear system matrix = precond matrix:<br>
        Matrix Object:<br>
          type=seqaij, rows=4096, cols=4096<br>
          total: nonzeros=79872, allocated nonzeros=79872<br>
            using I-node routines: found 1024 nodes, limit used is 5<br>
    Down solver (pre-smoother) ......<br>
<br>
</blockquote>
<br>
</div></div></blockquote></div><br>