<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Nov 13, 2014 at 6:24 AM, Massimiliano Leoni <span dir="ltr"><<a href="mailto:leoni.massimiliano1@gmail.com" target="_blank">leoni.massimiliano1@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">In data giovedì 13 novembre 2014 05:59:17, hai scritto:<br>
> On Thu, Nov 13, 2014 at 5:39 AM, Massimiliano Leoni <<br>
> We call Hypre directly to apply P^{-1}.<br>
><br>
> It should be. Send -ksp_view and we can verify.<br>
<br>
It looks like that's what's happening! :D<br>
Here it is<br>
<a href="http://pastebin.com/6J3Rj4x1" target="_blank">http://pastebin.com/6J3Rj4x1</a></blockquote><div><br></div><div>This is not exactly right because here you are just using additive fieldsplit. The mass matrix</div><div>is a good preconditioner for the Schur complement, not the zero matrix. Take a look at these</div><div>slides:</div><div><br></div><div>  <a href="http://www.mcs.anl.gov/petsc/documentation/tutorials/ParisTutorial.pdf">http://www.mcs.anl.gov/petsc/documentation/tutorials/ParisTutorial.pdf</a></div><div><br></div><div>at slide 115, or my talk from SIAM CS&E</div><div><br></div><div>  <a href="http://people.cs.uchicago.edu/~knepley/presentations/SIAMCSE13.pdf">http://people.cs.uchicago.edu/~knepley/presentations/SIAMCSE13.pdf</a></div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><br>
> You should use<br>
>   -pc_fieldsplit_schur_factorization_type diag<br>
Ok, This is not clear to me.<br>
According to what I studied, Schur decomposition is when I formally eliminate<br>
a variable [for Stokes, velocity u] and get a system in the others [sole<br>
pressure], which I then solve by some smart method and substitute back to find<br>
u. This is done on the full Stokes matrix A, not on the preconditioner...<br>
It should be not what I want to do, as I'm facing the whole system at once.<br>
<br>
What does PETSc mean by "Schur" here?<br>
<br>
>    Thanks,<br>
><br>
>       Matt<br>
<br>
Thank you very much!<br>
<span class=""><font color="#888888">Massimiliano<br>
<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">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></div>