This Python package provides Python bindings for the software GUPTRI and can be used with NumPy and, optionally, SageMath.

GUPTRI is a Fortran library by Jim Demmel and Bo Kågström for robust computation of generalized eigenvalues of singular matrix pencils. Standard tools like LAPACK do not reliably handle singular generalized eigenvalue problems.

GUPTRI solves this by computing a generalized block upper triangular form (generalized Schur staircase form) of a matrix pencil, revealing the Kronecker structure of the pencil. For details, see the documentation and the references therein.

See the examples and documentation at https://guptri-py.readthedocs.io.

Requirements: NumPy and, optionally, SageMath (tested with Sage 9.6 on Arch Linux and with earlier versions on macOS).

First, clone the repository from GitHub:

git clone https://github.com/mwageringel/guptri_py.git && cd guptri_py

To install with Python 3 and NumPy, run the following command:

pip3 install --upgrade --no-index -v .

Alternatively, for use with Sage, run this command:

sage -pip install --upgrade --no-index -v .

To install into the Python user install directory (no root access required), use:

sage -pip install --upgrade --no-index -v --user .

After successful installation, run the tests with Sage:

sage -t guptri_py

• With NumPy ≤ 1.17, it may be necessary to set:

  export NPY_DISTUTILS_APPEND_FLAGS=1

  to fix a linking problem. See numpy/numpy#12799.