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

Programming libraries for C linear algebra

  1. Apr 2, 2007 #1

    Can anyone recommend a freely available library in C or C++ that is able to calculate eigenvalues and eigenvectors of generalized eigenvalue problems. e.g. of the type:

    Au = LBu

    where A and B are matrices, L is the eigenvalue, and u is the eigenvector.

    N.B. the GNU Scientific library was my first thought, but it can only do standard eigenvalue problems, not those of the generalized type.
  2. jcsd
  3. Apr 2, 2007 #2
    blas doesn't have eigenvalue stuff?. BUt you can try CLAPACK and maybe look at lidia. you can tyr to find some at netlib.org.

    Or you can try numerical recipes code(usually for teaching purposes but still usable).
  4. Apr 3, 2007 #3
    CLAPACK does the job

    Thanks a lot! - CLAPACK does the job.

    For anyone else interested in doing matrix-type algebra in C here's some useful CLAPACK links:

    CLAPACK can be downloaded from here: http://www.netlib.org/clapack/
    Documentation can be found here: http://www.netlib.org/lapack/lug/index.html

    Explanations for function parameters are found in the code.

    It's a bit fiddly to install, and link to the libraries (all three of them) bundled in the package. The interface will feel quite archaic to anyone accustomed to the niceities of object orientated languages (such as myself). But it does do the job.

    Oh, and one last hint for anyone trying to use this: take note of the documentation explaining how matrices are stored in memory: particularly in what order array elements correspond to matrix elements. I didn't bother; made an assumption; got it wrong (damn! I had a 50:50 chance!), and spent a few hours wondering why things weren't working quite right!
  5. Dec 27, 2009 #4
    Hey, guys, I installed CLAPACK on my Mac without much trouble. However it failes on solving this simple equation Ax = b, where A is singular, yet the routine dgesv_() didn't detect it. Can you test the case on you machine? Here is the input.

    -10.0 -5.0
    6.0 3.0

    0 0.000000
    1 1.000000

    info = 0

    Solution from CLAPACK:
    0 1125899906842624.000000
    1 -2251799813685248.000000
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook