The symmetric matrix A could be represented as A=Q·T·Q ^T, where Q is an orthogonal matrix, and T is a tridiagonal matrix. We can say that matrix A is reduced to a tridiagonal matrix via a similarity transformation. Q ^T·A·Q = T.

Like other algorithms of orthogonal factorization (for example, QR and LQ decomposition algorithms), this algorithm uses a sequence of elementary reflections to transform the matrix A. The transformations are applied to the matrix both from the right and from the left, preserving the matrix symmetry on each stage and sequentially removing non-diagonal elements.

As a result of SMatrixTD subroutine, matrix A is replaced by the tridiagonal matrix T and a sequence of reflections transformations stored in a compact form. The format of the matrix and the subroutine parameters are described in detail in the subroutine comments; there we can note an analogy with QR-decomposition, that uses the lower triangular part of the matrix R to store the matrix Q and utilizes a very similar data storage format. As with QR decomposition, a subroutine for "unpacking" the matrix Q is presented: SMatrixTDUnpackQ.

This algorithm is transferred from the LAPACK library.

Report a bug

C#

C# 1.0 source.
tridiagonal.csharp.zip - Reduction of a symmetric matrix to tridiagonal form

C++

C++ source.
tridiagonal.cpp.zip - Reduction of a symmetric matrix to tridiagonal form
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.
tridiagonal.mpfr.zip - Reduction of a symmetric matrix to tridiagonal form
mpfr.zip - precompiled Win32 MPFR/GMP binaries

Delphi

Delphi source.
Can be compiled under FPC (in Delphi compatibility mode).
tridiagonal.delphi.zip - Reduction of a symmetric matrix to tridiagonal form

Visual Basic 6

Visual Basic 6 source.
tridiagonal.vb6.zip - Reduction of a symmetric matrix to tridiagonal form

Zonnon ^beta

Zonnon source.
Zonnon is an experimental language developed at ETH Zurich.
See www.zonnon.ethz.ch for more information.
tridiagonal.zonnon.zip - Reduction of a symmetric matrix to tridiagonal form