<div class="gmail_extra">If you set a null space, you should make the right hand side consistent.</div><div class="gmail_extra"><br></div><div class="gmail_extra">Have you checked the _accuracy_ of the solution in the case where a null space was not used? When the system is singular, the error can be large even with a small residual, but presumably the rest of your algorithm wants a minimum norm solution.<br>
<br><div class="gmail_quote">On Thu, Nov 8, 2012 at 2:48 AM, Thomas Witkowski <span dir="ltr"><<a href="mailto:thomas.witkowski@tu-dresden.de" target="_blank">thomas.witkowski@tu-dresden.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":2mz">Maybe the following information is also useful to find the problem or/and to understand what's going on: The inner solver solves with a Laplace matrix (discretized with free boundaries). So I set the constant nullspace to the KSP object. This inner solver is used inside the Schur complement solver of a PCFIELDSPLIT object. The Schur complement solver (KSPPREONLY and PCSHELL) is also set to have a constant null space. But the vector, on which the PCSHELL is applied, is not orthogonal to the constant null space. Can you help me to understand why this is still the case?<br>

<br>
To go back to the problem of the influence of the KSP monitor to the solution process: When I project out the constant null space before calling the KSPRICHARDSON with PCHYPRE, the monitor has no influence anymore on the solution process.</div>
</blockquote></div><br></div>