> -- can the SparseDirectUMFPACK class be made to work with a > BlockMatrixArray (or some other type of block matrix)?
I bet it can be made to work somehow. That said: - it would someone to do the work. The place to start is to take a look at the existing SparseDirectUMFPACK class and, in particular, the functions that copy matrix elements into the UMFPACK data structures - the situation is a bit awkward since the BlockMatrixArray is not really meant to allow element-by-element access. Rather, it is supposed to just be an operator. As a consequence, getting individual elements out to the direct solver is not something the current interface really supports. > I already have my matrix blocks assembled, I enter them into the > BlockMatrixArray, and I solve my system using SolverGMRES. My problem is > relatively small, so I'd like to test a direct sparse solver as well. Can i > use the blocks that I assembled, or do I need to think of another approach > to do If you already have the blocks in some way or other, you might be luckier putting it all into a BlockMatrix object for which SparseDirectUMFPACK is already implemented if I recall correctly. Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
