preconditioner for matrix-free linear system

Satish Balay balay at mcs.anl.gov
Tue Apr 14 13:44:51 CDT 2009


On Tue, 14 Apr 2009, Matthew Knepley wrote:

> I think all the PCs you want rely on AIJ storage. People in linear
> algebra don't think much of generic interfaces, and thus do not
> write their preconditioners to respect them (or you would only have
> to implement GetRow()).

If a user provides a MATSHELL with MatGetRow() interface - perhaps we
should be able to provide a generic MatConvert() interface as well [to
MATSHELL] - that gives a AIJ matrix - useable for PMat? [with access
to all PCs]

[To get efficient assembly we might have to call MatGetRow() twice?]

Satish


More information about the petsc-users mailing list