Programming libraries for C linear algebra

  • Thread starter odiakkoh
  • Start date
  • #1
8
0
Hi,

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.
 

Answers and Replies

  • #2
1,356
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).
www.nr.com.
 
  • #3
8
0
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!
 
  • #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.

Matrix:
-10.0 -5.0
6.0 3.0

Vector:
0 0.000000
1 1.000000

info = 0

Solution from CLAPACK:
0 1125899906842624.000000
1 -2251799813685248.000000
 
Top