[aspect-devel] Direct Solver in Aspect

[Timo Heister]

> I don't know what solver you intend to use, but the SparseDirectUMFPACK
> solver can work with BlockSparseMatrix and BlockVector. It would probably
> not be very difficult to extend it to use BlockMatrixArray, or to do the
> same with the interfaces to other direct solvers. This may turn out to be
> less intrusive and more extendible than rearranging the block structure.

No, I plan to use TrilinosWrappers::SolverDirect that can only deal
with a Trilinos SparseMatrix (not a Trilinos::BlockSparseMatrix).

>> 2. fix all the places that assume the current layout (say T is in block 3)
>
> Can this be done through appropriate modifications/additions to the
> Introspection class, rather than scattering if/elses throughout?

As much as possible yes. But there are many places where we assume
certain block indices (I am going through them right now).


-- 
Timo Heister
http://www.math.clemson.edu/~heister/