Numerical algorithms for the computation of the Smith normal form of integral matrices are described. More specifically, the compound matrix method, methods based on elementary row or column operations and methods using modular or p-adic arithmetic are presented. A variety of examples and numerical results are given illustrating the execution of the algorithms.
History
Citation
C. Koukouvinos, M. Mitrouli and Jennifer Seberry, Numerical algorithms for the computation of the Smith normal form of integral matrices, Congressus Numerantium, 133 (1998), 127-162.