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. On Mar 9, 2011 9:51 PM, "Matthew Knepley" <[email protected]> wrote: Each process calls FieldSplitSetIS() with the indices for THAT FIELD that are owned by that process. Each process call FieldSplitSetIS() n+1 times if you have n+1 fields. Identifying fields with processes is a mistake. Matt On Wed, Mar 9, 2011 at 1:50 PM, Thomas Witkowski < [email protected]> wrote:

As I ...

-- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener