The command MatrixInverse(A, rc, R) returns two lists.

The first list the command returns is a list of pairs  where  is a regular chain and  is the inverse of A modulo the saturated ideal of .

The second list the command returns is a list of triplets  where  is a regular chain and A is the input matrix such that A is not invertible modulo the saturated ideal of .

All the returned regular chains  form a triangular decomposition of rc (in the sense of Kalkbrener).

It is assumed that rc is strongly normalized.

The algorithm is an adaptation of the algorithm of Bareiss.

This command is part of the RegularChains[MatrixTools] package, so it can be used in the form MatrixInverse(..) only after executing the command with(RegularChains[MatrixTools]).  However, it can always be accessed through the long form of the command by using RegularChains[MatrixTools][MatrixInverse](..).