dolfin team mailing list archive

[Johan Jansson <johanjan@xxxxxxxxxxxxxxxx>]

On Wed, Jan 19, 2005 at 10:07:51AM -0600, Anders Logg wrote:
> Good to see that the linear algebra is coming along. Some comments:
> 
> - solve() should be static to enable NewGMRES::solve(), same as
>   GMRES::solve().

Ok.

> 
> - There should be a version of solve() that can handle VirtualMatrix.
>   I think the code will be the same, just call VirtualMatrix::mat().

Yes, I thought about several different solutions to enable
VirtualMatrix. This one is the most efficient I think.

> 
> - We need to add a class NewPreconditioner that can be passed as an
>   argument to NewGMRES, something like
> 

...

Is this really the job of DOLFIN though? PETSc already has an
infrastructure for preconditioners, shouldn't DOLFIN simply be a user
of the linear solver system of PETSc? If something is missing in the
way PETSc provides preconditioning, I think it makes more sense to add
it to PETSc, and then use it in DOLFIN. Or what do you think?

  Johan