<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Feb 8, 2018 at 5:29 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">Thanks Matt,<br>
I'm using KSPFGMRES and I'm sorry but what do you mean with initial residual?<br>
I also force a non-zero initial guess.<br>
</blockquote><div><br></div><div>If your initial residual has a component in the null space of the operator, it is likely to stay there.</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">Thanks again<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>
<br>
______________________________<wbr>__________<br>
Da: Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>><br>
Inviato: mercoledì 7 febbraio 2018 18:57:56<br>
A: Marco Cisternino<br>
Cc: petsc-users<br>
Oggetto: Re: [petsc-users] Elliptic operator with Neumann conditions<br>
<br>
On Thu, Feb 8, 2018 at 4:43 AM, Marco Cisternino <<a href="mailto:marco.cisternino@optimad.it">marco.cisternino@optimad.it</a><<wbr>mailto:<a href="mailto:marco.cisternino@optimad.it">marco.cisternino@<wbr>optimad.it</a>>> wrote:<br>
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>
<br>
It depends on the Krylov method you use and the initial residual. We do not do anything special.<br>
<br>
  Thanks,<br>
<br>
     Matt<br>
<br>
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><<wbr>mailto:<a href="mailto:marco.cisternino@optimad.it">marco.cisternino@<wbr>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><tel:%2B3901119719782><br>
<a href="http://www.optimad.it" rel="noreferrer" target="_blank">www.optimad.it</a><<a href="http://www.optimad.it" rel="noreferrer" target="_blank">http://www.<wbr>optimad.it</a>><br>
<br>
<br>
<br>
<br>
--<br>
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<br>
<br>
<a href="https://www.cse.buffalo.edu/~knepley/" rel="noreferrer" target="_blank">https://www.cse.buffalo.edu/~<wbr>knepley/</a><<a href="http://www.caam.rice.edu/~mk51/" rel="noreferrer" target="_blank">http://www.caam.rice.<wbr>edu/~mk51/</a>><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>