org.ejml.alg.dense.decomposition.bidiagonal

Class BidiagonalDecompositionTall_D64

java.lang.Object

org.ejml.alg.dense.decomposition.bidiagonal.BidiagonalDecompositionTall_D64

All Implemented Interfaces:

BidiagonalDecomposition<DenseMatrix64F>, DecompositionInterface<DenseMatrix64F>

public class BidiagonalDecompositionTall_D64
extends Object
implements BidiagonalDecomposition<DenseMatrix64F>

BidiagonalDecomposition specifically designed for tall matrices. First step is to perform QR decomposition on the input matrix. Then R is decomposed using a bidiagonal decomposition. By performing the bidiagonal decomposition on the smaller matrix computations can be saved if m/n > 5/3 and if U is NOT needed.

A = [Q₁ Q₂][U1 0; 0 I] [B1;0] V^T
U=[Q₁*U1 Q₂]
B=[B1;0]
A = U*B*V^T

A QRP decomposition is used internally. That decomposition relies an a fixed threshold for selecting singular values and is known to be less stable than SVD. There is the potential for a degregation of stability by using BidiagonalDecompositionTall instead of BidiagonalDecomposition. A few simple tests have shown that loss in stability to be insignificant.

See page 404 in "Fundamentals of Matrix Computations", 2nd by David S. Watkins.

Author:

Peter Abeles

Constructor Summary

Constructors

Constructor and Description
BidiagonalDecompositionTall_D64()

Method Summary

Methods

Modifier and Type	Method and Description
boolean	decompose(DenseMatrix64F orig) Computes the decomposition of the input matrix.
DenseMatrix64F	getB(DenseMatrix64F B, boolean compact) Returns the bidiagonal matrix.
void	getDiagonal(double[] diag, double[] off) Extracts the diagonal and off diagonal elements from the decomposition.
DenseMatrix64F	getU(DenseMatrix64F U, boolean transpose, boolean compact) Returns the orthogonal U matrix.
DenseMatrix64F	getV(DenseMatrix64F V, boolean transpose, boolean compact) Returns the orthogonal V matrix.
boolean	inputModified() Is the input matrix to DecompositionInterface.decompose(org.ejml.data.Matrix64F) is modified during the decomposition process.

Methods inherited from class java.lang.Object

clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Constructor Detail

BidiagonalDecompositionTall_D64

public BidiagonalDecompositionTall_D64()

Method Detail

getDiagonal

public void getDiagonal(double[] diag,
               double[] off)

Description copied from interface: BidiagonalDecomposition

Extracts the diagonal and off diagonal elements from the decomposition.

Specified by:

getDiagonal in interface BidiagonalDecomposition<DenseMatrix64F>

Parameters:

diag - diagonal elements from B.

off - off diagonal elements form B.

getB

public DenseMatrix64F getB(DenseMatrix64F B,
                  boolean compact)

Description copied from interface: BidiagonalDecomposition

Returns the bidiagonal matrix.

Specified by:

getB in interface BidiagonalDecomposition<DenseMatrix64F>

Parameters:

B - If not null the results are stored here, if null a new matrix is created.

Returns:

The bidiagonal matrix.

getU

public DenseMatrix64F getU(DenseMatrix64F U,
                  boolean transpose,
                  boolean compact)

Description copied from interface: BidiagonalDecomposition

Returns the orthogonal U matrix.

Specified by:

getU in interface BidiagonalDecomposition<DenseMatrix64F>

Parameters:

U - If not null then the results will be stored here. Otherwise a new matrix will be created.

Returns:

The extracted Q matrix.

getV

public DenseMatrix64F getV(DenseMatrix64F V,
                  boolean transpose,
                  boolean compact)

Description copied from interface: BidiagonalDecomposition

Returns the orthogonal V matrix.

Specified by:

getV in interface BidiagonalDecomposition<DenseMatrix64F>

Parameters:

V - If not null then the results will be stored here. Otherwise a new matrix will be created.

Returns:

The extracted Q matrix.

decompose

public boolean decompose(DenseMatrix64F orig)

Description copied from interface: DecompositionInterface

Computes the decomposition of the input matrix. Depending on the implementation the input matrix might be stored internally or modified. If it is modified then the function DecompositionInterface.inputModified() will return true and the matrix should not be modified until the decomposition is no longer needed.

Specified by:

decompose in interface DecompositionInterface<DenseMatrix64F>

Parameters:

orig - The matrix which is being decomposed. Modification is implementation dependent.

Returns:

Returns if it was able to decompose the matrix.

inputModified

public boolean inputModified()

Description copied from interface: DecompositionInterface

Is the input matrix to DecompositionInterface.decompose(org.ejml.data.Matrix64F) is modified during the decomposition process.

Specified by:

inputModified in interface DecompositionInterface<DenseMatrix64F>

Returns:

true if the input matrix to decompose() is modified.