On Efficient Sparse Integer Matrix Smith Normal Form Computations

Jean-Guillaume Dumas (1) , B. David Saunders (2) , Gilles Villard (3)

1 ID-IMAG - Informatique et Distribution
2 University of Delaware [Newark]
3 LIP - Laboratoire de l'Informatique du Parallélisme

We present a new algorithm to compute the Integer Smith normal form of large sparse matrices. We reduce the computation of the Smith form to independent, and therefore parallel, computations modulo powers of word-size primes. Consequently, the algorithm does not suffer from coefficient growth. We have implemented several variants of this algorithm (elimination and/or black box techniques) since practical performance depends strongly on the memory available. Our method has proven useful in algebraic topology for the computation of the homology of some large simplicial complexes.

Domaines: Calcul formel [cs.SC]

Jean-Guillaume Dumas, B. David Saunders, Gilles Villard. On Efficient Sparse Integer Matrix Smith Normal Form Computations. Journal of Symbolic Computation, 2001, 32 (1-2), pp.71-99. ⟨10.1006/jsco.2001.0451⟩. ⟨hal-02018782⟩