Ax=b Gauss elimination or?

  • Thread starter Milentije
  • Start date
  • #1
48
0

Main Question or Discussion Point

I have
Ax=b problem
where A 1484x1484 matrix,b 1484x1.
A is sparse(95% zeros) but if I go for LU decomposition what should be the next step?
Or is there any other method,I forgot algebra,learned it when I was undergrad long time ago.
 

Answers and Replies

  • #2
hilbert2
Science Advisor
Insights Author
Gold Member
1,409
445
Simple Gauss elimination is mainly a pedagogical tool used in teaching linear systems, and it's far from being the most effective method (especially when talking about specialized problems like sparse systems). If you're writing a code that solves your linear system, you can find free open source libraries/software packages for that. Just Google "sparse linear system open source".
 
  • #3
612
23
I have
Ax=b problem
where A 1484x1484 matrix,b 1484x1.
A is sparse(95% zeros) but if I go for LU decomposition what should be the next step?
Or is there any other method,I forgot algebra,learned it when I was undergrad long time ago.
1484? :bugeye:

May I ask what this is for, out of curiosity?
 
  • #4
SteamKing
Staff Emeritus
Science Advisor
Homework Helper
12,796
1,667
The best method will depend on several factors. Is A banded and symmetric? If is narrow banded and symmetric, Gauss, Choleski decomp, or other methods may be used. If A is large, sparse, and the zero entries are somewhat randomly distributed, with no symmetry or banded layout of the non-zero terms, then an iterative method might be more suitable.
 
  • #5
48
0
Yes,I am creating input file where velocities need to be calculated for every node.Total number of velocities in model is 1484,quite big.
Regarding software,I have problem to install SUPER Lu from LBNL,are there any links for simple code that to not require libraries(like BLAS( IN THIS CASE?
 
  • #6
SteamKing
Staff Emeritus
Science Advisor
Homework Helper
12,796
1,667
I would make my first stop at the netlib: http://www.netlib.org/

Depending on what method you use, you can always google (or dogpile) and find scads of code.
 
  • Like
Likes 1 person

Related Threads on Ax=b Gauss elimination or?

Replies
1
Views
2K
Replies
19
Views
404
  • Last Post
Replies
3
Views
2K
  • Last Post
Replies
1
Views
3K
  • Last Post
Replies
1
Views
4K
  • Last Post
Replies
4
Views
2K
Replies
3
Views
907
Replies
3
Views
2K
  • Last Post
Replies
2
Views
3K
Replies
1
Views
6K
Top