On Fri, Jun 22, 2012 at 12:15 PM, Anton Popov <span dir="ltr"><<a href="mailto:popov@uni-mainz.de" target="_blank">popov@uni-mainz.de</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi petsc team,<br>
I have a question about using nested matrices and vectors for coupled solution of Stokes system with filedsplit preconditioner.<br>
<br>
Is there a simple example explaining how I can create nested matrix and vector using<br>
"MatCreateNest" and "VecCreateNest" functions, provided that I already have all constituting blocks correctly assembled?<br>
I could not find these example online, and if I do it straightforwardly according to documentation, petsc complains that some vector blocks are not setup.<br></blockquote><div><br></div><div>You do not need MatNest to use FS. Some of us have a crazy MatNest attachment, so it seems like you do,</div>
<div>but it will work with any matrix.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Do I need to define velocity and pressure splits using "PCFieldSplitSetIS", because if I don't, petsc says "Unhandled case, must have at least two fields, not 1!"?<br></blockquote><div><br></div><div>
If you really have Stokes (a 0 pressure block), use -pc_fieldsplit_detect_saddle_point.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Do I need to explicitly specify the blocksize, if I want to use the filedsplit preconditioner for the velocity block?<br></blockquote><div><br></div><div>No.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

A small example would be extremely helpful<br></blockquote><div><br></div><div>SNES ex62 runs it, but you need my FEM assembly configured to do that. I will post instructions soon.</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 a lot,<br>
<br>
Anton<br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <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>