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
|
// File: BRepClass_Intersector3d.cdl
// Created: Fri Apr 1 09:14:28 1994
// Author: Laurent BUCHARD
// <lbr@fuegox>
//-Copyright: Matra Datavision 1994
// Modified by skv - Fri Mar 4 12:07:34 2005 OCC7966
#include <BRepClass3d_Intersector3d.ixx>
#include <IntCurveSurface_IntersectionPoint.hxx>
#include <gp_Lin.hxx>
#include <TopoDS_Face.hxx>
#include <TopAbs.hxx>
#include <IntCurveSurface_HInter.hxx>
#include <BRepAdaptor_HSurface.hxx>
#include <Geom_Line.hxx>
#include <gp_Pnt2d.hxx>
#include <BRepClass_FaceClassifier.hxx>
#include <GeomAdaptor_Curve.hxx>
#include <GeomAdaptor_HCurve.hxx>
#include <BRepAdaptor_HSurface.hxx>
//============================================================================
BRepClass3d_Intersector3d::BRepClass3d_Intersector3d()
: done(Standard_False),hasapoint(Standard_False)
{
}
//============================================================================
void BRepClass3d_Intersector3d::Perform(const gp_Lin& L,
const Standard_Real /*Prm*/,
const Standard_Real Tol,
const TopoDS_Face& Face) {
static IntCurveSurface_HInter HICS;
static BRepAdaptor_Surface surface;
static BRepClass_FaceClassifier classifier2d;
Handle(Geom_Line) geomline = new Geom_Line(L);
GeomAdaptor_Curve LL(geomline);
surface.Initialize(Face,Standard_True);
Standard_Boolean IsUPer, IsVPer;
Standard_Real uperiod=0, vperiod=0;
if ((IsUPer = surface.IsUPeriodic()) == Standard_True)
uperiod = surface.UPeriod();
if ((IsVPer = surface.IsVPeriodic()) == Standard_True)
vperiod = surface.VPeriod();
Standard_Real U1, U2, V1, V2;
U1 = surface.FirstUParameter();
U2 = surface.LastUParameter();
V1 = surface.FirstVParameter();
V2 = surface.LastVParameter();
//--
Handle(GeomAdaptor_HCurve) HLL = new GeomAdaptor_HCurve(LL);
Handle(BRepAdaptor_HSurface) Hsurface = new BRepAdaptor_HSurface(surface);
//--
HICS.Perform(HLL,Hsurface);
W=RealLast();
if(HICS.IsDone()) {
for(Standard_Integer index=HICS.NbPoints(); index>=1; index--) {
gp_Pnt2d Puv(HICS.Point(index).U(),HICS.Point(index).V());
Standard_Integer N1 = 0;
Standard_Integer N2 = 0;
Standard_Real X = Puv.X();
Standard_Real Y = Puv.Y();
if(IsUPer) {
if(X > U2) {
N1 = RealToInt( (X - U1) / uperiod );
}
if(X < U1) {
N1 = RealToInt( (X - U2) / uperiod );
}
Puv.SetX(X - uperiod * N1);
}
if(IsVPer) {
if(Y > V2) {
N2 = RealToInt ( (Y - V1) / vperiod );
}
if(Y < V1) {
N2 = RealToInt ( (Y - V2) / vperiod );
}
Puv.SetY(Y - vperiod * N2);
}
classifier2d.Perform(Face,Puv,Tol);
TopAbs_State currentstate = classifier2d.State();
if(currentstate==TopAbs_IN || currentstate==TopAbs_ON) {
const IntCurveSurface_IntersectionPoint& HICSPoint = HICS.Point(index);
Standard_Real HICSW = HICSPoint.W();
// Modified by skv - Fri Mar 4 12:07:34 2005 OCC7966 Begin
if((W > HICSW) && (HICSW>-Tol)) {
// if(W > HICSW) {
// Modified by skv - Fri Mar 4 12:07:34 2005 OCC7966 End
hasapoint = Standard_True;
U = HICSPoint.U();
V = HICSPoint.V();
W = HICSW;
transition = HICSPoint.Transition();
pnt = HICSPoint.Pnt();
state = currentstate;
face = Face;
if(Face.Orientation()==TopAbs_REVERSED) {
if(transition == IntCurveSurface_In)
transition = IntCurveSurface_Out;
else
transition = IntCurveSurface_In;
}
}
} //-- classifier state is IN or ON
done = Standard_True;
} //-- Loop on Intersection points.
} //-- HICS.IsDone()
}
|
