I wasn't sure if I should start a new thread as this is a very related and basic question. Before I got the impression that switching from an exact LU (UMFAPCK) to an ILU would be a good first try when UMFPACK started too much memory for my problem. However, it appears that the SparseILU class does not play nicely with BlockSparseMatrix.

That is correct. But if you want to just build a preconditioner for the entire matrix, then there is also no reason to substructure the matrix into individual blocks -- you might as well work with one matrix and one vector that correspond to all variables at once.

Block matrices/vectors are a tool to build preconditioners that work on specific blocks of a matrix. If that's not your goal, there is no need to use this tool.

Best
 W.

--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 [email protected]
                           www: http://www.math.colostate.edu/~bangerth/

--
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 [email protected].
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/8eda9083-5ffb-c4d7-8e50-b22724af6428%40colostate.edu.

Reply via email to