|July 3, 2011, 04:26||
Sparse solver for block-diagonal matrix
Join Date: Aug 2009
Posts: 3Rep Power: 7
I would like to solve a very large (and very sparse) linear system, whose sparsity pattern is reported below:
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.
|July 5, 2011, 06:37||
Join Date: May 2009
Location: St.Petersburg, Russia
Posts: 22Rep Power: 8
See Direct methods
|block, iterative, linear, sparse|
|Thread||Thread Starter||Forum||Replies||Last Post|
|How force solver to update variable with nearly diagonal matrix||olwi||OpenFOAM Running, Solving & CFD||0||February 29, 2008 04:52|
|matrix solver for vector machine||TakeharuM||Main CFD Forum||1||June 5, 2001 16:16|
|complex sparse matrix||Stephane Gasser||Main CFD Forum||1||May 14, 2001 13:15|
|matrix solver about block structured grid||young min, Sohn||Main CFD Forum||1||July 11, 2000 19:35|
|solver for linear system with large sparse matrix||Yangang Bao||Main CFD Forum||1||October 25, 1999 04:22|