<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Feb 8, 2018 at 4:43 AM, Marco Cisternino <span dir="ltr"><<a href="mailto:marco.cisternino@optimad.it" target="_blank">marco.cisternino@optimad.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi everybody,<br>
I would like to ask what solution is computed if I try to solve the linear system relative to the problem in subject without creating the null space.<br>
I tried with and without the call to<br>
MatNullSpaceCreate(m_<wbr>communicator, PETSC_TRUE, 0, NULL, &nullspace);<br>
and I get zero averaged solution with and the same solution plus a constant without.<br>
How does PETSc  work in the second case?<br>
</blockquote><div><br></div><div>It depends on the Krylov method you use and the initial residual. We do not do anything special.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Does it check the matrix singularity? And is it able to create the null space with the constant automatically?<br>
Thanks.<br>
<br>
<br>
Marco Cisternino, PhD<br>
<a href="mailto:marco.cisternino@optimad.it">marco.cisternino@optimad.it</a><br>
______________________________<wbr>_<br>
OPTIMAD Engineering srl<br>
Via Giacinto Collegno 18, Torino, Italia.<br>
<a href="tel:%2B3901119719782" value="+3901119719782">+3901119719782</a><br>
<a href="http://www.optimad.it" rel="noreferrer" target="_blank">www.optimad.it</a><br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>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</div><div><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>