<p>Thomas, I think Matt may have misunderstood what you wanted. PCFieldSplit does not solve multiple splits concurrently. You can emulate the first solve you need by defining a block diagonal system with one block per process. Then solving it with block Jacobi is the same as a direct solve. To expose more concurrency, you would need to create a new PC.</p>

<p><blockquote type="cite">On Mar 9, 2011 9:51 PM, "Matthew Knepley" <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br><br>Each process calls FieldSplitSetIS() with the indices for THAT FIELD that are owned by that process. Each<div>
process call FieldSplitSetIS() n+1 times if you have n+1 fields. Identifying fields with processes is a mistake.</div>
<div><br></div><div>   Matt<p><font color="#500050"><br><br>On Wed, Mar 9, 2011 at 1:50 PM, Thomas Witkowski <<a href="mailto:Thomas.Witkowski@tu-dresden.de">Thomas.Witkowski@tu-dresden.de</a>> wrote:<br>><br>> As I ...</font></p>
-- <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>
</blockquote></p>