Smith normal form

Let A be a non-zero m×n matrix over a principal ideal domain (PID). Then, there exist invertible m×m and n×n matrices S and T, so that

S·A·T = α10000α200000αr000

where the diagonal elements satisfy αi|αi+1. They are called invariant factors.

Syntax

smith_normal_form(Matrix)

Description

Given a matrix, returns a vector with the invariant factors.

linear_algebra.smith_normal_form.calc.png