Numeric and symbolic evaluation of the Pfaffian of general skew-symmetric matrices

DOI10.1016/J.CPC.2011.04.025MaRDI QIDQ2275564zbMATH OpenOpenAlexFDO

Authors C. González-Ballestero, George F. Bertsch, L. M. Robledo

Publication date 9 August 2011

Published in Computer Physics Communications (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1012.5022

zbMATH Keywords

\texttt{FORTRAN 77}\texttt{FORTRAN 90}Aitken's block diagonalizationCholesky's factorizationHouseholder transformationsLU-decompositionnumerical examplesPfaffianpivoting strategyskew symmetric matricestridiagonalization

Mathematics Subject Classification ID

Matrix exponential and similar functions of matrices (15A16) Pfaffian systems (58A17) Direct numerical methods for linear systems and matrix inversion (65F05)

Abstract: Evaluation of pfaffians arises in a number of physics applications, and for some of them a direct method is preferable to using the determinantal formula. We discuss two methods for the numerical evaluation of pfaffians. The first is tridiagonalization based on Householder transformations. The main advantage of this method is its numerical stability that makes unnecessary the implementation of a pivoting strategy. The second method considered is based on Aitken's block diagonalization formula. It yields to a kind of LU (similar to Cholesky's factorization) decomposition (under congruence) of arbitrary skew-symmetric matrices that is well suited both for the numeric and symbolic evaluations of the pfaffian. Fortran subroutines (FORTRAN 77 and 90) implementing both methods are given. We also provide simple implementations in Python and Mathematica for purpose of testing, or for exploratory studies of methods that make use of pfaffians.

Recommendations

• Algorithm 923: Efficient numerical computation of the Pfaffian for dense and banded skew-symmetric matrices
• A factorization algorithm to compute Pfaffians
• The Faddeev-LeVerrier algorithm and the Pfaffian
• On the Pfaffians and determinants of some skew-centrosymmetric matrices
• Properties of Pfaffians

Cites work

• A Note on the Stable Decomposition of Skew-Symmetric Matrices
• An extended set of FORTRAN basic linear algebra subprograms
• Coherent states of Fermi operators and the path integral.
• scientific article; zbMATH DE number 2094864 (Why is no real title available?)
• scientific article; zbMATH DE number 3243649 (Why is no real title available?)
• scientific article; zbMATH DE number 961607 (Why is no real title available?)
• On the reduction of a symmetric matrix to tridiagonal form
• Pfaffian
• Sequential Reduction to Tridiagonal Form
• Some remarks on the tensorial linearization of general and linearly separable boolean functions
• The Schur complement and its applications
• Unitary Triangularization of a Nonsymmetric Matrix

Cited in

(12)

• On the domain of attraction of a Tracy-Widom law with applications to testing multiple largest roots
• A factorization algorithm to compute Pfaffians
• Off-shell Bethe states and the six-vertex model
• Solution of the Skyrme-Hartree-Fock-Bogolyubov equations in the Cartesian deformed harmonic-oscillator basis. (VIII) \textsc{hfodd} (v2.73y): a new version of the program
• Optimized implementation for calculation and fast-update of Pfaffians installed to the open-source fermionic variational solver mVMC
• On the Pfaffians and determinants of some skew-centrosymmetric matrices
• Algorithm 923: Efficient numerical computation of the Pfaffian for dense and banded skew-symmetric matrices
• Zeilberger's holonomic ansatz for Pfaffians
• Universality in the tripartite information after global quenches: spin flip and semilocal charges
• Thermodynamics and maximal entropy of a monomer-dimer system on a ladder geometry
• The Faddeev-LeVerrier algorithm and the Pfaffian
• Determinant and Pfaffian of sum of skew symmetric matrices

Describes a project that uses

Uses Software

• Pfaffian
• BLAS

This page was built for publication: Numeric and symbolic evaluation of the Pfaffian of general skew-symmetric matrices

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2275564)