Faster Tensor Canonicalization. Benjamin E. Niehoff February, 2017.
Faster Tensor Canonicalization [link]Paper  abstract   bibtex   
The Butler-Portugal algorithm for obtaining the canonical form of a tensor expression with respect to slot symmetries and dummy-index renaming suffers, in certain cases with a high degree of symmetry, from O(n!) explosion in both computation time and memory. We present a modified algorithm which alleviates this problem in the most common cases—tensor expressions with subsets of indices which are totally symmetric or totally antisymmetric—in polynomial time. We also present an implementation of the label-renaming mechanism which improves upon that of the original Butler-Portugal algorithm, thus providing a significant speed increase for the average case as well as the highly-symmetric special case. The worst-case behavior remains O(n!), although it occurs in more limited situations unlikely to appear in actual computations. We comment on possible strategies to take if the nature of a computation should make these situations more likely.
@article{benjamin_e._niehoff_faster_2017,
	title = {Faster {Tensor} {Canonicalization}},
	url = {https://arxiv.org/abs/1702.08114#},
	abstract = {The Butler-Portugal algorithm for obtaining the canonical form of a tensor expression with respect to slot symmetries and dummy-index renaming suffers, in certain cases with a high degree of symmetry, from O(n!) explosion in both computation time and memory. We present a modified algorithm which alleviates this problem in the most common cases---tensor expressions with subsets of indices which are totally symmetric or totally antisymmetric---in polynomial time. We also present an implementation of the label-renaming mechanism which improves upon that of the original Butler-Portugal algorithm, thus providing a significant speed increase for the average case as well as the highly-symmetric special case. The worst-case behavior remains O(n!), although it occurs in more limited situations unlikely to appear in actual computations. We comment on possible strategies to take if the nature of a computation should make these situations more likely.},
	author = {{Benjamin E. Niehoff}},
	month = feb,
	year = {2017},
	keywords = {Multilinear algebra, Tensors, mentions sympy},
}

Downloads: 0