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
|
// File: Geom2dHatch_Intersector.cxx
// Created: Wed Mar 23 11:29:17 1994
// Author: Jean Marc LACHAUME
// <jml@phylox>
#include <Geom2dHatch_Intersector.ixx>
#include <ElCLib.hxx>
#include <Geom2d_Line.hxx>
#include <Precision.hxx>
#include <Geom2dLProp_CLProps2d.hxx>
//=======================================================================
//function : Geom2dHatch_Intersector
//purpose :
//=======================================================================
Geom2dHatch_Intersector::Geom2dHatch_Intersector() :
myConfusionTolerance(0.0),
myTangencyTolerance(0.0)
{
}
//=======================================================================
//function : Perform
//purpose :
//=======================================================================
void Geom2dHatch_Intersector::Perform(const gp_Lin2d& L,
const Standard_Real P,
const Standard_Real Tol,
const Geom2dAdaptor_Curve& C)
{
//Standard_Real pfbid,plbid;
IntRes2d_Domain DL;
if(P!=RealLast())
DL.SetValues(L.Location(),0.,Tol,ElCLib::Value(P,L),P,Tol);
else
DL.SetValues(L.Location(),0.,Tol,Standard_True);
IntRes2d_Domain DE(C.Value(C.FirstParameter()),
C.FirstParameter(),Precision::PIntersection(),
C.Value(C.LastParameter()),
C.LastParameter(),Precision::PIntersection());
Handle(Geom2d_Line) GL= new Geom2d_Line(L);
Geom2dAdaptor_Curve CGA(GL);
void *ptrpoureviterlesproblemesdeconst = (void *)(&C);
Geom2dInt_GInter Inter(CGA,
DL,
*((Geom2dAdaptor_Curve *)ptrpoureviterlesproblemesdeconst),
DE,
Precision::PConfusion(),
Precision::PIntersection());
this->SetValues(Inter);
}
//=======================================================================
//function : LocalGeometry
//purpose :
//=======================================================================
void Geom2dHatch_Intersector::LocalGeometry(const Geom2dAdaptor_Curve& E,
const Standard_Real U,
gp_Dir2d& Tang,
gp_Dir2d& Norm,
Standard_Real& C) const
{
//Standard_Real f,l;
Geom2dLProp_CLProps2d Prop(E.Curve(),U,2,Precision::PConfusion());
if(!Prop.IsTangentDefined()) return;
Prop.Tangent(Tang);
C = Prop.Curvature();
if (C > Precision::PConfusion() && C<RealLast())
Prop.Normal(Norm);
else
Norm.SetCoord(Tang.Y(),-Tang.X());
}
|
