Hi, I need to perform a Smith decomposition on (a priori not square) matrices to find a certain change of basis related to the Smith invariants.
The Smith normal form <https://en.wikipedia.org/wiki/Smith_normal_form> is already implemented in SymPy in `sympy.matrices.normalforms`, but it returns only the Smith normal form, not its decomposition itself. In general, the Smith decomposition of a matrix A is defined in terms of three matrices V, D, W so that A = V * D * W where V,W are both square, invertible, and integer-valued matrices. Is there any way of obtaining those matrices without reimplementing the algorithm myself? I assume that these matrices are computed at least implicitely, but I could not find a way of returning them. thanks! -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to sympy+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/sympy/6705806e-1277-4688-bc5b-9e31112e9c5bn%40googlegroups.com.
