0031007: Coding - eliminate warnings issued while compiling with -pedantic flag
[occt.git] / src / SWDRAW / SWDRAW_ShapeFix.cxx
1 // Created on: 1999-03-09
2 // Created by: data exchange team
3 // Copyright (c) 1999-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
5 //
6 // This file is part of Open CASCADE Technology software library.
7 //
8 // This library is free software; you can redistribute it and/or modify it under
9 // the terms of the GNU Lesser General Public License version 2.1 as published
10 // by the Free Software Foundation, with special exception defined in the file
11 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
12 // distribution for complete text of the license and disclaimer of any warranty.
13 //
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
16
17
18 #include <BRep_Builder.hxx>
19 #include <BRep_Tool.hxx>
20 #include <BRepBuilderAPI.hxx>
21 #include <BRepTopAdaptor_FClass2d.hxx>
22 #include <DBRep.hxx>
23 #include <Draw.hxx>
24 #include <Draw_ProgressIndicator.hxx>
25 #include <gp_Pnt2d.hxx>
26 #include <gp_XYZ.hxx>
27 #include <Message_ListIteratorOfListOfMsg.hxx>
28 #include <Message_Msg.hxx>
29 #include <Precision.hxx>
30 #include <ShapeAnalysis_Edge.hxx>
31 #include <ShapeAnalysis_FreeBounds.hxx>
32 #include <ShapeAnalysis_Wire.hxx>
33 #include <ShapeAnalysis_WireOrder.hxx>
34 #include <ShapeAnalysis_WireVertex.hxx>
35 #include <ShapeBuild_ReShape.hxx>
36 #include <ShapeExtend_DataMapIteratorOfDataMapOfShapeListOfMsg.hxx>
37 #include <ShapeExtend_DataMapOfShapeListOfMsg.hxx>
38 #include <ShapeExtend_MsgRegistrator.hxx>
39 #include <ShapeExtend_WireData.hxx>
40 #include <ShapeFix.hxx>
41 #include <ShapeFix_Face.hxx>
42 #include <ShapeFix_FixSmallFace.hxx>
43 #include <ShapeFix_Shape.hxx>
44 #include <ShapeFix_ShapeTolerance.hxx>
45 #include <ShapeFix_Wire.hxx>
46 #include <ShapeFix_Wireframe.hxx>
47 #include <ShapeFix_WireVertex.hxx>
48 #include <SWDRAW.hxx>
49 #include <SWDRAW_ShapeFix.hxx>
50 #include <TCollection_AsciiString.hxx>
51 #include <TColStd_DataMapIteratorOfDataMapOfAsciiStringInteger.hxx>
52 #include <TColStd_DataMapOfAsciiStringInteger.hxx>
53 #include <TopAbs_State.hxx>
54 #include <TopExp.hxx>
55 #include <TopExp_Explorer.hxx>
56 #include <TopoDS.hxx>
57 #include <TopoDS_Edge.hxx>
58 #include <TopoDS_Face.hxx>
59 #include <TopoDS_Iterator.hxx>
60 #include <TopoDS_Shape.hxx>
61 #include <TopoDS_Wire.hxx>
62 #include <TopTools_DataMapOfShapeListOfShape.hxx>
63 #include <TopTools_HSequenceOfShape.hxx>
64 #include <TopTools_IndexedMapOfShape.hxx>
65 #include <TopTools_MapOfShape.hxx>
66
67 #ifdef AIX
68 #include <strings.h>
69 #endif
70 #include <stdio.h>
71 #include <ShapeExtend_WireData.hxx>
72
73 //=======================================================================
74 //function : edgesameparam
75 //purpose  : 
76 //=======================================================================
77
78 static Standard_Integer edgesameparam (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
79 {
80 //  const Standard_CString arg1 = argv[1];
81   const Standard_CString arg2 (argc > 2 ? argv[2] : NULL);
82 //        ****    Edge:SameParameter         ****
83   if (argc < 2) { di<<"shapename ,  option f to force, else only Edges not-SameParameter are computed\n"; return 1 /* Error */; }
84   TopoDS_Shape Shape = DBRep::Get(argv[1]);
85
86   if (!ShapeFix::SameParameter(Shape, (argc > 2 && arg2[0] == 'f') , BRepBuilderAPI::Precision()) )
87     di<<"Some edges were not processed\n";
88   di<<"\n";
89   return 0; // Done
90 }
91
92 //=======================================================================
93 //function : settolerance
94 //purpose  : 
95 //=======================================================================
96
97 static Standard_Integer settolerance (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
98 {
99   if (argc < 3) {
100     di<<    "myshape val   : forces tolerances at <val>\n"
101       <<      "myshape < max : sets maximum tolerance to <max>\n"
102         <<    "myshape > min : sets minimum tolerance to <min>\n"
103           <<  "myshape min max : bounds tolerances between <min> and <max>\n"
104             <<"myshape mode=v-e-f other args : idem but only on vertex-edge-face\n";
105     return (argc < 2 ? 0 : 1 /* Error */);
106   }
107   Standard_CString arg1 = argv[1];
108   Standard_CString arg2 = argv[2];
109   TopoDS_Shape Shape = DBRep::Get(arg1);
110   if (Shape.IsNull()) { di<<"Shape unknown : "<<arg1<<"\n"; return 1 /* Error */; }
111   char mod2 = arg2[0];
112   Standard_Integer premarg = 2;
113   TopAbs_ShapeEnum styp = TopAbs_SHAPE;
114   if (mod2 == 'v') { styp = TopAbs_VERTEX; premarg = 3; }
115   if (mod2 == 'e') { styp = TopAbs_EDGE;   premarg = 3; }
116   if (mod2 == 'w') { styp = TopAbs_WIRE;   premarg = 3; }
117   if (mod2 == 'f') { styp = TopAbs_FACE;   premarg = 3; }
118   if (mod2 == 'a') { styp = TopAbs_SHAPE;  premarg = 3; }
119
120   Standard_Real tmin,tmax;
121   mod2 = argv[premarg][0];
122   if (mod2 == '=') tmin = tmax = Draw::Atof (argv[argc-1]);
123   else if (mod2 == '<') { tmin = 0; tmax = Draw::Atof (argv[argc-1]); }
124   else if (mod2 == '>') { tmin = Draw::Atof (argv[argc-1]); tmax = 0; }
125   else { tmin = Draw::Atof (argv[premarg]);  tmax = Draw::Atof (argv[argc-1]); }
126
127   if (argc == premarg + 1 || tmin == tmax) di<<"Setting Tolerance to "<<tmin<<"\n";
128   else if (tmax < tmin) di<<"Minimum Tolerance to "<<tmin<<"\n";
129   else if (tmin <= 0)   di<<"Maximum Tolerance to "<<tmax<<"\n";
130   else                  di<<"Tolerance Limited between "<<tmin<<" and "<<tmax<<"\n";
131   ShapeFix_ShapeTolerance sat;
132   sat.LimitTolerance (Shape,tmin,tmax,styp);
133   return 0; // Done
134 }
135
136 //=======================================================================
137 //function : stwire
138 //purpose  : 
139 //=======================================================================
140
141 static Standard_Integer stwire (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
142 {
143   if (argc < 2) { // help
144     di<<"Donner nom shape depart + nom shape resultat + option\n";
145     di<<"Options de chargement :  x add connected  (sinon add simple)\n"
146       <<  "Options de traitement :  l fix little/BRepBuilderAPI\n"
147       <<"Options de sortie : aucune make simple\n"
148       <<"  m MakeAPI r avec reorder  v vertex\n"
149       <<"Autres (se cumulent) :  q quid(stats)\n";
150     return 0;
151   }
152   if (argc < 4) { di<<"stwire tout court pour help\n"; return 1 /* Error */; }
153   Standard_CString arg1 = argv[1];
154   Standard_CString arg2 = argv[2];
155
156 //  Options
157   Standard_Integer i;  int ox,ol,om,orint,oq,ov;  ox=ol=om=orint=oq=ov=0;
158   for (i = 3; i < argc; i ++) {
159     int valopt = 1;
160     char opt = argv[i][0];
161     if (opt == '+')   opt = argv[i][1];
162     if (opt == '-') { opt = argv[i][1]; valopt = 0; }
163     switch (opt) {
164     case 'l' : ol = valopt; break;
165     case 'm' : om = valopt; break;
166     case 'q' : oq = valopt; break;
167     case 'r' : orint = valopt; break;
168     case 'v' : ov = valopt; break;
169     case 'x' : ox = valopt; break;
170     default : break;
171     }
172   }
173   TopoDS_Shape Shape = DBRep::Get(arg1);
174   if (Shape.IsNull()) { di<<arg1<<" inconnu\n"; return 1 /* Error */; }
175
176 //  On y va
177   Handle(ShapeExtend_WireData) sbwd = new ShapeExtend_WireData;
178   Handle(ShapeAnalysis_Wire) saw = new ShapeAnalysis_Wire;
179   saw->Load ( sbwd );
180   
181   TopoDS_Shape awire;  // en principe un Wire
182   if (Shape.ShapeType() == TopAbs_WIRE) {
183     awire = Shape;
184   } else if (Shape.ShapeType() == TopAbs_FACE) {
185     saw->SetFace (TopoDS::Face(Shape));
186     TopExp_Explorer expw(Shape,TopAbs_WIRE);
187     if (expw.More()) awire = expw.Current();
188     saw->SetPrecision (BRepBuilderAPI::Precision());
189   }
190   if (awire.IsNull()) {
191     di<<"Neither FACE nor WIRE : "<<arg1<<"\n";
192     di<<"Considering as list of edges ...\n";
193     awire = Shape;
194 //    return 1 /* Error */;
195   }
196
197 // Chargement : normal ou par connected(oriented)
198 //  if (ox) {
199     for (TopExp_Explorer exp(Shape,TopAbs_EDGE); exp.More(); exp.Next()) {
200       TopoDS_Edge E = TopoDS::Edge (exp.Current());
201       Standard_Integer orient = saw->CheckShapeConnect (E);
202       di<<"Orientation : "<<orient<<" LowerDist : "<< saw->MinDistance3d() << "\n";
203       if (ox) sbwd->AddOriented (E,orient);
204       else sbwd->Add (E);
205     }
206 //  }
207 //  else sbwd->Init (awire);
208   
209   Handle(ShapeFix_Wire) sfw = new ShapeFix_Wire;
210   sfw->Init ( saw );
211
212 //   Traitement en cours
213   if (ol) {
214     Standard_Integer nb = sfw->NbEdges();
215     for (i = 1; i <= nb; i ++) {
216       Standard_Boolean stat = sfw->FixSmall (i, Standard_True, 0.0);
217       //std::cout<<"FixSmall for"<<i<<(stat ? " done" : " not done"); //:sw <<" StatusFix="<<STW.StatusFix()<<"\n";
218       di<<"FixSmall for"<<i;
219       if (!stat) {
220         di<<" not done";
221       } else {
222         di<<" done";
223       }
224       if ( sfw->LastFixStatus ( ShapeExtend_FAIL ) ) di << " (failed)";
225       di << "\n";
226     }
227   }
228
229 //  Traitement
230   if (orint) {  // reorder ?
231     ShapeAnalysis_WireOrder WO ( (Shape.ShapeType() != TopAbs_FACE), BRepBuilderAPI::Precision());
232
233     Standard_Integer stwo = saw->CheckOrder (WO);
234     Standard_Integer nb = WO.NbEdges();
235     di<<"Reorder status : "<<stwo<<"  NbEdges="<<nb<<"\n";
236     for (i = 1; i <= nb; i ++) {
237       Standard_Integer iord = WO.Ordered(i);
238       di<<"Edge n0 "<<i;
239       if ( sbwd->Edge(iord).Orientation() == TopAbs_REVERSED) di<<" REV";
240       else di<<" FWD";
241       di<<" ordered to "<<iord<<" Gap="<<WO.Gap(i)<<"\n";
242     }
243     di<<"Reorder not yet done\n";
244     sfw->FixReorder (WO);
245 //   Mais on n execute pas
246   }
247
248 //  Statistiques
249   if (oq) {
250     ShapeAnalysis_Edge sae;
251     Standard_Integer nb = sbwd->NbEdges();
252     di<<"NbEdges : "<<nb<<"\n";
253     for (i = 1; i <= nb; i ++) {
254       TopoDS_Edge E = sbwd->Edge(i);
255       di<<"Edge "<<i;
256       if (E.Orientation() == TopAbs_REVERSED) di<<" REV";
257       else di<<" FWD";
258       if (BRep_Tool::Degenerated(E)) di<<" DGNR";
259       if (sbwd->IsSeam(i))        di<<" SEAM_WIRE";
260       if (Shape.ShapeType() == TopAbs_FACE &&
261           sae.IsSeam(E,TopoDS::Face(Shape)))  di<<" SEAM_FACE";
262       if (Shape.ShapeType() == TopAbs_FACE ) {
263         if (sae.HasPCurve(E,TopoDS::Face(Shape)))     di<<" PCU";
264         else di<<" NO_PCU";
265       }
266       if (sae.HasCurve3d(E))    di<<" C3D";
267       else di<<" NO_C3D";
268       if (sae.IsClosed3d(E))  di<<" CLOSED";
269       di<<"\n";
270     }
271   }
272
273 //  Resultat
274   TopoDS_Wire result;
275   if (ov) {
276     ShapeAnalysis_WireVertex sawv;
277     sawv.Init ( sbwd, saw->Precision() );
278     sawv.Analyze();
279     Standard_Integer nb = sbwd->NbEdges();
280     di<<"Nb(End)Vertex : "<<nb<<"\n";
281     for (i = 1; i <= nb; i ++) {
282       gp_XYZ pos; Standard_Real upre,ufol;
283       Standard_Integer stat = sawv.Data (i,pos,upre,ufol);
284       di<<i<<" : ";
285       switch (stat) {
286       case 0 : di<<"Same Vertex\n"; break;
287       case 1 : di<<"Same Coords with recorded precisions (but not Same Vertex)\n"; break;
288       case 2 : di<<"Close (with preci="<< saw->Precision()<<")\n"; break;
289       case 3 : di<<"End of "<<i<<" OK, Start of "<<(i == nb ? 1 : i+1)<<" at U="<<ufol; break;
290       case 4 : di<<"End of "<<i<<" at U="<<upre<<", Start of "<<(i == nb ? 1 : i+1)<<" OK"; break;
291       case 5 : di<<"Intersection, End of "<<i<<" at U="<<upre<<", Start of "<<(i == nb ? 1 : i+1)<<" at U="<<ufol; break;
292         default : di<<"Disjoined\n";
293       }
294       if (stat >= 3 && stat <= 5) di<<"\n   - Position : "<<pos.X()<<"  "<<pos.Y()<<"  "<<pos.Z()<<"\n";
295     }
296     ShapeFix_WireVertex sfwv;
297     sfwv.Init ( sawv );
298     di<<"Nb Fixed Vertex : "<< sfwv.Fix() <<"\n";
299   }
300   /*
301   if (oc) {
302     if ( Shape.ShapeType() == TopAbs_FACE ) {
303       Correct_Wire CW ( TopoDS::Face(Shape) );
304       Standard_Integer i, nb = sbwd->NbEdges();
305       Standard_Integer num = 1;
306       for (i = 1; i <= nb; i ++)  CW.Add (sbwd->Edge(i));
307       CW.Perform ( saw->Precision() );
308       nb = CW.NbWires();
309       if (nb != 1) {
310         //  On prend celui qui a le plus d edges
311         Standard_Integer nbe, maxe = 0;
312         for (i = 1; i <= nb; i ++) {
313           TopoDS_Wire wir = CW.Wire(i);
314           nbe = 0;
315           for (TopoDS_Iterator ite(wir); ite.More(); ite.Next()) nbe ++;
316           if (nbe > maxe) { num = i; maxe = nbe; }
317         }
318         di<<"Correct_Wire produced "<<nb<<" Wires, taken n0 "<<num<<"\n";
319       }
320       result = CW.Wire (num);
321     }
322     else di << "Cannot apply Correct_Wire: face not defined\n";
323   }
324   */
325   else if (om) result = sbwd->WireAPIMake();
326   else         result = sbwd->Wire();
327   if (result.IsNull()) {
328     di<<"Pas de resultat, desole\n";
329     return 1;  // Fail
330   }
331   DBRep::Set (arg2,result);
332   return 0; // Done
333 }
334
335 //=======================================================================
336 //function : reface
337 //purpose  : 
338 //=======================================================================
339
340 static Standard_Integer reface (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
341 {
342   if (argc < 3) { di<<"Donner un nom de SHAPE (SHELL ou FACE) + un nom de RESULTAT\n"; return 1 /* Error */; }
343   Standard_CString arg1 = argv[1];
344   Standard_CString arg2 = argv[2];
345   TopoDS_Shape Shape = DBRep::Get(arg1);
346   if (Shape.IsNull()) { di<<"Shape unknown : "<<arg1<<"\n"; return 1 /* Error */; }
347
348   Standard_Boolean rebuild = Standard_False;
349
350   Handle(ShapeFix_Face) STF = new ShapeFix_Face;
351 //  Options ?
352   Standard_Integer i; // svv Jan11 2000 : porting on DEC
353   for (i = 3; i < argc; i ++) {
354     Standard_Boolean valopt = Standard_True;
355     char opt = argv[i][0];
356     if (opt == '+')   opt = argv[i][1];
357     if (opt == '-') { opt = argv[i][1]; valopt = Standard_False; }
358     //std::cout<<(valopt ? ".." : ".. NO");
359     if (!valopt) {
360       di<<".. NO";
361     } else {
362       di<<"..";
363     }
364
365     if (opt == 'R') { di<<" REBUILD-ANYWAY .."; rebuild = valopt; }
366     if (opt == 'd') { di<<" fix-dgnr .."; STF->FixWireTool()->FixDegeneratedMode() = valopt; }
367     if (opt == 'r') { di<<" fix-reorder-wire .."; STF->FixWireTool()->FixReorderMode()    = valopt; }
368     if (opt == 'k') {
369     }
370   }
371
372   TopoDS_Face  face;
373   ShapeBuild_ReShape resh;
374
375   Standard_Integer nbf = 0, nbfc = 0;
376   for (TopExp_Explorer EF (Shape,TopAbs_FACE); EF.More(); EF.Next()) {
377     TopoDS_Face F = TopoDS::Face (EF.Current());    face = F;
378     nbf ++;
379     Standard_Boolean newface = Standard_False;
380 //    on va voir si ShapeTool_Face trouve qqchose a redire
381 //:sw    ShapeTool_Wire STW;
382 //:sw    STW.SetFace (F);
383
384     STF->Init (F);  // qui fait tout
385     STF->Perform();
386     face = STF->Face();
387     newface = STF->Status(ShapeExtend_DONE) || rebuild;
388     if (newface) { nbfc ++; resh.Replace (F,face); }
389   }
390   if (nbfc > 0) {
391     di<<"Faces reprises par ShapeFix_Face : "<<nbfc<<"\n";
392     DBRep::Set (arg2,resh.Apply (Shape,TopAbs_FACE,2));
393     return 0; // Done
394   }
395   else di<<"ShapeFix_Face n a rien trouve a redire\n";
396   return 0;
397 }
398
399
400 //=======================================================================
401 //function : fixshape
402 //purpose  : 
403 //=======================================================================
404
405 static Standard_Integer fixshape (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
406 {
407   Handle(ShapeExtend_MsgRegistrator) msg = new ShapeExtend_MsgRegistrator;
408   Handle(ShapeFix_Shape) sfs = new ShapeFix_Shape;
409   sfs->SetMsgRegistrator ( msg );
410   
411   Standard_CString res = 0;
412   Standard_Integer par = 0, mess=0;
413   for ( Standard_Integer i=1; i < argc; i++ )
414   {
415     if (strlen(argv[i]) == 2 &&
416       (argv[i][0] == '-' || argv[i][0] == '+' || argv[i][0] == '*'))
417     {
418       Standard_Integer val = ( argv[i][0] == '-' ? 0 : argv[i][0] == '+' ? 1 : -1 );
419       switch ( argv[i][1] ) {
420       case 'l': sfs->FixWireTool()->FixLackingMode()          = val; break;
421       case 'o': sfs->FixFaceTool()->FixOrientationMode()      = val; break;
422       case 'h': sfs->FixWireTool()->FixShiftedMode()          = val; break;
423       case 'm': sfs->FixFaceTool()->FixMissingSeamMode()      = val; break;
424       case 'd': sfs->FixWireTool()->FixDegeneratedMode()      = val; break;
425       case 's': sfs->FixWireTool()->FixSmallMode()            = val; break;
426       case 'i': sfs->FixWireTool()->FixSelfIntersectionMode() = val; break;
427       case 'n': sfs->FixWireTool()->FixNotchedEdgesMode()     = val; break;
428       case '?': mess                                          = val; break;
429       }
430       continue;
431     }
432     else if (!strcmp(argv[i], "-maxtaila"))
433     {
434       if (++i >= argc)
435       {
436         break;
437       }
438
439       sfs->FixWireTool()->SetMaxTailAngle(Draw::Atof(argv[i]) * (M_PI / 180));
440     }
441     else if (!strcmp(argv[i], "-maxtailw"))
442     {
443       if (++i >= argc)
444       {
445         break;
446       }
447
448       sfs->FixWireTool()->SetMaxTailWidth(Draw::Atof(argv[i]));
449       sfs->FixWireTool()->FixTailMode() = 1;
450     }
451     else
452     {
453       switch ( par ) {
454       case 0: res = argv[i];                       break;
455       case 1: {
456         TopoDS_Shape initShape =  DBRep::Get(argv[i]);
457         if(initShape.IsNull()) continue;
458         sfs->Init ( initShape );
459       } break;
460       case 2: sfs->SetPrecision   (Draw::Atof(argv[i])); break;
461       case 3: sfs->SetMaxTolerance(Draw::Atof(argv[i])); break;
462       }
463     }
464     par++;
465   }
466
467   if ( par <2 ) {
468     di << "Use: " << argv[0] << " result shape [tolerance [max_tolerance]] [switches]\n"
469       "[-maxtaila <degrees>] [-maxtailw <width>]\n";
470     di << "Switches allow to tune parameters of ShapeFix\n"; 
471     di << "The following syntax is used: <symbol><parameter>\n"; 
472     di << "- symbol may be - to set parameter off, + to set on or * to set default\n"; 
473     di << "- parameters are identified by letters:\n"; 
474     di << "  l - FixLackingMode\n"; 
475     di << "  o - FixOrientationMode\n"; 
476     di << "  h - FixShiftedMode\n"; 
477     di << "  m - FixMissingSeamMode\n"; 
478     di << "  d - FixDegeneratedMode\n"; 
479     di << "  s - FixSmallMode\n"; 
480     di << "  i - FixSelfIntersectionMode\n"; 
481     di << "  n - FixNotchedEdgesMode\n"; 
482     di << "For enhanced message output, use switch '+?'\n"; 
483     return 1;
484   }
485
486   Handle(Draw_ProgressIndicator) aProgress = new Draw_ProgressIndicator (di, 1);
487   sfs->Perform (aProgress);
488   DBRep::Set (res,sfs->Shape());
489
490   if ( mess ) 
491   {
492     TColStd_DataMapOfAsciiStringInteger aMapOfNumberOfFixes;
493     Standard_SStream aSStream;
494     TopoDS_Compound aCompound;
495     BRep_Builder aBuilder;
496     aBuilder.MakeCompound (aCompound);
497     const ShapeExtend_DataMapOfShapeListOfMsg &map = msg->MapShape();
498     // Counting the number of each type of fixes. If the switch '*?' store all modified shapes in compound.
499     for ( ShapeExtend_DataMapIteratorOfDataMapOfShapeListOfMsg it ( map ); it.More(); it.Next() )
500     {
501       for ( Message_ListIteratorOfListOfMsg iter ( it.Value() ); iter.More(); iter.Next() )
502       {
503         if ( aMapOfNumberOfFixes.IsBound ( iter.Value().Value() ) )
504         {
505           aMapOfNumberOfFixes ( iter.Value().Value() )++;
506         }
507         else
508         {
509           aMapOfNumberOfFixes.Bind ( iter.Value().Value(), 1 );
510         }
511       }
512       if ( mess < 0 )
513       {
514         aBuilder.Add ( aCompound, it.Key() );
515       }
516     }
517     
518     aSStream << " Fix" << std::setw (58) << "Count\n";
519     aSStream << " ------------------------------------------------------------\n";
520     for ( TColStd_DataMapIteratorOfDataMapOfAsciiStringInteger anIter ( aMapOfNumberOfFixes ); anIter.More(); anIter.Next() )
521     {
522       aSStream << " " << anIter.Key() << std::setw ( 60 - anIter.Key().Length() ) << anIter.Value() << "\n";
523     }
524     aSStream << " ------------------------------------------------------------\n";
525     di << aSStream;
526
527     if ( mess < 0 )
528     {
529       char buff[256];
530       Sprintf ( buff, "%s_%s", res, "m" );
531       di << " Modified shapes saved in compound: " << buff;
532       DBRep::Set (buff, aCompound);
533     }
534   }
535   
536   return 0; // Done
537 }
538
539 //=======================================================================
540 //function : fixgaps
541 //purpose  : 
542 //=======================================================================
543
544 Standard_Integer fixgaps(Draw_Interpretor& di, Standard_Integer n, const char** a)
545 {
546   if (n < 3) return 1;
547   
548   TopoDS_Shape S = DBRep::Get(a[2]);
549   if (S.IsNull()) {
550     di << " Shape is null\n";
551     return 1;
552   }
553
554   Handle(ShapeFix_Wireframe) SFWF = new ShapeFix_Wireframe(S);
555   Standard_Real prec = ( n >3 ? Draw::Atof(a[3]) : 0. );
556   SFWF->SetPrecision(prec);
557   if ( SFWF->FixWireGaps() ) {
558     DBRep::Set(a[1],SFWF->Shape());
559     di<<" Wireframe gaps fixed on shape\n";
560   }
561
562   return 0;
563 }
564
565 //=======================================================================
566 //function : fixsmall
567 //purpose  : 
568 //=======================================================================
569
570 Standard_Integer fixsmall(Draw_Interpretor& di, Standard_Integer n, const char** a)
571 {
572   if (n < 3) return 1;
573   
574   TopoDS_Shape S = DBRep::Get(a[2]);
575   if (S.IsNull()) {
576     di << " Shape is null\n";
577     return 1;
578   }
579
580   Standard_Real prec = (n == 4)? Draw::Atof(a[3]) : 1.;
581   ShapeFix_Wireframe SFWF(S);
582   SFWF.SetPrecision(prec);
583   
584   if (SFWF.FixSmallEdges()) {
585     DBRep::Set(a[1],SFWF.Shape());
586     di<<" Small edges fixed on shape\n";
587   }
588
589   return 0;
590 }
591
592 //=======================================================================
593 //function : fixsmalledges
594 //purpose  : 
595 //=======================================================================
596
597 static Standard_Integer fixsmalledges(Draw_Interpretor& di, Standard_Integer n, const char** a)
598 {
599   if( n < 3) {
600     di<<"Invalid number of arguments\n";
601     return 1;
602   }
603   TopoDS_Shape Sh = DBRep::Get(a[2]);
604  
605   Standard_Integer k = 3;
606   Standard_Real tol = 100000;
607   Standard_Integer mode = 2;
608   Standard_Real tolang = M_PI/2;
609   if(n > k) 
610     tol = Draw::Atof(a[k++]);
611   
612   if(n > k) 
613     mode=  Draw::Atoi(a[k++]);
614    
615   if(n > k) 
616     tolang = Draw::Atof(a[k++]);
617      
618   Handle(ShapeFix_Wireframe) aSfwr = new ShapeFix_Wireframe();
619   Handle(ShapeBuild_ReShape) aReShape = new ShapeBuild_ReShape;
620   aSfwr->SetContext(aReShape);
621   aSfwr->Load(Sh);
622   aSfwr->SetPrecision(tol);
623   Standard_Boolean aModeDrop = Standard_True;
624   if(mode == 2) 
625     aModeDrop = Standard_False;
626    
627   TopTools_MapOfShape theSmallEdges, theMultyEdges;
628   TopTools_DataMapOfShapeListOfShape theEdgeToFaces,theFaceWithSmall;
629   aSfwr->CheckSmallEdges(theSmallEdges,theEdgeToFaces,theFaceWithSmall, theMultyEdges);
630   aSfwr->MergeSmallEdges (theSmallEdges,theEdgeToFaces,theFaceWithSmall, theMultyEdges, aModeDrop,tolang);
631    //aSfwr->FixSmallEdges(); 
632   TopoDS_Shape resShape = aSfwr->Shape();
633   DBRep::Set ( a[1], resShape );
634   return 0;
635 }
636
637 //=======================================================================
638 //function : fixsmallfaces
639 //purpose  : 
640 //=======================================================================
641
642 Standard_Integer fixsmallfaces (Draw_Interpretor& theDI, Standard_Integer theArgc, const char** theArgv)
643 {
644   if (theArgc < 3)
645   {
646     std::cerr << "Use: " << theArgv[0] << " result shape [tolerance]" << std::endl;
647     return 1;
648   }
649   
650   TopoDS_Shape aShape = DBRep::Get(theArgv[2]);
651   if (aShape.IsNull()) {
652     theDI << "Error: Shape " << theArgv[2] << " is null\n";
653     return 1;
654   }
655
656   Standard_Real aPrec = (theArgc < 4 ? 1. : Draw::Atof(theArgv[3]));
657
658   ShapeFix_FixSmallFace aFixSmallFaces;
659   aFixSmallFaces.Init (aShape);
660   aFixSmallFaces.SetPrecision(aPrec);
661   aFixSmallFaces.Perform();
662   if (! aFixSmallFaces.Shape().IsSame (aShape))
663   {
664     theDI << "Small faces are removed";
665     aShape = aFixSmallFaces.Shape();
666   }
667
668   DBRep::Set (theArgv[1], aShape);
669
670   return 0;
671 }
672
673 //=======================================================================
674 //function : checkoverlapedges
675 //purpose  : 
676 //=======================================================================
677
678 static Standard_Integer checkoverlapedges(Draw_Interpretor& di, Standard_Integer n, const char** a)
679 {
680    if( n < 3) {
681      di<<"Invalid number of arguments\n";
682      return 1;
683    }
684    TopoDS_Shape Sh1 = DBRep::Get(a[1]);
685    TopoDS_Shape Sh2 = DBRep::Get(a[2]);
686    if(Sh1.IsNull() || Sh2.IsNull()) {
687      di<<"Invalid arguments\n";
688      return 1;
689    }
690    TopoDS_Edge e1 = TopoDS::Edge(Sh1);
691    TopoDS_Edge e2 = TopoDS::Edge(Sh2);
692    if(e1.IsNull() || e2.IsNull()) {
693      di<<"Invalid type of arguments\n";
694      return 1;
695    }
696
697    if (BRep_Tool::Degenerated(e1))
698    {
699      di << a[1] << " is degenerated\n";
700      return 1;
701    }
702
703    if (BRep_Tool::Degenerated(e2))
704    {
705      di << a[2] << " is degenerated\n";
706      return 1;
707    }
708
709    Standard_Real aTol = Precision::Confusion();
710    Standard_Real aDistDomain = 0.0;
711    Standard_Integer k = 3;
712    if(k < n)
713      aTol = Draw::Atof(a[k++]);
714    if(k < n)
715      aDistDomain = Draw::Atof(a[k++]);
716    
717    ShapeAnalysis_Edge sae;
718    if(sae.CheckOverlapping(e1,e2,aTol,aDistDomain)) {
719      if(aDistDomain ==0.0)
720        di<<"Edges is overlaping comletly\n";
721      else {
722         di<<"Edges is overlaped\n";
723         di<<"with tolerance = "<<aTol<<"\n";
724         di<<"on segment length = "<<aDistDomain<<"\n";
725      }
726    }
727    else di<<"Edges is not overlaped\n";
728    return 0;
729 }
730
731 //=======================================================================
732 //function : checkfclass2d
733 //purpose  : 
734 //=======================================================================
735
736 static Standard_Integer checkfclass2d(Draw_Interpretor& di, Standard_Integer n, const char** a)
737 {
738   if( n < 4) {
739     di<<"Invalid number of arguments\n";
740     return 1;
741   }
742   TopoDS_Shape Sh1 = DBRep::Get(a[1]);
743   Standard_Real ucoord = Draw::Atof(a[2]);
744   Standard_Real vcoord = Draw::Atof(a[3]);
745   if(Sh1.IsNull() || Sh1.ShapeType()!= TopAbs_FACE) {
746     di<<"Invalid arguments\n";
747     return 1;
748   }
749   TopoDS_Face aFace = TopoDS::Face(Sh1);
750   gp_Pnt2d p2d(ucoord,vcoord);
751   BRepTopAdaptor_FClass2d f2d(aFace,Precision::Confusion());
752   TopAbs_State stat = f2d.Perform(p2d);
753   if(stat == TopAbs_OUT)
754     di<<"Point is OUT\n";
755   else if(stat == TopAbs_IN)
756     di<<"Point is IN\n";
757   else if(stat == TopAbs_ON)
758     di<<"Point is ON\n";
759   else
760     di<<"Point is UNKNOWN\n";
761   return 0;
762 }
763
764 static Standard_Integer connectedges(Draw_Interpretor& di, Standard_Integer n, const char** a)
765 {
766   if( n < 3) {
767     di<<"Invalid number of arguments. Should be : result shape [toler shared]\n";
768     return 1;
769   }
770   TopoDS_Shape aSh1 = DBRep::Get(a[2]);
771   if(aSh1.IsNull()) {
772     di<<"Shape is null\n";
773     return 1;
774   }
775   Standard_Real aTol = Precision::Confusion();
776   if( n > 3)
777     aTol = Draw::Atof(a[3]);
778   
779   Standard_Boolean shared = Standard_True;
780   if( n > 4)
781     shared = (Draw::Atoi(a[4]) == 1);
782   TopExp_Explorer aExpE(aSh1,TopAbs_EDGE);
783   Handle(TopTools_HSequenceOfShape) aSeqEdges = new TopTools_HSequenceOfShape;
784   Handle(TopTools_HSequenceOfShape) aSeqWires = new TopTools_HSequenceOfShape;
785   TopTools_IndexedMapOfShape aMapEdges;
786   for( ; aExpE.More(); aExpE.Next())
787   {
788     aSeqEdges->Append(aExpE.Current());
789     aMapEdges.Add(aExpE.Current());
790   }
791   
792   ShapeAnalysis_FreeBounds::ConnectEdgesToWires(aSeqEdges,aTol,shared,aSeqWires );
793   TopoDS_Compound aComp;
794   BRep_Builder aB;
795   aB.MakeCompound(aComp);
796   Standard_Integer i = 1;
797   for( ; i <= aSeqWires->Length() ; i++)
798   {
799     TopoDS_Shape aW = aSeqWires->Value(i);
800     di<<"Wire - "<<i<<" : \n";
801     
802     TopExp_Explorer aExp1(aW, TopAbs_EDGE);
803     for( ; aExp1.More(); aExp1.Next())
804     {
805       if(shared)
806       {
807         Standard_Integer ind = aMapEdges.FindIndex(aExp1.Current());
808         di<<ind<<" ";
809       }
810        else
811       {
812         TopoDS_Vertex aV1, aV2;
813         TopExp::Vertices(TopoDS::Edge(aExp1.Current()), aV1,aV2);
814         gp_Pnt aP = BRep_Tool::Pnt(aV1);
815         di<<aP.X()<<" "<<aP.Y()<<" "<<aP.Z()<<"\n";
816       }
817     }
818    
819     di<<"\n";
820     aB.Add( aComp,aSeqWires->Value(i));
821   }
822   DBRep::Set(a[1],aComp);
823   return 0;
824   
825 }
826
827 //=======================================================================
828 //function : InitCommands
829 //purpose  : 
830 //=======================================================================
831
832  void SWDRAW_ShapeFix::InitCommands(Draw_Interpretor& theCommands) 
833 {
834   static Standard_Integer initactor = 0;
835   if (initactor)
836   {
837     return;
838   }
839   initactor = 1;
840   
841   Standard_CString g = SWDRAW::GroupName();
842   
843   theCommands.Add ("edgesameparam","nom shape draw ou * [+ option force]",
844                    __FILE__,edgesameparam,g);
845   theCommands.Add ("settolerance","shape [mode=v-e-f-a] val(fix value) or tolmin tolmax",
846                    __FILE__,settolerance,g);
847   theCommands.Add ("stwire","stwire tout court pour help complet",
848                    __FILE__,stwire,g);
849   theCommands.Add ("reface","shape result : controle sens wire",
850                    __FILE__,reface,g);
851   theCommands.Add ("fixshape",
852 "res shape [preci [maxpreci]] [{switches}]\n"
853 "  [-maxtaila <degrees>] [-maxtailw <width>]",
854                    __FILE__,fixshape,g);
855 //  theCommands.Add ("testfill","result edge1 edge2",
856 //                 __FILE__,XSHAPE_testfill,g);
857   theCommands.Add ("fixwgaps","result shape [toler=0]",
858                    __FILE__,fixgaps,g);
859   theCommands.Add ("fixsmall","result shape [toler=1.]",
860                    __FILE__,fixsmall,g);
861   theCommands.Add ("fixsmalledges","result shape [toler mode amxangle]",
862                    __FILE__,fixsmalledges,g);
863   theCommands.Add ("fixsmallfaces","result shape [toler=1.]",
864                    __FILE__,fixsmallfaces,g);
865   theCommands.Add ("checkoverlapedges","edge1 edge2 [toler domaindist]",
866                    __FILE__,checkoverlapedges,g);
867   theCommands.Add ("checkfclass2d","face ucoord vcoord",
868                    __FILE__,checkfclass2d,g);
869   theCommands.Add ("connectedges","res shape [toler shared]",
870                    __FILE__,connectedges,g);
871   
872 }
873