Hi Ian, you are correct, that that block matrices are not accepted by our TrilinosWrappers. One thing we typically do is that we use deal.ii solvers, which work fine with those. If you then write a preconditioner, you can use preconditioners from Trilinos for individual blocks. See for example step-32.
On Mon, Apr 23, 2012 at 2:42 PM, Ian Rose <[email protected]> wrote: > Hi all, > I am trying to experiment with a few different solvers with a system matrix > of type TrilinosWrappers::BlockSparseMatrix. As far as I can tell, no > version of the wrappers for the Trilinos solvers are able to take this as an > argument, as TrilinosWrappers::BlockSparseMatrix derives from > BlockMatrixBase rather than from TrilinosWrappers::SparseMatrix. Am I > understanding this correctly? In order to use the Trilinos solvers, do I > need to copy this into a TrilinosWrappers::SparseMatrix object? > Thanks a lot, > Ian Rose > > > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii > -- Timo Heister http://www.math.tamu.edu/~heister/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
