Efficient numerical diagonalization of hermitian 3x3 matrices

Physics – Computational Physics

Scientific paper

Rate now

  [ 0.00 ] – not rated yet Voters 0   Comments 0

Details

13 pages, no figures, new hybrid algorithm added, matches published version, typo in Eq. (39) corrected; software library avai

Scientific paper

10.1142/S0129183108012303

A very common problem in science is the numerical diagonalization of symmetric or hermitian 3x3 matrices. Since standard "black box" packages may be too inefficient if the number of matrices is large, we study several alternatives. We consider optimized implementations of the Jacobi, QL, and Cuppen algorithms and compare them with an analytical method relying on Cardano's formula for the eigenvalues and on vector cross products for the eigenvectors. Jacobi is the most accurate, but also the slowest method, while QL and Cuppen are good general purpose algorithms. The analytical algorithm outperforms the others by more than a factor of 2, but becomes inaccurate or may even fail completely if the matrix entries differ greatly in magnitude. This can mostly be circumvented by using a hybrid method, which falls back to QL if conditions are such that the analytical calculation might become too inaccurate. For all algorithms, we give an overview of the underlying mathematical ideas, and present detailed benchmark results. C and Fortran implementations of our code are available for download from http://www.mpi-hd.mpg.de/~globes/3x3/ .

No associations

LandOfFree

Say what you really think

Search LandOfFree.com for scientists and scientific papers. Rate them and share your experience with other people.

Rating

Efficient numerical diagonalization of hermitian 3x3 matrices does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Efficient numerical diagonalization of hermitian 3x3 matrices, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Efficient numerical diagonalization of hermitian 3x3 matrices will most certainly appreciate the feedback.

Rate now

     

Profile ID: LFWR-SCP-O-344411

  Search
All data on this website is collected from public sources. Our data reflects the most accurate information available at the time of publication.