1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
|
// file GccAna_Circ2dTanOnRad_5.cxx, REG 08/07/91
#include <GccAna_Circ2dTanOnRad.jxx>
#include <ElCLib.hxx>
#include <IntAna2d_AnaIntersection.hxx>
#include <IntAna2d_IntPoint.hxx>
#include <Standard_NegativeValue.hxx>
#include <Standard_OutOfRange.hxx>
#include <gp_Dir2d.hxx>
//=========================================================================
// Cercle tangent a un point Point1. +
// centre sur un cercle OnCirc. +
// de rayon Radius. +
// +
// On initialise le tableau de solutions cirsol ainsi que tous les +
// champs. +
// On elimine les cas ne presentant pas de solution. +
// On cree le cercle centre en Point1 de rayon Radius. +
// On intersecte ce cercle avec OnCirc et on obtient les points de +
// centre des solutions recherchees. +
// On cree ces solutions cirsol. +
//=========================================================================
GccAna_Circ2dTanOnRad::
GccAna_Circ2dTanOnRad (const gp_Pnt2d& Point1 ,
const gp_Circ2d& OnCirc ,
const Standard_Real Radius ,
const Standard_Real Tolerance ):
cirsol(1,2) ,
qualifier1(1,2) ,
TheSame1(1,2) ,
pnttg1sol(1,2),
pntcen3(1,2) ,
par1sol(1,2) ,
pararg1(1,2) ,
parcen3(1,2)
{
gp_Dir2d dirx(1.0,0.0);
Standard_Real Tol = Abs(Tolerance);
WellDone = Standard_False;
NbrSol = 0;
Standard_Real Roncirc = OnCirc.Radius();
Standard_Real dist1 = Point1.Distance(OnCirc.Location())-Roncirc;
Standard_Real dist2 = Point1.Distance(OnCirc.Location())+Roncirc;
if (Radius < 0.0) { Standard_NegativeValue::Raise(); }
else if ((dist1-Radius > Tol) || (Tol < Radius-dist2)) {
WellDone = Standard_True;
}
else {
Standard_Integer signe = 0;
if (Abs(dist1-Radius) < Tol) { signe = 1; }
else if (Abs(dist2-Radius) < Tol) { signe = -1; }
if (signe != 0) {
gp_Ax2d axe(gp_Pnt2d(OnCirc.Location().XY()-Roncirc*
gp_Dir2d(OnCirc.Location().X()-signe*Point1.X(),
OnCirc.Location().Y()-signe*Point1.Y()).XY()),dirx);
cirsol(1) = gp_Circ2d(axe,Radius);
// ================================
qualifier1(1) = GccEnt_noqualifier;
TheSame1(1) = 0;
pnttg1sol(1) = Point1;
pntcen3(1) = cirsol(1).Location();
pararg1(1) = 0.0;
par1sol(1)=ElCLib::Parameter(cirsol(1),pnttg1sol(1));
parcen3(1) = ElCLib::Parameter(OnCirc,pntcen3(1));
WellDone = Standard_True;
NbrSol = 1;
}
else {
IntAna2d_AnaIntersection Intp(OnCirc,gp_Circ2d(gp_Ax2d(Point1,dirx),
Radius));
if (Intp.IsDone()) {
if (!Intp.IsEmpty()) {
for (Standard_Integer i = 1 ; i <= Intp.NbPoints() ; i++) {
NbrSol++;
gp_Pnt2d Center(Intp.Point(i).Value());
cirsol(NbrSol) = gp_Circ2d(gp_Ax2d(Center,dirx),Radius);
// =======================================================
qualifier1(1) = GccEnt_noqualifier;
TheSame1(1) = 0;
pnttg1sol(1) = Point1;
pntcen3(1) = cirsol(1).Location();
par1sol(1)=ElCLib::Parameter(cirsol(1),pnttg1sol(1));
parcen3(1) = ElCLib::Parameter(OnCirc,pntcen3(1));
pararg1(1) = 0.0;
}
}
WellDone = Standard_True;
}
}
}
}
|
