Package org.ejml.dense.block.linsol.chol

Class CholeskyOuterSolver_FDRB

java.lang.Object

org.ejml.dense.block.linsol.chol.CholeskyOuterSolver_FDRB

All Implemented Interfaces:

LinearSolver<FMatrixRBlock,FMatrixRBlock>, LinearSolverDense<FMatrixRBlock>

@Generated("org.ejml.dense.block.linsol.chol.CholeskyOuterSolver_DDRB")
public class CholeskyOuterSolver_FDRB
extends Object
implements LinearSolverDense<FMatrixRBlock>

Linear solver that uses a block cholesky decomposition.

Solver works by using the standard Cholesky solving strategy:
A=L*L^T
A*x=b
L*L^T*x = b
L*y = b
L^T*x = y
x = L^-Ty

It is also possible to use the upper triangular cholesky decomposition.

Constructor Summary

Constructors

Constructor	Description
CholeskyOuterSolver_FDRB()

Method Summary

Modifier and Type	Method	Description
CholeskyDecomposition_F32<FMatrixRBlock>	getDecomposition()	If a decomposition class was used internally then this will return that class.
void	invert(FMatrixRBlock A_inv)	Computes the inverse of of the 'A' matrix passed into LinearSolver.setA(Matrix) and writes the results to the provided matrix.
boolean	modifiesA()	Returns true if the passed in matrix to LinearSolver.setA(Matrix) is modified.
boolean	modifiesB()	Returns true if the passed in 'B' matrix to LinearSolver.solve(Matrix, Matrix) is modified.
double	quality()	Returns a very quick to compute measure of how singular the system is.
boolean	setA(FMatrixRBlock A)	Decomposes and overwrites the input matrix.
void	solve(FMatrixRBlock B, @Nullable FMatrixRBlock X)	If X == null then the solution is written into B.

Methods inherited from class java.lang.Object

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

Constructor Details

CholeskyOuterSolver_FDRB

public CholeskyOuterSolver_FDRB()

Method Details

setA

public boolean setA(FMatrixRBlock A)

Decomposes and overwrites the input matrix.

Specified by:

setA in interface LinearSolver<FMatrixRBlock,FMatrixRBlock>

Parameters:

A - Semi-Positive Definite (SPD) system matrix. Modified. Reference saved.

Returns:

If the matrix can be decomposed. Will always return false of not SPD.

quality

public double quality()

Description copied from interface: LinearSolver

Returns a very quick to compute measure of how singular the system is. This measure will be invariant to the scale of the matrix and always be positive, with larger values indicating it is less singular. If not supported by the solver then the runtime exception IllegalArgumentException is thrown. This is NOT the matrix's condition.

How this function is implemented is not specified. One possible implementation is the following: In many decompositions a triangular matrix is extracted. The determinant of a triangular matrix is easily computed and once normalized to be scale invariant and its absolute value taken it will provide functionality described above.

Specified by:

quality in interface LinearSolver<FMatrixRBlock,FMatrixRBlock>

Returns:

The quality of the linear system.

solve

public void solve(FMatrixRBlock B,
 @Nullable
 @Nullable FMatrixRBlock X)

If X == null then the solution is written into B. Otherwise the solution is copied from B into X.

Specified by:

solve in interface LinearSolver<FMatrixRBlock,FMatrixRBlock>

Parameters:

B - A matrix ℜ ^{m × p}. Might be modified.

X - A matrix ℜ ^{n × p}, where the solution is written to. Modified.

invert

public void invert(FMatrixRBlock A_inv)

Description copied from interface: LinearSolverDense

Computes the inverse of of the 'A' matrix passed into LinearSolver.setA(Matrix) and writes the results to the provided matrix. If 'A_inv' needs to be different from 'A' is implementation dependent.

Specified by:

invert in interface LinearSolverDense<FMatrixRBlock>

Parameters:

A_inv - Where the inverted matrix saved. Modified.

modifiesA

public boolean modifiesA()

Description copied from interface: LinearSolver

Returns true if the passed in matrix to LinearSolver.setA(Matrix) is modified.

Specified by:

modifiesA in interface LinearSolver<FMatrixRBlock,FMatrixRBlock>

Returns:

true if A is modified in setA().

modifiesB

public boolean modifiesB()

Description copied from interface: LinearSolver

Returns true if the passed in 'B' matrix to LinearSolver.solve(Matrix, Matrix) is modified.

Specified by:

modifiesB in interface LinearSolver<FMatrixRBlock,FMatrixRBlock>

Returns:

true if B is modified in solve(B,X).

getDecomposition

public CholeskyDecomposition_F32<FMatrixRBlock> getDecomposition()

Description copied from interface: LinearSolver

If a decomposition class was used internally then this will return that class. Most linear solvers decompose the input matrix into a more simplistic form. However some solutions do not require decomposition, e.g. inverse by minor.

Specified by:

getDecomposition in interface LinearSolver<FMatrixRBlock,FMatrixRBlock>

Returns:

Internal decomposition class. If there is none then null.