Package org.ejml.dense.row.decomposition.chol

Class CholeskyDecompositionBlock_FDRM

java.lang.Object

org.ejml.dense.row.decomposition.chol.CholeskyDecompositionCommon_FDRM

org.ejml.dense.row.decomposition.chol.CholeskyDecompositionBlock_FDRM

All Implemented Interfaces:

org.ejml.interfaces.decomposition.CholeskyDecomposition<org.ejml.data.FMatrixRMaj>, org.ejml.interfaces.decomposition.CholeskyDecomposition_F32<org.ejml.data.FMatrixRMaj>, org.ejml.interfaces.decomposition.DecompositionInterface<org.ejml.data.FMatrixRMaj>

public class CholeskyDecompositionBlock_FDRM
extends CholeskyDecompositionCommon_FDRM

This is an implementation of Cholesky that processes internal submatrices as blocks. This is done to reduce the number of cache issues.

Field Summary

Fields inherited from class org.ejml.dense.row.decomposition.chol.CholeskyDecompositionCommon_FDRM

det, lower, maxWidth, n, t, T, vv

Constructor Summary

Constructors

Constructor	Description
CholeskyDecompositionBlock_FDRM(int blockWidth)	Creates a CholeksyDecomposition capable of decomposing a matrix that is n by n, where n is the width.

Method Summary

Modifier and Type	Method	Description
protected boolean	decomposeLower()	Performs Choleksy decomposition on the provided matrix.
protected boolean	decomposeUpper()	Performs an upper triangular decomposition.
void	setExpectedMaxSize(int numRows, int numCols)	Declares additional internal data structures.
static void	solveL_special(float[] L, org.ejml.data.FMatrixRMaj b_src, int indexSrc, int indexDst, org.ejml.data.FMatrixRMaj B)	This is a variation on the TriangularSolver_FDRM.solveL(float[], float[], int) function.
static void	symmRankTranA_sub(org.ejml.data.FMatrixRMaj a, org.ejml.data.FMatrixRMaj c, int startIndexC)	Performs this operation: c = c - aTa where c is a submatrix.

Methods inherited from class org.ejml.dense.row.decomposition.chol.CholeskyDecompositionCommon_FDRM

_getVV, computeDeterminant, decompose, getT, getT, inputModified, isLower

Methods inherited from class java.lang.Object

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

Constructor Detail

CholeskyDecompositionBlock_FDRM

public CholeskyDecompositionBlock_FDRM(int blockWidth)

Creates a CholeksyDecomposition capable of decomposing a matrix that is n by n, where n is the width.

Parameters:

blockWidth - The width of a block.

Method Detail

setExpectedMaxSize

public void setExpectedMaxSize(int numRows,
                               int numCols)

Declares additional internal data structures.

Overrides:

setExpectedMaxSize in class CholeskyDecompositionCommon_FDRM

decomposeLower

protected boolean decomposeLower()

Performs Choleksy decomposition on the provided matrix.

If the matrix is not positive definite then this function will return false since it can't complete its computations. Not all errors will be found.

Specified by:

decomposeLower in class CholeskyDecompositionCommon_FDRM

Returns:

True if it was able to finish the decomposition.

decomposeUpper

protected boolean decomposeUpper()

Description copied from class: CholeskyDecompositionCommon_FDRM

Performs an upper triangular decomposition.

Specified by:

decomposeUpper in class CholeskyDecompositionCommon_FDRM

Returns:

true if the matrix was decomposed.

solveL_special

public static void solveL_special(float[] L,
                                  org.ejml.data.FMatrixRMaj b_src,
                                  int indexSrc,
                                  int indexDst,
                                  org.ejml.data.FMatrixRMaj B)

This is a variation on the TriangularSolver_FDRM.solveL(float[], float[], int) function. It grabs the input from the top right row rectangle of the source matrix then writes the results to the lower bottom column rectangle. The rectangle matrices just matrices are submatrices of the matrix that is being decomposed. The results are also written to B.

Parameters:

L - A lower triangular matrix.

b_src - matrix with the vectors that are to be solved for

indexSrc - First index of the submatrix where the inputs are coming from.

indexDst - First index of the submatrix where the results are going to.

B -

symmRankTranA_sub

public static void symmRankTranA_sub(org.ejml.data.FMatrixRMaj a,
                                     org.ejml.data.FMatrixRMaj c,
                                     int startIndexC)

Performs this operation:

c = c - a^Ta
where c is a submatrix.

Only the upper triangle is updated.

Parameters:

a - A matrix.

c - A matrix.

startIndexC - start of the submatrix in c.