Local Smith Form at p : via elimination
Not efficient to compute in Z(p) ? elimination mod pe
| Gauß-Bareiss with pivot invertible mod pe-k+1
| When no more invertible element
| Return actual rank as rank mod pk
| Divide every remaining non-zero coefficient by p
If rank mod pe = integer rank then stop
- Try first with pe close to word size