[petsc-dev] FieldSplit restrictions

Matthew Knepley knepley at gmail.com
Sat May 15 22:19:09 CDT 2010


I agree with Barry.

By "making command line splits permanent" you mean they are not overridden
by API splits?

  Matt

On Sat, May 15, 2010 at 5:20 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:

>
> On May 15, 2010, at 5:09 PM, Jed Brown wrote:
>
> > Is there a case for putting a field in more than one split?
>
>     I believe absolutely. It is like overlapping Schwarz but not with
> overlapping grid points instead overlapping components.
>
>    Barry
>
> > The current
> > implementation allows this, but it makes it easier to accidentally have
> > more splits than you wanted.  The underlying challenge for Doing the
> > Right Thing is that the type is not normally set until
> > PCSetFromOptions().  Any calls to PCFieldSplitSet* will be ignored prior
> > to this.  The problem occurs when the user specifies
> > -pc_fieldsplit_*_fields AND the code adds splits with
> > PCFieldSplitSet(IS|Fields).  This used to create a bunch of splits which
> > is probably not what anyone asked for, I just pushed a some code that
> > makes command-line split definitions permanent.  This isn't entirely
> > satisfactory because there isn't currently a way to get rid of them, but
> > it's better than having lots of spurious splits.
> >
> > Lisandro, let me know if the current code behaves the way you expect.
> >
> > Jed
>
>


-- 
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/20100515/13c83836/attachment.html>


More information about the petsc-dev mailing list