<html>

  <head>

    <meta content="text/html; charset=ISO-8859-1"

      http-equiv="Content-Type">

  </head>

  <body text="#000000" bgcolor="#FFFFFF">

    <font face="Ubuntu">Thank you very much. I will try and let you know.<br>

      <br>

      Michele<br>

      <br>

    </font>

    <div class="moz-cite-prefix">On 05/17/2013 07:01 PM, Jed Brown

      wrote:<br>

    </div>

    <blockquote cite="mid:87ehd5m4gn.fsf@mcs.anl.gov" type="cite">

      <pre wrap="">Michele Rosso <a class="moz-txt-link-rfc2396E" href="mailto:mrosso@uci.edu"><mrosso@uci.edu></a> writes:

</pre>

      <blockquote type="cite">

        <pre wrap="">I noticed that the problem appears even if I use CG with the default 

preconditioner: commenting KSPSetDM() solves the problem.

</pre>

      </blockquote>

      <pre wrap="">

Okay, this issue can't show up if you use SNES, but it's a consequence

of making geometric multigrid work with a pure KSP interface.  You can

either use KSPSetComputeOperators() to put your assembly in a function

(which will also be called on coarse levels if you use -pc_type mg

without Galerkin coarse operators) or you can can provide the Jacobian

using KSPSetOperators() as usual, but also call KSPSetDMActive() so that

the DM is not used for computing/updating the Jacobian.

The logic is cleaner in petsc-3.4 and I think it just does the right

thing in your case.

</pre>

      <blockquote type="cite">

        <pre wrap="">So basically without a proper grid (it seems no grid with an even 

numbers of nodes qualifies) and with my own system matrix, I cannot use 

any type of multigrid

pre-conditioner?

</pre>

      </blockquote>

      <pre wrap="">

You can use all the AMG methods without setting a DM.

</pre>

    </blockquote>

    <br>

  </body>

</html>