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
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
|
// File: GCE2d_MakeArcOfCircle.cxx
// Created: Fri Oct 2 16:30:57 1992
// Author: Remi GILET
// <reg@topsn3>
#include <GCE2d_MakeArcOfCircle.ixx>
#include <gp_Vec2d.hxx>
#include <Geom2d_Circle.hxx>
#include <gce_MakeCirc2d.hxx>
#include <gce_MakeLin2d.hxx>
#include <ElCLib.hxx>
#include <StdFail_NotDone.hxx>
#include <IntAna2d_AnaIntersection.hxx>
GCE2d_MakeArcOfCircle::GCE2d_MakeArcOfCircle(const gp_Pnt2d& P1 ,
const gp_Pnt2d& P2 ,
const gp_Pnt2d& P3 )
{
gce_MakeCirc2d Cir = gce_MakeCirc2d(P1,P2,P3);
TheError = Cir.Status();
if (TheError == gce_Done) {
gp_Circ2d C(Cir.Value());
Standard_Real Alpha1 = ElCLib::Parameter(C,P1);
Standard_Real Alpha2 = ElCLib::Parameter(C,P3);
Handle(Geom2d_Circle) Circ = new Geom2d_Circle(C);
TheArc= new Geom2d_TrimmedCurve(Circ,Alpha1,Alpha2,Standard_True);
}
}
GCE2d_MakeArcOfCircle::GCE2d_MakeArcOfCircle(const gp_Pnt2d& P1 ,
const gp_Vec2d& V ,
const gp_Pnt2d& P2 )
{
Standard_Boolean Sense;
gp_Circ2d cir;
gp_Lin2d corde = gce_MakeLin2d(P1,P2);
gp_Dir2d dir(corde.Direction());
gp_Lin2d bis(gp_Pnt2d((P1.X()+P2.X())/2.,(P1.Y()+P2.Y())/2.),
gp_Dir2d(-dir.Y(),dir.X()));
gp_Lin2d norm(P1,gp_Dir2d(-V.Y(),V.X()));
TheError = gce_ConfusedPoints;
IntAna2d_AnaIntersection Intp(bis,norm);
if (Intp.IsDone())
{
if (!Intp.IsEmpty())
{
gp_Pnt2d center(Intp.Point(1).Value());
Standard_Real rad = (center.Distance(P1)+center.Distance(P2))/2.;
cir = gce_MakeCirc2d(center,rad);
TheError = gce_Done;
}
}
if (TheError == gce_Done) {
Standard_Real Alpha1 = ElCLib::Parameter(cir,P1);
Standard_Real Alpha2 = ElCLib::Parameter(cir,P2);
Handle(Geom2d_Circle) Circ = new Geom2d_Circle(cir);
gp_Vec2d vv(dir);
Standard_Real cross = V^vv;
Sense = cross > 0.;
TheArc= new Geom2d_TrimmedCurve(Circ,Alpha1,Alpha2,Sense);
}
}
GCE2d_MakeArcOfCircle::GCE2d_MakeArcOfCircle(const gp_Circ2d& Circ ,
const gp_Pnt2d& P1 ,
const gp_Pnt2d& P2 ,
const Standard_Boolean Sense )
{
Standard_Real Alpha1 = ElCLib::Parameter(Circ,P1);
Standard_Real Alpha2 = ElCLib::Parameter(Circ,P2);
Handle(Geom2d_Circle) C = new Geom2d_Circle(Circ);
TheArc= new Geom2d_TrimmedCurve(C,Alpha1,Alpha2,Sense);
TheError = gce_Done;
}
GCE2d_MakeArcOfCircle::GCE2d_MakeArcOfCircle(const gp_Circ2d& Circ ,
const gp_Pnt2d& P ,
const Standard_Real Alpha ,
const Standard_Boolean Sense )
{
Standard_Real Alphafirst = ElCLib::Parameter(Circ,P);
Handle(Geom2d_Circle) C = new Geom2d_Circle(Circ);
TheArc= new Geom2d_TrimmedCurve(C,Alphafirst,Alpha,Sense);
TheError = gce_Done;
}
GCE2d_MakeArcOfCircle::GCE2d_MakeArcOfCircle(const gp_Circ2d& Circ ,
const Standard_Real Alpha1 ,
const Standard_Real Alpha2 ,
const Standard_Boolean Sense )
{
Handle(Geom2d_Circle) C = new Geom2d_Circle(Circ);
TheArc= new Geom2d_TrimmedCurve(C,Alpha1,Alpha2,Sense);
TheError = gce_Done;
}
const Handle(Geom2d_TrimmedCurve)& GCE2d_MakeArcOfCircle::Value() const
{
StdFail_NotDone_Raise_if(!TheError == gce_Done,"");
return TheArc;
}
const Handle(Geom2d_TrimmedCurve)& GCE2d_MakeArcOfCircle::Operator() const
{
return Value();
}
GCE2d_MakeArcOfCircle::operator Handle(Geom2d_TrimmedCurve) () const
{
return Value();
}
|