Physics – Computational Physics
Scientific paper
2006-10-24
Int.J.Mod.Phys.C19:523-548,2008
Physics
Computational Physics
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
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.
Profile ID: LFWR-SCP-O-344411