# A fast algorithm for computing the Smith normal form with multipliers for a nonsingular integer matrix

A Las Vegas randomized algorithm is given to compute the Smith multipliers for a nonsingular integer matrix A, that is, unimodular matrices U and V such that AV=US, with S the Smith normal form of A. The expected running time of the algorithm is about the same as required to multiply together two matrices of the same dimension and size of entries as A. Explicit bounds are given for the size of the entries in both unimodular multipliers. The main tool used by the algorithm is the Smith massager, a relaxed version of V, the unimodular matrix specifying the column operations of the Smith computation. From the perspective of efficiency, the main tools used are fast linear solving and partial linearization of integer matrices. As an application of the Smith with multipliers algorithm, a fast algorithm is given to find the fractional part of the inverse of the input matrix.

• 2 publications
• 15 publications
• 4 publications
research
01/31/2021

### Linear Computation Coding

We introduce the new concept of computation coding. Similar to how rate-...
research
09/21/2022

### A cubic algorithm for computing the Hermite normal form of a nonsingular integer matrix

A Las Vegas randomized algorithm is given to compute the Hermite normal ...
research
02/05/2016

### A fast, deterministic algorithm for computing a Hermite Normal Form of a polynomial matrix

Given a square, nonsingular matrix of univariate polynomials F∈K[x]^n × ...
research
05/23/2011

### A Tool for Integer Homology Computation: Lambda-At Model

In this paper, we formalize the notion of lambda-AT-model (where λ is a ...
research
12/14/2022

### On Eigenvalue Gaps of Integer Matrices

Given an n× n matrix with integer entries in the range [-h,h], how close...
research
03/11/2018

### Banded Matrix Fraction Representation of Triangular Input Normal Pairs

An input pair (A,B) is triangular input normal if and only if A is trian...
research
03/17/2018

### Orthogonal Representations for Output System Pairs

A new class of canonical forms is given proposed in which (A, C) is in H...