src/BRepTest/BRepTest_CheckCommands.cxx

   1 // Created on: 1996-02-23
   2 // Created by: Jacques GOUSSARD
   3 // Copyright (c) 1996-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 #include <BRepTest.hxx>
  18 #include <BRepAlgo.hxx>
  19 #include <BRepCheck_Analyzer.hxx>
  20 #include <BRepCheck_Result.hxx>
  21 #include <BRepCheck_ListIteratorOfListOfStatus.hxx>
  22 #include <TopoDS_Iterator.hxx>
  23 #include <TopExp_Explorer.hxx>
  24 #include <TopTools_DataMapOfShapeListOfShape.hxx>
  25 #include <TopTools_DataMapIteratorOfDataMapOfShapeListOfShape.hxx>
  26 #include <TopTools_DataMapIteratorOfDataMapOfIntegerListOfShape.hxx>
  27 #include <TopTools_ListOfShape.hxx>
  28 #include <TopTools_ListIteratorOfListOfShape.hxx>
  29 #include <DBRep_DrawableShape.hxx>
  30 #include <Draw_SequenceOfDrawable3D.hxx>
  31 #include <BRepCheck.hxx>
  32 #include <BRepCheck_Edge.hxx>
  33 #include <Draw_Interpretor.hxx>
  34 #include <Draw_Appli.hxx>
  35 #include <Draw.hxx>
  36 #include <DBRep.hxx>
  37 #include <BRepTest.hxx>
  38 #include <GeometryTest.hxx>
  39 #include <Precision.hxx>
  40 #include <LocalAnalysis.hxx>
  41 #include <LocalAnalysis_SurfaceContinuity.hxx>
  42 #include <Geom_SphericalSurface.hxx>
  43 #include <Geom_Surface.hxx>
  44 #include <Geom_Curve.hxx>
  45 #include <Geom2d_TrimmedCurve.hxx>
  46 #include <Geom2d_Curve.hxx>
  47 #include <DrawTrSurf.hxx>
  48 #include <GeomAbs_Shape.hxx>
  49 #include <TCollection_AsciiString.hxx>
  50 #include <TopoDS.hxx>
  51 #include <TopExp.hxx>
  52 #include <TopTools_IndexedDataMapOfShapeListOfShape.hxx>
  53 #include <TopTools_MapOfShape.hxx>
  54 #include <TopTools_MapIteratorOfMapOfShape.hxx>
  55 #include <TopoDS_Shape.hxx>
  56 #include <TopoDS_Edge.hxx>
  57 #include <TopoDS_Face.hxx>
  58 #include <BRep_Tool.hxx>
  59
  60 #include <TopOpeBRepTool_PurgeInternalEdges.hxx>
  61 //#include <TopOpeBRepTool_FuseEdges.hxx>
  62 #include <BRepLib.hxx>
  63 #include <BRepLib_FuseEdges.hxx>
  64
  65 #include <TopTools_HSequenceOfShape.hxx>
  66 #include <BRep_Builder.hxx>
  67 #include <TopoDS_Compound.hxx>
  68 #include <TColStd_HArray1OfInteger.hxx>
  69
  70 #include <Standard_ErrorHandler.hxx>
  71 #include <Standard_Failure.hxx>
  72
  73 #include <stdio.h>
  74
  75 //Number of BRepCheck_Statuses in BRepCheck_Status.hxx file
  76 //(BRepCheck_NoError is not considered, i.e. general status
  77 //is smaller by one specified in file)
  78 static const Standard_Integer NumberOfStatus = 36;
  79
  80 static char* checkfaultyname = NULL;
  81 Standard_EXPORT void BRepTest_CheckCommands_SetFaultyName(const char* name)
  82 {
  83   if (checkfaultyname != NULL) {
  84     free(checkfaultyname);
  85     checkfaultyname = NULL;
  86   }
  87   if (name == NULL) {
  88     checkfaultyname = (char*)malloc(strlen("faulty_")+1);
  89     strcpy(checkfaultyname,"faulty_");
  90   }
  91   else {
  92     checkfaultyname = (char*)malloc(strlen(name)+1);
  93     strcpy(checkfaultyname,name);
  94   }
  95 }
  96
  97
  98 static TopTools_DataMapOfShapeListOfShape theMap;
  99 static Standard_Integer nbfaulty = 0;
 100 static Draw_SequenceOfDrawable3D lfaulty;
 101
 102 Standard_IMPORT Standard_Integer BRepCheck_Trace(const Standard_Integer phase);
 103
 104 //=======================================================================
 105 //function : FindNamed
 106 //=======================================================================
 107 static Standard_Boolean FindNamed(const TopoDS_Shape& S,
 108                                   char*& Name)
 109 {
 110   for (Standard_Integer i = 1 ;i <= lfaulty.Length(); i++) {
 111     Handle(DBRep_DrawableShape) DS =
 112       Handle(DBRep_DrawableShape)::DownCast(lfaulty(i));
 113     if (DS->Shape().IsSame(S)) {
 114       Name = (char*)DS->Name();
 115       return Standard_True;
 116     }
 117   }
 118   return Standard_False;
 119 }
 120
 121
 122 //=======================================================================
 123 //function : Contains
 124 //=======================================================================
 125 static Standard_Boolean Contains(const TopTools_ListOfShape& L,
 126                                  const TopoDS_Shape& S)
 127 {
 128   TopTools_ListIteratorOfListOfShape it;
 129   for (it.Initialize(L); it.More(); it.Next()) {
 130     if (it.Value().IsSame(S)) {
 131       return Standard_True;
 132     }
 133   }
 134   return Standard_False;
 135 }
 136
 137
 138
 139 //=======================================================================
 140 //function : PrintSub
 141 //=======================================================================
 142 static void PrintSub(Standard_OStream& OS,
 143                      const BRepCheck_Analyzer& Ana,
 144                      const TopoDS_Shape& S,
 145                      const TopAbs_ShapeEnum Subtype)
 146
 147 {
 148   char* Name;
 149   BRepCheck_ListIteratorOfListOfStatus itl;
 150   TopExp_Explorer exp;
 151   for (exp.Init(S,Subtype); exp.More(); exp.Next()) {
 152     const Handle(BRepCheck_Result)& res = Ana.Result(exp.Current());
 153     const TopoDS_Shape& sub = exp.Current();
 154     for (res->InitContextIterator();
 155          res->MoreShapeInContext();
 156          res->NextShapeInContext()) {
 157       if (res->ContextualShape().IsSame(S) &&
 158           !Contains(theMap(sub),S)) {
 159         theMap(sub).Append(S);
 160         itl.Initialize(res->StatusOnShape());
 161         if (itl.Value() != BRepCheck_NoError) {
 162           if (!FindNamed(sub,Name)) {
 163             nbfaulty++;
 164             Name = (char*)malloc(18*sizeof(char));
 165             Sprintf(Name,"%s%d",checkfaultyname,nbfaulty);
 166             DBRep::Set(Name,sub);
 167             lfaulty.Append(Draw::Get((Standard_CString&)Name));
 168           }
 169           OS << "Shape " << Name << " ";
 170           if (!FindNamed(S,Name)) {
 171             nbfaulty++;
 172             Name = (char*)malloc(18*sizeof(char));
 173             Sprintf(Name,"%s%d",checkfaultyname,nbfaulty);
 174             DBRep::Set(Name,S);
 175             lfaulty.Append(Draw::Get((Standard_CString&)Name));
 176           }
 177           OS << " on shape " << Name << " :\n";
 178           for (;itl.More(); itl.Next()) {
 179             BRepCheck::Print(itl.Value(),OS);
 180           }
 181         }
 182         break;
 183       }
 184     }
 185   }
 186 }
 187
 188
 189 //=======================================================================
 190 //function : Print
 191 //=======================================================================
 192 static void Print(Standard_OStream& OS,
 193                   const BRepCheck_Analyzer& Ana,
 194                   const TopoDS_Shape& S)
 195 {
 196   for (TopoDS_Iterator iter(S); iter.More(); iter.Next()) {
 197     Print(OS,Ana,iter.Value());
 198   }
 199
 200   char* Name;
 201   TopAbs_ShapeEnum styp = S.ShapeType();
 202   BRepCheck_ListIteratorOfListOfStatus itl;
 203   if (!Ana.Result(S).IsNull() && !theMap.IsBound(S)) {
 204     itl.Initialize(Ana.Result(S)->Status());
 205     if (itl.Value() != BRepCheck_NoError) {
 206       if (!FindNamed(S,Name)) {
 207         nbfaulty++;
 208         Name = (char*)malloc(18*sizeof(char));
 209         Sprintf(Name,"%s%d",checkfaultyname,nbfaulty);
 210         DBRep::Set(Name,S);
 211         lfaulty.Append(Draw::Get((Standard_CString&)Name));
 212       }
 213       OS << "On Shape " << Name << " :\n";
 214
 215       for (;itl.More(); itl.Next()) {
 216         if (itl.Value() != BRepCheck_NoError)
 217         BRepCheck::Print(itl.Value(),OS);
 218       }
 219     }
 220   }
 221   if (!theMap.IsBound(S)) {
 222     TopTools_ListOfShape thelist;
 223     theMap.Bind(S, thelist);
 224   }
 225
 226   switch (styp) {
 227   case TopAbs_EDGE:
 228     PrintSub(OS,Ana,S,TopAbs_VERTEX);
 229     break;
 230   case TopAbs_WIRE:
 231     PrintSub(OS,Ana,S,TopAbs_EDGE);
 232     PrintSub(OS,Ana,S,TopAbs_VERTEX);
 233     break;
 234   case TopAbs_FACE:
 235     PrintSub(OS,Ana,S,TopAbs_WIRE);
 236     PrintSub(OS,Ana,S,TopAbs_EDGE);
 237     PrintSub(OS,Ana,S,TopAbs_VERTEX);
 238     break;
 239   case TopAbs_SHELL:
 240 //    PrintSub(OS,Ana,S,TopAbs_FACE);
 241     break;
 242   case TopAbs_SOLID:
 243 //    PrintSub(OS,Ana,S,TopAbs_EDGE);
 244     PrintSub(OS,Ana,S,TopAbs_SHELL);
 245     break;
 246   default:
 247     break;
 248   }
 249
 250 }
 251
 252 //=======================================================================
 253 //function : computetolerance
 254 //purpose  :
 255 //=======================================================================
 256 static Standard_Integer computetolerance(Draw_Interpretor& di,
 257                                          Standard_Integer narg, const char** a)
 258 {
 259   if (narg < 2) {
 260     //std::cout << "Usage: computetolerance shape" << std::endl;
 261     di << "Usage: computetolerance shape\n";
 262     return 1;
 263   }
 264   TopoDS_Shape S = DBRep::Get(a[1]);
 265   Standard_Real tol;
 266   if (S.ShapeType() == TopAbs_EDGE) {
 267     BRepCheck_Edge bce(TopoDS::Edge(S));
 268     tol=bce.Tolerance();
 269     //std::cout<< "Tolerance de " << (void*) &(*S.TShape()) << " : " << tol << std::endl;
 270     Standard_SStream aSStream1;
 271     aSStream1<< "Tolerance de " << (void*) &(*S.TShape()) << " : " << tol << "\n";
 272     di << aSStream1;
 273   }
 274   else {
 275     TopTools_MapOfShape theEdges;
 276     TopExp_Explorer exp;
 277     for (exp.Init(S, TopAbs_EDGE); exp.More(); exp.Next()) {
 278       if (theEdges.Add(exp.Current())) {
 279         BRepCheck_Edge bce(TopoDS::Edge(exp.Current()));
 280         tol=bce.Tolerance();
 281         //std::cout<< "Tolerance de " << (void*) &(*exp.Current().TShape()) << " : " << tol << "\n";
 282         Standard_SStream aSStream2;
 283         aSStream2<< "Tolerance de " << (void*) &(*exp.Current().TShape()) << " : " << tol << "\n";
 284         di << aSStream2;
 285       }
 286     }
 287     //std::cout << std::endl;
 288     di << "\n";
 289   }
 290   return 0;
 291 }
 292
 293 //=======================================================================
 294 //function : checksection
 295 //purpose  : Checks the closure of a section line
 296 //=======================================================================
 297 static Standard_Integer checksection(Draw_Interpretor& di,
 298                                      Standard_Integer narg, const char** a)
 299 {
 300   if (narg < 2) {
 301     di << a[0] << " shape [-r <ref_val>]\n";
 302     return 1;
 303   }
 304
 305   Standard_Integer aCompareValue = -1;
 306   TopoDS_Shape S = DBRep::Get(a[1]);
 307
 308   for (Standard_Integer anAI = 2; anAI < narg; anAI++)
 309   {
 310     if (!strcmp(a[anAI], "-r"))
 311     {
 312       aCompareValue = Draw::Atoi(a[++anAI]);
 313     }
 314     else
 315     {
 316       di << "Error: Wrong option" << a[anAI] << "\n";
 317     }
 318   }
 319
 320   TopTools_MapOfShape theVertices;
 321   TopExp_Explorer exp;
 322   for (exp.Init(S, TopAbs_VERTEX); exp.More(); exp.Next()) {
 323     if (!theVertices.Add(exp.Current()))
 324       theVertices.Remove(exp.Current());
 325   }
 326   //std::cout << " nb alone Vertices : " << theVertices.Extent() << std::endl;
 327   di << " nb alone Vertices : " << theVertices.Extent() << "\n";
 328
 329   if (aCompareValue >= 0)
 330   {
 331     if (theVertices.Extent() == aCompareValue)
 332     {
 333       di << "Section is OK\n";
 334     }
 335     else
 336     {
 337       di << "Error: "<< aCompareValue << " vertices are expected but " <<
 338                         theVertices.Extent() << " are found.\n";
 339     }
 340   }
 341
 342   char Name[32];
 343   Standard_Integer ipp=0;
 344   TopTools_MapIteratorOfMapOfShape itvx;
 345   for (itvx.Initialize(theVertices); itvx.More(); itvx.Next()) {
 346     ipp++;
 347     Sprintf(Name,"alone_%d",ipp);
 348     DBRep::Set(Name, itvx.Key());
 349     //std::cout << Name << " " ;
 350     di << Name << " " ;
 351   }
 352   //std::cout << std::endl;
 353   di << "\n";
 354   return 0;
 355 }
 356 //=======================================================================
 357
 358 //=======================================================================
 359 //function : checkdiff
 360 //purpose  : Checks the differences beetween a result and his arguments
 361 //=======================================================================
 362 static Standard_Integer checkdiff(Draw_Interpretor& di,
 363                                   Standard_Integer narg, const char** a)
 364 {
 365   const char* syntaxe = "checkdiff arg1 [arg2..argn] result [closedSolid (0/1)] [geomCtrl (1/0)]";
 366   if (narg < 3) {
 367     if (narg==2) {
 368       Standard_Integer bcrtrace=Draw::Atoi(a[narg-1]);
 369       bcrtrace=BRepCheck_Trace(bcrtrace);
 370       //std::cout << "BRepCheck_Trace : " << bcrtrace << std::endl;
 371       di << "BRepCheck_Trace : " << bcrtrace << "\n";
 372       //std::cout << syntaxe << std::endl;
 373       di << syntaxe << "\n";
 374       return 0;
 375     }
 376     //std::cout << syntaxe << std::endl;
 377     di << syntaxe << "\n";
 378     return 1;
 379   }
 380
 381   Standard_Integer lastArg=narg-2;
 382   Standard_Boolean closedSolid = Standard_False;
 383   Standard_Boolean geomCtrl = Standard_True;
 384   TopoDS_Shape resu = DBRep::Get(a[narg-1]);
 385   if (resu.IsNull()) {
 386     if (narg < 4) {
 387       //std::cout << syntaxe << std::endl;
 388       di << syntaxe << "\n";
 389       return 1;
 390     }
 391     closedSolid = Draw::Atoi(a[narg-1]) != 0;
 392     resu = DBRep::Get(a[narg-2]);
 393     lastArg=narg-3;
 394     if (resu.IsNull()) {
 395       if (narg < 5) {
 396         //std::cout << syntaxe << std::endl;
 397         di << syntaxe << "\n";
 398         return 1;
 399       }
 400       geomCtrl=closedSolid;
 401       closedSolid = Draw::Atoi(a[narg-2]) != 0;
 402       resu = DBRep::Get(a[narg-3]);
 403       lastArg=narg-4;
 404       if (resu.IsNull()) {
 405         //std::cout << syntaxe << std::endl;
 406         di << syntaxe << "\n";
 407         return 1;
 408       }
 409     }
 410   }
 411
 412   TopTools_ListOfShape lesArgs;
 413   for (Standard_Integer id=1; id <=lastArg; id++) {
 414     lesArgs.Append(DBRep::Get(a[id]));
 415   }
 416
 417   if (BRepAlgo::IsValid(lesArgs, resu, closedSolid, geomCtrl)) {
 418     //std::cout << "Difference is Valid." << std::endl;
 419     di << "Difference is Valid.\n";
 420   } else {
 421     //std::cout << "error : Difference is Not Valid !" << std::endl;
 422     di << "error : Difference is Not Valid !\n";
 423   }
 424
 425   return 0;
 426 }
 427 //=======================================================================
 428
 429 //  Modified by skv - Tue Apr 27 13:38:44 2004 Begin
 430 //=======================================================================
 431 //function : CHK
 432 //purpose  : Checks a shape
 433 //=======================================================================
 434
 435 // static Standard_Integer CHK(Draw_Interpretor& theCommands,
 436 //                          Standard_Integer narg, const char** a)
 437 // {
 438 //   if (narg < 2) {
 439 //     return 1;
 440 //   }
 441
 442 //   Standard_Boolean doprint = Standard_True;
 443 //   if (narg == 3) { if (!strcmp(a[2],"-short")) doprint = Standard_False; }
 444
 445 //   TopoDS_Shape S = DBRep::Get(a[1]);
 446 //   if (S.IsNull()) {
 447 //     std::cout<<"not a topological shape"<<std::endl;
 448 //     return 1;
 449 //   }
 450
 451 //   Standard_Boolean GeomCtrl = Standard_True;
 452 //   if (!strcasecmp(a[0],"CHECKTOPSHAPE")) {
 453 //     GeomCtrl = Standard_False;
 454 //   }
 455
 456 //   BRepCheck_Analyzer ana(S,GeomCtrl);
 457 //   if (ana.IsValid()) {
 458 //     theCommands<<"This shape seems to be valid";
 459 //   }
 460 //   else {
 461 //     theMap.Clear();
 462 //     nbfaulty = 0;
 463 //     lfaulty.Clear();
 464 //     theMap.Clear();
 465 //     if (doprint) {
 466 //       Print(cout,ana,S);
 467 //       std::cout<<"\n";
 468 //       theMap.Clear();
 469 //       if (nbfaulty !=0)
 470 //      std::cout<<"Faulty shapes in variables "<<checkfaultyname<<"1 to "<<checkfaultyname<<nbfaulty<<" \n";
 471 //       std::cout<<std::endl;
 472 //     }
 473 //     else {
 474 //       theCommands<<"This shape has faulty shapes";
 475 //     }
 476 //   }
 477 //   return 0;
 478 // }
 479
 480 //=======================================================================
 481 //function : ContextualDump
 482 //purpose  : Contextual (modeling) style of output.
 483 //=======================================================================
 484
 485 //void ContextualDump(const BRepCheck_Analyzer &theAna,
 486 //                  const TopoDS_Shape       &theShape)
 487 void ContextualDump(Draw_Interpretor& theCommands,
 488                     const BRepCheck_Analyzer &theAna,
 489                     const TopoDS_Shape       &theShape)
 490 {
 491   theMap.Clear();
 492   nbfaulty = 0;
 493   lfaulty.Clear();
 494
 495   //Print(cout, theAna, theShape);
 496   Standard_SStream aSStream;
 497   Print(aSStream, theAna, theShape);
 498   theCommands << aSStream;
 499   //std::cout<<"\n";
 500   theCommands<<"\n";
 501   theMap.Clear();
 502
 503   if (nbfaulty !=0)
 504     theCommands<<"Faulty shapes in variables "<<checkfaultyname<<"1 to "<<checkfaultyname<<nbfaulty<<" \n";
 505     //std::cout<<"Faulty shapes in variables "<<checkfaultyname<<"1 to "<<checkfaultyname<<nbfaulty<<" \n";
 506
 507   //std::cout<<std::endl;
 508   theCommands<<"\n";
 509 }
 510
 511
 512 //=======================================================================
 513 //function : FillProblems
 514 // purpose : auxilary for StructuralDump
 515 //=======================================================================
 516 static void FillProblems(const BRepCheck_Status stat,
 517                          Handle(TColStd_HArray1OfInteger)& NbProblems)
 518 {
 519
 520   const Standard_Integer anID = static_cast<Standard_Integer> (stat);
 521
 522   if((NbProblems->Upper() < anID) || (NbProblems->Lower() > anID))
 523     return;
 524
 525   NbProblems->SetValue(anID, NbProblems->Value(anID)+1);
 526
 527 }
 528
 529
 530 //=======================================================================
 531 //function : GetProblemSub
 532 // purpose : auxilary for StructuralDump
 533 //=======================================================================
 534 static void GetProblemSub(const BRepCheck_Analyzer& Ana,
 535                           const TopoDS_Shape& Shape,
 536                           Handle(TopTools_HSequenceOfShape)& sl,
 537                           Handle(TColStd_HArray1OfInteger)& NbProblems,
 538                           const TopAbs_ShapeEnum Subtype)
 539 {
 540   BRepCheck_ListIteratorOfListOfStatus itl;
 541   TopExp_Explorer exp;
 542   for (exp.Init(Shape,Subtype); exp.More(); exp.Next()) {
 543     const Handle(BRepCheck_Result)& res = Ana.Result(exp.Current());
 544
 545     const TopoDS_Shape& sub = exp.Current();
 546     for (res->InitContextIterator();
 547          res->MoreShapeInContext();
 548          res->NextShapeInContext()) {
 549       if (res->ContextualShape().IsSame(Shape) &&
 550           !Contains(theMap(sub),Shape)) {
 551         theMap(sub).Append(Shape);
 552         itl.Initialize(res->StatusOnShape());
 553
 554         if (itl.Value() != BRepCheck_NoError) {
 555           Standard_Integer ii = 0;
 556
 557           for(ii=1; ii<=sl->Length(); ii++)
 558             if(sl->Value(ii).IsSame(sub)) break;
 559
 560           if(ii>sl->Length()) {
 561             sl->Append(sub);
 562             FillProblems(itl.Value(),NbProblems);
 563           }
 564           for(ii=1; ii<=sl->Length(); ii++)
 565             if(sl->Value(ii).IsSame(Shape)) break;
 566           if(ii>sl->Length()) {
 567             sl->Append(Shape);
 568             FillProblems(itl.Value(),NbProblems);
 569           }
 570         }
 571         break;
 572       }
 573     }
 574   }
 575 }
 576
 577
 578 //=======================================================================
 579 //function : GetProblemShapes
 580 // purpose : auxilary for StructuralDump
 581 //=======================================================================
 582 static void GetProblemShapes(const BRepCheck_Analyzer& Ana,
 583                              const TopoDS_Shape& Shape,
 584                              Handle(TopTools_HSequenceOfShape)& sl,
 585                              Handle(TColStd_HArray1OfInteger)& NbProblems)
 586 {
 587   for (TopoDS_Iterator iter(Shape); iter.More(); iter.Next()) {
 588     GetProblemShapes(Ana,iter.Value(),sl, NbProblems);
 589   }
 590   TopAbs_ShapeEnum styp = Shape.ShapeType();
 591   BRepCheck_ListIteratorOfListOfStatus itl;
 592   if (!Ana.Result(Shape).IsNull() && !theMap.IsBound(Shape)) {
 593     itl.Initialize(Ana.Result(Shape)->Status());
 594
 595     if (itl.Value() != BRepCheck_NoError) {
 596       sl->Append(Shape);
 597       FillProblems(itl.Value(),NbProblems);
 598     }
 599   }
 600   if (!theMap.IsBound(Shape)) {
 601     TopTools_ListOfShape thelist;
 602     theMap.Bind(Shape, thelist);
 603   }
 604
 605   switch (styp) {
 606   case TopAbs_EDGE:
 607     GetProblemSub(Ana, Shape, sl, NbProblems, TopAbs_VERTEX);
 608     break;
 609   case TopAbs_FACE:
 610     GetProblemSub(Ana, Shape, sl, NbProblems, TopAbs_WIRE);
 611     GetProblemSub(Ana, Shape, sl, NbProblems, TopAbs_EDGE);
 612     GetProblemSub(Ana, Shape, sl, NbProblems, TopAbs_VERTEX);
 613     break;
 614   case TopAbs_SHELL:
 615     break;
 616   case TopAbs_SOLID:
 617     GetProblemSub(Ana, Shape, sl, NbProblems, TopAbs_SHELL);
 618     break;
 619   default:
 620     break;
 621   }
 622
 623 }
 624
 625 //=======================================================================
 626 //function : StructuralDump
 627 //purpose  : Structural (data exchange) style of output.
 628 //=======================================================================
 629
 630 //void StructuralDump(const BRepCheck_Analyzer &theAna,
 631 //                    const Standard_CString   ShName,
 632 //                    const Standard_CString   Pref,
 633 //                  const TopoDS_Shape       &theShape)
 634 void StructuralDump(Draw_Interpretor& theCommands,
 635                     const BRepCheck_Analyzer &theAna,
 636                     const Standard_CString   ShName,
 637                     const Standard_CString   Pref,
 638                     const TopoDS_Shape       &theShape)
 639 {
 640   Standard_Integer i;
 641   //std::cout << "StructuralDump" << std::endl;
 642   //std::cout << " -- The Shape " << ShName << " has problems :"<<std::endl;
 643   //std::cout<<"  Check                                    Count"<<std::endl;
 644   //std::cout<<" ------------------------------------------------"<<std::endl;
 645   theCommands << " -- The Shape " << ShName << " has problems :\n";
 646   theCommands<<"  Check                                    Count\n";
 647   theCommands<<" ------------------------------------------------\n";
 648
 649   Handle(TColStd_HArray1OfInteger) NbProblems = new
 650                               TColStd_HArray1OfInteger(1,NumberOfStatus);
 651   for(i=1; i<=NumberOfStatus; i++) NbProblems->SetValue(i,0);
 652   Handle(TopTools_HSequenceOfShape) sl,slv,sle,slw,slf,sls,slo;
 653   sl = new TopTools_HSequenceOfShape();
 654   theMap.Clear();
 655   GetProblemShapes(theAna, theShape, sl, NbProblems);
 656   theMap.Clear();
 657
 658   Standard_Integer aProblemID = static_cast<Standard_Integer>(BRepCheck_InvalidPointOnCurve);
 659   if(NbProblems->Value(aProblemID) > 0)
 660     theCommands<<"  Invalid Point on Curve ................... "<<NbProblems->Value(aProblemID)<<"\n";
 661
 662   aProblemID = static_cast<Standard_Integer>(BRepCheck_InvalidPointOnCurveOnSurface);
 663   if(NbProblems->Value(aProblemID)>0)
 664     theCommands<<"  Invalid Point on CurveOnSurface .......... "<<NbProblems->Value(aProblemID)<<"\n";
 665
 666   aProblemID = static_cast<Standard_Integer>(BRepCheck_InvalidPointOnSurface);
 667   if(NbProblems->Value(aProblemID)>0)
 668     theCommands<<"  Invalid Point on Surface ................. "<<NbProblems->Value(aProblemID)<<"\n";
 669
 670   aProblemID = static_cast<Standard_Integer>(BRepCheck_No3DCurve);
 671   if(NbProblems->Value(aProblemID)>0)
 672     theCommands<<"  No 3D Curve .............................. "<<NbProblems->Value(aProblemID)<<"\n";
 673
 674   aProblemID = static_cast<Standard_Integer>(BRepCheck_Multiple3DCurve);
 675   if(NbProblems->Value(aProblemID)>0)
 676     theCommands<<"  Multiple 3D Curve ........................ "<<NbProblems->Value(aProblemID)<<"\n";
 677
 678   aProblemID = static_cast<Standard_Integer>(BRepCheck_Invalid3DCurve);
 679   if(NbProblems->Value(aProblemID)>0)
 680     theCommands<<"  Invalid 3D Curve ......................... "<<NbProblems->Value(aProblemID)<<"\n";
 681
 682   aProblemID = static_cast<Standard_Integer>(BRepCheck_NoCurveOnSurface);
 683   if(NbProblems->Value(aProblemID)>0)
 684     theCommands<<"  No Curve on Surface ...................... "<<NbProblems->Value(aProblemID)<<"\n";
 685
 686   aProblemID = static_cast<Standard_Integer>(BRepCheck_InvalidCurveOnSurface);
 687   if(NbProblems->Value(aProblemID)>0)
 688     theCommands<<"  Invalid Curve on Surface ................. "<<NbProblems->Value(aProblemID)<<"\n";
 689
 690   aProblemID = static_cast<Standard_Integer>(BRepCheck_InvalidCurveOnClosedSurface);
 691   if(NbProblems->Value(aProblemID)>0)
 692     theCommands<<"  Invalid Curve on closed Surface .......... "<<NbProblems->Value(aProblemID)<<"\n";
 693
 694   aProblemID = static_cast<Standard_Integer>(BRepCheck_InvalidSameRangeFlag);
 695   if(NbProblems->Value(aProblemID)>0)
 696     theCommands<<"  Invalid SameRange Flag ................... "<<NbProblems->Value(aProblemID)<<"\n";
 697
 698   aProblemID = static_cast<Standard_Integer>(BRepCheck_InvalidSameParameterFlag);
 699   if(NbProblems->Value(aProblemID)>0)
 700     theCommands<<"  Invalid SameParameter Flag ............... "<<NbProblems->Value(aProblemID)<<"\n";
 701
 702   aProblemID = static_cast<Standard_Integer>(BRepCheck_InvalidDegeneratedFlag);
 703   if(NbProblems->Value(aProblemID)>0)
 704     theCommands<<"  Invalid Degenerated Flag ................. "<<NbProblems->Value(aProblemID)<<"\n";
 705
 706   aProblemID = static_cast<Standard_Integer>(BRepCheck_FreeEdge);
 707   if(NbProblems->Value(aProblemID)>0)
 708     theCommands<<"  Free Edge ................................ "<<NbProblems->Value(aProblemID)<<"\n";
 709
 710   aProblemID = static_cast<Standard_Integer>(BRepCheck_InvalidMultiConnexity);
 711   if(NbProblems->Value(aProblemID)>0)
 712     theCommands<<"  Invalid MultiConnexity ................... "<<NbProblems->Value(aProblemID)<<"\n";
 713
 714   aProblemID = static_cast<Standard_Integer>(BRepCheck_InvalidRange);
 715   if(NbProblems->Value(aProblemID)>0)
 716     theCommands<<"  Invalid Range ............................ "<<NbProblems->Value(aProblemID)<<"\n";
 717
 718   aProblemID = static_cast<Standard_Integer>(BRepCheck_EmptyWire);
 719   if(NbProblems->Value(aProblemID)>0)
 720     theCommands<<"  Empty Wire ............................... "<<NbProblems->Value(aProblemID)<<"\n";
 721
 722   aProblemID = static_cast<Standard_Integer>(BRepCheck_RedundantEdge);
 723   if(NbProblems->Value(aProblemID)>0)
 724     theCommands<<"  Redundant Edge ........................... "<<NbProblems->Value(aProblemID)<<"\n";
 725
 726   aProblemID = static_cast<Standard_Integer>(BRepCheck_SelfIntersectingWire);
 727   if(NbProblems->Value(aProblemID)>0)
 728     theCommands<<"  Self Intersecting Wire ................... "<<NbProblems->Value(aProblemID)<<"\n";
 729
 730   aProblemID = static_cast<Standard_Integer>(BRepCheck_NoSurface);
 731   if(NbProblems->Value(aProblemID)>0)
 732     theCommands<<"  No Surface ............................... "<<NbProblems->Value(aProblemID)<<"\n";
 733
 734   aProblemID = static_cast<Standard_Integer>(BRepCheck_InvalidWire);
 735   if(NbProblems->Value(aProblemID)>0)
 736     theCommands<<"  Invalid Wire ............................. "<<NbProblems->Value(aProblemID)<<"\n";
 737
 738   aProblemID = static_cast<Standard_Integer>(BRepCheck_RedundantWire);
 739   if(NbProblems->Value(aProblemID)>0)
 740     theCommands<<"  Redundant Wire ........................... "<<NbProblems->Value(aProblemID)<<"\n";
 741
 742   aProblemID = static_cast<Standard_Integer>(BRepCheck_IntersectingWires);
 743   if(NbProblems->Value(aProblemID)>0)
 744     theCommands<<"  Intersecting Wires ....................... "<<NbProblems->Value(aProblemID)<<"\n";
 745
 746   aProblemID = static_cast<Standard_Integer>(BRepCheck_InvalidImbricationOfWires);
 747   if(NbProblems->Value(aProblemID)>0)
 748     theCommands<<"  Invalid Imbrication of Wires ............. "<<NbProblems->Value(aProblemID)<<"\n";
 749
 750   aProblemID = static_cast<Standard_Integer>(BRepCheck_EmptyShell);
 751   if(NbProblems->Value(aProblemID)>0)
 752     theCommands<<"  Empty Shell .............................. "<<NbProblems->Value(aProblemID)<<"\n";
 753
 754   aProblemID = static_cast<Standard_Integer>(BRepCheck_RedundantFace);
 755   if(NbProblems->Value(aProblemID)>0)
 756     theCommands<<"  Redundant Face ........................... "<<NbProblems->Value(aProblemID)<<"\n";
 757
 758   aProblemID = static_cast<Standard_Integer>(BRepCheck_UnorientableShape);
 759   if(NbProblems->Value(aProblemID)>0)
 760     theCommands<<"  Unorientable Shape ....................... "<<NbProblems->Value(aProblemID)<<"\n";
 761
 762   aProblemID = static_cast<Standard_Integer>(BRepCheck_NotClosed);
 763   if(NbProblems->Value(aProblemID)>0)
 764     theCommands<<"  Not Closed ............................... "<<NbProblems->Value(aProblemID)<<"\n";
 765
 766   aProblemID = static_cast<Standard_Integer>(BRepCheck_NotConnected);
 767   if(NbProblems->Value(aProblemID)>0)
 768     theCommands<<"  Not Connected ............................ "<<NbProblems->Value(aProblemID)<<"\n";
 769
 770   aProblemID = static_cast<Standard_Integer>(BRepCheck_SubshapeNotInShape);
 771   if(NbProblems->Value(aProblemID)>0)
 772     theCommands<<"  Subshape not in Shape .................... "<<NbProblems->Value(aProblemID)<<"\n";
 773
 774   aProblemID = static_cast<Standard_Integer>(BRepCheck_BadOrientation);
 775   if(NbProblems->Value(aProblemID)>0)
 776     theCommands<<"  Bad Orientation .......................... "<<NbProblems->Value(aProblemID)<<"\n";
 777
 778   aProblemID = static_cast<Standard_Integer>(BRepCheck_BadOrientationOfSubshape);
 779   if(NbProblems->Value(aProblemID)>0)
 780     theCommands<<"  Bad Orientation of Subshape .............. "<<NbProblems->Value(aProblemID)<<"\n";
 781
 782   aProblemID = static_cast<Standard_Integer>(BRepCheck_InvalidToleranceValue);
 783   if(NbProblems->Value(aProblemID)>0)
 784     theCommands<<"  Invalid tolerance value................... "<<NbProblems->Value(aProblemID)<<"\n";
 785
 786   aProblemID = static_cast<Standard_Integer>(BRepCheck_InvalidPolygonOnTriangulation);
 787   if(NbProblems->Value(aProblemID)>0)
 788     theCommands<<"  Invalid polygon on triangulation.......... "<<NbProblems->Value(aProblemID)<<"\n";
 789
 790   aProblemID = static_cast<Standard_Integer>(BRepCheck_InvalidImbricationOfShells);
 791   if(NbProblems->Value(aProblemID)>0)
 792     theCommands<<"  Invalid Imbrication of Shells............. "<<NbProblems->Value(aProblemID)<<"\n";
 793
 794  aProblemID = static_cast<Standard_Integer>(BRepCheck_EnclosedRegion);
 795   if(NbProblems->Value(aProblemID)>0)
 796     theCommands<<"  Enclosed Region........................... "<<NbProblems->Value(aProblemID)<<"\n";
 797
 798   aProblemID = static_cast<Standard_Integer>(BRepCheck_CheckFail);
 799   if(NbProblems->Value(aProblemID)>0)
 800     theCommands<<"  checkshape failure........................ "<<NbProblems->Value(aProblemID)<<"\n";
 801
 802
 803
 804   theCommands<<" ------------------------------------------------\n";
 805   theCommands<<"*** Shapes with problems : "<<sl->Length()<<"\n";
 806
 807   slv = new TopTools_HSequenceOfShape();
 808   sle = new TopTools_HSequenceOfShape();
 809   slw = new TopTools_HSequenceOfShape();
 810   slf = new TopTools_HSequenceOfShape();
 811   sls = new TopTools_HSequenceOfShape();
 812   slo = new TopTools_HSequenceOfShape();
 813
 814   for(i=1; i<=sl->Length(); i++) {
 815     TopoDS_Shape shi = sl->Value(i);
 816     TopAbs_ShapeEnum sti = shi.ShapeType();
 817     switch (sti) {
 818       case TopAbs_VERTEX : slv->Append (shi); break;
 819       case TopAbs_EDGE   : sle->Append (shi); break;
 820       case TopAbs_WIRE   : slw->Append (shi); break;
 821       case TopAbs_FACE   : slf->Append (shi); break;
 822       case TopAbs_SHELL  : sls->Append (shi); break;
 823       case TopAbs_SOLID  : slo->Append (shi); break;
 824       default            : break;
 825     }
 826   }
 827
 828   BRep_Builder B;
 829   if(slv->Length()>0) {
 830     TopoDS_Compound comp;
 831     B.MakeCompound(comp);
 832     Standard_Integer nb = slv->Length();
 833     for(i=1; i<=nb; i++)
 834       B.Add(comp,slv->Value(i));
 835     char aName[20];
 836     Sprintf(aName,"%s_v",Pref);
 837     DBRep::Set(aName,comp);
 838     //std::cout<<"VERTEX        : "<<(nb > 9 ? "" : " ")<<nb<<" Items -> compound named "<<aName<<std::endl;
 839     if (nb > 9)
 840       theCommands<<"VERTEX      : "<<nb<<" Items -> compound named "<<aName<<"\n";
 841     else
 842       theCommands<<"VERTEX      :  "<<nb<<" Items -> compound named "<<aName<<"\n";
 843   }
 844   if(sle->Length()>0) {
 845     TopoDS_Compound comp;
 846     B.MakeCompound(comp);
 847     Standard_Integer nb = sle->Length();
 848     for(i=1; i<=nb; i++)
 849       B.Add(comp,sle->Value(i));
 850     char aName[20];
 851     Sprintf(aName,"%s_e",Pref);
 852     DBRep::Set(aName,comp);
 853     //std::cout<<"EDGE  : "<<(nb > 9 ? "" : " ")<<nb<<" Items -> compound named "<<aName<<std::endl;
 854     if (nb > 9)
 855       theCommands<<"EDGE        : "<<nb<<" Items -> compound named "<<aName<<"\n";
 856     else
 857       theCommands<<"EDGE        :  "<<nb<<" Items -> compound named "<<aName<<"\n";
 858   }
 859   if(slw->Length()>0) {
 860     TopoDS_Compound comp;
 861     B.MakeCompound(comp);
 862     Standard_Integer nb = slw->Length();
 863     for(i=1; i<=nb; i++)
 864       B.Add(comp,slw->Value(i));
 865     char aName[20];
 866     Sprintf(aName,"%s_w",Pref);
 867     DBRep::Set(aName,comp);
 868     //std::cout<<"WIRE  : "<<(nb > 9 ? "" : " ")<<nb<<" Items -> compound named "<<aName<<std::endl;
 869     if (nb > 9)
 870       theCommands<<"WIRE        : "<<nb<<" Items -> compound named "<<aName<<"\n";
 871     else
 872       theCommands<<"WIRE        :  "<<nb<<" Items -> compound named "<<aName<<"\n";
 873   }
 874   if(slf->Length()>0) {
 875     TopoDS_Compound comp;
 876     B.MakeCompound(comp);
 877     Standard_Integer nb = slf->Length();
 878     for(i=1; i<=nb; i++)
 879       B.Add(comp,slf->Value(i));
 880     char aName[20];
 881     Sprintf(aName,"%s_f",Pref);
 882     DBRep::Set(aName,comp);
 883     //std::cout<<"FACE  : "<<(nb > 9 ? "" : " ")<<nb<<" Items -> compound named "<<aName<<std::endl;
 884     if (nb > 9)
 885       theCommands<<"FACE        : "<<nb<<" Items -> compound named "<<aName<<"\n";
 886     else
 887       theCommands<<"FACE        :  "<<nb<<" Items -> compound named "<<aName<<"\n";
 888   }
 889   if(sls->Length()>0) {
 890     TopoDS_Compound comp;
 891     B.MakeCompound(comp);
 892     Standard_Integer nb = sls->Length();
 893     for(i=1; i<=nb; i++)
 894       B.Add(comp,sls->Value(i));
 895     char aName[20];
 896     Sprintf(aName,"%s_s",Pref);
 897     DBRep::Set(aName,comp);
 898     //std::cout<<"SHELL : "<<(nb > 9 ? "" : " ")<<nb<<" Items -> compound named "<<aName<<std::endl;
 899     if (nb > 9)
 900       theCommands<<"SHELL       : "<<nb<<" Items -> compound named "<<aName<<"\n";
 901     else
 902       theCommands<<"SHELL       :  "<<nb<<" Items -> compound named "<<aName<<"\n";
 903   }
 904   if(slo->Length()>0) {
 905     TopoDS_Compound comp;
 906     B.MakeCompound(comp);
 907     Standard_Integer nb = slo->Length();
 908     for(i=1; i<=nb; i++)
 909       B.Add(comp,slo->Value(i));
 910     char aName[20];
 911     Sprintf(aName,"%s_o",Pref);
 912     DBRep::Set(aName,comp);
 913     //std::cout<<"SOLID : "<<(nb > 9 ? "" : " ")<<nb<<" Items -> compound named "<<aName<<std::endl;
 914     if (nb > 9)
 915       theCommands<<"SOLID       : "<<nb<<" Items -> compound named "<<aName<<"\n";
 916     else
 917       theCommands<<"SOLID       :  "<<nb<<" Items -> compound named "<<aName<<"\n";
 918   }
 919 }
 920
 921 //=======================================================================
 922 //function : checkshape
 923 //purpose  : Checks a shape
 924 //=======================================================================
 925
 926 static Standard_Integer checkshape(Draw_Interpretor& theCommands,
 927                                    Standard_Integer narg, const char** a)
 928 {
 929   if (narg == 1) {
 930     //std::cout << std::endl;
 931     //std::cout << "Usage : checkshape [-top] shape [result] [-short]" << std::endl;
 932     //std::cout << std::endl;
 933     //std::cout << "Where :" << std::endl;
 934     //std::cout << "   -top   - check topology only." << std::endl;
 935     //std::cout << "   shape  - the name of the shape to test." << std::endl;
 936     //std::cout << "   result - the prefix of the output shape names. If it is used, structural" << std::endl;
 937     //std::cout << "            output style will be used. Otherwise - contextual one." << std::endl;
 938     //std::cout << "   -short - short description of check." << std::endl;
 939     theCommands << "\n";
 940     theCommands << "Usage : checkshape [-top] shape [result] [-short]\n";
 941     theCommands << "\n";
 942     theCommands << "Where :\n";
 943     theCommands << "   -top   - check topology only.\n";
 944     theCommands << "   shape  - the name of the shape to test.\n";
 945     theCommands << "   result - the prefix of the output shape names. If it is used, structural\n";
 946     theCommands << "            output style will be used. Otherwise - contextual one.\n";
 947     theCommands << "   -short - short description of check.\n";
 948
 949     return 0;
 950   }
 951
 952   if (narg > 5) {
 953     //std::cout << "Invalid number of args!!!" << std::endl;
 954     //std::cout << "No args to have help." << std::endl;
 955     theCommands << "Invalid number of args!!!\n";
 956     theCommands << "No args to have help.\n";
 957
 958     return 1;
 959   }
 960
 961   Standard_Boolean aGeomCtrl = Standard_True;
 962   Standard_Integer aCurInd  = 1;
 963
 964   if (!strcmp(a[1],"-top")) {
 965     aGeomCtrl = Standard_False;
 966     aCurInd++;
 967   }
 968
 969   if (aCurInd > narg - 1) {
 970     //std::cout << "Invalid number of args!!!" << std::endl;
 971     //std::cout << "No args to have help." << std::endl;
 972     theCommands << "Invalid number of args!!!\n";
 973     theCommands << "No args to have help.\n";
 974
 975     return 1;
 976   }
 977
 978   Standard_CString aShapeName = a[aCurInd];
 979   TopoDS_Shape     aShape     = DBRep::Get(aShapeName);
 980
 981   if (aShape.IsNull()) {
 982     //std::cout << a[aCurInd] << " is not a topological shape!!!" << std::endl;
 983     theCommands << a[aCurInd] << " is not a topological shape!!!\n";
 984
 985     return 1;
 986   }
 987
 988   Standard_Boolean IsShortDump   = Standard_False;
 989   Standard_Boolean IsContextDump = Standard_True;
 990   Standard_Integer aBackInd      = narg - 1;
 991
 992   if (aCurInd < aBackInd) {
 993     if (!strcmp(a[aBackInd],"-short")) {
 994       IsShortDump = Standard_True;
 995       aBackInd--;
 996     }
 997   }
 998
 999   if (aCurInd < aBackInd - 1) {
1000     //std::cout << "Invalid number of args!!!" << std::endl;
1001     //std::cout << "No args to have help." << std::endl;
1002     theCommands << "Invalid number of args!!!\n";
1003     theCommands << "No args to have help.\n";
1004
1005     return 1;
1006   } else if (aCurInd < aBackInd) {
1007     IsContextDump = Standard_False;
1008   }
1009
1010   try {
1011     OCC_CATCH_SIGNALS
1012     BRepCheck_Analyzer anAna(aShape,aGeomCtrl);
1013     Standard_Boolean   isValid = anAna.IsValid();
1014
1015     if (isValid) {
1016       if (IsContextDump) {
1017         theCommands << "This shape seems to be valid";
1018       } else {
1019         theCommands << " -- The Shape " << aShapeName << " looks OK";
1020       }
1021     } else {
1022       if (IsShortDump) {
1023         theCommands<<"This shape has faulty shapes";
1024       } else {
1025         if (IsContextDump) {
1026           //ContextualDump(anAna, aShape);
1027           ContextualDump(theCommands, anAna, aShape);
1028         } else {
1029           Standard_CString aPref = a[aCurInd+1];
1030           //StructuralDump(anAna, aShapeName, aPref, aShape);
1031           StructuralDump(theCommands, anAna, aShapeName, aPref, aShape);
1032         }
1033       }
1034     }
1035   }
1036   catch(Standard_Failure const& anException) {
1037     theCommands<<"checkshape exception : ";
1038     theCommands << anException.GetMessageString();
1039     theCommands<<"\n";
1040     return 1;
1041   }
1042
1043   return 0;
1044 }
1045 //  Modified by skv - Tue Apr 27 13:38:24 2004 End
1046 /***************************************************************/
1047 static void InitEpsSurf(Standard_Real& epsnl,Standard_Real& epsdis, Standard_Real& epsangk1,
1048                         Standard_Real& epsangk2, Standard_Real& epsangn1,
1049                         Standard_Real& perce,Standard_Real& maxlen )
1050 {
1051   epsnl  = 0.001;
1052   epsdis = 0.001;
1053   epsangk1 = 0.001 ;
1054   epsangk2 = 0.001;
1055   epsangn1 = 0.001 ;
1056   perce = 0.01;
1057   maxlen = 10000 ;
1058 }
1059
1060 static Standard_Integer shapeG1continuity (Draw_Interpretor& di, Standard_Integer n, const char** a)
1061
1062 { Standard_Real epsnl,epsC0, epsC1, epsC2, epsG1,percent,maxlen;
1063   Standard_Integer nbeval;
1064   InitEpsSurf(epsnl,epsC0, epsC1, epsC2, epsG1,percent,maxlen);
1065   Standard_Boolean ISG1=Standard_True;
1066   if (n<4) return 1;
1067   TopoDS_Face face1,face2;
1068   Standard_Real f1,f2,l1,l2;
1069   TopoDS_Shape shape = DBRep::Get(a[1],TopAbs_SHAPE);
1070   if (shape.IsNull()) return 1;
1071   TopoDS_Shape  edge  = DBRep::Get(a[2],TopAbs_EDGE);
1072   if (edge.IsNull()) return 1;
1073 // calcul des deux faces
1074   TopTools_IndexedDataMapOfShapeListOfShape lface;
1075   TopExp::MapShapesAndAncestors(shape,TopAbs_EDGE,TopAbs_FACE,lface);
1076   const TopTools_ListOfShape& lfac = lface.FindFromKey(edge);
1077
1078   Standard_Integer nelem= lfac.Extent();
1079   if(nelem!=2) return 1;
1080   TopTools_ListIteratorOfListOfShape It;
1081   It.Initialize(lfac);
1082   face1=TopoDS::Face(It.Value());
1083   It.Next();
1084   face2=TopoDS::Face(It.Value());
1085
1086   Standard_Boolean IsSeam = face1.IsEqual(face2);
1087
1088 // calcul des deux pcurves
1089   const Handle(Geom2d_Curve) c1 = BRep_Tool::CurveOnSurface
1090                                (TopoDS::Edge(edge),face1,f1,l1);
1091   if (c1.IsNull()) return 1;
1092
1093   if (IsSeam)
1094     edge.Reverse();
1095   const  Handle(Geom2d_Curve) c2  = BRep_Tool::CurveOnSurface
1096                                (TopoDS::Edge(edge),face2,f2,l2);
1097   if (c2.IsNull()) return 1;
1098
1099    Handle(Geom2d_Curve) curv1= new Geom2d_TrimmedCurve(c1,f1,l1);
1100
1101    Handle(Geom2d_Curve) curv2= new Geom2d_TrimmedCurve(c2,f2,l2);
1102
1103 // calcul dees deux surfaces
1104   TopLoc_Location L1,L2;
1105   TopoDS_Face aLocalFace = face1;
1106   const Handle(Geom_Surface)& s1 = BRep_Tool::Surface(aLocalFace,L1);
1107 //  const Handle(Geom_Surface)& s1 =
1108 //    BRep_Tool::Surface(TopoDS::Face(face1),L1);
1109   if (s1.IsNull()) return 1;
1110   aLocalFace = face2;
1111   const Handle(Geom_Surface)& s2 = BRep_Tool::Surface(aLocalFace,L2);
1112 //  const Handle(Geom_Surface)& s2 =
1113 //    BRep_Tool::Surface(TopoDS::Face(face2),L2);
1114   if (s2.IsNull()) return 1;
1115
1116
1117    Handle(Geom_Surface) surf1 = Handle(Geom_Surface)::
1118                                  DownCast(s1->Transformed(L1.Transformation()));
1119    if (surf1.IsNull()) return 1;
1120    Handle(Geom_Surface) surf2 = Handle(Geom_Surface)::
1121                                 DownCast(s2->Transformed(L2.Transformation()));
1122    if (surf2.IsNull()) return 1;
1123
1124
1125
1126   nbeval = (Standard_Integer ) Draw::Atof( a[3]);
1127
1128   switch(n)
1129   {
1130   case 7  : epsG1 = Draw::Atof(a[6]);
1131     Standard_FALLTHROUGH
1132   case 6  : epsC0 = Draw::Atof(a[5]);
1133     Standard_FALLTHROUGH
1134   case 5  : epsnl = Draw::Atof(a[4]);
1135     Standard_FALLTHROUGH
1136   case 4  : break;
1137   default : return 1;
1138   }
1139
1140
1141   Standard_Real pard1, parf1, U, Uf, deltaU, nb = 0;
1142   Standard_Boolean isconti = Standard_True;
1143   Standard_Boolean isdone = Standard_True;
1144   pard1 = curv1->FirstParameter();
1145   parf1 = curv1->LastParameter();
1146   Standard_Real  MaxG0Value=0, MaxG1Angle=0;
1147   U = Min( pard1,parf1);
1148   Uf = Max (pard1,parf1);
1149
1150   deltaU =  Abs(parf1- pard1)/nbeval;
1151
1152   do
1153      { if ( nb == nbeval)
1154           { LocalAnalysis_SurfaceContinuity res(curv1, curv2, Uf,surf1, surf2,            GeomAbs_G1, epsnl,epsC0, epsC1, epsC2, epsG1,percent,maxlen );
1155             isdone = res.IsDone();
1156            if ( isdone) { isconti = res.IsG1();
1157                              if (isconti)
1158                               {if (res.C0Value()>MaxG0Value)  MaxG0Value = res.C0Value();
1159                                if (res.G1Angle ()>MaxG1Angle) MaxG1Angle = res.G1Angle();}}
1160             else isconti = Standard_False;}
1161
1162           else {LocalAnalysis_SurfaceContinuity res (curv1, curv2, (U+nb*deltaU                 ), surf1,surf2, GeomAbs_G1,epsnl,epsC0, epsC1, epsC2, epsG1,
1163                  percent,maxlen);
1164                  isdone = res.IsDone();
1165                  if ( isdone) { isconti = res.IsG1();
1166                                   if ( nb == 0) { MaxG0Value = res.C0Value();
1167                                                   MaxG1Angle = res.G1Angle();}
1168                                   if (res.C0Value()> MaxG0Value)  MaxG0Value = res.C0Value();
1169                                   if (res.G1Angle()> MaxG1Angle)  MaxG1Angle= res.G1Angle();;}
1170
1171                  else isconti = Standard_False;}
1172        if (!isconti) ISG1=Standard_False;
1173        nb++;
1174
1175      }
1176    while ((nb<nbeval)&& isdone );
1177
1178   //if (!isdone)  { std::cout<<" Problem in computation "<<std::endl; return 1;}
1179   //if (ISG1)
1180   //    {std::cout<<" the continuity is G1 "<<std::endl;}
1181   //else { std::cout<<" the continuity is not G1  "<<std::endl;}
1182   //std::cout<<"MaxG0Value :"<< MaxG0Value << std::endl;
1183   //std::cout<<"MaxG1Angle:"<<  MaxG1Angle << std::endl;
1184   if (!isdone)  { di<<" Problem in computation \n"; return 1;}
1185   if (ISG1)
1186       {di<<" the continuity is G1 \n";}
1187   else { di<<" the continuity is not G1  \n";}
1188   di<<"MaxG0Value :"<< MaxG0Value << "\n";
1189   di<<"MaxG1Angle:"<<  MaxG1Angle << "\n";
1190   return 0;
1191 }
1192 /*****************************************************************************/
1193 static Standard_Integer shapeG0continuity (Draw_Interpretor& di, Standard_Integer n, const char** a)
1194
1195 { Standard_Real epsnl,epsC0, epsC1, epsC2, epsG1,percent,maxlen;
1196   Standard_Integer nbeval;
1197   Standard_Boolean ISG0;
1198   InitEpsSurf(epsnl,epsC0, epsC1, epsC2, epsG1,percent,maxlen  );
1199
1200   if (n<4) return 1;
1201   TopoDS_Face face1,face2;
1202   Standard_Real f1,f2,l1,l2;
1203   TopoDS_Shape shape = DBRep::Get(a[1],TopAbs_SHAPE);
1204   if (shape.IsNull()) return 1;
1205   TopoDS_Shape  edge  = DBRep::Get(a[2],TopAbs_EDGE);
1206   if (edge.IsNull()) return 1;
1207 // calcul des deux faces
1208   TopTools_IndexedDataMapOfShapeListOfShape lface;
1209   TopExp::MapShapesAndAncestors(shape,TopAbs_EDGE,TopAbs_FACE,lface);
1210   const TopTools_ListOfShape& lfac = lface.FindFromKey(edge);
1211
1212   Standard_Integer nelem= lfac.Extent();
1213   if(nelem!=2) return 1;
1214   TopTools_ListIteratorOfListOfShape It;
1215   It.Initialize(lfac);
1216   face1=TopoDS::Face(It.Value());
1217   It.Next();
1218   face2=TopoDS::Face(It.Value());
1219
1220   Standard_Boolean IsSeam = face1.IsEqual(face2);
1221
1222 // calcul des deux pcurves
1223   const Handle(Geom2d_Curve) c1 = BRep_Tool::CurveOnSurface
1224                                (TopoDS::Edge(edge),face1,f1,l1);
1225   if (c1.IsNull()) return 1;
1226
1227   if (IsSeam)
1228     edge.Reverse();
1229   const  Handle(Geom2d_Curve) c2  = BRep_Tool::CurveOnSurface
1230                                (TopoDS::Edge(edge),face2,f2,l2);
1231   if (c2.IsNull()) return 1;
1232
1233    Handle(Geom2d_Curve) curv1= new Geom2d_TrimmedCurve(c1,f1,l1);
1234
1235    Handle(Geom2d_Curve) curv2= new Geom2d_TrimmedCurve(c2,f2,l2);
1236
1237 // calcul des deux surfaces
1238   TopLoc_Location L1,L2;
1239   TopoDS_Face aLocalFace = face1;
1240   const Handle(Geom_Surface)& s1 = BRep_Tool::Surface(aLocalFace,L1);
1241 //  const Handle(Geom_Surface)& s1 =
1242 //    BRep_Tool::Surface(TopoDS::Face(face1),L1);
1243   if (s1.IsNull()) return 1;
1244   aLocalFace = face2;
1245   const Handle(Geom_Surface)& s2 = BRep_Tool::Surface(aLocalFace,L2);
1246 //  const Handle(Geom_Surface)& s2 =
1247 //    BRep_Tool::Surface(TopoDS::Face(face2),L2);
1248   if (s2.IsNull()) return 1;
1249
1250
1251    Handle(Geom_Surface) surf1 = Handle(Geom_Surface)::
1252                                  DownCast(s1->Transformed(L1.Transformation()));
1253    if (surf1.IsNull()) return 1;
1254    Handle(Geom_Surface) surf2 = Handle(Geom_Surface)::
1255                                 DownCast(s2->Transformed(L2.Transformation()));
1256    if (surf2.IsNull()) return 1;
1257
1258
1259
1260   nbeval = (Standard_Integer ) Draw::Atof( a[3]);
1261
1262   switch(n)
1263   {
1264   case 6  : epsC0   = Draw::Atof(a[5]);
1265     Standard_FALLTHROUGH
1266   case 5  : epsnl    = Draw::Atof(a[4]);
1267     Standard_FALLTHROUGH
1268   case 4  : break;
1269   default : return 1;
1270   }
1271
1272   Standard_Real pard1, parf1, U, Uf, deltaU, nb = 0;
1273   Standard_Boolean isconti = Standard_True;
1274   Standard_Boolean isdone = Standard_True;
1275   pard1 = curv1->FirstParameter();
1276   parf1 = curv1->LastParameter();
1277   Standard_Real  MaxG0Value=0;
1278   U = Min( pard1,parf1);
1279   Uf = Max (pard1,parf1);
1280
1281   deltaU =  Abs(parf1- pard1)/nbeval;
1282   ISG0=Standard_True;
1283   do
1284      { if ( nb == nbeval)
1285           { LocalAnalysis_SurfaceContinuity res (curv1, curv2, Uf,surf1, surf2, GeomAbs_C0,epsnl,epsC0, epsC1, epsC2, epsG1,percent,maxlen  );
1286             isdone = res.IsDone();
1287             if ( isdone) { isconti = res.IsC0( );
1288                           if (isconti)
1289                           if (res.C0Value()>MaxG0Value)  MaxG0Value= res.C0Value();}
1290             else isconti = Standard_False;}
1291
1292        else {LocalAnalysis_SurfaceContinuity res (curv1, curv2, (U+nb*deltaU), surf1, surf2, GeomAbs_C0,epsnl,epsC0, epsC1, epsC2, epsG1, percent,maxlen );
1293              isdone = res.IsDone();
1294              if ( isdone) { isconti = res.IsC0() ;
1295                                   if ( nb == 0) { MaxG0Value = res.C0Value();}
1296                                   if (res.C0Value()> MaxG0Value)  MaxG0Value = res.C0Value();}
1297
1298                  else isconti = Standard_False;}
1299
1300        nb++;
1301        if (!isconti) ISG0=Standard_False;
1302
1303      }
1304    while ((nb<nbeval)&& isdone );
1305
1306   //f (!isdone)  { std::cout<<" Problem in computation "<<std::endl; return 1;}
1307   //if (ISG0)
1308   //    {std::cout<<" the continuity is G0 "<<std::endl;}
1309
1310   //else { std::cout<<" the continuity is not G0  "<<std::endl;}
1311   //std::cout<<"MaxG0Value :"<< MaxG0Value << std::endl;
1312   if (!isdone)  { di<<" Problem in computation \n"; return 1;}
1313   if (ISG0)
1314       {di<<" the continuity is G0 \n";}
1315
1316   else { di<<" the continuity is not G0  \n";}
1317   di<<"MaxG0Value :"<< MaxG0Value << "\n";
1318   return 0;
1319 }
1320 /*****************************************************************************************/
1321 static Standard_Integer shapeG2continuity (Draw_Interpretor& di, Standard_Integer n, const char** a)
1322
1323 { Standard_Real epsnl,epsC0, epsC1, epsC2, epsG1, percent,maxlen;
1324   Standard_Boolean ISG2=Standard_True;
1325   Standard_Integer nbeval;
1326   Standard_Real MaxG0Value=0,MaxG1Angle=0,MaxG2Curvature=0;
1327   InitEpsSurf(epsnl,epsC0, epsC1, epsC2, epsG1,percent,maxlen );
1328
1329   if (n<4) return 1;
1330   TopoDS_Face face1,face2;
1331   Standard_Real f1,f2,l1,l2;
1332   TopoDS_Shape shape = DBRep::Get(a[1],TopAbs_SHAPE);
1333   if (shape.IsNull()) return 1;
1334   TopoDS_Shape  edge  = DBRep::Get(a[2],TopAbs_EDGE);
1335   if (edge.IsNull()) return 1;
1336 // calcul des deux faces
1337   TopTools_IndexedDataMapOfShapeListOfShape lface;
1338   TopExp::MapShapesAndAncestors(shape,TopAbs_EDGE,TopAbs_FACE,lface);
1339   const TopTools_ListOfShape& lfac = lface.FindFromKey(edge);
1340
1341   Standard_Integer nelem= lfac.Extent();
1342   if(nelem!=2) return 1;
1343   TopTools_ListIteratorOfListOfShape It;
1344   It.Initialize(lfac);
1345   face1=TopoDS::Face(It.Value());
1346   It.Next();
1347   face2=TopoDS::Face(It.Value());
1348
1349   Standard_Boolean IsSeam = face1.IsEqual(face2);
1350
1351 // calcul des deux pcurves
1352   const Handle(Geom2d_Curve) c1 = BRep_Tool::CurveOnSurface
1353                                (TopoDS::Edge(edge),face1,f1,l1);
1354   if (c1.IsNull()) return 1;
1355
1356   if (IsSeam)
1357     edge.Reverse();
1358   const  Handle(Geom2d_Curve) c2  = BRep_Tool::CurveOnSurface
1359                                (TopoDS::Edge(edge),face2,f2,l2);
1360   if (c2.IsNull()) return 1;
1361
1362    Handle(Geom2d_Curve) curv1= new Geom2d_TrimmedCurve(c1,f1,l1);
1363
1364    Handle(Geom2d_Curve) curv2= new Geom2d_TrimmedCurve(c2,f2,l2);
1365
1366 // calcul des deux surfaces
1367   TopLoc_Location L1,L2;
1368   TopoDS_Face aLocalFace = face1;
1369   const Handle(Geom_Surface)& s1 = BRep_Tool::Surface(aLocalFace,L1);
1370 //  const Handle(Geom_Surface)& s1 =
1371 //    BRep_Tool::Surface(TopoDS::Face(face1),L1);
1372   if (s1.IsNull()) return 1;
1373   aLocalFace = face2;
1374   const Handle(Geom_Surface)& s2 = BRep_Tool::Surface(aLocalFace,L2);
1375 //  const Handle(Geom_Surface)& s2 =
1376 //    BRep_Tool::Surface(TopoDS::Face(face2),L2);
1377   if (s2.IsNull()) return 1;
1378
1379
1380    Handle(Geom_Surface) surf1 = Handle(Geom_Surface)::
1381                                  DownCast(s1->Transformed(L1.Transformation()));
1382    if (surf1.IsNull()) return 1;
1383    Handle(Geom_Surface) surf2 = Handle(Geom_Surface)::
1384                                 DownCast(s2->Transformed(L2.Transformation()));
1385    if (surf2.IsNull()) return 1;
1386
1387
1388
1389   nbeval = (Standard_Integer ) Draw::Atof( a[3]);
1390
1391   switch(n)
1392   {
1393   case 9  :  maxlen = Draw::Atof(a[8]);
1394     Standard_FALLTHROUGH
1395   case 8   : percent = Draw::Atof(a[7]);
1396     Standard_FALLTHROUGH
1397   case 7   : epsG1 = Draw::Atof(a[6]);
1398     Standard_FALLTHROUGH
1399   case 6  :  epsC0 = Draw::Atof(a[5]);
1400     Standard_FALLTHROUGH
1401   case 5  :  epsnl = Draw::Atof(a[4]);
1402     Standard_FALLTHROUGH
1403   case 4  : break;
1404   default : return 1;
1405   }
1406
1407
1408   Standard_Real pard1, parf1, U, Uf, deltaU, nb = 0;
1409   Standard_Boolean isconti = Standard_True;
1410   Standard_Boolean isdone = Standard_True;
1411   pard1 = curv1->FirstParameter();
1412   parf1 = curv1->LastParameter();
1413   U = Min( pard1,parf1);
1414   Uf = Max (pard1,parf1);
1415
1416   deltaU =  Abs(parf1- pard1)/nbeval;
1417
1418   do
1419      { if ( nb == nbeval)
1420           { LocalAnalysis_SurfaceContinuity res (curv1, curv2, Uf,surf1, surf2,             GeomAbs_G2,epsnl,epsC0, epsC1, epsC2, epsG1,percent,maxlen );
1421             isdone = res.IsDone();
1422             if (isdone){isconti = res.IsG2();
1423                         if (isconti)
1424                         {if(res.C0Value()>MaxG0Value) MaxG0Value=res.C0Value();
1425                          if(res.G1Angle()>MaxG1Angle) MaxG1Angle=res.G1Angle();
1426                          if(res.G2CurvatureGap()>MaxG2Curvature)
1427                            MaxG2Curvature=res.G2CurvatureGap();
1428                         }}
1429             else isconti = Standard_False;}
1430
1431        else { LocalAnalysis_SurfaceContinuity res (curv1, curv2, (U+nb*deltaU), surf1,surf2, GeomAbs_G2,epsnl,epsC0, epsC1, epsC2, epsG1,percent,maxlen );
1432              isdone = res.IsDone();
1433              if ( isdone) {isconti = res.IsG2();
1434                            if (nb==0){MaxG0Value=res.C0Value();
1435                                       MaxG1Angle=res.G1Angle();
1436                                       MaxG2Curvature=res.G2CurvatureGap();}
1437                            if(res.C0Value()>MaxG0Value)
1438                                 MaxG0Value=res.C0Value();
1439                            if(res.G1Angle()>MaxG1Angle)
1440                                 MaxG1Angle=res.G1Angle();
1441                            if(res.G2CurvatureGap()>MaxG2Curvature)
1442                               MaxG2Curvature=res.G2CurvatureGap();}
1443                  else isconti = Standard_False;}
1444
1445        nb++;
1446        if (!isconti) ISG2=Standard_False;
1447
1448      }
1449    while ((nb<nbeval)&& isdone );
1450
1451   //if (!isdone)  { std::cout<<" Problem in computation "<<std::endl; return 1;}
1452   //if (ISG2)
1453   //std::cout<<" the continuity is G2 "<<std::endl;
1454   //else std::cout<<" the continuity is not G2  "<<std::endl;
1455   //std::cout<<"MaxG0Value :"<< MaxG0Value << std::endl;
1456   //std::cout<<"MaxG1Angle:"<<  MaxG1Angle << std::endl;
1457   //std::cout<<"MaxG2Curvature:"<<MaxG2Curvature<<std::endl;
1458   if (!isdone)  { di<<" Problem in computation \n"; return 1;}
1459   if (ISG2)
1460   di<<" the continuity is G2 \n";
1461   else di<<" the continuity is not G2  \n";
1462   di<<"MaxG0Value :"<< MaxG0Value << "\n";
1463   di<<"MaxG1Angle:"<<  MaxG1Angle << "\n";
1464   di<<"MaxG2Curvature:"<<MaxG2Curvature<<"\n";
1465   return 0;
1466 }
1467
1468
1469 //=======================================================================
1470 //function : clintedge
1471 //purpose  :
1472 //=======================================================================
1473 static Standard_Integer clintedge(Draw_Interpretor& di,
1474                                          Standard_Integer narg, const char** a)
1475 {
1476   char newname[255];
1477
1478   if (narg < 2) {
1479     //std::cout << "Usage: clintedge shape" << std::endl;
1480     di << "Usage: clintedge shape\n";
1481     return 1;
1482   }
1483   TopoDS_Shape S = DBRep::Get(a[1]);
1484
1485   TopTools_DataMapOfShapeListOfShape mymap;
1486   TopOpeBRepTool_PurgeInternalEdges mypurgealgo(S);
1487   Standard_Integer nbedges = mypurgealgo.NbEdges();
1488   if (nbedges > 0)
1489   {
1490     //std::cout<<nbedges<<" internal (or external) edges to be removed"<<std::endl;
1491     di<<nbedges<<" internal (or external) edges to be removed\n";
1492
1493     Standard_Integer i = 1;
1494     char* temp = newname;
1495
1496     Sprintf(newname,"%s_%d",a[1],i);
1497     DBRep::Set(temp,mypurgealgo.Shape());
1498     //std::cout<<newname<<" ";
1499     di<<newname<<" ";
1500
1501     //std::cout<<std::endl;
1502     di<<"\n";
1503   }
1504   else
1505     di << "no internal (or external) edges\n";
1506     //std::cout << "no internal (or external) edges"<<std::endl;
1507
1508   return 0;
1509 }
1510
1511
1512 //=======================================================================
1513 //function : facintedge
1514 //purpose  :
1515 //=======================================================================
1516 static Standard_Integer facintedge(Draw_Interpretor& di,
1517                                          Standard_Integer narg, const char** a)
1518 {
1519   char newname[255];
1520
1521   if (narg < 2) {
1522     //std::cout << "Usage: facintedge shape" << std::endl;
1523     di << "Usage: facintedge shape\n";
1524     return 1;
1525   }
1526   TopoDS_Shape S = DBRep::Get(a[1]);
1527
1528   TopTools_DataMapOfShapeListOfShape mymap;
1529   TopOpeBRepTool_PurgeInternalEdges mypurgealgo(S);
1530   mypurgealgo.Faces(mymap);
1531
1532   Standard_Integer i = 1;
1533   char* temp = newname;
1534
1535   TopTools_DataMapIteratorOfDataMapOfShapeListOfShape itFacEdg;
1536   for (itFacEdg.Initialize(mymap); itFacEdg.More(); itFacEdg.Next()) {
1537       Sprintf(newname,"%s_%d",a[1],i);
1538       DBRep::Set(temp,itFacEdg.Key());
1539       //std::cout<<newname<<" ";
1540       di<<newname<<" ";
1541       i++;
1542     }
1543
1544   //std::cout<<std::endl;
1545   di<<"\n";
1546
1547   return 0;
1548 }
1549
1550 //=======================================================================
1551 //function : fuseedge
1552 //purpose  :
1553 //=======================================================================
1554 static Standard_Integer fuseedge(Draw_Interpretor& di,
1555                                          Standard_Integer narg, const char** a)
1556 {
1557   char newname[255];
1558
1559   if (narg < 2) {
1560     //std::cout << "Usage: fuseedge shape" << std::endl;
1561     di << "Usage: fuseedge shape\n";
1562     return 1;
1563   }
1564   TopoDS_Shape S = DBRep::Get(a[1]);
1565
1566   TopTools_DataMapOfIntegerListOfShape mymap;
1567   //TopOpeBRepTool_FuseEdges myfusealgo(S);
1568   BRepLib_FuseEdges myfusealgo(S);
1569   myfusealgo.SetConcatBSpl();
1570   Standard_Integer nbvertices;
1571   nbvertices = myfusealgo.NbVertices();
1572
1573   if (nbvertices > 0) {
1574
1575     //std::cout<<nbvertices<<" vertices to be removed"<<std::endl;
1576     di<<nbvertices<<" vertices to be removed\n";
1577
1578     Standard_Integer i = 1;
1579     char* temp = newname;
1580
1581     Sprintf(newname,"%s_%d",a[1],i);
1582     DBRep::Set(temp,myfusealgo.Shape());
1583     //std::cout<<newname<<" ";
1584     di<<newname<<" ";
1585
1586     //std::cout<<std::endl;
1587     di<<"\n";
1588   }
1589   else
1590     di << "no vertices to remove\n";
1591     //std::cout << "no vertices to remove"<<std::endl;
1592
1593   return 0;
1594 }
1595
1596
1597 //=======================================================================
1598 //function : listfuseedge
1599 //purpose  :
1600 //=======================================================================
1601 static Standard_Integer listfuseedge(Draw_Interpretor& di,
1602                                          Standard_Integer narg, const char** a)
1603 {
1604   char newname[255];
1605
1606   if (narg < 2) {
1607     //std::cout << "Usage: listfuseedge shape" << std::endl;
1608     di << "Usage: listfuseedge shape\n";
1609     return 1;
1610   }
1611   TopoDS_Shape S = DBRep::Get(a[1]);
1612
1613   TopTools_DataMapOfIntegerListOfShape mymap;
1614   BRepLib_FuseEdges myfusealgo(S);
1615   myfusealgo.Edges(mymap);
1616
1617   Standard_Integer i;
1618   char* temp = newname;
1619
1620   TopTools_DataMapIteratorOfDataMapOfIntegerListOfShape itLstEdg;
1621   for (itLstEdg.Initialize(mymap); itLstEdg.More(); itLstEdg.Next()) {
1622       const Standard_Integer& iLst = itLstEdg.Key();
1623       const TopTools_ListOfShape& LmapEdg = mymap.Find(iLst);
1624       TopTools_ListIteratorOfListOfShape itEdg;
1625       i = 1;
1626       for (itEdg.Initialize(LmapEdg); itEdg.More(); itEdg.Next()) {
1627         Sprintf(newname,"%s_%d_%d",a[1],iLst,i);
1628         DBRep::Set(temp,itEdg.Value());
1629         //std::cout<<newname<<" ";
1630         di<<newname<<" ";
1631         i++;
1632       }
1633
1634     }
1635
1636   //std::cout<<std::endl;
1637   di<<"\n";
1638
1639   return 0;
1640 }
1641
1642 //=======================================================================
1643 //function : tolsphere
1644 //purpose  :
1645 //=======================================================================
1646 static Standard_Integer tolsphere(Draw_Interpretor& di, Standard_Integer n, const char** a)
1647 {
1648   if (n != 2)
1649   {
1650     di << "use toolsphere shape\n";
1651     return 1;
1652   }
1653
1654   TopoDS_Shape aS = DBRep::Get(a[1]);
1655   if (aS.IsNull())
1656   {
1657     di << "No such shape " << a[1] << "\n";
1658     return 1;
1659   }
1660
1661   TopTools_IndexedMapOfShape aMapV;
1662   TopExp::MapShapes(aS, TopAbs_VERTEX, aMapV);
1663   for (Standard_Integer i = 1; i <= aMapV.Extent(); i++)
1664   {
1665     const TopoDS_Vertex& aV = TopoDS::Vertex(aMapV.FindKey(i));
1666     Standard_Real aRadius = BRep_Tool::Tolerance(aV);
1667     gp_Pnt aCenter = BRep_Tool::Pnt(aV);
1668     Handle(Geom_Surface) aSph = new Geom_SphericalSurface(gp_Ax2(aCenter,gp::DZ()), aRadius);
1669     TCollection_AsciiString aName(a[1]);
1670     aName = aName + "_v" + i;
1671     DrawTrSurf::Set(aName.ToCString(), aSph);
1672     di << aName << " ";
1673   }
1674   return 0;
1675 }
1676
1677 //=======================================================================
1678 //function : validrange
1679 //purpose  :
1680 //=======================================================================
1681 static Standard_Integer validrange(Draw_Interpretor& di,
1682   Standard_Integer narg, const char** a)
1683 {
1684   if (narg < 2)
1685   {
1686     di << "usage: validrange edge [(out) u1 u2]";
1687     return 1;
1688   }
1689
1690   TopoDS_Edge aE = TopoDS::Edge(DBRep::Get(a[1],TopAbs_EDGE, true));
1691   if (aE.IsNull())
1692     return 1;
1693
1694   Standard_Real u1, u2;
1695   if (BRepLib::FindValidRange(aE, u1, u2))
1696   {
1697     if (narg > 3)
1698     {
1699       Draw::Set(a[2], u1);
1700       Draw::Set(a[3], u2);
1701     }
1702     else
1703     {
1704       di << u1 << " " << u2;
1705     }
1706   }
1707   else
1708     di << "edge has no valid range";
1709   return 0;
1710 }
1711
1712 //=======================================================================
1713 //function : CheckCommands
1714 //purpose  :
1715 //=======================================================================
1716
1717 void BRepTest::CheckCommands(Draw_Interpretor& theCommands)
1718 {
1719   static Standard_Boolean done = Standard_False;
1720   if (done) return;
1721   done = Standard_True;
1722
1723   BRepTest_CheckCommands_SetFaultyName("faulty_");
1724   DBRep::BasicCommands(theCommands);
1725
1726   const char* g = "TOPOLOGY Check commands";
1727
1728 //  Modified by skv - Tue Apr 27 13:35:35 2004 Begin
1729   theCommands.Add("checkshape",
1730                   "checkshape : no args to have help",
1731                   __FILE__,
1732                   checkshape,
1733                   g);
1734 //   theCommands.Add("checkshape",
1735 //                "checks the validity of a shape : checkshape name,\n                      short description of check : checkshape name -short",
1736 //                __FILE__,
1737 //                CHK,
1738 //                g);
1739 //   theCommands.Add("checktopshape",
1740 //                "checks the topological validity of a shape : checktopshape name",
1741 //                __FILE__,
1742 //                CHK,
1743 //                g);
1744 //  Modified by skv - Tue Apr 27 13:35:39 2004 End
1745
1746   theCommands.Add("checksection",
1747                   "checks the closure of a section : checksection name [-r <RefVal>]\n"
1748                   "\"-r\" - allowed number of allone vertices.",
1749                   __FILE__,
1750                   checksection,
1751                   g);
1752
1753   theCommands.Add("checkdiff",
1754                   "checks the validity of the diff beetween the shapes arg1..argn and result :\n checkdiff arg1 [arg2..argn] result [closedSolid (1/0)] [geomCtrl (1/0)]",
1755                   __FILE__,
1756                   checkdiff,
1757                   g);
1758
1759 g = "TOPOLOGY Analysis of shapes ";
1760
1761 theCommands.Add("shapeG0continuity",
1762                  "shapeG0continuity  shape  edge nbeval [epsnul [epsG0]]",
1763                   __FILE__,
1764                   shapeG0continuity, g);
1765
1766 theCommands.Add("shapeG1continuity",
1767                 "shapeG1continuity  shape  edge nbeval [epsnul [epsG0 [epsG1]]]",
1768                    __FILE__,
1769                 shapeG1continuity ,g);
1770 theCommands.Add("shapeG2continuity",
1771                   "shapeG2continuity shape  edge  nbeval [epsnul [epsG0 [epsG1 [maxlen [perce]]]]]",
1772                   __FILE__,
1773                   shapeG2continuity,g);
1774
1775 theCommands.Add("computetolerance",
1776                   "computetolerance shape",
1777                   __FILE__,
1778                   computetolerance,g);
1779
1780 theCommands.Add("clintedge",
1781                   "clintedge shape",
1782                   __FILE__,
1783                   clintedge,g);
1784
1785 theCommands.Add("facintedge",
1786                   "facintedge shape",
1787                   __FILE__,
1788                   facintedge,g);
1789
1790 theCommands.Add("fuseedge",
1791                   "fuseedge shape",
1792                   __FILE__,
1793                   fuseedge,g);
1794
1795 theCommands.Add("listfuseedge",
1796                   "listfuseedge shape",
1797                   __FILE__,
1798                   listfuseedge,g);
1799 theCommands.Add("tolsphere", "toolsphere shape\n"
1800                 "\t\tshows vertex tolerances by drawing spheres",
1801                 __FILE__, tolsphere, g);
1802 theCommands.Add("validrange",
1803                 "validrange edge [(out) u1 u2]\n"
1804                 "\t\tcomputes valid range of the edge, and\n"
1805                 "\t\tprints first and last values or sets the variables u1 and u2",
1806                 __FILE__, validrange, g);
1807 }
1808