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
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
|
// 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
#ifndef _Standard_HeaderFile
#include <Standard.hxx>
#endif
#ifndef _Standard_Macro_HeaderFile
#include <Standard_Macro.hxx>
#endif
#ifndef _math_Matrix_HeaderFile
#include <math_Matrix.hxx>
#endif
#ifndef _Standard_Boolean_HeaderFile
#include <Standard_Boolean.hxx>
#endif
#ifndef _Standard_Integer_HeaderFile
#include <Standard_Integer.hxx>
#endif
#ifndef _Standard_Real_HeaderFile
#include <Standard_Real.hxx>
#endif
#ifndef _Standard_OStream_HeaderFile
#include <Standard_OStream.hxx>
#endif
class StdFail_NotDone;
class Standard_OutOfRange;
class Standard_DimensionError;
class Standard_ConstructionError;
class math_Matrix;
class math_Vector;
//! This class implements the least square solution of a set of <br>
//! linear equations of m unknowns (n >= m) using the Householder <br>
//! method. It solves A.X = B. <br>
//! This algorithm has more numerical stability than <br>
//! GaussLeastSquare but is longer. <br>
//! It must be used if the matrix is singular or nearly singular. <br>
//! It is about 16% longer than GaussLeastSquare if there is only <br>
//! one member B to solve. <br>
//! It is about 30% longer if there are twenty B members to solve. <br>
class math_Householder {
public:
void* operator new(size_t,void* anAddress)
{
return anAddress;
}
void* operator new(size_t size)
{
return Standard::Allocate(size);
}
void operator delete(void *anAddress)
{
if (anAddress) Standard::Free((Standard_Address&)anAddress);
}
//! Given an input matrix A with n>= m, given an input matrix B <br>
//! this constructor performs the least square resolution of <br>
//! the set of linear equations A.X = B for each column of B. <br>
//! If a column norm is less than EPS, the resolution can't <br>
//! be done. <br>
//! Exception DimensionError is raised if the row number of B <br>
//! is different from the A row number. <br>
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 <br>
//! this constructor performs the least square resolution of <br>
//! the set of linear equations A.X = B for each column of B. <br>
//! If a column norm is less than EPS, the resolution can't <br>
//! be done. <br>
//! Exception DimensionError is raised if the row number of B <br>
//! is different from the A row number. <br>
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 <br>
//! this constructor performs the least square resolution of <br>
//! the set of linear equations A.X = B. <br>
//! If a column norm is less than EPS, the resolution can't <br>
//! be done. <br>
//! Exception DimensionError is raised if the length of B <br>
//! is different from the A row number. <br>
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. <br>
Standard_Boolean IsDone() const;
//! Given the integer Index, this routine returns the <br>
//! corresponding least square solution sol. <br>
//! Exception NotDone is raised if the resolution has not be <br>
//! done. <br>
//! Exception OutOfRange is raised if Index <=0 or <br>
//! Index is more than the number of columns of B. <br>
void Value(math_Vector& sol,const Standard_Integer Index = 1) const;
//! Returns the matrix sol of all the solutions of the system <br>
//! A.X = B. <br>
//! Exception NotDone is raised is the resolution has not be <br>
//! done. <br>
const math_Matrix& AllValues() const;
//! Prints informations on the current state of the object. <br>
Standard_EXPORT void Dump(Standard_OStream& o) const;
protected:
//! This method is used internally for each constructor <br>
//! above and can't be used directly. <br>
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>
// other Inline functions and methods (like "C++: function call" methods)
#endif
|
