[petsc-dev] Using PCFieldSplitSetIS

[Jed Brown]

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" <knepley at gmail.com> 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 <
Thomas.Witkowski at tu-dresden.de> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20110310/c5ee64fc/attachment.html>