> I got an idea: you can try to wrap SparseMatrix in a class that
> applies tvmult instead of vmult. Code (not tested, not sure it even
> compiles):
> template<class MAT>
> class TMatrix
> {
> public:
> TMatrix(const MAT &m)
> : mat(m)
> {}
>
...and under the name TransposeMatrix, we could even add it to the
library. But, base it on the implementation of ScaledMatrix, using
MatrixPointerBase and only the vector type as template argument.
> Obviously you don't have access to preconditioners like this...
>
You can use the same TransposeMatrix for your preconditioner. Many of
them have a Tvmult already.
Guido
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
