Engee documentation

LDL Factorization

Decomposition of a square Hermite positively definite matrix into lower, upper and diagonal components.

ldl factorizations

Description

Block LDL Factorization uniquely decomposes the Hermite positively defined input matrix as

where:

  • - lower triangular matrix with unit diagonal elements (unitriangular matrix).

  • - diagonal matrix.

  • - Hermite (complex-conjugate) transposed matrix L.

Only the diagonal and lower triangle of the input matrix are used. Any imaginary component of the diagonal elements is ignored.

The LDL decomposition requires half the computation of the LU decomposition and is always stable. It is more efficient than the Choletsky decomposition because it avoids the computation of square roots of diagonal elements.

Ports

Input

Port_1 - initial matrix S
`matrix M by M

The input signal is in the form of a square matrix.

The algorithm requires the input matrix to be square and Hermite positively defined. When the input is not positively defined, the block responds as specified by the input parameter Non-positive Definity.

Data types: Float64

Support for complex numbers: yes

Output

Port_1 - output matrix
`matrix M by M

The output is a composite matrix with elements below the main diagonal of , diagonal elements of and elements above the main diagonal of ,

where:

  • - lower triangular matrix with unit diagonal elements (unitriangular matrix).

  • - diagonal matrix.

  • - Hermite (complex-conjugate) transposed matrix L.

The output format is shown below for a 5 by 5 matrix.

ldl factorizations 1

Data types: Float64

Parameters

Main

Non-positive definite input - response to non-positive definite input matrix
Ignore (by default) | Warning | Error

Sets the action when non-positive definite input matrix parameters occur:

  • Ignore - the simulation continues and no warning is issued. The result is not a valid decomposition. A partial decomposition will be present in the upper left corner of the output matrix.

  • Warning - the simulation continues, but a warning message is displayed in the Engee command window. The result is not a valid decomposition. A partial decomposition will be present in the upper left corner of the output matrix.

  • Error - an error dialogue box is displayed and the simulation stops.