port SmithNormalForm (CoCoA-4)package to CoCoA-5
#4 Updated by John Abbott 6 months ago
- Status changed from New to In Progress
- Target version changed from CoCoA-5.?.? to CoCoA-5.2.4
- % Done changed from 0 to 10
Anna says she has "cleaned up" the original CoCoA-4 code so that CoCoA-5 accepts it.
The code is in
MatNormalForm.cpkg5 (original file was
The package exports two names
There is no entry in the CoCoA-5 documentation.
I would prefer
SmithNormalForm rather than just
The code does not run properly. I got what I believe to be an infinite loop with the following input (given to a fresh CoCoA-5):
M := mat(ZZ,[[random(-9,9) | j in 1..4] | i in 1..4]); Smith(M);
matrix(ZZ, [[2, -8, -2, -1], [-1, 6, -1, 4], [-2, -1, -7, 0], [-2, -3, -7, -3]])