On Wed, Aug 8, 2012 at 10:51 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</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">
<div class="im">On Wed, Aug 8, 2012 at 10:23 PM, Colin McAuliffe <span dir="ltr"><<a href="mailto:cjm2176@columbia.edu" target="_blank">cjm2176@columbia.edu</a>></span> wrote:<br></div><div class="gmail_quote"><div class="im">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Thanks all, regarding use of DM in 3.3, is the procedure now to create the fields with PCFieldSplitSetIS and then use DMCreateFieldDecompositionDM to create a new DM based from the new fields and the DM for the original problem?<br>

</blockquote><div><br></div></div><div>1) Use petsc-dev</div><div><br></div><div>2) PCFieldSplitSetIS() is independent. This allows you to define splits however you want, but then recursive gets harder</div><div><br></div>
<div>
3) In 3.3., it uses DMCreateFieldDecompositionDM() to split all fields apart at once</div><div><br></div><div>4) In petsc-dev, it uses DMCreateSubDM() which can split off any combination of fields, which from the command line is something like</div>

<div><br></div><div>     -pc_fieldsplit_0_fields 2,0 -pc_fieldsplit_1_fields 1,3</div></div></blockquote><div><br></div><div>I should have shown recursive:</div><div><br></div><div>  -fieldsplit_0_pc_type fieldsplit</div>
<div><br></div><div>will split 2,0 into two blocks.</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"><div class="gmail_quote">
<div>     Matt</div><div class="im"><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

Colin<br><br><span><font color="#888888">
-- <br>
Colin McAuliffe<br>
PhD Candidate<br>
Columbia University<br>
Department of Civil Engineering and Engineering Mechanics<br>
</font></span></blockquote></div></div><div class="HOEnZb"><div class="h5"><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>
</div></div></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>