The complex triangular matrix could be inverted by using both the general algorithm and the specialized algorithm which consider the form of the matrix. We can pick out the following types of triangular matrices:

• upper triangular
• upper triangular with unit diagonal
• lower triangular lower triangular with unit diagonal

When inverting an upper triangular matrix we get an upper triangular matrix. When inverting an upper triangular matrix unit diagonal we get an upper triangular matrix unit diagonal and so on. Therefore, the algorithm returns an inverse matrix in the same format. The elements below the main diagonal (for the upper triangular matrix, and similarly in other cases) are not modified by the algorithm.

It should be noted that the algorithm can invert matrices with a unit diagonal (which isn't stored). Similarly, the inverse matrix will be the matrix with the unit diagonal stored in the same format.

The subroutine CMatrixTRInverse performs the inversion of matrices of all listed types. The parameters IsUpper and IsUnitTriangular define the type of matrix. The subroutine parameters are described in more detail in the program comment.

This algorithm is transferred from the LAPACK library.

Report a bug

C#

C# 1.0 source.
ctrinverse.csharp.zip - Complex triangular matrix inverse

C++

C++ source.
ctrinverse.cpp.zip - Complex triangular matrix inverse
ablas.zip - optimized basic linear algebra subroutines with SSE2 support (for C++ sources only)

C++, multiple precision arithmetic

C++ source. MPFR/GMP is used.
GMP source is available from gmplib.org. MPFR source is available from www.mpfr.org.
ctrinverse.mpfr.zip - Complex triangular matrix inverse
mpfr.zip - precompiled Win32 MPFR/GMP binaries

Delphi

Delphi source.
Can be compiled under FPC (in Delphi compatibility mode).
ctrinverse.delphi.zip - Complex triangular matrix inverse

Visual Basic 6

Visual Basic 6 source.
ctrinverse.vb6.zip - Complex triangular matrix inverse

Zonnon ^beta

Zonnon source.
Zonnon is an experimental language developed at ETH Zurich.
See www.zonnon.ethz.ch for more information.
ctrinverse.zonnon.zip - Complex triangular matrix inverse