Dear Forum, Is there a possibility to find in GAP the Hermite normal form (and the transformation matrix) for a large matrix over the integers with few nonzero entries (i.e. a sparse matrix)? The package "Gauss" seems to offer this functionality for fields and finite residue rings only. Anvita