# Sparse solver for block-diagonal matrix

 Register Blogs Members List Search Today's Posts Mark Forums Read

 July 3, 2011, 04:26 Sparse solver for block-diagonal matrix #1 New Member   Alberto Cuoci Join Date: Aug 2009 Posts: 3 Rep Power: 16 Hello everybody, I would like to solve a very large (and very sparse) linear system, whose sparsity pattern is reported below: Immagine.bmp There are N dense blocks (whose local dimension is NB) along the main diagonal. The off-diagonal blocks are not dense, but only the elements along the local diagonal are present (as reported in the Figure). The sparsity pattern is symmetric, but in terms of values, the matrix is NOT symmetric. Usually the block dimension NB is equal to 50 and the number of blocks N is equal to 100,000. The max number of off-diagonal blocks is 4 (of course per each main block, which means that the maximum total number of off-diagonal blocks is 4*N). I was looking for a good iterative solver to manage and exploit this kind of structure, or the most appropriate technique available for solving this system in the most efficient way. Up to now I used a generic sparse iterative solver for linear systems, without exploiting the particular pattern here described. Do you have any suggestions or comments? Thanks a lot in advance. Alberto

 July 5, 2011, 06:37 #2 New Member   Alexey Join Date: May 2009 Location: St.Petersburg, Russia Posts: 23 Rep Power: 17 See Direct methods http://www.cise.ufl.edu/research/sparse/codes/ Iterative methods http://www-users.cs.umn.edu/~saad/so...KIT/index.html

 Tags block, iterative, linear, sparse