The Generalized Eigenvalue Decomposition

All variations of eigenvalue decompositions (real symmetric and non-symmetric, complex Hermitian and non-Hermitian) are represented by the GeneralizedEigenvalueDecompositionT class. It has no constructors. Instead, it is created by calling the GetEigenvalueDecomposition method on the primary matrix (A). This method has two overloads. The first overload has one argument: the secondary matrix (B). Both matrices must be square. The second overload takes an additional Boolean argument that specifies whether the contents of the matrices may be destroyed during the calculation of the eigenvalue decomposition.

var A = Matrix.CreateSymmetric(3, new double[]
    { 1,2,3, 2,4,9, 3,9,15 },
    MatrixTriangle.Upper, MatrixElementOrder.ColumnMajor);
var B = Matrix.CreateSymmetric(3, new double[]
    { 3,-2,3, -2,6,-9, 3,-9,17 },
    MatrixTriangle.Upper, MatrixElementOrder.ColumnMajor);
var geig = A.GetEigenvalueDecomposition(B);

Dim A = Matrix.CreateSymmetric(3, New Double() {
    1, 2, 3, 2, 4, 9, 3, 9, 15},
    MatrixTriangle.Upper, MatrixElementOrder.ColumnMajor)
Dim B = Matrix.CreateSymmetric(3, New Double() {
    3, -2, 3, -2, 6, -9, 3, -9, 17},
    MatrixTriangle.Upper, MatrixElementOrder.ColumnMajor)
Dim geig = A.GetEigenvalueDecomposition(B)

let A = Matrix.CreateSymmetric(3, 
            [| 1;2;3; 2;4;9; 3;9;15 |],
            MatrixTriangle.Upper, MatrixElementOrder.ColumnMajor)
let B = Matrix.CreateSymmetric(3, 
            [| 3;-2;3; -2;6;-9; 3;-9;17 |],
            MatrixTriangle.Upper, MatrixElementOrder.ColumnMajor)
let geig = A.GetEigenvalueDecomposition(B)

The Decompose method performs the actual decomposition. This method makes a copy of the matrix if necessary. It then calls the appropriate LAPACK routine to perform the actual decomposition. This method is called by other methods as needed. You will rarely need to call it explicitly.

Because real matrices may have complex eigenvalues and eigenvectors, the eigenvalue decomposition of a real matrix deserves some special attention.

For symmetric matrices where the secondary matrix is positive definite, all eigenvalues and eigenvectors are always real. The Eigenvalues property returns a vector containing the eigenvalues. The Eigenvectors property returns a matrix that has the corresponding eigenvectors as its columns:

var L = geig.Eigenvalues;
var X = geig.Eigenvectors;

Dim L = geig.Eigenvalues
Dim X = geig.Eigenvectors

let L = geig.Eigenvalues
let X = geig.Eigenvectors

Otherwise, there may be complex eigenvalues, which occur in complex conjugate pairs. The corresponding eigenvectors are also complex and also come in complex conjugate pairs. The Eigenvalues and Eigenvectors return only the real eigenvalues and eigenvectors. To get all eigenvalues and eigenvectors, use the ComplexEigenvalues and ComplexEigenvectors properties.

The RawEigenvectors property returns a matrix whose columns contain all the information about the eigenvectors. For real eigenvalues, the corresponding column contains the corresponding eigenvector. For complex eigenvalues, which always occur in complex conjugate pairs, two adjoining columns contain the real and imaginary elements of the corresponding eigenvectors.

var A = Matrix.Create(4, 4, new double[] {
        3.9,   4.3,   4.3,  4.4,
        12.5,  21.5,  21.5, 26.0,
        -34.5, -47.5, -43.5,-46.0,
        -0.5,  7.5,   3.5,  6.0 }, 
    MatrixElementOrder.ColumnMajor);
var B = Matrix.Create(4, 4, new double[] {
    1.0, 1.0, 1.0, 1.0,
    2.0, 3.0, 3.0, 3.0,
    -3.0, -5.0, -4.0, -4.0,
    1.0, 4.0, 3.0, 4.0 },
    MatrixElementOrder.ColumnMajor);
var geig = A.GetEigenvalueDecomposition(B);

var L = geig.Eigenvalues;
var X = geig.Eigenvectors;
var complexL = geig.ComplexEigenvalues;
var complexX = geig.ComplexEigenvectors;
var rawX = geig.RawEigenvectors;

Dim A = Matrix.Create(4, 4, New Double() {
        3.9, 4.3, 4.3, 4.4,
        12.5, 21.5, 21.5, 26.0,
        -34.5, -47.5, -43.5, -46.0,
        -0.5, 7.5, 3.5, 6.0},
    MatrixElementOrder.ColumnMajor)
Dim B = Matrix.Create(4, 4, New Double() {
        1.0, 1.0, 1.0, 1.0,
        2.0, 3.0, 3.0, 3.0,
        -3.0, -5.0, -4.0, -4.0,
        1.0, 4.0, 3.0, 4.0},
        MatrixElementOrder.ColumnMajor)
Dim geig = A.GetEigenvalueDecomposition(B)

Dim L = geig.Eigenvalues
Dim X = geig.Eigenvectors
Dim complexL = geig.ComplexEigenvalues
Dim complexX = geig.ComplexEigenvectors
Dim rawX = geig.RawEigenvectors

let A = Matrix.Create(4, 4, 
          [|
            3.9;   4.3;   4.3;  4.4;
            12.5;  21.5;  21.5; 26.0;
            -34.5; -47.5; -43.5;-46.0;
            -0.5;  7.5;   3.5;  6.0 
          |], MatrixElementOrder.ColumnMajor)
let B = Matrix.Create(4, 4,
          [|
            1.0; 1.0; 1.0; 1.0;
            2.0; 3.0; 3.0; 3.0;
            -3.0; -5.0; -4.0; -4.0;
            1.0; 4.0; 3.0; 4.0 
          |], MatrixElementOrder.ColumnMajor)
let geig = A.GetEigenvalueDecomposition(B)

let L = geig.Eigenvalues
let X = geig.Eigenvectors
let complexL = geig.ComplexEigenvalues
let complexX = geig.ComplexEigenvectors
let rawX = geig.RawEigenvectors

The generalized eigenvalues and eigenvalues of a pair of complex matrices are also complex. The Eigenvalues property returns a vector containing the eigenvalues. The Eigenvectors property returns a matrix that has the corresponding eigenvectors as its columns.

The generalized eigenvalues of a pair of Hermitian matrices where the secondary matrix is also positive definite, are also real. This simply means that the imaginary component of the eigenvalue vector is identically zero.

The ComplexEigenvalues, ComplexEigenvectors, and RawEigenvectors properties are meaningless for complex eigenvalue decompositions.

If the secondary matrix B is singular, then one or more of the generalized eigenvalues of the pair is infinite. The generalized eigenvalues can also be returned as pairs whose quotient is the generalized eigenvalue. When the denominator of the quotient is zero, the generalized eigenvalue is infinite.

Three properties give access to the quotient form of the generalized eigenvalues. The EigenvalueDenominators property returns a vector containing the denominators of the quotients. For real matrices, it is always real. The EigenvalueNumerators property returns a vector containing the numerators of the quotients. As with the generalized eigenvalues themselves, this property only returns real numerators for real matrices. To access the the full set of numerators, use the ComplexEigenvalueNumerators property.