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 class="gmail_quote"><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>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><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">

Colin<br><br><span class="HOEnZb"><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><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>