b311480e |
1 | // Copyright (c) 1995-1999 Matra Datavision |
973c2be1 |
2 | // Copyright (c) 1999-2014 OPEN CASCADE SAS |
b311480e |
3 | // |
973c2be1 |
4 | // This file is part of Open CASCADE Technology software library. |
b311480e |
5 | // |
d5f74e42 |
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 |
973c2be1 |
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. |
b311480e |
11 | // |
973c2be1 |
12 | // Alternatively, this file may be used under the terms of Open CASCADE |
13 | // commercial license or contractual agreement. |
b311480e |
14 | |
7fd59977 |
15 | //Modif JCV 10/01/91 |
16 | |
d5477f8c |
17 | #include <gp_Elips2d.hxx> |
18 | |
42cf5bc1 |
19 | #include <gp_Ax2d.hxx> |
20 | #include <gp_Ax22d.hxx> |
42cf5bc1 |
21 | #include <gp_Pnt2d.hxx> |
22 | #include <gp_Trsf2d.hxx> |
23 | #include <gp_Vec2d.hxx> |
24 | #include <Standard_ConstructionError.hxx> |
7fd59977 |
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 | |