Unimodular Equivalence of Integral Simplices
Feihu Liu, Sihao Tao, Guoce Xin
TL;DR
This work addresses the problem of detecting unimodular equivalence between full-dimensional integral simplices by reducing it to unimodular permutation (UP) equivalence of vertex matrices. It introduces the HEM algorithm, a Hermite-normal-form–based approach that uses permuted Hermite normal forms and pattern groups to obtain canonical representatives, achieving average-case quasi-polynomial time and provable polynomial-time performance with failure probability below $2.5\times 10^{-7}$. A dimension-free resolution to a key open problem is provided: two $d$-dimensional simplices are unimodularly equivalent if and only if their $n$-dimensional pyramids are unimodularly equivalent, for all $d$ and $n\ge d$. The results are complemented by an average-case analysis, SNF-based acceleration, and practical considerations that link to Ehrhart theory and related matrix invariants, with potential implications for broader equivalence testing problems including graph isomorphism.
Abstract
Testing the unimodular equivalence of two full-dimensional integral simplices can be reduced to testing unimodular permutation (UP) equivalence of two nonsingular matrices. We conduct a systematic study of UP-equivalence, which leads to the first average-case quasi-polynomial time algorithm, called \texttt{HEM}, for deciding the unimodular equivalence of $d$-dimensional integral simplices, as well as achieving a polynomial-time complexity with a failure probability less than $2.5 \times 10^{-7}$. A key ingredient is the introduction of the \emph{permuted Hermite normal form} and its associated \emph{pattern group}, which streamlines the UP-equivalence test by comparing canonical forms derived from induced coset representatives. We also present an acceleration strategy based on Smith normal forms. As a theoretical by-product, we prove that two full-dimensional integral simplices are unimodularly equivalent if and only if their $n$-dimensional pyramids are unimodularly equivalent. This resolves an open question posed by Abney-McPeek et al.
