Physics – Condensed Matter – Mesoscale and Nanoscale Physics
Scientific paper
2011-02-16
Physics
Condensed Matter
Mesoscale and Nanoscale Physics
New algorithm (Parlett-Reid), overall speed-ups, now also Matlab implementation. Program code available as ancillary files, or
Scientific paper
Computing the Pfaffian of a skew-symmetric matrix is a problem that arises in various fields of physics. Both computing the Pfaffian and a related problem, computing the canonical form of a skew-symmetric matrix under unitary congruence, can be solved easily once the skew-symmetric matrix has been reduced to skew-symmetric tridiagonal form. We develop efficient numerical methods for computing this tridiagonal form based on Gauss transformations, using a skew-symmetric, blocked form of the Parlett-Reid algorithm, or based on unitary transformations, using block Householder transformations and Givens rotations, that are applicable to dense and banded matrices, respectively. We also give a complete and fully optimized implementation of these algorithms in Fortran, and also provide Python, Matlab and Mathematica implementations for convenience. Finally, we apply these methods to compute the topological charge of a class D nanowire, and show numerically the equivalence of definitions based on the Hamiltonian and the scattering matrix.
No associations
LandOfFree
Efficient numerical computation of the Pfaffian for dense and banded skew-symmetric 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 computation of the Pfaffian for dense and banded skew-symmetric matrices, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Efficient numerical computation of the Pfaffian for dense and banded skew-symmetric matrices will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-379643