I can't seem to find a way to use the BlockwiseDirect preconditioner on TrilinosWrappers::BlockSparseMatrix objects? The signature for BlockwiseDirect only takes the base SparseMatrix object.

Would it be possible to temporarily forget the block structure and recover the underlying SparseMatrix to construct the preconditioner?

There is no underlying SparseMatrix for block matrices. Each block stores its own SparseMatrix object. They know nothing of each other.

But you can of course create a preconditioner that applies a BlockwiseDirect method to each of the diagonal blocks of your block matrix, assuming that each of these blocks is invertible.

Best
 W.

--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/4c50901b-cd70-48ce-9874-ae1ad8aebba4%40colostate.edu.

Reply via email to