src/ShapeFix/ShapeFix_Solid.cxx

   1 // Copyright (c) 1999-2014 OPEN CASCADE SAS
   2 //
   3 // This file is part of Open CASCADE Technology software library.
   4 //
   5 // This library is free software; you can redistribute it and/or modify it under
   6 // the terms of the GNU Lesser General Public License version 2.1 as published
   7 // by the Free Software Foundation, with special exception defined in the file
   8 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
   9 // distribution for complete text of the license and disclaimer of any warranty.
  10 //
  11 // Alternatively, this file may be used under the terms of Open CASCADE
  12 // commercial license or contractual agreement.
  13
  14
  15 #include <Bnd_Box2d.hxx>
  16 #include <BRep_Builder.hxx>
  17 #include <BRep_Tool.hxx>
  18 #include <BRepClass3d_SolidClassifier.hxx>
  19 #include <Geom_Curve.hxx>
  20 #include <Geom_Surface.hxx>
  21 #include <gp_Pnt.hxx>
  22 #include <Message_Msg.hxx>
  23 #include <Message_ProgressIndicator.hxx>
  24 #include <Message_ProgressSentry.hxx>
  25 #include <Precision.hxx>
  26 #include <ShapeAnalysis.hxx>
  27 #include <ShapeAnalysis_Curve.hxx>
  28 #include <ShapeAnalysis_Edge.hxx>
  29 #include <ShapeAnalysis_FreeBounds.hxx>
  30 #include <ShapeBuild_ReShape.hxx>
  31 #include <ShapeExtend.hxx>
  32 #include <ShapeExtend_BasicMsgRegistrator.hxx>
  33 #include <ShapeExtend_WireData.hxx>
  34 #include <ShapeFix_Shell.hxx>
  35 #include <ShapeFix_Solid.hxx>
  36 #include <Standard_ErrorHandler.hxx>
  37 #include <Standard_Failure.hxx>
  38 #include <Standard_Type.hxx>
  39 #include <TopAbs.hxx>
  40 #include <TopExp.hxx>
  41 #include <TopoDS.hxx>
  42 #include <TopoDS_CompSolid.hxx>
  43 #include <TopoDS_Edge.hxx>
  44 #include <TopoDS_Iterator.hxx>
  45 #include <TopoDS_Shape.hxx>
  46 #include <TopoDS_Shell.hxx>
  47 #include <TopoDS_Solid.hxx>
  48 #include <TopoDS_Vertex.hxx>
  49 #include <TopoDS_Wire.hxx>
  50 #include <TopTools_DataMapIteratorOfDataMapOfShapeListOfShape.hxx>
  51 #include <TopTools_DataMapIteratorOfDataMapOfShapeShape.hxx>
  52 #include <TopTools_DataMapOfShapeInteger.hxx>
  53 #include <TopTools_DataMapOfShapeListOfShape.hxx>
  54 #include <TopTools_IndexedDataMapOfShapeListOfShape.hxx>
  55 #include <TopTools_IndexedDataMapOfShapeShape.hxx>
  56 #include <TopTools_IndexedMapOfShape.hxx>
  57 #include <TopTools_ListIteratorOfListOfShape.hxx>
  58 #include <TopTools_ListOfShape.hxx>
  59 #include <TopTools_MapIteratorOfMapOfShape.hxx>
  60 #include <TopTools_MapOfShape.hxx>
  61 #include <TopTools_SequenceOfShape.hxx>
  62
  63 IMPLEMENT_STANDARD_RTTIEXT(ShapeFix_Solid,ShapeFix_Root)
  64
  65 //======================================================
  66 //function : ShapeFix_Solid
  67 //purpose  :
  68 //=======================================================================
  69 ShapeFix_Solid::ShapeFix_Solid()
  70 {
  71   myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK);
  72   myFixShellMode = -1;
  73   myFixShell = new ShapeFix_Shell;
  74   myCreateOpenSolidMode = Standard_False;
  75 }
  76
  77 //=======================================================================
  78 //function : ShapeFix_Solid
  79 //purpose  :
  80 //=======================================================================
  81
  82 ShapeFix_Solid::ShapeFix_Solid(const TopoDS_Solid& solid)
  83 {
  84   myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK);
  85   myFixShellMode = -1;
  86   myFixShell = new ShapeFix_Shell;
  87   myCreateOpenSolidMode = Standard_False;
  88   Init(solid);
  89 }
  90
  91 //=======================================================================
  92 //function : Init
  93 //purpose  :
  94 //=======================================================================
  95
  96  void ShapeFix_Solid::Init(const TopoDS_Solid& solid)
  97 {
  98   mySolid = solid;
  99   //mySolid = TopoDS::Solid(shape.EmptyCopied());
 100   //BRep_Builder B;
 101  // for( TopoDS_Iterator iter(solid); iter.More(); iter.Next())
 102  //  B.Add(mySolid,TopoDS::Shell(iter.Value()));
 103   myShape = solid;
 104 }
 105 #ifdef OCCT_DEBUG_GET_MIDDLE_POINT
 106 //=======================================================================
 107 //function : GetMiddlePoint
 108 //purpose  :
 109 //=======================================================================
 110 static void GetMiddlePoint(const TopoDS_Shape& aShape, gp_Pnt& pmid)
 111 {
 112   TopExp_Explorer aExp(aShape,TopAbs_EDGE);
 113   gp_XYZ center(0.0,0.0,0.0);
 114   Standard_Integer numpoints =0;
 115   for( ; aExp.More(); aExp.Next()) {
 116     TopoDS_Edge e1 = TopoDS::Edge(aExp.Current());
 117     Standard_Real f,l;
 118     Handle(Geom_Curve) c3d = BRep_Tool::Curve(e1,f,l);
 119     if(!c3d.IsNull()) {
 120       for(Standard_Integer i =1 ; i <=5; i++) {
 121         Standard_Real param = f+(l-f)/4*(i-1);
 122         gp_Pnt pt;
 123         numpoints++;
 124         c3d->D0(param,pt);
 125         center+=pt.XYZ();
 126
 127       }
 128     }
 129   }
 130   center /= numpoints;
 131   pmid.SetXYZ(center);
 132 }
 133 #endif
 134 //=======================================================================
 135 //function : CollectSolids
 136 //purpose  :
 137 //=======================================================================
 138 static void CollectSolids(const TopTools_SequenceOfShape& aSeqShells ,
 139                           TopTools_IndexedDataMapOfShapeListOfShape& anIndexedMapShellHoles,
 140                           TopTools_DataMapOfShapeInteger& theMapStatus)
 141 {
 142   TopTools_MapOfShape aMapHoles;
 143   TopTools_DataMapOfShapeListOfShape aMapShellHoles;
 144   for ( Standard_Integer i1 = 1; i1 <= aSeqShells.Length(); i1++ ) {
 145     TopoDS_Shell aShell1 = TopoDS::Shell(aSeqShells.Value(i1));
 146     TopTools_ListOfShape lshells;
 147     aMapShellHoles.Bind(aShell1,lshells);
 148   }
 149   //Finds roots shells and hole shells.
 150   for ( Standard_Integer i = 1; i <= aSeqShells.Length(); i++ ) {
 151     TopoDS_Shell aShell1 =  TopoDS::Shell(aSeqShells.Value(i));
 152     TopExp_Explorer aExpEdges(aShell1,TopAbs_EDGE);
 153     if(!BRep_Tool::IsClosed(aShell1) || !aExpEdges.More()) continue;
 154     TopoDS_Solid solid;
 155     BRep_Builder B;
 156     B.MakeSolid (solid);
 157     B.Add (solid,aShell1);
 158     try {
 159       OCC_CATCH_SIGNALS
 160      TopAbs_State infinstatus = TopAbs_UNKNOWN;
 161      BRepClass3d_SolidClassifier bsc3d (solid);
 162      Standard_Integer st = 0;
 163      if(!theMapStatus.IsBound(aShell1)) {
 164
 165          bsc3d.PerformInfinitePoint(Precision::Confusion());
 166          infinstatus = bsc3d.State();
 167
 168          if(infinstatus != TopAbs_UNKNOWN &&  infinstatus !=TopAbs_ON)
 169            st = (infinstatus == TopAbs_IN ? 1 :2);
 170          theMapStatus.Bind(aShell1,st);
 171
 172
 173        }
 174       else {
 175         st = theMapStatus.Find(aShell1);
 176         if(st)
 177           infinstatus =  (theMapStatus.Find(aShell1) == 1 ? TopAbs_IN : TopAbs_OUT);
 178       }
 179       if(!st) continue;
 180       for ( Standard_Integer j = 1; j <= aSeqShells.Length(); j++ ) {
 181         if(i==j) continue;
 182         TopoDS_Shape aShell2 = aSeqShells.Value(j);
 183         if(!BRep_Tool::IsClosed(aShell2)) continue;
 184         if(aMapHoles.Contains(aShell2)) continue;
 185         if(aMapShellHoles.IsBound(aShell2)) {
 186           Standard_Boolean isAnalysis = Standard_False;
 187           const TopTools_ListOfShape& ls = aMapShellHoles.Find(aShell2);
 188           for(TopTools_ListIteratorOfListOfShape li(ls); li.More() && !isAnalysis; li.Next())
 189             isAnalysis = li.Value().IsSame(aShell1);
 190           if(isAnalysis) continue;
 191         }
 192         TopAbs_State pointstatus = TopAbs_UNKNOWN;
 193         Standard_Integer numon =0;
 194         TopTools_IndexedMapOfShape amapVert;
 195         for(TopExp_Explorer aExpVert(aShell2,TopAbs_VERTEX); aExpVert.More() && amapVert.Extent() < 10; aExpVert.Next())
 196           amapVert.Add(aExpVert.Current());
 197         for(Standard_Integer k = 1; k <= amapVert.Extent() &&
 198             (pointstatus ==TopAbs_UNKNOWN || (pointstatus ==TopAbs_ON && numon < 3)); k++) {
 199           TopoDS_Vertex aV = TopoDS::Vertex(amapVert.FindKey(k));
 200           gp_Pnt aPf =  BRep_Tool::Pnt(aV);
 201           bsc3d.Perform(aPf,Precision::Confusion());
 202           pointstatus =bsc3d.State();
 203           if(pointstatus ==TopAbs_ON) numon++;
 204         }
 205
 206         if(numon == 3 && pointstatus ==TopAbs_ON) {
 207 #ifdef OCCT_DEBUG_GET_MIDDLE_POINT
 208           gp_Pnt pmid;
 209           GetMiddlePoint(aShell2,pmid);
 210           bsc3d.Perform(pmid,Precision::Confusion());
 211 #endif
 212           pointstatus = /*(bsc3d.State() == TopAbs_IN ? TopAbs_IN :*/TopAbs_OUT;
 213         }
 214         if(pointstatus != infinstatus) {
 215           aMapShellHoles.ChangeFind(aShell1).Append(aShell2);
 216           if( aMapHoles.Contains(aShell2))
 217             aMapHoles.Remove(aShell2);
 218           else aMapHoles.Add(aShell2);
 219         }
 220       }
 221     }
 222     catch(Standard_Failure const& anException) {
 223 #ifdef OCCT_DEBUG
 224       cout << "Warning: ShapeFix_Solid::SolidFromShell: Exception: ";
 225       anException.Print(cout); cout << endl;
 226 #endif
 227       (void)anException;
 228       continue;
 229     }
 230   }
 231   TopTools_DataMapIteratorOfDataMapOfShapeListOfShape aItShellHoles( aMapShellHoles);
 232   for(; aItShellHoles.More();aItShellHoles.Next()) {
 233     if(aMapHoles.Contains(aItShellHoles.Key())) continue;
 234     const TopTools_ListOfShape& lHoles =aItShellHoles.Value();
 235     if(lHoles.IsEmpty()) continue;
 236     for(TopTools_ListIteratorOfListOfShape lItHoles(lHoles);lItHoles.More(); lItHoles.Next()) {
 237       if(aMapHoles.Contains(lItHoles.Value())) {
 238         const TopTools_ListOfShape& lUnHoles = aMapShellHoles.Find(lItHoles.Value());
 239         for(TopTools_ListIteratorOfListOfShape lItUnHoles(lUnHoles);lItUnHoles.More(); lItUnHoles.Next())
 240           aMapHoles.Remove(lItUnHoles.Value());
 241       }
 242     }
 243   }
 244   for(TopTools_MapIteratorOfMapOfShape aIterHoles(aMapHoles);aIterHoles.More(); aIterHoles.Next())
 245     aMapShellHoles.UnBind (aIterHoles.Key());
 246
 247   for (Standard_Integer i = 1; i <= aSeqShells.Length(); ++i) {
 248     const TopoDS_Shape& aShell1 = aSeqShells.Value (i);
 249     if (aMapShellHoles.IsBound (aShell1)) {
 250       const TopTools_ListOfShape& ls = aMapShellHoles.Find (aShell1);
 251       anIndexedMapShellHoles.Add (aShell1, ls);
 252     }
 253   }
 254 }
 255 //=======================================================================
 256 //function : CreateSolids
 257 //purpose  :
 258 //=======================================================================
 259
 260 static Standard_Boolean CreateSolids(const TopoDS_Shape theShape,TopTools_IndexedMapOfShape& aMapSolids)
 261 {
 262   TopTools_SequenceOfShape aSeqShells;
 263   Standard_Boolean isDone = Standard_False;
 264
 265   for(TopExp_Explorer aExpShell(theShape,TopAbs_SHELL); aExpShell.More(); aExpShell.Next()) {
 266     aSeqShells.Append(aExpShell.Current());
 267   }
 268   TopTools_IndexedDataMapOfShapeListOfShape aMapShellHoles;
 269   TopTools_DataMapOfShapeInteger aMapStatus;
 270   CollectSolids(aSeqShells,aMapShellHoles,aMapStatus);
 271   TopTools_IndexedDataMapOfShapeShape ShellSolid;
 272   //Defines correct orientation of shells
 273   for (Standard_Integer i = 1; i <= aMapShellHoles.Extent(); ++i) {
 274     TopoDS_Shell aShell = TopoDS::Shell(aMapShellHoles.FindKey(i));
 275     TopExp_Explorer aExpEdges(aShell,TopAbs_EDGE);
 276     if(!BRep_Tool::IsClosed(aShell) || !aExpEdges.More()) {
 277       ShellSolid.Add(aShell,aShell);
 278       isDone = Standard_True;
 279       continue;
 280     }
 281     BRep_Builder B;
 282     TopAbs_State infinstatus = TopAbs_UNKNOWN;
 283     TopoDS_Solid aSolid;
 284     B.MakeSolid (aSolid);
 285     B.Add (aSolid,aShell);
 286     if(aMapStatus.IsBound(aShell)) {
 287       Standard_Integer st = aMapStatus.Find(aShell);
 288       if(st)
 289         infinstatus =  (aMapStatus.Find(aShell) == 1 ? TopAbs_IN : TopAbs_OUT);
 290     }
 291
 292     else {
 293
 294     try {
 295       OCC_CATCH_SIGNALS
 296       BRepClass3d_SolidClassifier bsc3d (aSolid);
 297       bsc3d.PerformInfinitePoint(Precision::Confusion());
 298       infinstatus = bsc3d.State();
 299       }
 300     catch(Standard_Failure const& anException) {
 301 #ifdef OCCT_DEBUG
 302       cout << "Warning: ShapeFix_Solid::SolidFromShell: Exception: ";
 303       anException.Print(cout); cout << endl;
 304 #endif
 305       (void)anException;
 306       ShellSolid.Add(aShell,aSolid);
 307       continue;
 308     }
 309   }
 310   if (infinstatus == TopAbs_IN) {
 311     isDone = Standard_True;
 312     aShell.Reverse();
 313     TopoDS_Solid aTmpSolid;
 314     B.MakeSolid (aTmpSolid);
 315     B.Add (aTmpSolid,aShell);
 316     aSolid = aTmpSolid;
 317   }
 318
 319     const TopTools_ListOfShape& lHoles = aMapShellHoles (i);
 320     for(TopTools_ListIteratorOfListOfShape lItHoles(lHoles); lItHoles.More();lItHoles.Next()) {
 321       TopoDS_Shell aShellHole = TopoDS::Shell(lItHoles.Value());
 322       if(aMapStatus.IsBound(aShellHole)) {
 323         infinstatus = (aMapStatus.Find(aShellHole) == 1 ? TopAbs_IN : TopAbs_OUT);
 324       }
 325       else {
 326         TopoDS_Solid solid;
 327         B.MakeSolid (solid);
 328         B.Add (solid,aShellHole);
 329         BRepClass3d_SolidClassifier bsc3dHol (solid);
 330         bsc3dHol.PerformInfinitePoint(Precision::Confusion());
 331         infinstatus = bsc3dHol.State();
 332       }
 333       if (infinstatus == TopAbs_OUT) {
 334         aShellHole.Reverse();
 335         isDone = Standard_True;
 336       }
 337       B.Add(aSolid,aShellHole);
 338     }
 339     ShellSolid.Add(aShell,aSolid);
 340   }
 341   //Creation of compsolid from shells containing shared faces.
 342   TopTools_IndexedDataMapOfShapeListOfShape aMapFaceShells;
 343   TopExp::MapShapesAndAncestors(theShape,TopAbs_FACE,TopAbs_SHELL,aMapFaceShells);
 344   for(Standard_Integer i =1; i <= aMapFaceShells.Extent(); i++) {
 345     const TopTools_ListOfShape& lshells = aMapFaceShells.FindFromIndex(i);
 346     if(lshells.Extent() <2) continue;
 347     TopoDS_CompSolid aCompSolid;
 348     BRep_Builder aB;
 349     aB.MakeCompSolid(aCompSolid);
 350     isDone = (theShape.ShapeType() != TopAbs_COMPSOLID || isDone);
 351     Standard_Integer nbSol = 0;
 352
 353     for(TopTools_ListIteratorOfListOfShape lItSh(lshells);lItSh.More(); lItSh.Next()) {
 354       if(ShellSolid.Contains(lItSh.Value())) {
 355         const TopoDS_Shape& aShape = ShellSolid.FindFromKey(lItSh.Value());
 356         TopExp_Explorer aExpSol(aShape, TopAbs_SOLID);
 357
 358         for(;aExpSol.More(); aExpSol.Next())
 359         {
 360           aB.Add(aCompSolid,aExpSol.Current());
 361           nbSol++;
 362         }
 363
 364       }
 365     }
 366     if(nbSol >1)
 367     {
 368       for(TopTools_ListIteratorOfListOfShape lItSh1(lshells);lItSh1.More(); lItSh1.Next())
 369       {
 370         if(ShellSolid.Contains(lItSh1.Value()))
 371           ShellSolid.ChangeFromKey(lItSh1.Value()) = aCompSolid;
 372       }
 373     }
 374
 375   }
 376   for(Standard_Integer kk =1 ; kk <= ShellSolid.Extent();kk++)
 377     if(!aMapSolids.Contains(ShellSolid.FindFromIndex(kk)))
 378        aMapSolids.Add(ShellSolid.FindFromIndex(kk));
 379   isDone = (aMapSolids.Extent() >1 || isDone);
 380   return isDone;
 381 }
 382 //=======================================================================
 383 //function : Perform
 384 //purpose  :
 385 //=======================================================================
 386
 387 Standard_Boolean ShapeFix_Solid::Perform(const Handle(Message_ProgressIndicator)& theProgress)
 388 {
 389
 390   Standard_Boolean status = Standard_False;
 391   if ( Context().IsNull() )
 392     SetContext ( new ShapeBuild_ReShape );
 393   myFixShell->SetContext(Context());
 394
 395   Standard_Integer NbShells = 0;
 396   TopoDS_Shape S = Context()->Apply ( myShape );
 397
 398   // Calculate number of underlying shells
 399   Standard_Integer aNbShells = 0;
 400   for ( TopExp_Explorer aExpSh(S, TopAbs_SHELL); aExpSh.More(); aExpSh.Next() )
 401     aNbShells++;
 402
 403   // Start progress scope (no need to check if progress exists -- it is safe)
 404   Message_ProgressSentry aPSentry(theProgress, "Fixing solid stage", 0, 2, 1);
 405
 406   if ( NeedFix(myFixShellMode) )
 407   {
 408     // Start progress scope (no need to check if progress exists -- it is safe)
 409     Message_ProgressSentry aPSentryFixShell(theProgress, "Fixing shell", 0, aNbShells, 1);
 410
 411     // Fix shell by shell using ShapeFix_Shell tool
 412     for ( TopExp_Explorer aExpSh(S, TopAbs_SHELL); aExpSh.More() && aPSentryFixShell.More(); aExpSh.Next(), aPSentryFixShell.Next() )
 413     {
 414       TopoDS_Shape sh = aExpSh.Current();
 415
 416       myFixShell->Init( TopoDS::Shell(sh) );
 417       if ( myFixShell->Perform(theProgress) )
 418       {
 419         status = Standard_True;
 420         myStatus |= ShapeExtend::EncodeStatus ( ShapeExtend_DONE1 );
 421       }
 422       NbShells+= myFixShell->NbShells();
 423     }
 424
 425     // Halt algorithm in case of user's abort
 426     if ( !aPSentryFixShell.More() )
 427       return Standard_False;
 428   }
 429   else
 430   {
 431     NbShells = aNbShells;
 432   }
 433
 434   // Switch to the second stage
 435   aPSentry.Next();
 436
 437   if(NbShells ==1) {
 438     TopoDS_Shape tmpShape = Context()->Apply(myShape);
 439     TopExp_Explorer aExp(tmpShape,TopAbs_SHELL);
 440     Standard_Boolean isClosed = Standard_False;
 441     if(aExp.More()) {
 442       TopoDS_Shell  aShtmp = TopoDS::Shell(aExp.Current());
 443       ShapeAnalysis_FreeBounds sfb(aShtmp);
 444       TopoDS_Compound aC1 = sfb.GetClosedWires();
 445       TopoDS_Compound aC2 = sfb.GetOpenWires();
 446       Standard_Integer numedge =0;
 447       TopExp_Explorer aExp1(aC1,TopAbs_EDGE);
 448       for( ; aExp1.More(); aExp1.Next())
 449         numedge++;
 450       for(aExp1.Init(aC2,TopAbs_EDGE) ; aExp1.More(); aExp1.Next())
 451         numedge++;
 452       isClosed = (!numedge);
 453       aShtmp.Closed(isClosed);
 454     }
 455
 456     if(isClosed || myCreateOpenSolidMode) {
 457       TopoDS_Iterator itersh(tmpShape);
 458       TopoDS_Shell aShell;
 459       if(itersh.More() && itersh.Value().ShapeType() == TopAbs_SHELL)
 460         aShell = TopoDS::Shell(itersh.Value());
 461       if(!aShell.IsNull()) {
 462         TopoDS_Solid aSol = SolidFromShell(aShell);
 463         if(ShapeExtend::DecodeStatus(myStatus,ShapeExtend_DONE2)) {
 464           SendWarning (Message_Msg ("FixAdvSolid.FixOrientation.MSG20"));// Orientaion of shell was corrected.
 465           Context()->Replace(tmpShape,aSol);
 466           tmpShape = aSol;
 467         }
 468       }
 469       mySolid  = TopoDS::Solid(tmpShape);
 470     }
 471     else {
 472       status = Standard_True;
 473       myStatus |= ShapeExtend::EncodeStatus ( ShapeExtend_DONE3 );
 474       TopoDS_Iterator aIt(tmpShape,Standard_False);
 475       Context()->Replace(tmpShape,aIt.Value());
 476       SendFail (Message_Msg ("FixAdvSolid.FixShell.MSG10")); // Solid can not be created from open shell.
 477     }
 478   }
 479   else {
 480     TopoDS_Shape aResShape = Context()->Apply(myShape);
 481     TopTools_SequenceOfShape aSeqShells;
 482     TopTools_IndexedMapOfShape aMapSolids;
 483     if(CreateSolids(aResShape,aMapSolids)) {
 484       SendWarning (Message_Msg ("FixAdvSolid.FixOrientation.MSG20"));// Orientaion of shell was corrected..
 485       if(aMapSolids.Extent() ==1) {
 486         TopoDS_Shape aResSol = aMapSolids.FindKey(1);
 487         if(aResShape.ShapeType() == TopAbs_SHELL && myCreateOpenSolidMode) {
 488           TopoDS_Solid solid;
 489           BRep_Builder B;
 490           B.MakeSolid (solid);
 491           B.Add (solid,aResSol);
 492           mySolid = solid;
 493         }
 494         else {
 495           mySolid =aResSol;
 496           if(aResSol.ShapeType() == TopAbs_SHELL)
 497             SendFail (Message_Msg ("FixAdvSolid.FixShell.MSG10")); // Solid can not be created from open shell.
 498         }
 499         Context()->Replace(aResShape,mySolid);
 500       }
 501       else if(aMapSolids.Extent() >1) {
 502         SendWarning (Message_Msg ("FixAdvSolid.FixOrientation.MSG30"));// Bad connected solid a few solids were created.
 503         BRep_Builder aB;
 504         TopoDS_Compound aComp;
 505         aB.MakeCompound(aComp);
 506         Message_ProgressSentry aPSentryCreatingSolid(theProgress, "Creating solid",
 507                                         0, aMapSolids.Extent(), 1);
 508         for(Standard_Integer i =1; (i <= aMapSolids.Extent()) && (aPSentryCreatingSolid.More());
 509             i++, aPSentryCreatingSolid.Next())
 510         {
 511           TopoDS_Shape aResSh =aMapSolids.FindKey(i);
 512           if(aResShape.ShapeType() == TopAbs_SHELL && myCreateOpenSolidMode) {
 513             aResSh.Closed(Standard_False);
 514             TopoDS_Solid solid;
 515             BRep_Builder B;
 516             B.MakeSolid (solid);
 517             B.Add (solid,aResSh);
 518             aResSh = solid;
 519           }
 520           else if (aResShape.ShapeType() == TopAbs_SHELL)
 521             SendFail(Message_Msg ("FixAdvSolid.FixShell.MSG10")); // Solid can not be created from open shell.
 522           aB.Add(aComp,aResSh);
 523
 524         }
 525         if ( !aPSentryCreatingSolid.More() )
 526           return Standard_False; // aborted execution
 527         Context()->Replace(aResShape,aComp);
 528       }
 529     }
 530   }
 531   myShape = Context()->Apply(myShape);
 532   return status;
 533 }
 534 //=======================================================================
 535 //function : Shape
 536 //purpose  :
 537 //=======================================================================
 538
 539 TopoDS_Shape ShapeFix_Solid::Shape()
 540 {
 541  return myShape;
 542 }
 543 //=======================================================================
 544 //function : SolidFromShell
 545 //purpose  :
 546 //=======================================================================
 547
 548 TopoDS_Solid ShapeFix_Solid::SolidFromShell (const TopoDS_Shell& shell)
 549 {
 550   TopoDS_Shell sh = shell;
 551   if (!sh.Free ()) sh.Free(Standard_True);
 552
 553   TopoDS_Solid solid;
 554   BRep_Builder B;
 555   B.MakeSolid (solid);
 556   B.Add (solid,sh);
 557 //   Pas encore fini : il faut une bonne orientation
 558   try {
 559     OCC_CATCH_SIGNALS
 560     BRepClass3d_SolidClassifier bsc3d (solid);
 561     Standard_Real t = Precision::Confusion();    // tolerance moyenne
 562     bsc3d.PerformInfinitePoint(t);
 563
 564     if (bsc3d.State() == TopAbs_IN) {
 565       //         Ensuite, inverser C-A-D REPRENDRE LES SHELLS
 566       //         (l inversion du solide n est pas bien prise en compte)
 567       sh = shell;
 568       if (!sh.Free ()) sh.Free(Standard_True);
 569       TopoDS_Solid soli2;
 570       B.MakeSolid (soli2);    // on recommence
 571       sh.Reverse();
 572       B.Add (soli2,sh);
 573       solid = soli2;
 574       myStatus |= ShapeExtend::EncodeStatus ( ShapeExtend_DONE2 );
 575     }
 576   }
 577   catch(Standard_Failure const& anException) {
 578 #ifdef OCCT_DEBUG
 579     cout << "Warning: ShapeFix_Solid::SolidFromShell: Exception: ";
 580     anException.Print(cout); cout << endl;
 581 #endif
 582     (void)anException;
 583     return solid;
 584   }
 585   return solid;
 586 }
 587
 588 //=======================================================================
 589 //function : Status
 590 //purpose  :
 591 //=======================================================================
 592
 593 Standard_Boolean ShapeFix_Solid::Status (const ShapeExtend_Status theStatus) const
 594 {
 595   return ShapeExtend::DecodeStatus (myStatus, theStatus);
 596 }
 597
 598 //=======================================================================
 599 //function : Solid
 600 //purpose  :
 601 //=======================================================================
 602
 603  TopoDS_Shape ShapeFix_Solid::Solid() const
 604 {
 605  return mySolid;
 606 }
 607
 608 //=======================================================================
 609 //function : SetMsgRegistrator
 610 //purpose  :
 611 //=======================================================================
 612
 613 void ShapeFix_Solid::SetMsgRegistrator(const Handle(ShapeExtend_BasicMsgRegistrator)& msgreg)
 614 {
 615   ShapeFix_Root::SetMsgRegistrator ( msgreg );
 616   myFixShell->SetMsgRegistrator ( msgreg );
 617 }
 618
 619 //=======================================================================
 620 //function : SetPrecision
 621 //purpose  :
 622 //=======================================================================
 623
 624 void ShapeFix_Solid::SetPrecision (const Standard_Real preci)
 625 {
 626   ShapeFix_Root::SetPrecision ( preci );
 627   myFixShell->SetPrecision ( preci );
 628 }
 629
 630 //=======================================================================
 631 //function : SetMinTolerance
 632 //purpose  :
 633 //=======================================================================
 634
 635 void ShapeFix_Solid::SetMinTolerance (const Standard_Real mintol)
 636 {
 637   ShapeFix_Root::SetMinTolerance ( mintol );
 638   myFixShell->SetMinTolerance ( mintol );
 639 }
 640
 641 //=======================================================================
 642 //function : SetMaxTolerance
 643 //purpose  :
 644 //=======================================================================
 645
 646 void ShapeFix_Solid::SetMaxTolerance (const Standard_Real maxtol)
 647 {
 648   ShapeFix_Root::SetMaxTolerance ( maxtol );
 649   myFixShell->SetMaxTolerance ( maxtol );
 650 }