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
|
// File: IntCurveSurface_Intersection.cxx
// Created: Fri Apr 16 10:21:36 1993
// Author: Laurent BUCHARD
// <lbr@phobox>
#include <IntCurveSurface_Intersection.ixx>
#include <IntCurveSurface_IntersectionPoint.hxx>
#include <IntCurveSurface_IntersectionSegment.hxx>
#include <IntCurveSurface_TransitionOnCurve.hxx>
#define PARAMEQUAL(a,b) (Abs((a)-(b))< (1e-8))
//================================================================================
IntCurveSurface_Intersection::IntCurveSurface_Intersection(): done(Standard_False) {
}
//================================================================================
Standard_Boolean IntCurveSurface_Intersection::IsDone() const { return(done); }
//================================================================================
Standard_Integer IntCurveSurface_Intersection::NbPoints() const {
if (!done) {StdFail_NotDone::Raise();}
return lpnt.Length();
}
//================================================================================
Standard_Integer IntCurveSurface_Intersection::NbSegments() const {
if (!done) {StdFail_NotDone::Raise();}
return lseg.Length();
}
//================================================================================
const IntCurveSurface_IntersectionPoint& IntCurveSurface_Intersection::Point( const Standard_Integer N) const {
if (!done) {StdFail_NotDone::Raise();}
return lpnt.Value(N);
}
//================================================================================
const IntCurveSurface_IntersectionSegment& IntCurveSurface_Intersection::Segment( const Standard_Integer N) const {
if (!done) {StdFail_NotDone::Raise();}
return lseg.Value(N);
}
//================================================================================
void IntCurveSurface_Intersection::SetValues(const IntCurveSurface_Intersection& Other) {
if(Other.done) {
lseg.Clear();
lpnt.Clear();
Standard_Integer N = Other.lpnt.Length();
Standard_Integer i ;
for( i=1; i<= N ; i++) {
lpnt.Append(Other.lpnt.Value(i));
}
N = Other.lseg.Length();
for(i=1; i<= N ; i++) {
lseg.Append(Other.lseg.Value(i));
}
done=Standard_True;
}
else {
done=Standard_False;
}
}
//================================================================================
void IntCurveSurface_Intersection::Append(const IntCurveSurface_Intersection& Other,
// const Standard_Real a,
const Standard_Real ,
// const Standard_Real b)
const Standard_Real )
{
Standard_Integer i,ni;
if(Other.done) {
ni = Other.lpnt.Length();
for(i=1;i<=ni;i++) { Append(Other.Point(i)); }
ni = Other.lseg.Length();
for(i=1;i<=ni;i++) { Append(Other.Segment(i)); }
}
}
//================================================================================
void IntCurveSurface_Intersection::Append(const IntCurveSurface_IntersectionPoint& OtherPoint) {
Standard_Integer i,ni;
Standard_Real anu,anv,anw,u,v,w;
IntCurveSurface_TransitionOnCurve TrOnCurve,anTrOnCurve;
gp_Pnt P,anP;
ni = lpnt.Length();
for(i=1;i<=ni;i++) {
OtherPoint.Values(P,u,v,w,TrOnCurve);
lpnt(i).Values(anP,anu,anv,anw,anTrOnCurve);
if(PARAMEQUAL(u,anu)) {
if(PARAMEQUAL(v,anv)) {
if(PARAMEQUAL(w,anw)) {
if(anTrOnCurve==TrOnCurve) {
return;
}
}
}
}
}
lpnt.Append(OtherPoint);
}
//================================================================================
void IntCurveSurface_Intersection::Append(const IntCurveSurface_IntersectionSegment& OtherSegment) {
lseg.Append(OtherSegment);
}
//================================================================================
void IntCurveSurface_Intersection::ResetFields() {
if(done) {
lseg.Clear();
lpnt.Clear();
done=Standard_False;
}
}
//================================================================================
void IntCurveSurface_Intersection::Dump() const {
if(done) {
Standard_Integer i,ni;
ni = lpnt.Length();
for(i=1;i<=ni;i++) { Point(i).Dump(); }
ni = lseg.Length();
for(i=1;i<=ni;i++) { Segment(i).Dump(); }
}
else {
cout<<" Intersection NotDone"<<endl;
}
}
|
