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
|
// File : BRepBlend_BlendTool.cxx
// Created : Mon Dec 6 09:49:19 1993
// Author : Jacques GOUSSARD
// Copyright : Matra Datavision 1993
#include <BRepBlend_BlendTool.ixx>
#include <BRepClass_FaceClassifier.hxx>
#include <Extrema_EPCOfExtPC2d.hxx>
#include <Extrema_POnCurv2d.hxx>
#include <TopoDS_Vertex.hxx>
#include <TopoDS_Edge.hxx>
#include <BRep_Tool.hxx>
#include <Geom2d_Line.hxx>
#include <Geom2dAdaptor_Curve.hxx>
#include <gp_Dir2d.hxx>
#include <gp_Vec2d.hxx>
#include <Geom2dInt_GInter.hxx>
#include <IntRes2d_IntersectionPoint.hxx>
#include <BRepBlend_HCurve2dTool.hxx>
//=======================================================================
//function : Project
//purpose : Projection orthogonal d'un point sur une courbe
// pmn 8/10/98 : On retourne toujours une distance. (BUC60360)
//=======================================================================
Standard_Boolean BRepBlend_BlendTool::Project(const gp_Pnt2d& P,
const Handle(Adaptor3d_HSurface)&,
const Handle(Adaptor2d_HCurve2d)& C,
Standard_Real& Paramproj,
Standard_Real& Dist)
{
Paramproj = BRepBlend_HCurve2dTool::FirstParameter(C);
gp_Pnt2d P2d;
BRepBlend_HCurve2dTool::D0(C, Paramproj, P2d);
Dist = P2d.Distance(P);
const Standard_Real t = BRepBlend_HCurve2dTool::LastParameter(C);
BRepBlend_HCurve2dTool::D0(C, t, P2d);
if (P2d.Distance(P) < Dist) {
Paramproj = t;
Dist = P2d.Distance(P);
}
const Standard_Real epsX = 1.e-8;
const Standard_Integer Nbu = 20;
const Standard_Real Tol = 1.e-5;
Extrema_EPCOfExtPC2d extrema(P, C->Curve2d(), Nbu, epsX, Tol);
if (!extrema.IsDone())
return Standard_True;
const Standard_Integer Nbext = extrema.NbExt();
Standard_Real aDist2 = Dist * Dist;
for (Standard_Integer i=1; i<=Nbext; i++) {
if (extrema.SquareDistance(i) < aDist2) {
aDist2 = extrema.SquareDistance(i);
Paramproj = extrema.Point(i).Parameter();
}
}
Dist = sqrt (aDist2);
return Standard_True;
}
//=======================================================================
//function : Inters
//purpose : Intersection d'un segment avec une courbe
//=======================================================================
Standard_Boolean BRepBlend_BlendTool::Inters(const gp_Pnt2d& P1,
const gp_Pnt2d& P2,
const Handle(Adaptor3d_HSurface)&,
const Handle(Adaptor2d_HCurve2d)& C,
Standard_Real& Param,
Standard_Real& Dist)
{
const Standard_Real Tol = 1.e-8;
const gp_Vec2d v(P1,P2);
const Standard_Real mag = v.Magnitude();
if(mag < Tol) return Standard_False;
gp_Dir2d d(v);
Handle(Geom2d_Line) bid = new Geom2d_Line(P1,d);
Geom2dAdaptor_Curve seg(bid,-0.01*mag,1.01*mag);
Geom2dInt_GInter inter(seg,C->Curve2d(),Tol,Tol);
if (!inter.IsDone())
return Standard_False;
const Standard_Integer Nbint = inter.NbPoints();
if (Nbint == 0)
return Standard_False;
IntRes2d_IntersectionPoint ip = inter.Point(1);
Param = ip.ParamOnSecond();
Dist = P1.Distance(ip.Value());
return Standard_True;
}
Standard_Integer BRepBlend_BlendTool::NbSamplesV
(const Handle(Adaptor3d_HSurface)& S,
const Standard_Real u1,
const Standard_Real u2)
{
return 10;
}
Standard_Integer BRepBlend_BlendTool::NbSamplesU
(const Handle(Adaptor3d_HSurface)& S,
const Standard_Real u1,
const Standard_Real u2)
{
return 10;
}
void BRepBlend_BlendTool::Bounds(const Handle(Adaptor2d_HCurve2d)& A,
Standard_Real& Ufirst,
Standard_Real& Ulast)
{
Ufirst = BRepBlend_HCurve2dTool::FirstParameter(A);
Ulast = BRepBlend_HCurve2dTool::LastParameter(A);
}
|
