Salome HOME
Join modifications from BR_Dev_For_4_0 tag V4_1_1.
[modules/geom.git] / src / GEOMImpl / GEOMImpl_PartitionDriver.cxx
index fc91fd19cb281debee98be7cefc3530d04b55481..b681770981346109afa803f8a847d5c356b30871 100644 (file)
@@ -1,18 +1,18 @@
 // Copyright (C) 2005  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
-// 
+//
 // This library is free software; you can redistribute it and/or
 // modify it under the terms of the GNU Lesser General Public
-// License as published by the Free Software Foundation; either 
+// License as published by the Free Software Foundation; either
 // version 2.1 of the License.
-// 
-// This library is distributed in the hope that it will be useful 
-// but WITHOUT ANY WARRANTY; without even the implied warranty of 
-// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU 
+//
+// This library is distributed in the hope that it will be useful
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 // Lesser General Public License for more details.
 //
-// You should have received a copy of the GNU Lesser General Public  
-// License along with this library; if not, write to the Free Software 
+// You should have received a copy of the GNU Lesser General Public
+// License along with this library; if not, write to the Free Software
 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
 //
 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
 #include <GEOM_Object.hxx>
 #include <GEOM_Function.hxx>
 
-#include <NMTAlgo_Splitter1.hxx>
+//#include <NMTAlgo_Splitter1.hxx>
+#include <GEOMAlgo_Splitter.hxx>
+#include <TopTools_IndexedMapOfShape.hxx>
 
 #include <TDataStd_IntegerArray.hxx>
 
 #include <BRep_Tool.hxx>
 #include <BRepAlgo.hxx>
 
+#include <TopAbs.hxx>
+#include <TopExp.hxx>
 #include <TopoDS.hxx>
 #include <TopoDS_Shape.hxx>
 #include <TopoDS_Vertex.hxx>
 #include <TopoDS_Wire.hxx>
-#include <TopAbs.hxx>
-#include <TopExp.hxx>
+#include <TopoDS_Iterator.hxx>
 #include <TopTools_MapOfShape.hxx>
 #include <TopTools_ListIteratorOfListOfShape.hxx>
 
 //=======================================================================
 //function : GetID
 //purpose  :
-//======================================================================= 
+//=======================================================================
 const Standard_GUID& GEOMImpl_PartitionDriver::GetID()
 {
   static Standard_GUID aPartitionDriver("FF1BBB22-5D14-4df2-980B-3A668264EA16");
-  return aPartitionDriver; 
+  return aPartitionDriver;
 }
 
 
 //=======================================================================
 //function : GEOMImpl_PartitionDriver
-//purpose  : 
+//purpose  :
 //=======================================================================
-GEOMImpl_PartitionDriver::GEOMImpl_PartitionDriver() 
+GEOMImpl_PartitionDriver::GEOMImpl_PartitionDriver()
 {
 }
 
+//=======================================================================
+//function : SimplifyCompound
+//purpose  :
+//=======================================================================
+static void PrepareShapes (const TopoDS_Shape&   theShape,
+                           Standard_Integer      theType,
+                           TopTools_ListOfShape& theSimpleList)
+{
+  if (theType == PARTITION_NO_SELF_INTERSECTIONS ||
+      theShape.ShapeType() != TopAbs_COMPOUND) {
+    theSimpleList.Append(theShape);
+    return;
+  }
+
+  // explode compound on simple shapes to allow their intersections
+  TopoDS_Iterator It (theShape, Standard_True, Standard_True);
+  TopTools_MapOfShape mapShape;
+  for (; It.More(); It.Next()) {
+    if (mapShape.Add(It.Value())) {
+      TopoDS_Shape curSh = It.Value();
+      PrepareShapes(curSh, theType, theSimpleList);
+    }
+  }
+}
+
 //=======================================================================
 //function : Execute
 //purpose  :
-//======================================================================= 
+//=======================================================================
 Standard_Integer GEOMImpl_PartitionDriver::Execute(TFunction_Logbook& log) const
 {
-  if (Label().IsNull()) return 0;    
+  if (Label().IsNull()) return 0;
   Handle(GEOM_Function) aFunction = GEOM_Function::GetFunction(Label());
 
   GEOMImpl_IPartition aCI (aFunction);
   Standard_Integer aType = aFunction->GetType();
 
   TopoDS_Shape aShape;
-  NMTAlgo_Splitter1 PS;
+  //sklNMTAlgo_Splitter1 PS;
+  GEOMAlgo_Splitter PS;
 
-  if (aType == PARTITION_PARTITION) {
+  if (aType == PARTITION_PARTITION || aType == PARTITION_NO_SELF_INTERSECTIONS)
+  {
     Handle(TColStd_HSequenceOfTransient) aShapes  = aCI.GetShapes();
     Handle(TColStd_HSequenceOfTransient) aTools   = aCI.GetTools();
     Handle(TColStd_HSequenceOfTransient) aKeepIns = aCI.GetKeepIns();
     Handle(TColStd_HSequenceOfTransient) aRemIns  = aCI.GetRemoveIns();
     Handle(TColStd_HArray1OfInteger) aMaterials   = aCI.GetMaterials();
-    Standard_Boolean DoRemoveWebs = !aMaterials.IsNull();
-
-    unsigned int ind, nbshapes = 0;
-    nbshapes += aShapes->Length() + aTools->Length();
-    nbshapes += aKeepIns->Length() + aRemIns->Length();
+    //skl Standard_Boolean DoRemoveWebs = !aMaterials.IsNull();
 
-    TopTools_MapOfShape ShapesMap(nbshapes), ToolsMap(nbshapes);
+    unsigned int ind;
+    //unsigned int ind, nbshapes = 0;
+    //nbshapes += aShapes->Length() + aTools->Length();
+    //nbshapes += aKeepIns->Length() + aRemIns->Length();
+    //TopTools_MapOfShape ShapesMap(nbshapes), ToolsMap(nbshapes);
+    TopTools_MapOfShape ShapesMap, ToolsMap;
 
     // add object shapes that are in ListShapes;
     for (ind = 1; ind <= aShapes->Length(); ind++) {
@@ -104,11 +135,18 @@ Standard_Integer GEOMImpl_PartitionDriver::Execute(TFunction_Logbook& log) const
       if (aShape_i.IsNull()) {
         Standard_NullObject::Raise("In Partition a shape is null");
       }
-      if (ShapesMap.Add(aShape_i)) {
-        PS.AddShape(aShape_i);
-        if (DoRemoveWebs) {
-          if (aMaterials->Length() >= ind)
-            PS.SetMaterial(aShape_i, aMaterials->Value(ind));
+      //
+      TopTools_ListOfShape aSimpleShapes;
+      PrepareShapes(aShape_i, aType, aSimpleShapes);
+      TopTools_ListIteratorOfListOfShape aSimpleIter (aSimpleShapes);
+      for (; aSimpleIter.More(); aSimpleIter.Next()) {
+        const TopoDS_Shape& aSimpleSh = aSimpleIter.Value();
+        if (ShapesMap.Add(aSimpleSh)) {
+          PS.AddShape(aSimpleSh);
+          //skl if (DoRemoveWebs) {
+          //skl if (aMaterials->Length() >= ind)
+          //skl PS.SetMaterial(aSimpleSh, aMaterials->Value(ind));
+          //skl }
         }
       }
     }
@@ -120,8 +158,16 @@ Standard_Integer GEOMImpl_PartitionDriver::Execute(TFunction_Logbook& log) const
       if (aShape_i.IsNull()) {
         Standard_NullObject::Raise("In Partition a tool shape is null");
       }
-      if (!ShapesMap.Contains(aShape_i) && ToolsMap.Add(aShape_i))
-        PS.AddTool(aShape_i);
+      //
+      TopTools_ListOfShape aSimpleShapes;
+      PrepareShapes(aShape_i, aType, aSimpleShapes);
+      TopTools_ListIteratorOfListOfShape aSimpleIter (aSimpleShapes);
+      for (; aSimpleIter.More(); aSimpleIter.Next()) {
+        const TopoDS_Shape& aSimpleSh = aSimpleIter.Value();
+        if (!ShapesMap.Contains(aSimpleSh) && ToolsMap.Add(aSimpleSh)) {
+          PS.AddTool(aSimpleSh);
+        }
+      }
     }
 
     // add shapes that are in ListKeepInside, as object shapes;
@@ -131,8 +177,15 @@ Standard_Integer GEOMImpl_PartitionDriver::Execute(TFunction_Logbook& log) const
       if (aShape_i.IsNull()) {
         Standard_NullObject::Raise("In Partition a Keep Inside shape is null");
       }
-      if (!ToolsMap.Contains(aShape_i) && ShapesMap.Add(aShape_i))
-        PS.AddShape(aShape_i);
+      //
+      TopTools_ListOfShape aSimpleShapes;
+      PrepareShapes(aShape_i, aType, aSimpleShapes);
+      TopTools_ListIteratorOfListOfShape aSimpleIter (aSimpleShapes);
+      for (; aSimpleIter.More(); aSimpleIter.Next()) {
+        const TopoDS_Shape& aSimpleSh = aSimpleIter.Value();
+        if (!ToolsMap.Contains(aSimpleSh) && ShapesMap.Add(aSimpleSh))
+          PS.AddShape(aSimpleSh);
+      }
     }
 
     // add shapes that are in ListRemoveInside, as object shapes;
@@ -142,14 +195,25 @@ Standard_Integer GEOMImpl_PartitionDriver::Execute(TFunction_Logbook& log) const
       if (aShape_i.IsNull()) {
         Standard_NullObject::Raise("In Partition a Remove Inside shape is null");
       }
-      if (!ToolsMap.Contains(aShape_i) && ShapesMap.Add(aShape_i))
-        PS.AddShape(aShape_i);
+      //
+      TopTools_ListOfShape aSimpleShapes;
+      PrepareShapes(aShape_i, aType, aSimpleShapes);
+      TopTools_ListIteratorOfListOfShape aSimpleIter (aSimpleShapes);
+      for (; aSimpleIter.More(); aSimpleIter.Next()) {
+        const TopoDS_Shape& aSimpleSh = aSimpleIter.Value();
+        if (!ToolsMap.Contains(aSimpleSh) && ShapesMap.Add(aSimpleSh))
+          PS.AddShape(aSimpleSh);
+      }
     }
 
-    PS.Compute();
-    PS.SetRemoveWebs(DoRemoveWebs);
-    PS.Build((TopAbs_ShapeEnum) aCI.GetLimit());
+    PS.SetLimitMode(aCI.GetKeepNonlimitShapes());
+    PS.SetLimit( (TopAbs_ShapeEnum)aCI.GetLimit() );
+    PS.Perform();
 
+    //skl PS.Compute();
+    //skl PS.SetRemoveWebs(!DoRemoveWebs);
+    //skl PS.Build((TopAbs_ShapeEnum) aCI.GetLimit());
+    /*skl
     // suppress result outside of shapes in KInsideMap
     for (ind = 1; ind <= aKeepIns->Length(); ind++) {
       Handle(GEOM_Function) aRefShape = Handle(GEOM_Function)::DownCast(aKeepIns->Value(ind));
@@ -163,8 +227,10 @@ Standard_Integer GEOMImpl_PartitionDriver::Execute(TFunction_Logbook& log) const
       TopoDS_Shape aShape_i = aRefShape->GetValue();
       PS.RemoveShapesInside(aShape_i);
     }
-
-  } else if (aType == PARTITION_HALF) {
+    */
+  }
+  else if (aType == PARTITION_HALF)
+  {
     Handle(GEOM_Function) aRefShape = aCI.GetShape();
     Handle(GEOM_Function) aRefPlane = aCI.GetPlane();
     TopoDS_Shape aShapeArg = aRefShape->GetValue();
@@ -180,9 +246,10 @@ Standard_Integer GEOMImpl_PartitionDriver::Execute(TFunction_Logbook& log) const
     // add tool shapes that are in ListTools and not in ListShapes;
     PS.AddTool(aPlaneArg);
 
-    PS.Compute();
-    PS.SetRemoveWebs(Standard_False);
-    PS.Build(aShapeArg.ShapeType());
+    //skl PS.Compute();
+    PS.Perform();
+    //PS.SetRemoveWebs(Standard_False);
+    //PS.Build(aShapeArg.ShapeType());
 
   } else {
   }
@@ -200,6 +267,9 @@ Standard_Integer GEOMImpl_PartitionDriver::Execute(TFunction_Logbook& log) const
   TopTools_IndexedMapOfShape aResIndices;
   TopExp::MapShapes(aShape, aResIndices);
 
+  // Map: source_shape/images of source_shape in Result
+  const TopTools_IndexedDataMapOfShapeListOfShape& aMR = PS.ImagesResult();
+
   // history for all argument shapes
   TDF_LabelSequence aLabelSeq;
   aFunction->GetDependency(aLabelSeq);
@@ -222,7 +292,9 @@ Standard_Integer GEOMImpl_PartitionDriver::Execute(TFunction_Logbook& log) const
 
     for (Standard_Integer ie = 1; ie <= nbArgumentEntities; ie++) {
       TopoDS_Shape anEntity = anArgumentIndices.FindKey(ie);
-      const TopTools_ListOfShape& aModified = PS.Modified(anEntity);
+      if (!aMR.Contains(anEntity)) continue;
+
+      const TopTools_ListOfShape& aModified = aMR.FindFromKey(anEntity);
       Standard_Integer nbModified = aModified.Extent();
 
       if (nbModified > 0) {
@@ -230,10 +302,11 @@ Standard_Integer GEOMImpl_PartitionDriver::Execute(TFunction_Logbook& log) const
         Handle(TDataStd_IntegerArray) anAttr =
           TDataStd_IntegerArray::Set(aWhatHistoryLabel, 1, nbModified);
 
+        int ih = 1;
         TopTools_ListIteratorOfListOfShape itM (aModified);
-        for (int im = 1; itM.More(); itM.Next(), ++im) {
+        for (; itM.More(); itM.Next(), ++ih) {
           int id = aResIndices.FindIndex(itM.Value());
-          anAttr->SetValue(im, id);
+          anAttr->SetValue(ih, id);
         }
       }
     }
@@ -248,17 +321,17 @@ Standard_Integer GEOMImpl_PartitionDriver::Execute(TFunction_Logbook& log) const
 //=======================================================================
 //function :  GEOMImpl_PartitionDriver_Type_
 //purpose  :
-//======================================================================= 
+//=======================================================================
 Standard_EXPORT Handle_Standard_Type& GEOMImpl_PartitionDriver_Type_()
 {
 
   static Handle_Standard_Type aType1 = STANDARD_TYPE(TFunction_Driver);
   if ( aType1.IsNull()) aType1 = STANDARD_TYPE(TFunction_Driver);
   static Handle_Standard_Type aType2 = STANDARD_TYPE(MMgt_TShared);
-  if ( aType2.IsNull()) aType2 = STANDARD_TYPE(MMgt_TShared); 
+  if ( aType2.IsNull()) aType2 = STANDARD_TYPE(MMgt_TShared);
   static Handle_Standard_Type aType3 = STANDARD_TYPE(Standard_Transient);
   if ( aType3.IsNull()) aType3 = STANDARD_TYPE(Standard_Transient);
+
 
   static Handle_Standard_Transient _Ancestors[]= {aType1,aType2,aType3,NULL};
   static Handle_Standard_Type _aType = new Standard_Type("GEOMImpl_PartitionDriver",
@@ -273,7 +346,7 @@ Standard_EXPORT Handle_Standard_Type& GEOMImpl_PartitionDriver_Type_()
 //=======================================================================
 //function : DownCast
 //purpose  :
-//======================================================================= 
+//=======================================================================
 const Handle(GEOMImpl_PartitionDriver) Handle(GEOMImpl_PartitionDriver)::DownCast(const Handle(Standard_Transient)& AnObject)
 {
   Handle(GEOMImpl_PartitionDriver) _anOtherObject;