3be08b94da07fd2c8e28380bb89d7e8d276b45ef
[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_Ax2d.hxx>
18 #include <gp_Ax22d.hxx>
19 #include <gp_Elips2d.hxx>
20 #include <gp_Pnt2d.hxx>
21 #include <gp_Trsf2d.hxx>
22 #include <gp_Vec2d.hxx>
23 #include <Standard_ConstructionError.hxx>
24
25 void gp_Elips2d::Coefficients (Standard_Real& A, 
26                                Standard_Real& B, 
27                                Standard_Real& C, 
28                                Standard_Real& D, 
29                                Standard_Real& E, 
30                                Standard_Real& F) const 
31 {
32   Standard_Real DMin = minorRadius * minorRadius;
33   Standard_Real DMaj = majorRadius * majorRadius;
34   if (DMin <= gp::Resolution() && DMaj <= gp::Resolution()) {
35     A = B = C = D = E = F = 0.0;
36   }
37   else {
38     gp_Trsf2d T;
39     T.SetTransformation (pos.XAxis());
40     Standard_Real T11 = T.Value (1, 1);
41     Standard_Real T12 = T.Value (1, 2);
42     Standard_Real T13 = T.Value (1, 3);
43     if (DMin <= gp::Resolution()) {
44       A = T11 * T11;    B = T12 * T12;   C = T11 * T12;
45       D = T11 * T13;    E = T12 * T13;   F = T13 * T13 - DMaj;
46     }
47     else {
48       Standard_Real T21 = T.Value (2, 1);
49       Standard_Real T22 = T.Value (2, 2);
50       Standard_Real T23 = T.Value (2, 3);
51       A = (T11 * T11 / DMaj) + (T21 * T21 / DMin);
52       B = (T12 * T12 / DMaj) + (T22 * T22 / DMin);
53       C = (T11 * T12 / DMaj) + (T21 * T22 / DMin);
54       D = (T11 * T13 / DMaj) + (T21 * T23 / DMin);
55       E = (T12 * T13 / DMaj) + (T22 * T23 / DMin);
56       F = (T13 * T13 / DMaj) + (T23 * T23 / DMin) - 1.0;
57     }
58   }
59 }
60
61 void gp_Elips2d::Mirror (const gp_Pnt2d& P)
62 { pos.Mirror(P); }
63
64 gp_Elips2d gp_Elips2d::Mirrored (const gp_Pnt2d& P) const  
65 {
66   gp_Elips2d E = *this;
67   E.pos.Mirror (P);
68   return E; 
69 }
70
71 void gp_Elips2d::Mirror (const gp_Ax2d& A)
72 { pos.Mirror(A); }
73
74 gp_Elips2d gp_Elips2d::Mirrored (const gp_Ax2d& A) const  
75 {
76   gp_Elips2d E = *this;
77   E.pos.Mirror (A);
78   return E; 
79 }
80