<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Feb 7, 2018 at 1:45 PM, 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">I'm sorry Matt but I cannot understand what flexible gmres computes when no null space is created.<br>
Could you give me some hints, please? Even in very simple cases...<br>
</blockquote><div><br></div><div>I don't think I can clear it up any further.</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">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>
<br>
______________________________<wbr>__________<br>
Da: Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>><br>
Inviato: mercoledì 7 febbraio 2018 19:38:08<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 5:29 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>
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>
<br>
If your initial residual has a component in the null space of the operator, it is likely to stay there.<br>
<br>
   Matt<br>
<br>
Thanks again<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>
______________________________<wbr>__________<br>
Da: Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a><mailto:<a href="mailto:knepley@gmail.com">knep<wbr>ley@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>><mailto:<a href="mailto:marco.cisternino@optimad.it">marco.<wbr>cisternino@optimad.it</a><mailto:<a href="mailto:marco.cisternino@optimad.it">m<wbr>arco.cisternino@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>><mailto:<a href="mailto:marco.cisternino@optimad.it">marco.<wbr>cisternino@optimad.it</a><mailto:<a href="mailto:marco.cisternino@optimad.it">m<wbr>arco.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><tel:%2B3901119719782><tel:%<wbr>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>><<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>
<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>