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
|
// File: BRepTest_ExtremaCommands.cxx
// Created: Fri Sep 8 11:37:46 1995
// Author: Modelistation
// <model@ecolox>
//
// modified by mps (juillet 96) : ajout de la commande distmini
#include <DBRep.hxx>
#include <BRepTest.hxx>
#include <gp_Pnt.hxx>
#include <TopoDS_Shape.hxx>
#include <TopoDS_Edge.hxx>
#include <BRepExtrema_Poly.hxx>
#include <BRepExtrema_DistShapeShape.hxx>
#include <BRepLib_MakeEdge.hxx>
#include <BRepLib_MakeEdge.hxx>
#include <BRepLib_MakeVertex.hxx>
#include <Draw.hxx>
#include <Draw_Interpretor.hxx>
//#ifdef WNT
#include <stdio.h>
//#endif
//=======================================================================
//function : distance
//purpose :
//=======================================================================
static Standard_Integer distance (Draw_Interpretor& di,
Standard_Integer n,
const char** a)
{
if (n < 3) return 1;
const char *name1 = a[1];
const char *name2 = a[2];
TopoDS_Shape S1 = DBRep::Get(name1);
TopoDS_Shape S2 = DBRep::Get(name2);
if (S1.IsNull() || S2.IsNull()) return 1;
gp_Pnt P1,P2;
Standard_Real D;
if (!BRepExtrema_Poly::Distance(S1,S2,P1,P2,D)) return 1;
//cout << " distance : " << D << endl;
di << " distance : " << D << "\n";
TopoDS_Edge E = (TopoDS_Edge) BRepLib_MakeEdge(P1,P2);
DBRep::Set("distance",E);
return 0;
}
static Standard_Integer distmini(Draw_Interpretor& di, Standard_Integer n, const char** a)
{
Standard_Integer i1;
// gp_Pnt P;
if (n != 4) return 1;
const char *ns1 = (a[2]), *ns2 = (a[3]), *ns0 = (a[1]);
TopoDS_Shape S1(DBRep::Get(ns1)), S2(DBRep::Get(ns2)) ;
BRepExtrema_DistShapeShape dst(S1 ,S2 );
if (dst.IsDone())
{
#ifdef DEB
//dst.Dump(cout);
Standard_SStream aSStream;
dst.Dump(aSStream);
di << aSStream;
#endif
char named[100];
sprintf(named, "%s%s" ,ns0,"_val");
char* tempd = named;
Draw::Set(tempd,dst.Value());
di << named << " ";
for (i1=1; i1<= dst.NbSolution(); i1++)
{
gp_Pnt P1,P2;
P1 = (dst.PointOnShape1(i1));
P2 = (dst.PointOnShape2(i1));
if (dst.Value()<=1.e-9)
{
TopoDS_Vertex V = (TopoDS_Vertex) BRepLib_MakeVertex(P1);
char namev[100];
if (i1==1)
sprintf(namev, "%s" ,ns0);
else sprintf(namev, "%s%d" ,ns0,i1);
char* tempv = namev;
DBRep::Set(tempv,V);
di << namev << " ";
}
else
{char name[100];
TopoDS_Edge E = (TopoDS_Edge) BRepLib_MakeEdge (P1, P2);
if (i1==1)
{sprintf(name,"%s",ns0);}
else {sprintf(name,"%s%d",ns0,i1);}
char* temp = name;
DBRep::Set(temp,E);
di << name << " " ;
}
}
}
else di << "probleme"<< "\n";
//else cout << "probleme"<< endl;
return 0;
}
//=======================================================================
//function : ExtremaCommands
//purpose :
//=======================================================================
void BRepTest::ExtremaCommands(Draw_Interpretor& theCommands)
{
static Standard_Boolean done = Standard_False;
if (done) return;
done = Standard_True;
const char* g = "TOPOLOGY Extrema commands";
theCommands.Add("dist","dist Shape1 Shape2" ,__FILE__,distance,g);
theCommands.Add("distmini",
"distmini name Shape1 Shape2",
__FILE__,
distmini,g);
}
|
