On Wed, Aug 8, 2012 at 8:48 AM, 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">
Hello all,<br>
<br>
When using PCFieldSplitSetIS to define splits within splits, should the IS contain indices in the original matrix or in the next highest split level?<br></blockquote><div><br></div><div>The next highest split level. You might consider using a DM if you have some structure (like discretization on a mesh).</div>
<div>This should handle all the index-wrangling for you. It is new, but intended for just this kind of thing.</div><div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

Also, my matrix is in aij format but one of the sub fields has a block diagonal structure. Is it still possible use block jacobi on this field?<br></blockquote><div><br></div><div>Yes.</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 and all the best<span class="HOEnZb"><font color="#888888"><br>
Colin<br>
<br>
<br>
<br>
<br>
-- <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>