This does not make sense. See Wikipedia
sage: R = Zmod(6)
sage: R in Domains()
False
sage: M=matrix(R,[[4,3],[2,1]])
sage: M.smith_form()
---------------------------------------------------------------------------
Traceback (most recent call last)
TypeError: Smith form only defined over Noetherian integral domains
In fact, the error is misleading: it should say that the Smith form makes sense only over Principal ideal rings which are domains (I wrote it like that for emphasis: of course, I mean the Principal Ideal Domains). There is a similar structure theorem for modules over Dedekind domains, but now the torsion-free = projective which is free + fractional ideal. This probably is not implemented. At any rate, this is not Smith normal form (probably, this must be called the Steinitz normal form...).
Asked: 2014-08-08 08:38:16 +0100
Seen: 3,213 times
Last updated: Aug 08 '14
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.
