dolfin team mailing list archive

[Robert Kirby <robert.c.kirby@xxxxxxxxx>]

Hi all,Does DOLFIN support the notion of a "matrix action", what would be an
"Operator" in Trilinos or a "matrix shell" in PETSc?
The basic interface would tell the dimension of the domain and range and
support an "apply" method.  Indexing, whether
getting or setting values, is not allowed.  Of course, GenericMatrix should
implement this interface.

I ask because sometimes it is advantageous to form the action of a matrix
without forming it explicitly, especially for high
order methods or tensor-product domains.  I have some recent results on fast
matrix-free methods for the Bernstein basis,
and would like to try them using DOLFIN meshes and solvers.

Thanks,
Rob