0026937: Eliminate NO_CXX_EXCEPTION macro support
[occt.git] / src / GccAna / GccAna_Circ2dTanOnRad_5.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
16 #include <ElCLib.hxx>
17 #include <GccAna_Circ2dTanOnRad.hxx>
18 #include <GccEnt_BadQualifier.hxx>
19 #include <GccEnt_QualifiedCirc.hxx>
20 #include <GccEnt_QualifiedLin.hxx>
21 #include <gp_Circ2d.hxx>
22 #include <gp_Dir2d.hxx>
23 #include <gp_Lin2d.hxx>
24 #include <gp_Pnt2d.hxx>
25 #include <IntAna2d_AnaIntersection.hxx>
26 #include <IntAna2d_IntPoint.hxx>
27 #include <Standard_NegativeValue.hxx>
28 #include <Standard_OutOfRange.hxx>
29 #include <StdFail_NotDone.hxx>
30
31 //=========================================================================
32 //   Circle tangent to a point   Point1.                                  +
33 //          center on circle     OnCirc.                                  +
34 //          radius               Radius.                                  +
35 //                                                                        +
36 //  Initialize the table of solutions cirsol and all fields.              +
37 //  Eliminate cases not being the solution.                               +
38 //  Create the circle with center in Point1 of radius Radius.             +
39 //  Intersect this circle with OnCirc and obtain the center points        +
40 //  of found solutions.                                                   +
41 //  Create solutions cirsol.                                              +
42 //=========================================================================
43 GccAna_Circ2dTanOnRad::
44    GccAna_Circ2dTanOnRad (const gp_Pnt2d&     Point1    ,
45                           const gp_Circ2d&    OnCirc    ,
46                           const Standard_Real Radius    ,
47                           const Standard_Real Tolerance ):
48    cirsol(1,2)   ,
49    qualifier1(1,2) ,
50    TheSame1(1,2) ,
51    pnttg1sol(1,2),
52    pntcen3(1,2)  ,
53    par1sol(1,2)  ,
54    pararg1(1,2)  ,
55    parcen3(1,2)  
56 {
57
58    gp_Dir2d dirx(1.0,0.0);
59    Standard_Real Tol = Abs(Tolerance);
60    WellDone = Standard_False;
61    NbrSol = 0;
62    Standard_Real Roncirc = OnCirc.Radius();
63    Standard_Real dist1 = Point1.Distance(OnCirc.Location())-Roncirc;
64    Standard_Real dist2 = Point1.Distance(OnCirc.Location())+Roncirc;
65
66    if (Radius < 0.0) { throw Standard_NegativeValue(); }
67    else if ((dist1-Radius > Tol) || (Tol < Radius-dist2)) { 
68      WellDone = Standard_True; 
69    }
70    else {
71      Standard_Integer signe = 0;
72      if (Abs(dist1-Radius) < Tol) { signe = 1; }
73      else if (Abs(dist2-Radius) < Tol) { signe = -1; }
74      if (signe != 0) {
75        gp_Ax2d axe(gp_Pnt2d(OnCirc.Location().XY()-Roncirc*
76                          gp_Dir2d(OnCirc.Location().X()-signe*Point1.X(),
77                            OnCirc.Location().Y()-signe*Point1.Y()).XY()),dirx);
78        cirsol(1) = gp_Circ2d(axe,Radius);
79 //      ================================
80        qualifier1(1) = GccEnt_noqualifier;
81        TheSame1(1) = 0;
82        pnttg1sol(1) = Point1;
83        pntcen3(1) = cirsol(1).Location();
84        pararg1(1) = 0.0;
85        par1sol(1)=ElCLib::Parameter(cirsol(1),pnttg1sol(1));
86        parcen3(1) = ElCLib::Parameter(OnCirc,pntcen3(1));
87        WellDone = Standard_True;
88        NbrSol = 1;
89      }
90      else {
91        IntAna2d_AnaIntersection Intp(OnCirc,gp_Circ2d(gp_Ax2d(Point1,dirx),
92                                                       Radius));
93        if (Intp.IsDone()) {
94          if (!Intp.IsEmpty()) {
95            for (Standard_Integer i = 1 ; i <= Intp.NbPoints() ; i++) {
96              NbrSol++;
97              gp_Pnt2d Center(Intp.Point(i).Value());
98              cirsol(NbrSol) = gp_Circ2d(gp_Ax2d(Center,dirx),Radius);
99 //           =======================================================
100        qualifier1(1) = GccEnt_noqualifier;
101              TheSame1(1) = 0;
102              pnttg1sol(1) = Point1;
103              pntcen3(1) = cirsol(1).Location();
104              par1sol(1)=ElCLib::Parameter(cirsol(1),pnttg1sol(1));
105              parcen3(1) = ElCLib::Parameter(OnCirc,pntcen3(1));
106              pararg1(1) = 0.0;
107            }
108          }
109          WellDone = Standard_True;
110        }
111      }
112    }
113  }