/usr/include/oce/math_Householder.hxx is in liboce-foundation-dev 0.18.2-2build1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 | // This file is generated by WOK (CPPExt).
// Please do not edit this file; modify original file instead.
// The copyright and license terms as defined for the original file apply to
// this header file considered to be the "object code" form of the original source.
#ifndef _math_Householder_HeaderFile
#define _math_Householder_HeaderFile
#include <Standard.hxx>
#include <Standard_DefineAlloc.hxx>
#include <Standard_Macro.hxx>
#include <math_Matrix.hxx>
#include <Standard_Boolean.hxx>
#include <Standard_Integer.hxx>
#include <Standard_Real.hxx>
#include <math_Vector.hxx>
#include <Standard_OStream.hxx>
class StdFail_NotDone;
class Standard_OutOfRange;
class Standard_DimensionError;
class Standard_ConstructionError;
class math_Matrix;
//! 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.
class math_Householder
{
public:
DEFINE_STANDARD_ALLOC
//! 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.
Standard_EXPORT 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.
Standard_EXPORT 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 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.
Standard_EXPORT math_Householder(const math_Matrix& A, const math_Vector& B, const Standard_Real EPS = 1.0e-20);
//! Returns true if the computations are successful, otherwise returns false.
Standard_Boolean IsDone() 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.
void Value (math_Vector& sol, const Standard_Integer Index = 1) 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.
const math_Matrix& AllValues() const;
//! Prints informations on the current state of the object.
Standard_EXPORT void Dump (Standard_OStream& o) const;
protected:
//! This method is used internally for each constructor
//! above and can't be used directly.
Standard_EXPORT void Perform (const math_Matrix& A, const math_Matrix& B, const Standard_Real EPS) ;
private:
math_Matrix Sol;
math_Matrix Q;
Standard_Boolean Done;
Standard_Integer mylowerArow;
Standard_Integer myupperArow;
Standard_Integer mylowerAcol;
Standard_Integer myupperAcol;
};
#include <math_Householder.lxx>
#endif // _math_Householder_HeaderFile
|