Compute the Hermite normal form with torsion removed.
The input matrix.
The defactored Hermite normal form.
Compute the Hermite normal form with torsion removed.