dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Thu, Jan 13, 2005 at 04:53:38AM +0100, Johan Jansson wrote:

> By the way, I found out how to use the PETSc Krylov solvers by
> providing only the action of the operator:
> 
> http://www.epcc.ed.ac.uk/computing/services/cray_service/documents/local_libs/PETSC/manual.html#Node26
> 
> I don't think it's meaningful to wrap it, but that's what we need to
> use for the Newton solver in the ODE solver.

I've added a template for PETSc matrix-free matrices (shell) called
VirtualMatrix.

I started to implement the multi-adaptive Newton-solver without the
wrapper, but the number of PETSc-calls just to initialize the matrix
is so large I think a wrapper is needed.

/Anders