1 // Copyright (c) 1997-1999 Matra Datavision
2 // Copyright (c) 1999-2014 OPEN CASCADE SAS
4 // This file is part of Open CASCADE Technology software library.
6 // This library is free software; you can redistribute it and/or modify it under
7 // the terms of the GNU Lesser General Public License version 2.1 as published
8 // by the Free Software Foundation, with special exception defined in the file
9 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
10 // distribution for complete text of the license and disclaimer of any warranty.
12 // Alternatively, this file may be used under the terms of Open CASCADE
13 // commercial license or contractual agreement.
16 #define No_Standard_RangeError
17 #define No_Standard_OutOfRange
18 #define No_Standard_DimensionError
22 #include <math_Jacobi.hxx>
23 #include <math_Matrix.hxx>
24 #include <math_NotSquare.hxx>
25 #include <math_Recipes.hxx>
26 #include <StdFail_NotDone.hxx>
28 math_Jacobi::math_Jacobi(const math_Matrix& A) : AA(1, A.RowNumber(),
30 EigenValues(1, A.RowNumber()),
31 EigenVectors(1, A.RowNumber(),
34 math_NotSquare_Raise_if(A.RowNumber() != A.ColNumber(), " ");
37 Standard_Integer Error = Jacobi(AA, EigenValues, EigenVectors, NbRotations);
42 Done = Standard_False;
46 void math_Jacobi::Dump(Standard_OStream& o) const {
50 o << " Status = Done \n";
51 o << " The eigenvalues vector is: " << EigenValues << std::endl;
54 o << "Status = not Done \n";