src/Model/Model_BodyBuilder.cpp

   1 // Copyright (C) 2014-2017  CEA/DEN, EDF R&D
   2 //
   3 // This library is free software; you can redistribute it and/or
   4 // modify it under the terms of the GNU Lesser General Public
   5 // License as published by the Free Software Foundation; either
   6 // version 2.1 of the License, or (at your option) any later version.
   7 //
   8 // This library is distributed in the hope that it will be useful,
   9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 // Lesser General Public License for more details.
  12 //
  13 // You should have received a copy of the GNU Lesser General Public
  14 // License along with this library; if not, write to the Free Software
  15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
  16 //
  17 // See http://www.salome-platform.org/ or
  18 // email : webmaster.salome@opencascade.com<mailto:webmaster.salome@opencascade.com>
  19 //
  20
  21 #include <Model_BodyBuilder.h>
  22
  23 #include <Model_Data.h>
  24 #include <Model_Document.h>
  25 #include <TNaming_Builder.hxx>
  26 #include <TNaming_NamedShape.hxx>
  27 #include <TNaming_Iterator.hxx>
  28 #include <TNaming_Tool.hxx>
  29 #include <TDataStd_Name.hxx>
  30 #include <TDataStd_Integer.hxx>
  31 #include <TopoDS.hxx>
  32 #include <TopoDS_Face.hxx>
  33 #include <TDF_ChildIterator.hxx>
  34 #include <TDF_ChildIDIterator.hxx>
  35 #include <TDF_Reference.hxx>
  36 #include <TopTools_MapOfShape.hxx>
  37 #include <TopExp_Explorer.hxx>
  38 #include <TopTools_ListOfShape.hxx>
  39 #include <TopTools_ListIteratorOfListOfShape.hxx>
  40 #include <TopTools_DataMapOfShapeListOfShape.hxx>
  41 #include <TopTools_DataMapIteratorOfDataMapOfShapeListOfShape.hxx>
  42 #include <TopTools_DataMapIteratorOfDataMapOfShapeShape.hxx>
  43 #include <TopTools_MapIteratorOfMapOfShape.hxx>
  44 #include <TopTools_IndexedDataMapOfShapeListOfShape.hxx>
  45 #include <TopTools_IndexedMapOfShape.hxx>
  46 #include <TopTools_DataMapOfShapeShape.hxx>
  47 #include <TopExp.hxx>
  48 #include <BRepTools.hxx>
  49 #include <BRep_Tool.hxx>
  50 #include <GeomAPI_Shape.h>
  51 #include <GeomAlgoAPI_MakeShape.h>
  52 #include <GeomAlgoAPI_SortListOfShapes.h>
  53 #include <Config_PropManager.h>
  54 // DEB
  55 //#include <TCollection_AsciiString.hxx>
  56 //#include <TDF_Tool.hxx>
  57 //#define DEB_IMPORT 1
  58
  59 Model_BodyBuilder::Model_BodyBuilder(ModelAPI_Object* theOwner)
  60 : ModelAPI_BodyBuilder(theOwner)
  61 {
  62 }
  63
  64 // Converts evolution of naming shape to selection evelution and back to avoid
  65 // naming support on the disabled results. Deeply in the labels tree, recursively.
  66 static void evolutionToSelectionRec(TDF_Label theLab, const bool theFlag) {
  67   std::list<std::pair<TopoDS_Shape, TopoDS_Shape> > aShapePairs; // to store old and new shapes
  68   Handle(TNaming_NamedShape) aName;
  69   int anEvolution = -1;
  70   if (theLab.FindAttribute(TNaming_NamedShape::GetID(), aName)) {
  71     TNaming_Evolution aNSEvol = aName->Evolution();
  72     if ((aNSEvol == TNaming_SELECTED && theFlag) ||
  73         (aNSEvol != TNaming_SELECTED && !theFlag)) { // nothing to do, it is already correct
  74       return;
  75     }
  76     anEvolution = (int)(aNSEvol);
  77     if (!theFlag) {
  78       Handle(TDataStd_Integer) anAttrEvol;
  79       if (theLab.FindAttribute(TDataStd_Integer::GetID(), anAttrEvol)) {
  80         anEvolution = anAttrEvol->Get();
  81       }
  82     } else {
  83       TDataStd_Integer::Set(theLab, anEvolution);
  84     }
  85
  86     for(TNaming_Iterator anIter(aName); anIter.More(); anIter.Next()) {
  87       // iterator goes in reversed order relatively to the Builder, to, make the list reversed
  88       aShapePairs.push_front(std::pair<TopoDS_Shape, TopoDS_Shape>
  89         (anIter.OldShape(), anIter.NewShape()));
  90     }
  91
  92     // create new
  93     TNaming_Builder aBuilder(theLab);
  94     TNaming_Evolution anEvol = (TNaming_Evolution)(anEvolution);
  95     std::list<std::pair<TopoDS_Shape, TopoDS_Shape> >::iterator aPairsIter = aShapePairs.begin();
  96     for(; aPairsIter != aShapePairs.end(); aPairsIter++) {
  97       if (theFlag) { // disabled => make selection
  98         if (anEvolution == TNaming_DELETE) // issue 2274 : don't put too many same null shapes
  99           aBuilder.Select(aPairsIter->first, aPairsIter->first);
 100         else if (anEvolution == TNaming_PRIMITIVE)
 101           aBuilder.Select(aPairsIter->second, aPairsIter->second);
 102         else
 103           aBuilder.Select(aPairsIter->second, aPairsIter->first);
 104       } else if (anEvol == TNaming_GENERATED) {
 105         aBuilder.Generated(aPairsIter->first, aPairsIter->second);
 106       } else if (anEvol == TNaming_MODIFY) {
 107         aBuilder.Modify(aPairsIter->first, aPairsIter->second);
 108       } else if (anEvol == TNaming_DELETE) {
 109         aBuilder.Delete(aPairsIter->first);
 110       } else if (anEvol == TNaming_PRIMITIVE) {
 111         aBuilder.Generated(aPairsIter->second);
 112       } else if (anEvol == TNaming_SELECTED) {
 113         aBuilder.Select(aPairsIter->second, aPairsIter->first);
 114       }
 115     }
 116   }
 117   // recursive call for all sub-labels
 118   TDF_ChildIterator anIter(theLab, Standard_False);
 119   for(; anIter.More(); anIter.Next()) {
 120     evolutionToSelectionRec(anIter.Value(), theFlag);
 121   }
 122 }
 123
 124 void Model_BodyBuilder::evolutionToSelection(const bool theFlag)
 125 {
 126   std::shared_ptr<Model_Data> aData = std::dynamic_pointer_cast<Model_Data>(data());
 127   if (!aData) // unknown case
 128     return;
 129   TDF_Label& aShapeLab = aData->shapeLab();
 130   evolutionToSelectionRec(aShapeLab, theFlag);
 131 }
 132
 133 void Model_BodyBuilder::store(const std::shared_ptr<GeomAPI_Shape>& theShape,
 134                               const bool theIsStoreSameShapes)
 135 {
 136   std::shared_ptr<Model_Data> aData = std::dynamic_pointer_cast<Model_Data>(data());
 137   if (aData) {
 138     TDF_Label& aShapeLab = aData->shapeLab();
 139     // clean builders
 140     clean();
 141     // store the new shape as primitive
 142     TNaming_Builder aBuilder(aShapeLab);
 143     if (!theShape)
 144       return;  // bad shape
 145     TopoDS_Shape aShape = theShape->impl<TopoDS_Shape>();
 146     if (aShape.IsNull())
 147       return;  // null shape inside
 148
 149     if(!theIsStoreSameShapes) {
 150       Handle(TNaming_NamedShape) aNS = TNaming_Tool::NamedShape(aShape, aShapeLab);
 151       if(!aNS.IsNull() && !aNS->IsEmpty()) {
 152         // This shape is already in document, store reference instead of shape;
 153         const TDF_Label aFoundLabel = aNS->Label();
 154         TDF_Reference::Set(aShapeLab, aFoundLabel);
 155         aShapeLab.ForgetAttribute(TNaming_NamedShape::GetID());
 156         return;
 157       }
 158     }
 159
 160     aBuilder.Generated(aShape);
 161     // register name
 162     aShapeLab.ForgetAttribute(TDF_Reference::GetID());
 163     if(!aBuilder.NamedShape()->IsEmpty()) {
 164       Handle(TDataStd_Name) anAttr;
 165       if(aBuilder.NamedShape()->Label().FindAttribute(TDataStd_Name::GetID(),anAttr)) {
 166         std::string aName (TCollection_AsciiString(anAttr->Get()).ToCString());
 167         if(!aName.empty()) {
 168           std::shared_ptr<Model_Document> aDoc =
 169             std::dynamic_pointer_cast<Model_Document>(document());
 170           aDoc->addNamingName(aBuilder.NamedShape()->Label(), aName);
 171         }
 172       }
 173     }
 174   }
 175 }
 176
 177 void Model_BodyBuilder::storeGenerated(const std::shared_ptr<GeomAPI_Shape>& theFromShape,
 178   const std::shared_ptr<GeomAPI_Shape>& theToShape)
 179 {
 180   std::shared_ptr<Model_Data> aData = std::dynamic_pointer_cast<Model_Data>(data());
 181   if (aData) {
 182     TDF_Label& aShapeLab = aData->shapeLab();
 183     // clean builders
 184     clean();
 185     // store the new shape as primitive
 186     TNaming_Builder aBuilder(aShapeLab);
 187     if (!theFromShape || !theToShape)
 188       return;  // bad shape
 189     TopoDS_Shape aShapeBasis = theFromShape->impl<TopoDS_Shape>();
 190     if (aShapeBasis.IsNull())
 191       return;  // null shape inside
 192     TopoDS_Shape aShapeNew = theToShape->impl<TopoDS_Shape>();
 193     if (aShapeNew.IsNull())
 194       return;  // null shape inside
 195     aBuilder.Generated(aShapeBasis, aShapeNew);
 196     // register name
 197     if(!aBuilder.NamedShape()->IsEmpty()) {
 198       Handle(TDataStd_Name) anAttr;
 199       if(aBuilder.NamedShape()->Label().FindAttribute(TDataStd_Name::GetID(),anAttr)) {
 200         std::string aName (TCollection_AsciiString(anAttr->Get()).ToCString());
 201         if(!aName.empty()) {
 202           std::shared_ptr<Model_Document> aDoc =
 203             std::dynamic_pointer_cast<Model_Document>(document());
 204           aDoc->addNamingName(aBuilder.NamedShape()->Label(), aName);
 205         }
 206       }
 207     }
 208   }
 209 }
 210
 211 void Model_BodyBuilder::storeModified(const std::shared_ptr<GeomAPI_Shape>& theOldShape,
 212   const std::shared_ptr<GeomAPI_Shape>& theNewShape, const int theDecomposeSolidsTag)
 213 {
 214   std::shared_ptr<Model_Data> aData = std::dynamic_pointer_cast<Model_Data>(data());
 215   if (aData) {
 216     TDF_Label& aShapeLab = aData->shapeLab();
 217     // clean builders
 218     clean();
 219     // store the new shape as primitive
 220     TNaming_Builder aBuilder(aShapeLab);
 221     if (!theOldShape || !theNewShape)
 222       return;  // bad shape
 223     TopoDS_Shape aShapeOld = theOldShape->impl<TopoDS_Shape>();
 224     if (aShapeOld.IsNull())
 225       return;  // null shape inside
 226     TopoDS_Shape aShapeNew = theNewShape->impl<TopoDS_Shape>();
 227     if (aShapeNew.IsNull())
 228       return;  // null shape inside
 229     aBuilder.Modify(aShapeOld, aShapeNew);
 230     if(!aBuilder.NamedShape()->IsEmpty()) {
 231       Handle(TDataStd_Name) anAttr;
 232       if(aBuilder.NamedShape()->Label().FindAttribute(TDataStd_Name::GetID(),anAttr)) {
 233         std::string aName (TCollection_AsciiString(anAttr->Get()).ToCString());
 234         if(!aName.empty()) {
 235           std::shared_ptr<Model_Document> aDoc =
 236             std::dynamic_pointer_cast<Model_Document>(document());
 237           aDoc->addNamingName(aBuilder.NamedShape()->Label(), aName);
 238         }
 239       }
 240     }
 241   }
 242 }
 243
 244 void  Model_BodyBuilder::storeWithoutNaming(const std::shared_ptr<GeomAPI_Shape>& theShape)
 245 {
 246   std::shared_ptr<Model_Data> aData = std::dynamic_pointer_cast<Model_Data>(data());
 247   if (aData) {
 248     clean();
 249     if (!theShape.get())
 250       return; // bad shape
 251     TopoDS_Shape aShape = theShape->impl<TopoDS_Shape>();
 252     if (aShape.IsNull())
 253       return;  // null shape inside
 254     TNaming_Builder aBuilder(aData->shapeLab());
 255     aBuilder.Select(aShape, aShape);
 256   }
 257 }
 258
 259 void Model_BodyBuilder::clean()
 260 {
 261   TDF_Label aLab = std::dynamic_pointer_cast<Model_Data>(data())->shapeLab();
 262   if (aLab.IsNull())
 263     return;
 264   std::map<int, TNaming_Builder*>::iterator aBuilder = myBuilders.begin();
 265   for(; aBuilder != myBuilders.end(); aBuilder++) {
 266     delete aBuilder->second;
 267     // clear also shapes on cleaned sub-labels (#2241)
 268     Handle(TNaming_NamedShape) aNS;
 269     if (aLab.FindChild(aBuilder->first).FindAttribute(TNaming_NamedShape::GetID(), aNS)) {
 270       aNS->Clear();
 271     }
 272   }
 273   myBuilders.clear();
 274   // remove the old reference (if any)
 275   aLab.ForgetAttribute(TDF_Reference::GetID());
 276 }
 277
 278 Model_BodyBuilder::~Model_BodyBuilder()
 279 {
 280   clean();
 281 }
 282
 283 TNaming_Builder* Model_BodyBuilder::builder(const int theTag)
 284 {
 285   std::map<int, TNaming_Builder*>::iterator aFind = myBuilders.find(theTag);
 286   if (aFind == myBuilders.end()) {
 287     std::shared_ptr<Model_Data> aData = std::dynamic_pointer_cast<Model_Data>(data());
 288     myBuilders[theTag] = new TNaming_Builder(
 289       theTag == 0 ? aData->shapeLab() : aData->shapeLab().FindChild(theTag));
 290     aFind = myBuilders.find(theTag);
 291   }
 292   return aFind->second;
 293 }
 294
 295 void Model_BodyBuilder::buildName(const int theTag, const std::string& theName)
 296 {
 297   std::shared_ptr<Model_Document> aDoc = std::dynamic_pointer_cast<Model_Document>(document());
 298   //aDoc->addNamingName(builder(theTag)->NamedShape()->Label(), theName);
 299   TDataStd_Name::Set(builder(theTag)->NamedShape()->Label(), theName.c_str());
 300 }
 301 void Model_BodyBuilder::generated(
 302   const std::shared_ptr<GeomAPI_Shape>& theNewShape, const std::string& theName, const int theTag)
 303 {
 304   TopoDS_Shape aShape = theNewShape->impl<TopoDS_Shape>();
 305   builder(theTag)->Generated(aShape);
 306   if(!theName.empty())
 307     buildName(theTag, theName);
 308 }
 309
 310 void Model_BodyBuilder::generated(const std::shared_ptr<GeomAPI_Shape>& theOldShape,
 311   const std::shared_ptr<GeomAPI_Shape>& theNewShape, const std::string& theName, const int theTag)
 312 {
 313   TopoDS_Shape anOldShape = theOldShape->impl<TopoDS_Shape>();
 314   TopoDS_Shape aNewShape = theNewShape->impl<TopoDS_Shape>();
 315   builder(theTag)->Generated(anOldShape, aNewShape);
 316   if(!theName.empty())
 317     buildName(theTag, theName);
 318   TopAbs_ShapeEnum aGenShapeType = aNewShape.ShapeType();
 319   if(aGenShapeType == TopAbs_WIRE || aGenShapeType == TopAbs_SHELL) {
 320     TopAbs_ShapeEnum anExplodeShapeType = aGenShapeType == TopAbs_WIRE ? TopAbs_EDGE : TopAbs_FACE;
 321     const TDF_Label aLabel = builder(theTag)->NamedShape()->Label();
 322     int aTag = 1;
 323     std::shared_ptr<Model_Document> aDoc = std::dynamic_pointer_cast<Model_Document>(document());
 324     for(TopExp_Explorer anExp(aNewShape, anExplodeShapeType); anExp.More(); anExp.Next()) {
 325       TDF_Label aChildLabel = aLabel.FindChild(aTag);
 326       TNaming_Builder aBuilder(aChildLabel);
 327       aBuilder.Generated(anOldShape, anExp.Current());
 328       TCollection_AsciiString aChildName = TCollection_AsciiString((theName + "_").c_str()) + aTag;
 329       //aDoc->addNamingName(aChildLabel, aChildName.ToCString());
 330       TDataStd_Name::Set(aChildLabel, aChildName.ToCString());
 331       aTag++;
 332     }
 333   }
 334 }
 335
 336
 337 void Model_BodyBuilder::modified(const std::shared_ptr<GeomAPI_Shape>& theOldShape,
 338   const std::shared_ptr<GeomAPI_Shape>& theNewShape, const std::string& theName, const int theTag)
 339 {
 340   TopoDS_Shape anOldShape = theOldShape->impl<TopoDS_Shape>();
 341   TopoDS_Shape aNewShape = theNewShape->impl<TopoDS_Shape>();
 342   builder(theTag)->Modify(anOldShape, aNewShape);
 343   if(!theName.empty())
 344     buildName(theTag, theName);
 345 }
 346
 347 void Model_BodyBuilder::deleted(const std::shared_ptr<GeomAPI_Shape>& theOldShape,
 348   const int theTag)
 349 {
 350   TopoDS_Shape aShape = theOldShape->impl<TopoDS_Shape>();
 351   builder(theTag)->Delete(aShape);
 352 }
 353
 354 void Model_BodyBuilder::loadDeletedShapes (GeomAlgoAPI_MakeShape* theMS,
 355   std::shared_ptr<GeomAPI_Shape>  theShapeIn,
 356   const int  theKindOfShape,
 357   const int  theTag)
 358 {
 359   TopoDS_Shape aShapeIn = theShapeIn->impl<TopoDS_Shape>();
 360   TopTools_MapOfShape aView;
 361   TopExp_Explorer ShapeExplorer (aShapeIn, (TopAbs_ShapeEnum)theKindOfShape);
 362   for (; ShapeExplorer.More(); ShapeExplorer.Next ()) {
 363     const TopoDS_Shape& aRoot = ShapeExplorer.Current ();
 364     if (!aView.Add(aRoot)) continue;
 365     std::shared_ptr<GeomAPI_Shape> aRShape(new GeomAPI_Shape());
 366     aRShape->setImpl((new TopoDS_Shape(aRoot)));
 367     if (theMS->isDeleted (aRShape)) {
 368       builder(theTag)->Delete(aRoot);
 369     }
 370   }
 371 }
 372
 373 // Keep only the shapes with minimal shape type
 374 static void keepTopLevelShapes(ListOfShape& theShapes)
 375 {
 376   GeomAPI_Shape::ShapeType aKeepShapeType = GeomAPI_Shape::SHAPE;
 377   ListOfShape::iterator anIt = theShapes.begin();
 378   while (anIt != theShapes.end()) {
 379     GeomAPI_Shape::ShapeType aType = (*anIt)->shapeType();
 380     if (aType < aKeepShapeType) {
 381       // found a shape with lesser shape type => remove all previous shapes
 382       aKeepShapeType = aType;
 383       theShapes.erase(theShapes.begin(), anIt);
 384       ++anIt;
 385     } else if (aType > aKeepShapeType) {
 386       // shapes with greater shape type should be removed from the list
 387       ListOfShape::iterator aRemoveIt = anIt++;
 388       theShapes.erase(aRemoveIt);
 389     } else
 390       ++anIt;
 391   }
 392 }
 393
 394 void Model_BodyBuilder::loadAndOrientModifiedShapes (
 395   GeomAlgoAPI_MakeShape* theMS,
 396   std::shared_ptr<GeomAPI_Shape>  theShapeIn,
 397   const int  theKindOfShape,
 398   const int  theTag,
 399   const std::string& theName,
 400   GeomAPI_DataMapOfShapeShape& theSubShapes,
 401   const bool theIsStoreSeparate,
 402   const bool theIsStoreAsGenerated)
 403 {
 404   int anIndex = 1;
 405   int aTag = theTag;
 406   bool isBuilt = !theName.empty();
 407   std::string aName = theName;
 408   std::ostringstream aStream;
 409   GeomShapePtr aResultShape = shape();
 410   TopoDS_Shape aShapeIn = theShapeIn->impl<TopoDS_Shape>();
 411   TopTools_MapOfShape aView;
 412   TopExp_Explorer aShapeExplorer (aShapeIn, (TopAbs_ShapeEnum)theKindOfShape);
 413   for (; aShapeExplorer.More(); aShapeExplorer.Next ()) {
 414     const TopoDS_Shape& aRoot = aShapeExplorer.Current ();
 415     if (!aView.Add(aRoot)) continue;
 416     if (TNaming_Tool::NamedShape(aRoot, builder(theTag)->NamedShape()->Label()).IsNull())
 417       continue; // there is no sense to write history if old shape does not exist in the document
 418     ListOfShape aList;
 419     std::shared_ptr<GeomAPI_Shape> aRShape(new GeomAPI_Shape());
 420     aRShape->setImpl((new TopoDS_Shape(aRoot)));
 421     theMS->modified(aRShape, aList);
 422     keepTopLevelShapes(aList);
 423     // sort the list of images before naming
 424     GeomAlgoAPI_SortListOfShapes::sort(aList);
 425
 426     // to trace situation where several objects are produced by one parent (#2317)
 427     int aSameParentShapes = -1;
 428     std::list<std::shared_ptr<GeomAPI_Shape> >::const_iterator
 429       anIt = aList.begin(), aLast = aList.end();
 430     for (; anIt != aLast; anIt++) {
 431       TopoDS_Shape aNewShape = (*anIt)->impl<TopoDS_Shape>();
 432       if (theSubShapes.isBound(*anIt)) {
 433         std::shared_ptr<GeomAPI_Shape> aMapShape(theSubShapes.find(*anIt));
 434         aNewShape.Orientation(aMapShape->impl<TopoDS_Shape>().Orientation());
 435       }
 436       GeomShapePtr aGeomNewShape(new GeomAPI_Shape());
 437       aGeomNewShape->setImpl(new TopoDS_Shape(aNewShape));
 438       if(!aRoot.IsSame(aNewShape) && aResultShape->isSubShape(aGeomNewShape, false) &&
 439          !aResultShape->isSame(*anIt)) { // to avoid put of same shape on main label and sub
 440         int aBuilderTag = aTag;
 441         if (!theIsStoreSeparate)
 442           aSameParentShapes++;
 443         if (aSameParentShapes > 0) { // store in other label
 444           aBuilderTag = 100000 - aSameParentShapes * 10 - aTag;
 445         }
 446         if(theIsStoreAsGenerated) {
 447           // Here we store shapes as generated, to avoid problem when one parent shape produce
 448           // several child shapes. In this case naming could not determine which shape to select.
 449           builder(aBuilderTag)->Generated(aRoot, aNewShape);
 450         } else {
 451           builder(aBuilderTag)->Modify(aRoot, aNewShape);
 452         }
 453         if(isBuilt) {
 454           if(theIsStoreSeparate) {
 455             aStream.str(std::string());
 456             aStream.clear();
 457             aStream << theName << "_" << anIndex++;
 458             aName = aStream.str();
 459           }
 460           if (aSameParentShapes > 0) {
 461             aStream.str(std::string());
 462             aStream.clear();
 463             aStream << aName << "_" << aSameParentShapes << "divided";
 464             std::string aNameDiv = aStream.str();
 465             buildName(aBuilderTag, aNameDiv);
 466           } else {
 467             buildName(aBuilderTag, aName);
 468           }
 469         }
 470         if(theIsStoreSeparate) {
 471           aTag++;
 472         }
 473       } else if (aResultShape->isSame(*anIt)) {
 474         // keep the modification evolution on the root level (2241 - history propagation issue)
 475         if(theIsStoreAsGenerated) {
 476           builder(0)->Generated(aRoot, aNewShape);
 477         } else {
 478           builder(0)->Modify(aRoot, aNewShape);
 479         }
 480       }
 481     }
 482   }
 483 }
 484
 485 void Model_BodyBuilder::loadAndOrientGeneratedShapes (
 486   GeomAlgoAPI_MakeShape* theMS,
 487   std::shared_ptr<GeomAPI_Shape>  theShapeIn,
 488   const int  theKindOfShape,
 489   const int  theTag,
 490   const std::string& theName,
 491   GeomAPI_DataMapOfShapeShape& theSubShapes)
 492 {
 493   TopoDS_Shape aShapeIn = theShapeIn->impl<TopoDS_Shape>();
 494   TopTools_MapOfShape aView;
 495   bool isBuilt = !theName.empty();
 496   TopExp_Explorer aShapeExplorer (aShapeIn, (TopAbs_ShapeEnum)theKindOfShape);
 497   for (; aShapeExplorer.More(); aShapeExplorer.Next ()) {
 498     const TopoDS_Shape& aRoot = aShapeExplorer.Current ();
 499     if (!aView.Add(aRoot)) continue;
 500     //if (TNaming_Tool::NamedShape(aRoot, builder(theTag)->NamedShape()->Label()).IsNull())
 501     //  continue; // there is no sense to write history if old shape does not exist in the document
 502     ListOfShape aList;
 503     std::shared_ptr<GeomAPI_Shape> aRShape(new GeomAPI_Shape());
 504     aRShape->setImpl((new TopoDS_Shape(aRoot)));
 505     theMS->generated(aRShape, aList);
 506     keepTopLevelShapes(aList);
 507     std::list<std::shared_ptr<GeomAPI_Shape> >::const_iterator
 508       anIt = aList.begin(), aLast = aList.end();
 509     for (; anIt != aLast; anIt++) {
 510       TopoDS_Shape aNewShape = (*anIt)->impl<TopoDS_Shape>();
 511       if (theSubShapes.isBound(*anIt)) {
 512         std::shared_ptr<GeomAPI_Shape> aMapShape(theSubShapes.find(*anIt));
 513         aNewShape.Orientation(aMapShape->impl<TopoDS_Shape>().Orientation());
 514       }
 515       if (!aRoot.IsSame (aNewShape)) {
 516         builder(theTag)->Generated(aRoot,aNewShape);
 517         if(isBuilt)
 518           buildName(theTag, theName);
 519       }
 520       TopAbs_ShapeEnum aGenShapeType = aNewShape.ShapeType();
 521       if(aGenShapeType == TopAbs_WIRE || aGenShapeType == TopAbs_SHELL) {
 522         TopAbs_ShapeEnum anExplodeShapeType =
 523           aGenShapeType == TopAbs_WIRE ? TopAbs_EDGE : TopAbs_FACE;
 524         const TDF_Label aLabel = builder(theTag)->NamedShape()->Label();
 525         int aTag = 1;
 526         std::shared_ptr<Model_Document> aDoc =
 527           std::dynamic_pointer_cast<Model_Document>(document());
 528         for(TopExp_Explorer anExp(aNewShape, anExplodeShapeType); anExp.More(); anExp.Next()) {
 529           TDF_Label aChildLabel = aLabel.FindChild(aTag);
 530           TNaming_Builder aBuilder(aChildLabel);
 531           aBuilder.Generated(aRoot, anExp.Current());
 532           TCollection_AsciiString aChildName =
 533             TCollection_AsciiString((theName + "_").c_str()) + aTag;
 534           TDataStd_Name::Set(aChildLabel, aChildName.ToCString());
 535           aTag++;
 536         }
 537       }
 538     }
 539   }
 540 }
 541
 542 //=======================================================================
 543 int getDangleShapes(const TopoDS_Shape&           theShapeIn,
 544   const TopAbs_ShapeEnum        theGeneratedFrom,
 545   TopTools_DataMapOfShapeShape& theDangles)
 546 {
 547   theDangles.Clear();
 548   TopTools_IndexedDataMapOfShapeListOfShape subShapeAndAncestors;
 549   TopAbs_ShapeEnum GeneratedTo;
 550   if (theGeneratedFrom == TopAbs_FACE) GeneratedTo = TopAbs_EDGE;
 551   else if (theGeneratedFrom == TopAbs_EDGE) GeneratedTo = TopAbs_VERTEX;
 552   else return Standard_False;
 553   TopExp::MapShapesAndAncestors(theShapeIn, GeneratedTo, theGeneratedFrom, subShapeAndAncestors);
 554   for (Standard_Integer i = 1; i <= subShapeAndAncestors.Extent(); i++) {
 555     const TopoDS_Shape& mayBeDangle = subShapeAndAncestors.FindKey(i);
 556     const TopTools_ListOfShape& ancestors = subShapeAndAncestors.FindFromIndex(i);
 557     if (ancestors.Extent() == 1) theDangles.Bind(ancestors.First(), mayBeDangle);
 558   }
 559   return theDangles.Extent();
 560 }
 561
 562 //=======================================================================
 563 void loadGeneratedDangleShapes(
 564   const TopoDS_Shape&      theShapeIn,
 565   const TopAbs_ShapeEnum   theGeneratedFrom,
 566   TNaming_Builder *        theBuilder)
 567 {
 568   TopTools_DataMapOfShapeShape dangles;
 569   if (!getDangleShapes(theShapeIn, theGeneratedFrom, dangles)) return;
 570   TopTools_DataMapIteratorOfDataMapOfShapeShape itr(dangles);
 571   for (; itr.More(); itr.Next())
 572     theBuilder->Generated(itr.Key(), itr.Value());
 573 }
 574
 575 //=======================================================================
 576 void Model_BodyBuilder::loadNextLevels(std::shared_ptr<GeomAPI_Shape> theShape,
 577   const std::string& theName, int&  theTag)
 578 {
 579   if(theShape->isNull()) return;
 580   TopoDS_Shape aShape = theShape->impl<TopoDS_Shape>();
 581   std::string aName;
 582   if (aShape.ShapeType() == TopAbs_SOLID) {
 583     TopExp_Explorer expl(aShape, TopAbs_FACE);
 584     for (; expl.More(); expl.Next()) {
 585       builder(theTag)->Generated(expl.Current());
 586       TCollection_AsciiString aStr(theTag);
 587       aName = theName + aStr.ToCString();
 588       buildName(theTag, aName);
 589       theTag++;
 590     }
 591   }
 592   else if (aShape.ShapeType() == TopAbs_SHELL || aShape.ShapeType() == TopAbs_FACE) {
 593     // load faces and all the free edges
 594     TopTools_IndexedMapOfShape Faces;
 595     TopExp::MapShapes(aShape, TopAbs_FACE, Faces);
 596     if (Faces.Extent() > 1 || (aShape.ShapeType() == TopAbs_SHELL && Faces.Extent() == 1)) {
 597       TopExp_Explorer expl(aShape, TopAbs_FACE);
 598       for (; expl.More(); expl.Next()) {
 599         builder(theTag)->Generated(expl.Current());
 600         TCollection_AsciiString aStr(theTag);
 601         aName = theName + aStr.ToCString();
 602         buildName(theTag, aName);
 603         theTag++;
 604       }
 605     }
 606     TopTools_IndexedDataMapOfShapeListOfShape anEdgeAndNeighbourFaces;
 607     TopExp::MapShapesAndAncestors(aShape, TopAbs_EDGE, TopAbs_FACE, anEdgeAndNeighbourFaces);
 608     for (Standard_Integer i = 1; i <= anEdgeAndNeighbourFaces.Extent(); i++)
 609     {
 610       const TopTools_ListOfShape& aLL = anEdgeAndNeighbourFaces.FindFromIndex(i);
 611       if (aLL.Extent() < 2) {
 612         if (BRep_Tool::Degenerated(TopoDS::Edge(anEdgeAndNeighbourFaces.FindKey(i))))
 613           continue;
 614         builder(theTag)->Generated(anEdgeAndNeighbourFaces.FindKey(i));
 615         TCollection_AsciiString aStr(theTag);
 616         aName = theName + aStr.ToCString();
 617         buildName(theTag, aName);
 618         theTag++;
 619       } else {
 620         TopTools_ListIteratorOfListOfShape anIter(aLL);
 621         const TopoDS_Face& aFace = TopoDS::Face(anIter.Value());
 622         anIter.Next();
 623         if(aFace.IsEqual(anIter.Value())) {
 624           builder(theTag)->Generated(anEdgeAndNeighbourFaces.FindKey(i));
 625           TCollection_AsciiString aStr(theTag);
 626           aName = theName + aStr.ToCString();
 627           buildName(theTag, aName);
 628           theTag++;
 629         }
 630       }
 631     }
 632   } else if (aShape.ShapeType() == TopAbs_WIRE) {
 633     TopTools_IndexedMapOfShape Edges;
 634     BRepTools::Map3DEdges(aShape, Edges);
 635     if (Edges.Extent() == 1) {
 636       builder(theTag++)->Generated(Edges.FindKey(1));
 637       TopExp_Explorer expl(aShape, TopAbs_VERTEX);
 638       for (; expl.More(); expl.Next()) {
 639         builder(theTag)->Generated(expl.Current());
 640         TCollection_AsciiString aStr(theTag);
 641         aName = theName + aStr.ToCString();
 642         buildName(theTag, aName);
 643         theTag++;
 644       }
 645     } else {
 646       TopExp_Explorer expl(aShape, TopAbs_EDGE);
 647       for (; expl.More(); expl.Next()) {
 648         builder(theTag)->Generated(expl.Current());
 649         TCollection_AsciiString aStr(theTag);
 650         aName = theName + aStr.ToCString();
 651         buildName(theTag, aName);
 652         theTag++;
 653       }
 654       // and load generated vertices.
 655       TopTools_DataMapOfShapeShape generated;
 656       if (getDangleShapes(aShape, TopAbs_EDGE, generated))
 657       {
 658         TNaming_Builder* pBuilder = builder(theTag++);
 659         loadGeneratedDangleShapes(aShape, TopAbs_EDGE, pBuilder);
 660       }
 661     }
 662   } else if (aShape.ShapeType() == TopAbs_EDGE) {
 663     TopExp_Explorer expl(aShape, TopAbs_VERTEX);
 664     for (; expl.More(); expl.Next()) {
 665       builder(theTag)->Generated(expl.Current());
 666       TCollection_AsciiString aStr(theTag);
 667       aName = theName + aStr.ToCString();
 668       buildName(theTag, aName);
 669       theTag++;
 670     }
 671   }
 672 }
 673
 674 //=======================================================================
 675 int findAmbiguities(const TopoDS_Shape&           theShapeIn,
 676   TopTools_ListOfShape&   theList)
 677 {
 678   theList.Clear();
 679   // edges -> ancestor faces list
 680   TopTools_IndexedDataMapOfShapeListOfShape aSubShapeAndAncestors;
 681   TopExp::MapShapesAndAncestors(theShapeIn, TopAbs_EDGE, TopAbs_FACE, aSubShapeAndAncestors);
 682   // keeps the shapes which are already in the resulting list
 683   TopTools_MapOfShape alreadyThere;
 684   // map from faces identifier (combination of hash-codes) to list of edges produced such ID
 685   NCollection_DataMap<int, NCollection_List<TopoDS_Shape> > aFacesIDs;
 686
 687   TopTools_IndexedDataMapOfShapeListOfShape::Iterator anAncestorsIter(aSubShapeAndAncestors);
 688   for (; anAncestorsIter.More(); anAncestorsIter.Next()) {
 689     const TopTools_ListOfShape& ancestors = anAncestorsIter.Value();
 690     if (ancestors.Extent() < 2)
 691       continue;
 692     Standard_Integer anID = 0;
 693     for(TopTools_ListIteratorOfListOfShape aFaceIt(ancestors); aFaceIt.More(); aFaceIt.Next()) {
 694       anID ^= HashCode(aFaceIt.ChangeValue(), 1990657); // Pierpont prime
 695     }
 696     if (aFacesIDs.IsBound(anID)) { // there found same edge, check they really have same faces
 697       const NCollection_List<TopoDS_Shape>& aSameFaces1 =
 698         aSubShapeAndAncestors.FindFromKey(anAncestorsIter.Key());
 699       NCollection_List<TopoDS_Shape>::Iterator aSameEdge(aFacesIDs.ChangeFind(anID));
 700       for(; aSameEdge.More(); aSameEdge.Next()) {
 701         const NCollection_List<TopoDS_Shape>& aSameFaces2 =
 702           aSubShapeAndAncestors.FindFromKey(aSameEdge.Value());
 703         if (aSameFaces2.Extent() != aSameFaces1.Extent()) // the number of faces is different
 704           break;
 705
 706         NCollection_List<TopoDS_Shape>::Iterator aFaceIter1(aSameFaces1);
 707         for(; aFaceIter1.More(); aFaceIter1.Next()) {
 708           NCollection_List<TopoDS_Shape>::Iterator aFaceIter2(aSameFaces2);
 709           for(; aFaceIter2.More(); aFaceIter2.Next()) {
 710             if (aFaceIter1.Value().IsSame(aFaceIter2.Value()))
 711               break;
 712           }
 713           if (!aFaceIter2.More()) // aFaceIter1 contains a face, which is not in aFaceIter2
 714             break;
 715         }
 716         if (!aFaceIter1.More()) { // all the faces are same => put to the result
 717           if (alreadyThere.Add(aSameEdge.Value()))
 718             theList.Append(aSameEdge.Value());
 719           if (alreadyThere.Add(anAncestorsIter.Key()))
 720             theList.Append(anAncestorsIter.Key());
 721         }
 722       }
 723     } else { // ID is unique, just add this edge
 724       aFacesIDs.Bind(anID, NCollection_List<TopoDS_Shape>());
 725     }
 726     aFacesIDs.ChangeFind(anID).Append(anAncestorsIter.Key()); // add to the list anyway
 727   }
 728   return theList.Extent();
 729 }
 730
 731 //=======================================================================
 732 void Model_BodyBuilder::loadFirstLevel(
 733   std::shared_ptr<GeomAPI_Shape> theShape, const std::string& theName, int&  theTag)
 734 {
 735   if(theShape->isNull()) return;
 736   TopoDS_Shape aShape = theShape->impl<TopoDS_Shape>();
 737   std::string aName;
 738   if (aShape.ShapeType() == TopAbs_COMPOUND || aShape.ShapeType() == TopAbs_COMPSOLID) {
 739     TopoDS_Iterator itr(aShape);
 740     for (; itr.More(); itr.Next()) {
 741       builder(theTag)->Generated(itr.Value());
 742       TCollection_AsciiString aStr(theTag);
 743       aName = theName + aStr.ToCString();
 744       buildName(theTag, aName);
 745       if(!theName.empty()) buildName(theTag, aName);
 746       theTag++;
 747       if (itr.Value().ShapeType() == TopAbs_COMPOUND ||
 748         itr.Value().ShapeType() == TopAbs_COMPSOLID)
 749       {
 750         std::shared_ptr<GeomAPI_Shape> itrShape(new GeomAPI_Shape());
 751         itrShape->setImpl(new TopoDS_Shape(itr.Value()));
 752         loadFirstLevel(itrShape, theName, theTag);
 753       } else {
 754         std::shared_ptr<GeomAPI_Shape> itrShape(new GeomAPI_Shape());
 755         itrShape->setImpl(new TopoDS_Shape(itr.Value()));
 756         loadNextLevels(itrShape, theName, theTag);
 757       }
 758     }
 759   } else {
 760     std::shared_ptr<GeomAPI_Shape> itrShape(new GeomAPI_Shape());
 761     itrShape->setImpl(new TopoDS_Shape(aShape));
 762     loadNextLevels(itrShape, theName, theTag);
 763   }
 764   TopTools_ListOfShape   aList;
 765   if(findAmbiguities(aShape, aList)) {
 766     TopTools_ListIteratorOfListOfShape it(aList);
 767     for (; it.More(); it.Next(),theTag++) {
 768       builder(theTag)->Generated(it.Value());
 769       TCollection_AsciiString aStr(theTag);
 770       aName = theName + aStr.ToCString();
 771       buildName(theTag, aName);
 772     }
 773   }
 774 }
 775
 776 //=======================================================================
 777 void Model_BodyBuilder::loadDisconnectedEdges(
 778   std::shared_ptr<GeomAPI_Shape> theShape, const std::string& theName, int&  theTag)
 779 {
 780   if(theShape->isNull()) return;
 781   TopoDS_Shape aShape = theShape->impl<TopoDS_Shape>();
 782   TopTools_DataMapOfShapeListOfShape edgeNaborFaces;
 783   TopTools_ListOfShape empty;
 784   TopExp_Explorer explF(aShape, TopAbs_FACE);
 785   for (; explF.More(); explF.Next()) {
 786     const TopoDS_Shape& aFace = explF.Current();
 787     TopExp_Explorer explV(aFace, TopAbs_EDGE);
 788     for (; explV.More(); explV.Next()) {
 789       const TopoDS_Shape& anEdge = explV.Current();
 790       if (!edgeNaborFaces.IsBound(anEdge)) edgeNaborFaces.Bind(anEdge, empty);
 791       Standard_Boolean faceIsNew = Standard_True;
 792       TopTools_ListIteratorOfListOfShape itrF(edgeNaborFaces.Find(anEdge));
 793       for (; itrF.More(); itrF.Next()) {
 794         if (itrF.Value().IsSame(aFace)) {
 795           faceIsNew = Standard_False;
 796           break;
 797         }
 798       }
 799       if (faceIsNew)
 800         edgeNaborFaces.ChangeFind(anEdge).Append(aFace);
 801     }
 802   }
 803
 804   TopTools_MapOfShape anEdgesToDelete;
 805   TopExp_Explorer anEx(aShape,TopAbs_EDGE);
 806   std::string aName;
 807   for(;anEx.More();anEx.Next()) {
 808     Standard_Boolean aC0 = Standard_False;
 809     TopoDS_Shape anEdge1 = anEx.Current();
 810     if (edgeNaborFaces.IsBound(anEdge1)) {
 811       const TopTools_ListOfShape& aList1 = edgeNaborFaces.Find(anEdge1);
 812       if (aList1.Extent()<2) continue;
 813       TopTools_DataMapIteratorOfDataMapOfShapeListOfShape itr(edgeNaborFaces);
 814       for (; itr.More(); itr.Next()) {
 815         TopoDS_Shape anEdge2 = itr.Key();
 816         if(anEdgesToDelete.Contains(anEdge2)) continue;
 817         if (anEdge1.IsSame(anEdge2)) continue;
 818         const TopTools_ListOfShape& aList2 = itr.Value();
 819         // compare lists of the neighbour faces of edge1 and edge2
 820         if (aList1.Extent() == aList2.Extent()) {
 821           Standard_Integer aMatches = 0;
 822           for(TopTools_ListIteratorOfListOfShape aLIter1(aList1);aLIter1.More();aLIter1.Next())
 823             for(TopTools_ListIteratorOfListOfShape aLIter2(aList2);aLIter2.More();aLIter2.Next())
 824               if (aLIter1.Value().IsSame(aLIter2.Value())) aMatches++;
 825           if (aMatches == aList1.Extent()) {
 826             aC0=Standard_True;
 827             builder(theTag)->Generated(anEdge2);
 828             anEdgesToDelete.Add(anEdge2);
 829             TCollection_AsciiString aStr(theTag);
 830             aName = theName + aStr.ToCString();
 831             buildName(theTag, aName);
 832             theTag++;
 833           }
 834         }
 835       }
 836       TopTools_MapIteratorOfMapOfShape itDelete(anEdgesToDelete);
 837       for(;itDelete.More();itDelete.Next())
 838         edgeNaborFaces.UnBind(itDelete.Key());
 839       edgeNaborFaces.UnBind(anEdge1);
 840     }
 841     if (aC0) {
 842       builder(theTag)->Generated(anEdge1);
 843       TCollection_AsciiString aStr(theTag);
 844       aName = theName + aStr.ToCString();
 845       buildName(theTag, aName);
 846       theTag++;
 847     }
 848   }
 849 }
 850
 851 void Model_BodyBuilder::loadDisconnectedVertexes(std::shared_ptr<GeomAPI_Shape> theShape,
 852                                                  const std::string& theName, int&  theTag)
 853 {
 854   if(theShape->isNull()) return;
 855   TopoDS_Shape aShape = theShape->impl<TopoDS_Shape>();
 856   TopTools_DataMapOfShapeListOfShape vertexNaborEdges;
 857   TopTools_ListOfShape empty;
 858   TopExp_Explorer explF(aShape, TopAbs_EDGE);
 859   for (; explF.More(); explF.Next()) {
 860     const TopoDS_Shape& anEdge = explF.Current();
 861     TopExp_Explorer explV(anEdge, TopAbs_VERTEX);
 862     for (; explV.More(); explV.Next()) {
 863       const TopoDS_Shape& aVertex = explV.Current();
 864       if (!vertexNaborEdges.IsBound(aVertex)) vertexNaborEdges.Bind(aVertex, empty);
 865       Standard_Boolean faceIsNew = Standard_True;
 866       TopTools_ListIteratorOfListOfShape itrF(vertexNaborEdges.Find(aVertex));
 867       for (; itrF.More(); itrF.Next()) {
 868         if (itrF.Value().IsSame(anEdge)) {
 869           faceIsNew = Standard_False;
 870           break;
 871         }
 872       }
 873       if (faceIsNew) {
 874         vertexNaborEdges.ChangeFind(aVertex).Append(anEdge);
 875       }
 876     }
 877   }
 878   std::string aName;
 879   TopTools_DataMapIteratorOfDataMapOfShapeListOfShape itr(vertexNaborEdges);
 880   for (; itr.More(); itr.Next()) {
 881     const TopTools_ListOfShape& naborEdges = itr.Value();
 882     if (naborEdges.Extent() < 2) {
 883       builder(theTag)->Generated(itr.Key());
 884       TCollection_AsciiString aStr(theTag);
 885       aName = theName + aStr.ToCString();
 886       buildName(theTag, aName);
 887       theTag++;
 888     }
 889   }
 890 }
 891
 892 std::shared_ptr<GeomAPI_Shape> Model_BodyBuilder::shape()
 893 {
 894   std::shared_ptr<Model_Data> aData = std::dynamic_pointer_cast<Model_Data>(data());
 895   if (aData) {
 896     TDF_Label aShapeLab = aData->shapeLab();
 897     Handle(TDF_Reference) aRef;
 898     if (aShapeLab.FindAttribute(TDF_Reference::GetID(), aRef)) {
 899       aShapeLab = aRef->Get();
 900     }
 901     Handle(TNaming_NamedShape) aName;
 902     if (aShapeLab.FindAttribute(TNaming_NamedShape::GetID(), aName)) {
 903       TopoDS_Shape aShape = aName->Get();
 904       if (!aShape.IsNull()) {
 905         std::shared_ptr<GeomAPI_Shape> aRes(new GeomAPI_Shape);
 906         aRes->setImpl(new TopoDS_Shape(aShape));
 907         return aRes;
 908       }
 909     }
 910   }
 911   return std::shared_ptr<GeomAPI_Shape>();
 912 }
 913
 914 bool Model_BodyBuilder::isLatestEqual(const std::shared_ptr<GeomAPI_Shape>& theShape)
 915 {
 916   if (theShape.get()) {
 917     TopoDS_Shape aShape = theShape->impl<TopoDS_Shape>();
 918     std::shared_ptr<Model_Data> aData = std::dynamic_pointer_cast<Model_Data>(data());
 919     if (aData) {
 920       TDF_Label& aShapeLab = aData->shapeLab();
 921       Handle(TNaming_NamedShape) aName;
 922       if (aShapeLab.FindAttribute(TNaming_NamedShape::GetID(), aName)) {
 923         TopoDS_Shape aLatest = TNaming_Tool::CurrentShape(aName);
 924         if (aLatest.IsNull())
 925           return false;
 926         if (aLatest.IsEqual(aShape))
 927           return true;
 928         // check sub-shapes for comp-solids:
 929         for (TopExp_Explorer anExp(aShape, aLatest.ShapeType()); anExp.More(); anExp.Next()) {
 930           if (aLatest.IsEqual(anExp.Current()))
 931             return true;
 932         }
 933       }
 934     }
 935   }
 936   return false;
 937 }