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
_______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
