This class implements the least square solution of a set of linear equations of m unknowns (n >= m) using the Householder method. It solves A.X = B. This algorithm has more numerical stability than GaussLeastSquare but is longer. It must be used if the matrix is singular or nearly singular. It is about 16% longer than GaussLeastSquare if there is only one member B to solve. It is about 30% longer if there are twenty B members to solve. More...

#include <math_Householder.hxx>

Public Member Functions
	math_Householder (const math_Matrix &A, const math_Matrix &B, const Standard_Real EPS=1.0e-20)
	Given an input matrix A with n>= m, given an input matrix B this constructor performs the least square resolution of the set of linear equations A.X = B for each column of B. If a column norm is less than EPS, the resolution can't be done. Exception DimensionError is raised if the row number of B is different from the A row number. More...

	math_Householder (const math_Matrix &A, const math_Matrix &B, const Standard_Integer lowerArow, const Standard_Integer upperArow, const Standard_Integer lowerAcol, const Standard_Integer upperAcol, const Standard_Real EPS=1.0e-20)
	Given an input matrix A with n>= m, given an input matrix B this constructor performs the least square resolution of the set of linear equations A.X = B for each column of B. If a column norm is less than EPS, the resolution can't be done. Exception DimensionError is raised if the row number of B is different from the A row number. More...

	math_Householder (const math_Matrix &A, const math_Vector &B, const Standard_Real EPS=1.0e-20)
	Given an input matrix A with n>= m, given an input vector B this constructor performs the least square resolution of the set of linear equations A.X = B. If a column norm is less than EPS, the resolution can't be done. Exception DimensionError is raised if the length of B is different from the A row number. More...

Standard_Boolean	IsDone () const
	Returns true if the computations are successful, otherwise returns false. More...

void	Value (math_Vector &sol, const Standard_Integer Index=1) const
	Given the integer Index, this routine returns the corresponding least square solution sol. Exception NotDone is raised if the resolution has not be done. Exception OutOfRange is raised if Index <=0 or Index is more than the number of columns of B. More...

const math_Matrix &	AllValues () const
	Returns the matrix sol of all the solutions of the system A.X = B. Exception NotDone is raised is the resolution has not be done. More...

void	Dump (Standard_OStream &o) const
	Prints information on the current state of the object. More...

Protected Member Functions
void	Perform (const math_Matrix &A, const math_Matrix &B, const Standard_Real EPS)
	This method is used internally for each constructor above and can't be used directly. More...

Detailed Description

This class implements the least square solution of a set of linear equations of m unknowns (n >= m) using the Householder method. It solves A.X = B. This algorithm has more numerical stability than GaussLeastSquare but is longer. It must be used if the matrix is singular or nearly singular. It is about 16% longer than GaussLeastSquare if there is only one member B to solve. It is about 30% longer if there are twenty B members to solve.

Constructor & Destructor Documentation

◆ math_Householder() [1/3]

math_Householder::math_Householder	(	const math_Matrix &	A,
		const math_Matrix &	B,
		const Standard_Real	EPS = `1.0e-20`
	)

Given an input matrix A with n>= m, given an input matrix B this constructor performs the least square resolution of the set of linear equations A.X = B for each column of B. If a column norm is less than EPS, the resolution can't be done. Exception DimensionError is raised if the row number of B is different from the A row number.

◆ math_Householder() [2/3]

math_Householder::math_Householder	(	const math_Matrix &	A,
		const math_Matrix &	B,
		const Standard_Integer	lowerArow,
		const Standard_Integer	upperArow,
		const Standard_Integer	lowerAcol,
		const Standard_Integer	upperAcol,
		const Standard_Real	EPS = `1.0e-20`
	)

Given an input matrix A with n>= m, given an input matrix B this constructor performs the least square resolution of the set of linear equations A.X = B for each column of B. If a column norm is less than EPS, the resolution can't be done. Exception DimensionError is raised if the row number of B is different from the A row number.

◆ math_Householder() [3/3]

math_Householder::math_Householder	(	const math_Matrix &	A,
		const math_Vector &	B,
		const Standard_Real	EPS = `1.0e-20`
	)

Given an input matrix A with n>= m, given an input vector B this constructor performs the least square resolution of the set of linear equations A.X = B. If a column norm is less than EPS, the resolution can't be done. Exception DimensionError is raised if the length of B is different from the A row number.

Member Function Documentation

◆ AllValues()

const math_Matrix& math_Householder::AllValues ( ) const

Returns the matrix sol of all the solutions of the system A.X = B. Exception NotDone is raised is the resolution has not be done.

◆ Dump()

void math_Householder::Dump ( Standard_OStream & o ) const

Prints information on the current state of the object.

◆ IsDone()

Standard_Boolean math_Householder::IsDone ( ) const

Returns true if the computations are successful, otherwise returns false.

◆ Perform()

void math_Householder::Perform	(	const math_Matrix &	A,
		const math_Matrix &	B,
		const Standard_Real	EPS
	)

protected

This method is used internally for each constructor above and can't be used directly.

◆ Value()

void math_Householder::Value	(	math_Vector &	sol,
		const Standard_Integer	Index = `1`
	)		const

Given the integer Index, this routine returns the corresponding least square solution sol. Exception NotDone is raised if the resolution has not be done. Exception OutOfRange is raised if Index <=0 or Index is more than the number of columns of B.

The documentation for this class was generated from the following file:

math_Householder.hxx

Public Member Functions

Protected Member Functions