Package org.ejml.dense.row.decomposition.bidiagonal

Class BidiagonalDecompositionTall_MT_FDRM

java.lang.Object

org.ejml.dense.row.decomposition.bidiagonal.BidiagonalDecompositionTall_MT_FDRM

All Implemented Interfaces:

BidiagonalDecomposition<FMatrixRMaj>, BidiagonalDecomposition_F32<FMatrixRMaj>, DecompositionInterface<FMatrixRMaj>

@Generated("org.ejml.dense.row.decomposition.bidiagonal.BidiagonalDecompositionTall_FDRM")
public class BidiagonalDecompositionTall_MT_FDRM
extends Object
implements BidiagonalDecomposition_F32<FMatrixRMaj>

BidiagonalDecomposition_F32 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_F32. 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.

Constructor Summary

Constructors

Constructor	Description
BidiagonalDecompositionTall_MT_FDRM()

Method Summary

Modifier and Type	Method	Description
boolean	decompose(FMatrixRMaj orig)	Computes the decomposition of the input matrix.
FMatrixRMaj	getB(@Nullable FMatrixRMaj B, boolean compact)	Returns the bidiagonal matrix.
void	getDiagonal(float[] diag, float[] off)	Extracts the diagonal and off diagonal elements from the decomposition.
FMatrixRMaj	getU(@Nullable FMatrixRMaj U, boolean transpose, boolean compact)	Returns the orthogonal U matrix.
FMatrixRMaj	getV(@Nullable FMatrixRMaj V, boolean transpose, boolean compact)	Returns the orthogonal V matrix.
boolean	inputModified()	Checks if the input matrix to DecompositionInterface.decompose(org.ejml.data.Matrix) 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 Details

BidiagonalDecompositionTall_MT_FDRM

public BidiagonalDecompositionTall_MT_FDRM()

Method Details

getDiagonal

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

Description copied from interface: BidiagonalDecomposition_F32

Extracts the diagonal and off diagonal elements from the decomposition.

Specified by:

getDiagonal in interface BidiagonalDecomposition_F32<FMatrixRMaj>

Parameters:

diag - diagonal elements from B.

off - off diagonal elements form B.

getB

public FMatrixRMaj getB(@Nullable
 @Nullable FMatrixRMaj B,
 boolean compact)

Description copied from interface: BidiagonalDecomposition

Returns the bidiagonal matrix.

Specified by:

getB in interface BidiagonalDecomposition<FMatrixRMaj>

Parameters:

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

Returns:

The bidiagonal matrix.

getU

public FMatrixRMaj getU(@Nullable
 @Nullable FMatrixRMaj U,
 boolean transpose,
 boolean compact)

Description copied from interface: BidiagonalDecomposition

Returns the orthogonal U matrix.

Specified by:

getU in interface BidiagonalDecomposition<FMatrixRMaj>

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 FMatrixRMaj getV(@Nullable
 @Nullable FMatrixRMaj V,
 boolean transpose,
 boolean compact)

Description copied from interface: BidiagonalDecomposition

Returns the orthogonal V matrix.

Specified by:

getV in interface BidiagonalDecomposition<FMatrixRMaj>

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(FMatrixRMaj 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<FMatrixRMaj>

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

Checks if the input matrix to DecompositionInterface.decompose(org.ejml.data.Matrix) is modified during the decomposition process.

Specified by:

inputModified in interface DecompositionInterface<FMatrixRMaj>

Returns:

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