Hallo!

Ich bin auf der Suche nach Bibliotheken für C++ bzw. C, die das folgende Problem lösen:

Eingabe Eine mxn-Matrix A mit ganzzahligen (oder rationalen) Einträgen.

Ausgabe Eine invertierbare mxm-Matrix S und eine invertierbare nxn-Matrix T, sodass das Produkt SAT eine Diagonalmatrix ist.

Eine Möglichkeit wäre, die Smithsche Normalform von A (en.wikipedia.org/wiki/Smithnormalform) zu berechnen, was wohl mit linbox (linalg.org) möglich wäre.

Habt Ihr andere / bessere Vorschläge?

Vielen Dank für Eure Antworten!