] Shell preconditioner within a fieldsplit
Damn - the last part of my email is wrong. You want to set the PCType to "mat".
KSPType preonly is fine
On Mon, 14 Nov 2016 at 07:04, Dave May
mailto:dave.mayhe...@gmail.com>> wrote:
Looks like you want the contents of your mat shell, spe
Damn - the last part of my email is wrong. You want to set the PCType to
"mat". KSPType preonly is fine
On Mon, 14 Nov 2016 at 07:04, Dave May wrote:
> Looks like you want the contents of your mat shell, specifically the op
> Ax, to define the action of the preconditioner.
>
> You need to either
Looks like you want the contents of your mat shell, specifically the op Ax,
to define the action of the preconditioner.
You need to either create a PCShell (rather than a MatShell), and define
the operation called by PCApply(), or keep you current shell but change
"preonly" to "mat".
http://www.m
From the manual page for PCFieldSplitGetSubKSP()
You must call KSPSetUp() before calling PCFieldSplitGetSubKSP().
Please try this and let us know if it still doesn't call your
PressureBlock().
PETSc developers: we should figure out how to have it generate a very
useful error mess
Hello,
What is the proper way to set up a shell preconditioner within a fielsplit? I
have tried it on my own, but do not get the proper behavior. The relevant
portion looks like this:
__
// Global System
KSPSetOperators(ksp, A,