[petsc-users] Nullspaces for schur complement PCs

Jed Brown jed at jedbrown.org
Wed Nov 12 20:39:13 CST 2014


Dmitry Karpeyev <karpeev at mcs.anl.gov> writes:

> I think the problem-specific solver configuration can be carried out by an
> additional user-specified callback (attached, e.g., using PCSetUserSetUp())
> called from PCSetUp().
> The problem of propagating it down the hierarchy can be solved by using DMs
> (even DMShell would do, as it carries an ApplicationCtx) and a consistent
> use of DMFieldSplitRestrictHook (by analogy with other restriction hooks,
> like DMSubDomainRestrictHook).  DMFieldSplitRestrict() would be called when
> creating each split and one of the restriction hooks would transfer the
> necessary contexts.
>
> This is a bit heavyweight, but flexible and logically consistent, it seems
> to me, and no worse that what has to be done to propagate things like
> DMSNES/DMKSP through the MG hierarchy. (N)ASM uses something similar, too.

Previously discussed.

http://lists.mcs.anl.gov/pipermail/petsc-dev/2012-August/009291.html

It makes things more dependent on DM, but that may be less bad than the
alternatives.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20141113/837afb41/attachment-0001.pgp>


More information about the petsc-users mailing list