Anyway, what you need is a way for the matrix_mult function (whether it's global or a static member function) to establish what the proper MatrixFree object is that belongs to the matrix object it is called with. Here is a crude approach:
Actually, there is a better way: MatCreateShell allows you to give a "context" pointer. You currently pass NULL in this place. Pass *this here, and then in the static matrix_mult function get back the context using MatShellGetContext.
Best W. _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
