A faster hafnian formula for complex matrices and its benchmarking on a supercomputer

被引:21
作者
Björklund A. [1 ]
Gupt B. [2 ]
Quesada N. [2 ]
机构
[1] Department of Computer Science, Lund University
[2] Xanadu, Toronto
来源
ACM Journal of Experimental Algorithmics | 2019年 / 24卷 / 01期
关键词
Gaussian boson sampling; Hafnian; Perfect matching;
D O I
10.1145/3325111
中图分类号
学科分类号
摘要
We introduce new and simple algorithms for the calculation of the number of perfect matchings of complex weighted, undirected graphs with and without loops. Our compact formulas for the hafnian and loop hafnian of n × n complex matrices run in O(n32n/2) time, are embarrassingly parallelizable and, to the best of our knowledge, are the fastest exact algorithms to compute these quantities. Despite our highly optimized algorithm, numerical benchmarks on the Titan supercomputer with matrices up to size 56 × 56 indicate that one would require the 288,000 CPUs of this machine for about 6 weeks to compute the hafnian of a 100 × 100 matrix. © 2019 Association for Computing Machinery.
引用
收藏
相关论文
共 44 条
[1]  
Quesada N., Gupt B., Izaac J., Hafnian, (2018)
[2]  
Aaronson S., Arkhipov A., The computational complexity of linear optics, Proceedings of the 43rd Annual ACM Symposium on Theory of Computing, pp. 333-342, (2011)
[3]  
Abramowitz M., Stegun I.A., Handbook of Mathematical Functions: With Formulas, Graphs, and Mathematical Tables., 55, (1972)
[4]  
Anderson E., Bai Z., Bischof C., Susan Blackford L., Demmel J., Dongarra J., Du Croz J., Greenbaum A., Hammarling S., McKenney A., Et al., LAPACK Users' Guide., (1999)
[5]  
Balasubramanian K., Combinatorics and Diagonals of Matrices, (1980)
[6]  
Barvinok A., Polynomial time algorithms to approximate permanents and mixed discriminants within a simply exponential factor, Rand. Struct. Algor., 14, 1, pp. 29-61
[7]  
Barvinok A., Approximating permanents and hafnians, Discrete Analysis, 2, (2016)
[8]  
Barvinok A., Combinatorics and Complexity of Partition Functions., 274, (2016)
[9]  
Barvinok A.I., Two algorithmic results for the traveling salesman problem, Math. Operat. Res., 21, 1, pp. 65-84, (1996)
[10]  
Bax E., Franklin J., A finite-difference sieve to compute the permanent, Caltech-CS-TR-96-04, (1996)