0032137: Coding Rules - merge redundant .lxx files into header files within Package gp
[occt.git] / src / gp / gp_Elips2d.cxx
1 // Copyright (c) 1995-1999 Matra Datavision
2 // Copyright (c) 1999-2014 OPEN CASCADE SAS
3 //
4 // This file is part of Open CASCADE Technology software library.
5 //
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.
11 //
12 // Alternatively, this file may be used under the terms of Open CASCADE
13 // commercial license or contractual agreement.
14
15 //Modif JCV 10/01/91
16
17 #include <gp_Elips2d.hxx>
18
19 #include <gp_Ax2d.hxx>
20 #include <gp_Ax22d.hxx>
21 #include <gp_Pnt2d.hxx>
22 #include <gp_Trsf2d.hxx>
23 #include <gp_Vec2d.hxx>
24 #include <Standard_ConstructionError.hxx>
25
26 void gp_Elips2d::Coefficients (Standard_Real& A, 
27                                Standard_Real& B, 
28                                Standard_Real& C, 
29                                Standard_Real& D, 
30                                Standard_Real& E, 
31                                Standard_Real& F) const 
32 {
33   Standard_Real DMin = minorRadius * minorRadius;
34   Standard_Real DMaj = majorRadius * majorRadius;
35   if (DMin <= gp::Resolution() && DMaj <= gp::Resolution()) {
36     A = B = C = D = E = F = 0.0;
37   }
38   else {
39     gp_Trsf2d T;
40     T.SetTransformation (pos.XAxis());
41     Standard_Real T11 = T.Value (1, 1);
42     Standard_Real T12 = T.Value (1, 2);
43     Standard_Real T13 = T.Value (1, 3);
44     if (DMin <= gp::Resolution()) {
45       A = T11 * T11;    B = T12 * T12;   C = T11 * T12;
46       D = T11 * T13;    E = T12 * T13;   F = T13 * T13 - DMaj;
47     }
48     else {
49       Standard_Real T21 = T.Value (2, 1);
50       Standard_Real T22 = T.Value (2, 2);
51       Standard_Real T23 = T.Value (2, 3);
52       A = (T11 * T11 / DMaj) + (T21 * T21 / DMin);
53       B = (T12 * T12 / DMaj) + (T22 * T22 / DMin);
54       C = (T11 * T12 / DMaj) + (T21 * T22 / DMin);
55       D = (T11 * T13 / DMaj) + (T21 * T23 / DMin);
56       E = (T12 * T13 / DMaj) + (T22 * T23 / DMin);
57       F = (T13 * T13 / DMaj) + (T23 * T23 / DMin) - 1.0;
58     }
59   }
60 }
61
62 void gp_Elips2d::Mirror (const gp_Pnt2d& P)
63 { pos.Mirror(P); }
64
65 gp_Elips2d gp_Elips2d::Mirrored (const gp_Pnt2d& P) const  
66 {
67   gp_Elips2d E = *this;
68   E.pos.Mirror (P);
69   return E; 
70 }
71
72 void gp_Elips2d::Mirror (const gp_Ax2d& A)
73 { pos.Mirror(A); }
74
75 gp_Elips2d gp_Elips2d::Mirrored (const gp_Ax2d& A) const  
76 {
77   gp_Elips2d E = *this;
78   E.pos.Mirror (A);
79   return E; 
80 }
81