Main       Download       Commercial support       FAQ       Forum       About Us

Reduction of a symmetric matrix to tridiagonal form

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.

Manual entries

C++ ortfac subpackage   
C# ortfac subpackage   

This article is intended for personal use only.

Download ALGLIB

C#

C# source.

Downloads page

 

C++

C++ source.

Downloads page

 

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.

Downloads page

 

FreePascal

FreePascal version.

Downloads page

 

Delphi

Delphi version.

Downloads page

 

VB.NET

VB.NET version.

Downloads page

 

VBA

VBA version.

Downloads page

 

Python

Python version (CPython and IronPython are supported).

Downloads page

 

 

ALGLIB® - numerical analysis library, 1999-2012.
ALGLIB is a registered trademark of the ALGLIB Project.