CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

Sparse solver for block-diagonal matrix

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

LinkBack Thread Tools Display Modes
Old   July 3, 2011, 04:26
Default Sparse solver for block-diagonal matrix
New Member
Alberto Cuoci
Join Date: Aug 2009
Posts: 3
Rep Power: 8
alberto_cuoci is on a distinguished road
Hello everybody,

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.
alberto_cuoci is offline   Reply With Quote

Old   July 5, 2011, 06:37
New Member
Join Date: May 2009
Location: St.Petersburg, Russia
Posts: 22
Rep Power: 8
ignat is on a distinguished road
See Direct methods

Iterative methods
ignat is offline   Reply With Quote


block, iterative, linear, sparse

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
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

All times are GMT -4. The time now is 02:35.