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
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
|
// File: ChFiKPart_ComputeData_Fcts.cxx
// Created: Mon May 22 13:50:25 1995
// Author: Laurent BOURESCHE
// <lbo@phylox>
#include <ChFiKPart_ComputeData_Fcts.hxx>
#include <TColStd_Array1OfReal.hxx>
#include <TColStd_Array1OfInteger.hxx>
#include <ProjLib_ProjectedCurve.hxx>
#include <Geom2d_BezierCurve.hxx>
#include <Geom2d_BSplineCurve.hxx>
#include <Geom2d_Line.hxx>
#include <TColgp_Array1OfPnt2d.hxx>
#include <TopOpeBRepDS_Curve.hxx>
#include <TopOpeBRepDS_Surface.hxx>
#include <Standard_NotImplemented.hxx>
#include <GeomAdaptor_HCurve.hxx>
#include <GeomAdaptor_HSurface.hxx>
//=======================================================================
//function : InPeriod
//purpose :
//=======================================================================
Standard_Real ChFiKPart_InPeriod(const Standard_Real U,
const Standard_Real UFirst,
const Standard_Real ULast,
const Standard_Real Eps)
{
Standard_Real u = U, period = ULast - UFirst;
while (Eps < (UFirst-u)) u += period;
while (Eps > (ULast -u)) u -= period;
if ( u < UFirst) u = UFirst;
return u;
}
//=======================================================================
//function : PCurve
//purpose : Calculate a straight line in form of BSpline to guarantee
// the parameters.
//=======================================================================
Handle(Geom2d_BSplineCurve) ChFiKPart_PCurve(const gp_Pnt2d& UV1,
const gp_Pnt2d& UV2,
const Standard_Real Pardeb,
const Standard_Real Parfin)
{
TColgp_Array1OfPnt2d p(1,2);
TColStd_Array1OfReal k(1,2);
TColStd_Array1OfInteger m(1,2);
m.Init(2);
k(1) = Pardeb;
k(2) = Parfin;
p(1) = UV1;
p(2) = UV2;
Handle(Geom2d_BSplineCurve) Pcurv = new Geom2d_BSplineCurve(p,k,m,1);
return Pcurv;
}
//=======================================================================
//function : ProjPC
//purpose : For spherical corners the contours which of are not
// isos the circle is projected.
//=======================================================================
void ChFiKPart_ProjPC(const GeomAdaptor_Curve& Cg,
const GeomAdaptor_Surface& Sg,
Handle(Geom2d_Curve)& Pcurv)
{
if (Sg.GetType() < GeomAbs_BezierSurface) {
Handle(GeomAdaptor_HCurve) HCg = new GeomAdaptor_HCurve(Cg);
Handle(GeomAdaptor_HSurface) HSg = new GeomAdaptor_HSurface(Sg);
ProjLib_ProjectedCurve Projc (HSg,HCg);
switch (Projc.GetType()) {
case GeomAbs_Line :
{
Pcurv = new Geom2d_Line(Projc.Line());
}
break;
case GeomAbs_BezierCurve :
{
Handle(Geom2d_BezierCurve) BezProjc = Projc.Bezier();
TColgp_Array1OfPnt2d TP(1,BezProjc->NbPoles());
if (BezProjc->IsRational()) {
TColStd_Array1OfReal TW(1,BezProjc->NbPoles());
BezProjc->Poles(TP);
BezProjc->Weights(TW);
Pcurv = new Geom2d_BezierCurve(TP,TW);
}
else {
BezProjc->Poles(TP);
Pcurv = new Geom2d_BezierCurve(TP);
}
}
break;
#if 0
TColgp_Array1OfPnt2d TP(1,Projc.NbPoles());
if (Projc.IsRational()) {
TColStd_Array1OfReal TW(1,Projc.NbPoles());
Projc.PolesAndWeights(TP,TW);
Pcurv = new Geom2d_BezierCurve(TP,TW);
}
else {
Projc.Poles(TP);
Pcurv = new Geom2d_BezierCurve(TP);
}
}
break;
#endif
case GeomAbs_BSplineCurve :
{
Handle(Geom2d_BSplineCurve) BspProjc = Projc.BSpline();
TColgp_Array1OfPnt2d TP(1,BspProjc->NbPoles());
TColStd_Array1OfReal TK(1,BspProjc->NbKnots());
TColStd_Array1OfInteger TM(1,BspProjc->NbKnots());
BspProjc->Knots(TK);
BspProjc->Multiplicities(TM);
if (BspProjc->IsRational()) {
TColStd_Array1OfReal TW(1,BspProjc->NbPoles());
BspProjc->Poles(TP);
BspProjc->Weights(TW);
Pcurv = new Geom2d_BSplineCurve(TP,TW,TK,TM,BspProjc->Degree());
}
else {
BspProjc->Poles(TP);
Pcurv = new Geom2d_BSplineCurve(TP,TK,TM,BspProjc->Degree());
}
}
break;
#if 0
TColgp_Array1OfPnt2d TP(1,Projc.NbPoles());
TColStd_Array1OfReal TK(1,Projc.NbKnots());
TColStd_Array1OfInteger TM(1,Projc.NbKnots());
Projc.KnotsAndMultiplicities(TK,TM);
if (Projc.IsRational()) {
TColStd_Array1OfReal TW(1,Projc.NbPoles());
Projc.PolesAndWeights(TP,TW);
Pcurv = new Geom2d_BSplineCurve(TP,TW,TK,TM,Projc.Degree());
}
else {
Projc.Poles(TP);
Pcurv = new Geom2d_BSplineCurve(TP,TK,TM,Projc.Degree());
}
}
break;
#endif
default :
Standard_NotImplemented::Raise("failed approximation of the pcurve ");
}
}
else {
Standard_NotImplemented::Raise("approximate pcurve on the left surface");
}
}
//=======================================================================
//function : IndexCurveInDS
//purpose : Place a Curve in the DS and return its index.
//=======================================================================
Standard_Integer ChFiKPart_IndexCurveInDS(const Handle(Geom_Curve)& C,
TopOpeBRepDS_DataStructure& DStr)
{
return DStr.AddCurve(TopOpeBRepDS_Curve(C,0.));
}
//=======================================================================
//function : IndexSurfaceInDS
//purpose : Place a Surface in the DS and return its index.
//=======================================================================
Standard_Integer ChFiKPart_IndexSurfaceInDS(const Handle(Geom_Surface)& S,
TopOpeBRepDS_DataStructure& DStr)
{
return DStr.AddSurface(TopOpeBRepDS_Surface(S,0.));
}
|
