CholeskyMethod

class CholeskyMethod(*args)

Least squares solver using Cholesky decomposition.

Available constructors:

CholeskyMethod(proxy, weight, indices)

CholeskyMethod(proxy, indices)

CholeskyMethod(design)

Parameters
proxyDesignProxy

Input sample

weightsequence of float

Output weights

indicessequence of int

Indices allowed in the basis

design2-d sequence of float

A priori known design matrix

See also

LeastSquaresMethod, SVDMethod, QRMethod

Examples

Solves a linear least squares problem with Cholesky method:

>>> import openturns as ot
>>> A = ot.Matrix([[1,1],[1,2],[1,3],[1,4]])
>>> y = ot.Point([6,5,7,10])
>>> method = ot.CholeskyMethod(A)
>>> x = method.solve(y)
>>> print(x)
[3.5,1.4]

Methods

computeWeightedDesign(self[, whole])

Build the design matrix.

getBasis(self)

Accessor to the basis.

getClassName(self)

Accessor to the object’s name.

getCurrentIndices(self)

Current indices accessor.

getGramInverse(self)

Get the inverse Gram matrix of input sample.

getGramInverseDiag(self)

Get the diagonal of the inverse Gram matrix.

getGramInverseTrace(self)

Get the trace of the inverse Gram matrix.

getH(self)

Get the projection matrix H.

getHDiag(self)

Get the diagonal of the projection matrix H.

getId(self)

Accessor to the object’s id.

getInitialIndices(self)

Initial indices accessor.

getInputSample(self)

Input sample accessor.

getName(self)

Accessor to the object’s name.

getShadowedId(self)

Accessor to the object’s shadowed id.

getVisibility(self)

Accessor to the object’s visibility state.

getWeight(self)

Accessor to the weights.

hasName(self)

Test if the object is named.

hasVisibleName(self)

Test if the object has a distinguishable name.

setName(self, name)

Accessor to the object’s name.

setShadowedId(self, id)

Accessor to the object’s shadowed id.

setVisibility(self, visible)

Accessor to the object’s visibility state.

solve(self, rhs)

Solve the least-squares problem.

solveNormal(self, rhs)

Solve the least-squares problem using normal equation.

trashDecomposition(self)

Drop the current decomposition.

update(self, addedIndices, conservedIndices, …)

Update the current decomposition.
__init__(self, \*args)

Initialize self. See help(type(self)) for accurate signature.

computeWeightedDesign(self, whole=False)

Build the design matrix.

Parameters
wholebool, defaults to False

Whether to use the initial indices instead of the current indices

Returns
psiAkMatrix

The design matrix

getBasis(self)

Accessor to the basis.

Returns
basiscollection of Function

Basis.

getClassName(self)

Accessor to the object’s name.

Returns
class_namestr

The object class name (object.__class__.__name__).

getCurrentIndices(self)

Current indices accessor.

Returns
indicesIndices

Indices of the current decomposition in the global basis.

getGramInverse(self)

Get the inverse Gram matrix of input sample.

G^{-1} = (X^T * X)^{-1}

Returns
cCovarianceMatrix

The inverse Gram matrix.

getGramInverseDiag(self)

Get the diagonal of the inverse Gram matrix.

diag(G^{-1}) = diag((X^T * X)^{-1})

Returns
dPoint

The diagonal of the inverse Gram matrix.

getGramInverseTrace(self)

Get the trace of the inverse Gram matrix.

Tr(G^{-1}) = Tr(x^T * x)^{-1}

Returns
xScalar

The trace of inverse Gram matrix.

getH(self)

Get the projection matrix H.

H = X * (X^T * X)^{-1} * X^T

Returns
hSymmetricMatrix

The projection matrix H.

getHDiag(self)

Get the diagonal of the projection matrix H.

H = X * (X^T * X)^{-1} * X^T

Returns
dPoint

The diagonal of H.

getId(self)

Accessor to the object’s id.

Returns
idint

Internal unique identifier.

getInitialIndices(self)

Initial indices accessor.

Returns
indicesIndices

Initial indices of the terms in the global basis.

getInputSample(self)

Input sample accessor.

Returns
inputSampleSample

Input sample.

getName(self)

Accessor to the object’s name.

Returns
namestr

The name of the object.

getShadowedId(self)

Accessor to the object’s shadowed id.

Returns
idint

Internal unique identifier.

getVisibility(self)

Accessor to the object’s visibility state.

Returns
visiblebool

Visibility flag.

getWeight(self)

Accessor to the weights.

Returns
weightPoint

Weights.

hasName(self)

Test if the object is named.

Returns
hasNamebool

True if the name is not empty.

hasVisibleName(self)

Test if the object has a distinguishable name.

Returns
hasVisibleNamebool

True if the name is not empty and not the default one.

setName(self, name)

Accessor to the object’s name.

Parameters
namestr

The name of the object.

setShadowedId(self, id)

Accessor to the object’s shadowed id.

Parameters
idint

Internal unique identifier.

setVisibility(self, visible)

Accessor to the object’s visibility state.

Parameters
visiblebool

Visibility flag.

solve(self, rhs)

Solve the least-squares problem.

\vect{a} = \argmin_{\vect{x} \in \Rset^P} ||M\vect{x}-\vect{b}||^2

Parameters
bsequence of float

Second term of the equation

Returns
aPoint

The solution.

solveNormal(self, rhs)

Solve the least-squares problem using normal equation.

M^T*M*x=M^T*b

Parameters
bsequence of float

Second term of the equation

Returns
xPoint

The solution.

trashDecomposition(self)

Drop the current decomposition.

update(self, addedIndices, conservedIndices, removedIndices, row=False)

Update the current decomposition.

Parameters
addedIndicessequence of int

Indices of added basis terms.

conservedIndicessequence of int

Indices of conserved basis terms.

removedIndicessequence of int

Indices of removed basis terms.