projects / occt.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
0022898: IGES import fails in german environment
[occt.git] / src / BRepTest / BRepTest_OtherCommands.cxx
CommitLineData
b311480e 1// Created on: 1995-04-13
2// Created by: Jean Yves LEBEY
3// Copyright (c) 1995-1999 Matra Datavision
4// Copyright (c) 1999-2012 OPEN CASCADE SAS
5//
6// The content of this file is subject to the Open CASCADE Technology Public
7// License Version 6.5 (the "License"). You may not use the content of this file
8// except in compliance with the License. Please obtain a copy of the License
9// at http://www.opencascade.org and read it completely before using this file.
10//
11// The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
12// main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
13//
14// The Original Code and all software distributed under the License is
15// distributed on an "AS IS" basis, without warranty of any kind, and the
16// Initial Developer hereby disclaims all such warranties, including without
17// limitation, any warranties of merchantability, fitness for a particular
18// purpose or non-infringement. Please see the License for the specific terms
19// and conditions governing the rights and limitations under the License.
20
7fd59977 21
22#include <BRepTest.hxx>
23
91322f44 24#include <Draw.hxx>
7fd59977 25#include <Draw_Interpretor.hxx>
26#include <DBRep.hxx>
a15d6ace 27#include <DrawTrSurf.hxx>
28
29#include <string.h>
30#include <stdio.h>
31
32#include <Standard_ErrorHandler.hxx>
33#include <Precision.hxx>
34#include <TCollection_AsciiString.hxx>
35#include <gp_Lin.hxx>
36#include <gp_Pnt.hxx>
37#include <gp_Dir.hxx>
38#include <gp_Ax1.hxx>
39
40#include <ElCLib.hxx>
41
42#include <TColgp_SequenceOfPnt.hxx>
43
44#include <GeomAbs_JoinType.hxx>
45#include <Geom_Line.hxx>
46
47#include <IntCurvesFace_Intersector.hxx>
48
49#include <TopAbs.hxx>
50#include <TopAbs_Orientation.hxx>
51
52#include <TopoDS.hxx>
7fd59977 53#include <TopoDS_Shape.hxx>
54#include <TopoDS_Compound.hxx>
55#include <TopoDS_CompSolid.hxx>
56#include <TopoDS_Solid.hxx>
57#include <TopoDS_Shell.hxx>
58#include <TopoDS_Face.hxx>
59#include <TopoDS_Wire.hxx>
60#include <TopoDS_Edge.hxx>
61#include <TopoDS_Vertex.hxx>
a15d6ace 62#include <TopoDS_Iterator.hxx>
7fd59977 63
a15d6ace 64#include <BRep_Builder.hxx>
65#include <BRep_Tool.hxx>
7fd59977 66
7fd59977 67#include <TopExp_Explorer.hxx>
a15d6ace 68
69#include <TopTools_ListOfShape.hxx>
70#include <TopTools_ListIteratorOfListOfShape.hxx>
7fd59977 71#include <TopTools_MapOfShape.hxx>
72
a15d6ace 73#include <LocOpe_CSIntersector.hxx>
74#include <LocOpe_SequenceOfLin.hxx>
75#include <LocOpe_PntFace.hxx>
76#include <BRepFeat_MakeDPrism.hxx>
7fd59977 77
a15d6ace 78#include <BRepTools.hxx>
79#include <BRepIntCurveSurface_Inter.hxx>
80#include <BRepOffset.hxx>
81#include <BRepOffset_MakeOffset.hxx>
82#include <BRepClass3d_SolidClassifier.hxx>
83
84static
85 void SampleEdges (const TopoDS_Shape& theShape,
86 TColgp_SequenceOfPnt& theSeq);
87static
88 TopoDS_Face NextFaceForPrism (const TopoDS_Shape& shape,
89 const TopoDS_Shape& basis,
90 const gp_Ax1& ax1);
91static
92 void PrintState (Draw_Interpretor& aDI,
93 const TopAbs_State& aState);
94//
95static Standard_Integer emptyshape(Draw_Interpretor&, Standard_Integer, const char** );
96static Standard_Integer subshape (Draw_Interpretor&, Standard_Integer, const char** );
97static Standard_Integer brepintcs (Draw_Interpretor&, Standard_Integer, const char** );
98static Standard_Integer MakeBoss (Draw_Interpretor&, Standard_Integer, const char** );
99static Standard_Integer MakeShell (Draw_Interpretor&, Standard_Integer, const char** );
100static Standard_Integer xbounds (Draw_Interpretor&, Standard_Integer, const char** );
101static Standard_Integer xclassify (Draw_Interpretor&, Standard_Integer, const char** );
7fd59977 102
103//=======================================================================
a15d6ace 104//function : OtherCommands
105//purpose :
7fd59977 106//=======================================================================
a15d6ace 107void BRepTest::OtherCommands(Draw_Interpretor& theCommands)
108{
109 static Standard_Boolean done = Standard_False;
110 if (done) return;
111 done = Standard_True;
112
113 const char* g = "TOPOLOGY other commands";
114
115 theCommands.Add("shape",
116 "shape name V/E/W/F/Sh/So/CS/C; make a empty shape",__FILE__,emptyshape,g);
7fd59977 117
a15d6ace 118 theCommands.Add("subshape",
119 "subshape name V/E/W/F/Sh/So/CS/C index; get subsshape <index> of given type"
120 ,__FILE__,subshape,g);
121
122 theCommands.Add("BRepIntCS",
123 "Calcul d'intersection entre face et curve : BRepIntCS curve shape"
124 ,__FILE__,brepintcs,g);
125
126 theCommands.Add("makeboss", "create a boss on the shape myS", __FILE__, MakeBoss, g);
127 theCommands.Add("mksh", "create a shell on Shape", __FILE__, MakeShell, g);
128 theCommands.Add("xbounds", "xbounds face", __FILE__, xbounds, g);
129 theCommands.Add("xclassify", "use xclassify Solid [Tolerance=1.e-7]", __FILE__, xclassify, g);
130
131
132}
133//=======================================================================
134//function : emptyshape
135//purpose : shape : shape name V/E/W/F/SH/SO/CS/C
136//=======================================================================
137Standard_Integer emptyshape(Draw_Interpretor& , Standard_Integer n, const char** a)
7fd59977 138{
139 if (n <= 1) return 1;
140
141 BRep_Builder B;
142 TopoDS_Shape S;
143
144 if (n == 3) {
145 TCollection_AsciiString as(a[2]); as.LowerCase();
146 const char* a2 = as.ToCString();
147
148 if ( ! strcmp(a2,"c") ) {
149 TopoDS_Compound SS; B.MakeCompound(SS); S = SS;
150 }
151 else if ( ! strcmp(a2,"cs") ) {
152 TopoDS_CompSolid SS; B.MakeCompSolid(SS); S = SS;
153 }
154 else if ( ! strcmp(a2,"so") ) {
155 TopoDS_Solid SS; B.MakeSolid(SS); S = SS;
156 }
157 else if ( ! strcmp(a2,"sh") ) {
158 TopoDS_Shell SS; B.MakeShell(SS); S = SS;
159 }
160 else if ( ! strcmp(a2,"f") ) {
161 TopoDS_Face SS; B.MakeFace(SS); S = SS;
162 }
163 else if ( ! strcmp(a2,"w") ) {
164 TopoDS_Wire SS; B.MakeWire(SS); S = SS;
165 }
166 else if ( ! strcmp(a2,"e") ) {
167 TopoDS_Edge SS; B.MakeEdge(SS); S = SS;
168 }
169 else if ( ! strcmp(a2,"v") ) {
170 TopoDS_Vertex SS; B.MakeVertex(SS); S = SS;
171 }
172 else {
173 return 1;
174 }
175 }
176
177 const char *shapename = a[1];
178 DBRep::Set(shapename,S);
179 return 0;
180}
181
182//=======================================================================
a15d6ace 183//function : subshape
184//purpose :
7fd59977 185//=======================================================================
a15d6ace 186Standard_Integer subshape(Draw_Interpretor& di, Standard_Integer n, const char** a)
7fd59977 187{
188 if (n <= 2) return 1;
189
190
191 TopoDS_Shape S = DBRep::Get(a[1]);
192 if (S.IsNull()) return 0;
193 char newname[1024];
194 strcpy(newname,a[1]);
195 char* p = newname;
196 while (*p != '\0') p++;
197 *p = '_';
198 p++;
199 Standard_Integer i = 0;
200 if (n == 3) {
91322f44 201 Standard_Integer isub = Draw::Atoi(a[2]);
7fd59977 202 TopoDS_Iterator itr(S);
203 while (itr.More()) {
204 i++;
205 if ( i == isub ) {
91322f44 206 Sprintf(p,"%d",i);
7fd59977 207 DBRep::Set(newname,itr.Value());
208 di.AppendElement(newname);
209 break;
210 }
211 itr.Next();
212 }
213 }
214 else {
215 // explode a type
216 TopAbs_ShapeEnum typ;
217 switch (a[2][0]) {
218
219 case 'C' :
220 case 'c' :
221 typ = TopAbs_COMPSOLID;
222 break;
223
224 case 'S' :
225 case 's' :
226 if ((a[2][1] == 'O')||(a[2][1] == 'o'))
227 typ = TopAbs_SOLID;
228 else if ((a[2][1] == 'H')||(a[2][1] == 'h'))
229 typ = TopAbs_SHELL;
230 else
231 return 1;
232 break;
233
234 case 'F' :
235 case 'f' :
236 typ = TopAbs_FACE;
237 break;
238
239 case 'W' :
240 case 'w' :
241 typ = TopAbs_WIRE;
242 break;
243
244 case 'E' :
245 case 'e' :
246 typ = TopAbs_EDGE;
247 break;
248
249 case 'V' :
250 case 'v' :
251 typ = TopAbs_VERTEX;
252 break;
253
254 default :
255 return 1;
256 }
257
91322f44 258 Standard_Integer isub = Draw::Atoi(a[3]);
7fd59977 259 TopTools_MapOfShape M;
260 M.Add(S);
261 TopExp_Explorer ex(S,typ);
262 while (ex.More()) {
263 if (M.Add(ex.Current())) {
264 i++;
265 if ( i == isub ) {
91322f44 266 Sprintf(p,"%d",i);
7fd59977 267 DBRep::Set(newname,ex.Current());
268 di.AppendElement(newname);
269 break;
270 }
271 }
272 ex.Next();
273 }
274 }
275 return 0;
276}
7fd59977 277//=======================================================================
a15d6ace 278//function : brepintcs
279//purpose :
7fd59977 280//=======================================================================
a15d6ace 281Standard_Integer brepintcs(Draw_Interpretor& , Standard_Integer n, const char** a)
7fd59977 282{
283 if (n <= 2) return 1;
284 TopoDS_Shape S = DBRep::Get(a[n-1]);
285 if (S.IsNull()) return 3;
286
287 static BRepIntCurveSurface_Inter theAlg;
288 static double tol=1e-6;
289 static int nbpi=0;
290 static gp_Pnt curp;
291
292 if (n==3) {
293 Handle(Geom_Curve) C= DrawTrSurf::GetCurve(a[1]);
294 if (C.IsNull()) return 2;
295 GeomAdaptor_Curve acur(C);
296 theAlg.Init(S, acur, tol);
297 for (; theAlg.More(); theAlg.Next()) {
298 curp=theAlg.Pnt();
299 nbpi++;
300 char name[64];
301 char* temp = name; // pour portage WNT
91322f44 302 Sprintf(temp, "%s_%d", "brics", nbpi);
7fd59977 303 DrawTrSurf::Set(temp, curp);
304 }
305 }
306 else {
307 Handle(Geom_Line) hl;
308 gp_Lin thel;
309 for (Standard_Integer il = 1; il<n ; il++) {
310 hl= Handle(Geom_Line)::DownCast(DrawTrSurf::GetCurve(a[il]));
311 if (!hl.IsNull()) {
312 thel=hl->Lin();
313 theAlg.Init(S, thel, tol);
314 for (; theAlg.More(); theAlg.Next()) {
315 curp=theAlg.Pnt();
316 nbpi++;
317 char name[64];
318 char* temp = name; // pour portage WNT
91322f44 319 Sprintf(temp, "%s_%d", "brics", nbpi);
7fd59977 320 DrawTrSurf::Set(temp, curp);
321 }
322 }
323 }
324 }
325 //POP pour NT
326 return 0;
327}
7fd59977 328//=======================================================================
a15d6ace 329//function : MakeBoss
330//purpose :
7fd59977 331//=======================================================================
a15d6ace 332Standard_Integer MakeBoss(Draw_Interpretor& , Standard_Integer , const char** a)
7fd59977 333{
7fd59977 334 TopoDS_Shape myS = DBRep::Get( a[2] );
335
336 TopoDS_Shape myBasis = DBRep::Get( a[3] ) ;
337
338 Standard_Real ang = -0.05235987901687622;
339
340
341 TopoDS_Face basis = TopoDS::Face(myBasis);
342
343 BRepFeat_MakeDPrism DPRISM(myS, basis, basis, ang, 1, Standard_True);
344
345 TopoDS_Shape myFaceOnShape;
346 gp_Pnt Pnt(0.0, 0.0, 50.0);
347 gp_Dir Dir(-0.0, -0.0, -1.0);
348 gp_Ax1 ax(Pnt, Dir);
349
350 myFaceOnShape = NextFaceForPrism(myS, myBasis, ax);
351
352 DPRISM.Perform (myFaceOnShape);
353 DPRISM.Build();
354
355 if( DPRISM.IsDone() ) DBRep::Set( a[1], DPRISM.Shape() );
356
357 return 0;
358}
7fd59977 359//=======================================================================
a15d6ace 360//function : MakeShell
361//purpose :
7fd59977 362//=======================================================================
a15d6ace 363Standard_Integer MakeShell(Draw_Interpretor& , Standard_Integer , const char** a)
7fd59977 364{
365
366 TopoDS_Shape aShape = DBRep::Get( a[1] );
367 TopTools_ListOfShape Lst;
368 TopExp_Explorer Exp(aShape, TopAbs_FACE);
369 TopoDS_Shape InputShape(DBRep::Get( a[2] ));
370 TopoDS_Face F = TopoDS::Face(InputShape);
371// TopoDS_Face F = TopoDS::Face(DBRep::Get( a[2] ));
372
91322f44 373 Standard_Real Off = -Draw::Atof( a[3] );
7fd59977 374
375 BRepOffset_MakeOffset Offset;
376
377 Offset.Initialize( aShape, Off, 1.0e-3, BRepOffset_Skin,
378 Standard_True , Standard_False , GeomAbs_Arc );
379 Offset.AddFace( F );
380 Offset.MakeThickSolid();
381
382 if( Offset.IsDone() ) {
383 // SaveShape::Save(Offset.Shape(), "ss");
384 DBRep::Set( a[1], Offset.Shape() );
385 }
386 return 0;
387}
bd05fabf
S		 388//=======================================================================
389//function : xbounds
390//purpose :
391//=======================================================================
392Standard_Integer xbounds(Draw_Interpretor& di, Standard_Integer n, const char** a)
393{
394 if (n<2) {
395 di << "Usage : " << a[0] << " face" << "\n";
396 return 0;
397 }
398 //
399
400 Standard_Real aUMin, aUMax, aVMin, aVMax;
401 TopoDS_Shape aS;
402 TopoDS_Face aF;
403 //
404 aS=DBRep::Get(a[1]);
405 if (aS.IsNull()) {
406 di << " null shapes is not allowed here\n";
407 return 0;
408 }
409 if (aS.ShapeType()!=TopAbs_FACE) {
410 di << " shape" << a[1] <<" must be a face\n";
411 return 0;
412 }
413 //
414 aF=*((TopoDS_Face*)&aS);
415 //
416 BRepTools::UVBounds(aF, aUMin, aUMax, aVMin, aVMax);
417 //
418 TCollection_AsciiString aStr;
419 TCollection_AsciiString sUMin(aUMin);
420 TCollection_AsciiString sUMax(aUMax);
421 TCollection_AsciiString sVMin(aVMin);
422 TCollection_AsciiString sVMax(aVMax);
423 //
424 aStr=aStr+sUMin + "\n";
425 aStr=aStr+sUMax + "\n";
426 aStr=aStr+sVMin + "\n";
427 aStr=aStr+sVMax + "\n";
428 di <<aStr.ToCString();
429 //
430 return 0;
431}
7fd59977 432//=======================================================================
a15d6ace 433//function : xclassify
7fd59977 434//purpose :
435//=======================================================================
a15d6ace 436Standard_Integer xclassify (Draw_Interpretor& aDI, Standard_Integer n, const char** a)
7fd59977 437{
a15d6ace 438 if (n < 2) {
439 aDI<<" use xclassify Solid [Tolerance=1.e-7]\n";
440 return 1;
441 }
442
443 TopoDS_Shape aS = DBRep::Get(a[1]);
444 if (aS.IsNull()) {
445 aDI<<" Null Shape is not allowed here\n";
446 return 0;
447 }
448
449 if (aS.ShapeType()!=TopAbs_SOLID) {
450 aDI<< " Shape type must be SOLID\n";
451 return 0;
452 }
453 //
454 Standard_Real aTol=1.e-7;
455 TopAbs_State aState = TopAbs_UNKNOWN;
456 //
457 aTol=1.e-7;
458 if (n==3) {
91322f44 459 aTol=Draw::Atof(a[2]);
a15d6ace 460 }
461 //
462 BRepClass3d_SolidClassifier aSC(aS);
463 aSC.PerformInfinitePoint(aTol);
464
465 aState = aSC.State();
466 PrintState(aDI, aState);
467 //
468 return 0;
469}
470//=======================================================================
471//function : PrintState
472//purpose :
473//=======================================================================
474void PrintState (Draw_Interpretor& aDI,
475 const TopAbs_State& aState)
476{
477 aDI<<"state is: ";
478 switch (aState) {
479 case TopAbs_IN:
480 aDI<<"IN\n";
481 break;
482 case TopAbs_OUT:
483 aDI<<"OUT\n";
484 break;
485 case TopAbs_ON:
486 aDI<<"ON\n";
487 break;
488 case TopAbs_UNKNOWN:
489 default:
490 aDI<<"UNKNOWN\n";
491 break;
492 }
7fd59977 493}
7fd59977 494//=======================================================================
495//function : NextFaceForPrism
496//purpose : Search a face from <shape> which intersects with a line of
497// direction <ax1> and location a point of <basis>.
498//=======================================================================
a15d6ace 499TopoDS_Face NextFaceForPrism (const TopoDS_Shape& shape,
500 const TopoDS_Shape& basis,
501 const gp_Ax1& ax1)
7fd59977 502{
503 TopoDS_Face nextFace;
504
505 TColgp_SequenceOfPnt seqPnts;
506 SampleEdges(basis, seqPnts);
507
508 for (Standard_Integer i=1; i<=seqPnts.Length(); i++) {
509 const gp_Pnt& pt = seqPnts(i);
510 // find a axis through a face
511 gp_Dir dir = ax1.Direction();
512 gp_Ax1 ax1b(pt, dir);
513
514 LocOpe_CSIntersector ASI(shape);
515 LocOpe_SequenceOfLin slin;
516 slin.Append(ax1b);
517 ASI.Perform(slin);
518
519 if (ASI.IsDone()) {
520 Standard_Integer no=1, IndFrom, IndTo;
521 TopAbs_Orientation theOr;
522 Standard_Real min = 1.e-04, Tol = -Precision::Confusion();
523 if (ASI.LocalizeAfter (no, min, Tol, theOr, IndFrom, IndTo)) {
524 nextFace = ASI.Point(no, IndFrom).Face();
525 break;
526 }
527 }
528 }
529
530 return nextFace;
531}
532
7fd59977 533
534//=======================================================================
535//function : SampleEdges
536//purpose : Sampling of <theShape>.
537//design : Collect the vertices and points on the edges
538//=======================================================================
a15d6ace 539void SampleEdges (const TopoDS_Shape& theShape, TColgp_SequenceOfPnt& theSeq)
7fd59977 540{
541
542 theSeq.Clear();
543
544
545 TopTools_MapOfShape theMap;
546 TopExp_Explorer exp;
547
548 // Adds all vertices/pnt
549 for (exp.Init(theShape,TopAbs_VERTEX); exp.More(); exp.Next()) {
550 if (theMap.Add(exp.Current())) {
551 theSeq.Append (BRep_Tool::Pnt(TopoDS::Vertex(exp.Current())));
552 }
553 }
554
555 // Computes points on edge, but does not take the extremities into account
556 Standard_Integer NECHANT = 5;
557 Handle(Geom_Curve) C;
558 Standard_Real f,l,prm;
559 for (exp.Init (theShape,TopAbs_EDGE); exp.More(); exp.Next()) {
560 const TopoDS_Edge& edg = TopoDS::Edge(exp.Current());
561 if (theMap.Add(edg)) {
562 if (!BRep_Tool::Degenerated(edg)) {
563 C = BRep_Tool::Curve(edg,f,l);
564 for (Standard_Integer i=1; i < NECHANT; i++) {
565 prm = ((NECHANT-i)*f+i*l)/NECHANT;
566 theSeq.Append (C->Value(prm));
567 }
568 }
569 }
570 }
571}
572
Open CASCADE Technology repository