ffc team mailing list archive

[Jed Brown <jed@xxxxxxxx>]

Anders Logg wrote:
> On Fri, Jun 05, 2009 at 06:00:36PM -0500, Brandon Chabaud wrote:
>> Is there some way to inform the Krylov solver that constant pressure is in the
>> null space of the operator from Dolfin?  The only options that I could find for
>> the KrylovSolver command were setting the solver type and setting the
>> preconditioner type.
> 
> No, that's not included in the interface. You will have to extract the
> PETSc Mat pointer using mat() from PETScMatrix and use PETSc calls.

Almost.  You need to get the KSP object from the solver and call
KSPSetNullSpace().

Jed

Attachment: signature.asc
Description: OpenPGP digital signature