Amesos2 support for Tpetra Block CRS matrices
Created by: cihanuq
I'd like to flag interest for direct solver support for (experimental) block matrices in Tpetra. This would be very useful for geophysical inversion applications where iterative solvers do not work (well), and where each (FEM) mesh node has multiple degrees of freedom. There was a comment in the mailing list which suggested that there is a way to 'unroll' the block matrix and pass the resulting regular matrix to Amesos2. I could not find a way to do so, however.