Matrix diagonalisation

  • #1
Hi,

I'm looking for a matrix diagonalisation subroutine in fortran. Does anyone know of any online sources for such things? The matrices are small (The largest being 100x100) so LAPACK/BLAS seems like overkill.
 

Answers and Replies

  • #2
AlephZero
Science Advisor
Homework Helper
7,002
293
Why do you think LAPACK is "overkill"? The reason the code may look complicated is because it is carefully written to handle any numbers that are thrown at them, but it doesn't add anything significant to the computing time.

You don't need to understand what every line of the code is for. Just compile the libraries and use them.
 
  • #3
SteamKing
Staff Emeritus
Science Advisor
Homework Helper
12,798
1,670
LAPACK is large because it is designed to provide many different kinds of matrix solutions. You say your largest problem will be the diagonalization of a 100 x 100 matrix. How are these matrices obtained? For example, finite elements typically produce banded matrices, where zeros appear in a large number of the cells, and the remaining elements are symmetric about the main diagonal. Boundary elements typically produce fully populated matrices where zero elements are rare. Least squares regression analyses produce matrices which can be difficult to diagonalize unless you have access to special routines or have a machine capable of doing extra-long floating point precision. What's yer poison?
 

Related Threads on Matrix diagonalisation

  • Last Post
Replies
22
Views
1K
  • Last Post
Replies
1
Views
2K
Replies
2
Views
683
Replies
2
Views
1K
  • Last Post
Replies
4
Views
2K
Replies
5
Views
721
Replies
4
Views
963
  • Last Post
Replies
6
Views
3K
Replies
2
Views
715
Replies
7
Views
5K
Top