/usr/include/oce/math_Gauss.hxx is in liboce-foundation-dev 0.17.2-2.
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 | // 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_Gauss_HeaderFile
#define _math_Gauss_HeaderFile
#include <Standard.hxx>
#include <Standard_DefineAlloc.hxx>
#include <Standard_Macro.hxx>
#include <Standard_Boolean.hxx>
#include <math_Matrix.hxx>
#include <math_IntegerVector.hxx>
#include <Standard_Real.hxx>
#include <math_Vector.hxx>
#include <Standard_OStream.hxx>
class math_NotSquare;
class Standard_DimensionError;
class StdFail_NotDone;
class math_Matrix;
//! This class implements the Gauss LU decomposition (Crout algorithm)
//! with partial pivoting (rows interchange) of a square matrix and
//! the different possible derived calculation :
//! - solution of a set of linear equations.
//! - inverse of a matrix.
//! - determinant of a matrix.
class math_Gauss
{
public:
DEFINE_STANDARD_ALLOC
//! Given an input n X n matrix A this constructor performs its LU
//! decomposition with partial pivoting (interchange of rows).
//! This LU decomposition is stored internally and may be used to
//! do subsequent calculation.
//! If the largest pivot found is less than MinPivot the matrix A is
//! considered as singular.
//! Exception NotSquare is raised if A is not a square matrix.
Standard_EXPORT math_Gauss(const math_Matrix& A, const Standard_Real MinPivot = 1.0e-20);
//! Returns true if the computations are successful, otherwise returns false
Standard_Boolean IsDone() const;
//! Given the input Vector B this routine returns the solution X of the set
//! of linear equations A . X = B.
//! Exception NotDone is raised if the decomposition of A was not done
//! successfully.
//! Exception DimensionError is raised if the range of B is not
//! equal to the number of rows of A.
Standard_EXPORT void Solve (const math_Vector& B, math_Vector& X) const;
//! Given the input Vector B this routine solves the set of linear
//! equations A . X = B. B is replaced by the vector solution X.
//! Exception NotDone is raised if the decomposition of A was not done
//! successfully.
//! Exception DimensionError is raised if the range of B is not
//! equal to the number of rows of A.
Standard_EXPORT void Solve (math_Vector& B) const;
//! This routine returns the value of the determinant of the previously LU
//! decomposed matrix A.
//! Exception NotDone may be raised if the decomposition of A was not done
//! successfully, zero is returned if the matrix A was considered as singular.
Standard_EXPORT Standard_Real Determinant() const;
//! This routine outputs Inv the inverse of the previously LU decomposed
//! matrix A.
//! Exception DimensionError is raised if the ranges of B are not
//! equal to the ranges of A.
Standard_EXPORT void Invert (math_Matrix& Inv) const;
//! Prints on the stream o information on the current state
//! of the object.
//! Is used to redefine the operator <<.
Standard_EXPORT void Dump (Standard_OStream& o) const;
protected:
Standard_Boolean Singular;
math_Matrix LU;
math_IntegerVector Index;
Standard_Real D;
private:
Standard_Boolean Done;
};
#include <math_Gauss.lxx>
#endif // _math_Gauss_HeaderFile
|