pfapack: 高效计算稠密和带状斜对称矩阵的Pfaffian

代码和算法来自arXiv:1102.3440，由Michael Wimmer撰写。

安装

推荐方式（因为它包括更快的C/FORTRAN绑定）

conda install pfapack

或者使用

pip install pfapack

用法

from pfapack import pfaffian as pf
import numpy.matlib

A = numpy.matlib.rand(100, 100)
A = A - A.T
pfa1 = pf.pfaffian(A)
pfa2 = pf.pfaffian(A, method="H")
pfa3 = pf.pfaffian_schur(A)

print(pfa1, pfa2, pfa3)

如果使用conda安装，C/FORTRAN代码包含Python绑定，使用方法如下

from pfapack.ctypes import pfaffian as cpf

pfa1 = cpf(A)
pfa2 = cpf(A, method="H")

print(pfa1, pfa2)

引用

如果您在研究中使用了pfapack，请使用以下bib条目进行引用

@article{wimmer2012algorithm,
  title={Efficient numerical computation of the pfaffian for dense and banded skew-symmetric matrices},
  author={Michael Wimmer},
  journal={ACM Transactions on Mathematical Software (TOMS)},
  volume={38},
  number={4},
  pages={1--17},
  year={2012},
  publisher={ACM New York, NY, USA}
}

许可协议

MIT许可

贡献者

• Bas Nijholt
• Michael Wimmer（算法作者）