Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Ax=b Gauss elimination or?

  1. Jul 26, 2013 #1
    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.
     
  2. jcsd
  3. Jul 26, 2013 #2

    hilbert2

    User Avatar
    Science Advisor
    Gold Member

    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".
     
  4. Jul 26, 2013 #3
    1484? :bugeye:

    May I ask what this is for, out of curiosity?
     
  5. Jul 26, 2013 #4

    SteamKing

    User Avatar
    Staff Emeritus
    Science Advisor
    Homework Helper

    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.
     
  6. Jul 26, 2013 #5
    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?
     
  7. Jul 26, 2013 #6

    SteamKing

    User Avatar
    Staff Emeritus
    Science Advisor
    Homework Helper

    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.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: Ax=b Gauss elimination or?
  1. AX=B versus Ax=b (Replies: 2)

Loading...