src/HYDROData/HYDROData_SplitToZonesTool.cxx

   1 // Copyright (C) 2007-2015  CEA/DEN, EDF R&D, OPEN CASCADE
   2 //
   3 // Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   5 //
   6 // This library is free software; you can redistribute it and/or
   7 // modify it under the terms of the GNU Lesser General Public
   8 // License as published by the Free Software Foundation; either
   9 // version 2.1 of the License, or (at your option) any later version.
  10 //
  11 // This library is distributed in the hope that it will be useful,
  12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14 // Lesser General Public License for more details.
  15 //
  16 // You should have received a copy of the GNU Lesser General Public
  17 // License along with this library; if not, write to the Free Software
  18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  19 //
  20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  21 //
  22
  23 #include "HYDROData_SplitToZonesTool.h"
  24
  25 #include "HYDROData_PolylineXY.h"
  26 #include "HYDROData_ShapesGroup.h"
  27 #include <HYDROData_Transform.h>
  28 #include <BRepAlgoAPI_Cut.hxx>
  29 #include <BRepAlgoAPI_Common.hxx>
  30 #include <BRepBuilderAPI_MakeFace.hxx>
  31 #include <BRep_Builder.hxx>
  32 #include <TopExp_Explorer.hxx>
  33 #include <TopoDS.hxx>
  34 #include <TopoDS_Shape.hxx>
  35 #include <TopoDS_Compound.hxx>
  36 #include <TopoDS_Wire.hxx>
  37 #include <TopoDS_Edge.hxx>
  38 #include <BRepCheck_Analyzer.hxx>
  39 #include <TopTools_ListOfShape.hxx>
  40 #include <TopTools_ListIteratorOfListOfShape.hxx>
  41 #include <gp_Pln.hxx>
  42 #include <BRepGProp.hxx>
  43 #include <GProp_GProps.hxx>
  44 #include <Geom_Plane.hxx>
  45 #include <BRepBuilderAPI_FindPlane.hxx>
  46
  47 //#define DEB_SPLIT_TO_ZONES 1
  48 //#define DEB_SPLIT_TO_ZONES_CHECK_PARTITION 1
  49 #if (defined (DEB_SPLIT_TO_ZONES) || defined(DEB_SPLIT_TO_ZONES_CHECK_PARTITION))
  50 #include <BRepTools.hxx>
  51 static TCollection_AsciiString fileNameBefore("BeforeTranslation");
  52 #endif
  53
  54 TopoDS_Face HYDROData_SplitToZonesTool::SplitData::Face() const
  55 {
  56   TopoDS_Face aResFace;
  57
  58   if( !Shape.IsNull() )
  59   {
  60     if ( Shape.ShapeType() == TopAbs_FACE )
  61     {
  62       aResFace = TopoDS::Face( Shape );
  63     }
  64     else if ( Shape.ShapeType() == TopAbs_WIRE )
  65     {
  66       BRepBuilderAPI_MakeFace aMakeFace( TopoDS::Wire( Shape ), Standard_True );
  67       aMakeFace.Build();
  68       if( aMakeFace.IsDone() )
  69         aResFace = aMakeFace.Face();
  70     }
  71   }
  72
  73   return aResFace;
  74 }
  75
  76 Standard_Boolean buildLimFace(const TopoDS_Wire& theBndWire, TopoDS_Face& outFace)
  77 {
  78   GProp_GProps G;
  79   BRepGProp::LinearProperties(theBndWire,G);
  80   const gp_Pnt& aCPnt = G.CentreOfMass();
  81   gp_Pln aPln;
  82   BRepBuilderAPI_FindPlane fndPlane (theBndWire, Precision::Confusion());
  83   if(fndPlane.Found())
  84     aPln = fndPlane.Plane()->Pln();
  85   else
  86     aPln = gp_Pln(aCPnt, gp::OZ().Direction());
  87   BRepBuilderAPI_MakeFace  aMkFace(aPln, theBndWire);
  88   if(aMkFace.IsDone()) {
  89     outFace = aMkFace.Face();
  90     if(!outFace.IsNull()) {
  91 #ifdef DEB_SPLIT_TO_ZONES
  92 //      BRepTools::Write(limFace,"FL.brep");
  93 #endif
  94       return Standard_True;
  95     }
  96   }
  97   return Standard_False;
  98 }
  99 //======================================================================================
 100
 101 void HYDROData_SplitToZonesTool::SetFileNames(const QString& theNameBefore, const QString& theNameAfter)
 102 {
 103 #ifdef DEB_SPLIT_TO_ZONES_CHECK_PARTITION
 104   if(!theNameBefore.isNull() && !theNameBefore.isEmpty())
 105   {
 106     TCollection_AsciiString aName (theNameBefore.toStdString().data());
 107     fileNameBefore = aName;
 108   }
 109   if(!theNameAfter.isNull() && !theNameAfter.isEmpty())
 110   {
 111     TCollection_AsciiString aName (theNameAfter.toStdString().data());
 112     HYDROData_Transform::SetFileName (aName);
 113   }
 114 #endif
 115 }
 116 //======================================================================================
 117 Standard_Integer HYDROData_SplitToZonesTool::SplitFaces(const TopoDS_Compound& theComp,
 118                                                         HYDROData_Transform& theTool)
 119 {
 120   theTool.Clear();
 121   theTool.SetArgument(theComp);
 122   theTool.Detect();
 123   const Standard_Boolean bToTransform = theTool.ToTransform();
 124   theTool.SetToTransform(bToTransform);
 125   theTool.Perform();
 126   Standard_Integer anErr = theTool.ErrorStatus();
 127   return anErr;
 128 }
 129
 130 //======================================================================================
 131 HYDROData_SplitToZonesTool::SplitDataList
 132   HYDROData_SplitToZonesTool::Split( const HYDROData_SequenceOfObjects&  theObjectList,
 133                                      const HYDROData_SequenceOfObjects&  theGroupsList,
 134                                      const Handle(HYDROData_PolylineXY)& thePolyline )
 135 {
 136   SplitDataList anOutputSplitDataList;
 137   if(theObjectList.IsEmpty()) return anOutputSplitDataList;
 138   // Preparation.
 139   // Collect the object shapes to split. InputDataList will contain elements which will hold shape & name_of_shape.
 140   SplitDataList anInputSplitDataList;
 141 #ifdef DEB_SPLIT_TO_ZONES
 142   TCollection_AsciiString aNam("Faces_");
 143 #endif
 144   for( int anIndex = 1, aLength = theObjectList.Length(); anIndex <= aLength; anIndex++ )
 145   {
 146     Handle(HYDROData_Object) aGeomObj =
 147       Handle(HYDROData_Object)::DownCast( theObjectList.Value( anIndex ) );
 148     if( aGeomObj.IsNull() )
 149       continue;
 150
 151     TopoDS_Shape aShape = aGeomObj->GetTopShape();
 152     if ( aShape.IsNull() )
 153       continue;
 154 #ifdef DEB_SPLIT_TO_ZONES
 155    TCollection_AsciiString aName = aNam + anIndex + ".brep";
 156    BRepTools::Write(aShape, aName.ToCString());
 157 #endif
 158     if ( aShape.ShapeType() == TopAbs_COMPOUND ) {
 159       // Create split data for each face contained in the compound
 160       TopExp_Explorer anExp( aShape, TopAbs_FACE );
 161       for ( ; anExp.More(); anExp.Next() ) {
 162         const TopoDS_Face& aFace = TopoDS::Face( anExp.Current() );
 163         if ( !aFace.IsNull() ) {
 164           SplitData aSplitData( SplitData::Data_Zone, aFace, aGeomObj->GetName() );
 165           anInputSplitDataList.append( aSplitData );
 166         }
 167       }
 168     } else {
 169       SplitData aSplitData( SplitData::Data_Zone, aShape, aGeomObj->GetName() );
 170       anInputSplitDataList.append( aSplitData );
 171     }
 172   }
 173   //
 174   SplitDataList anInputGroupList;
 175   for( int anIndex = 1, aLength = theGroupsList.Length(); anIndex <=aLength; anIndex++ )
 176   {
 177     Handle(HYDROData_ShapesGroup) aGeomGroup =
 178       Handle(HYDROData_ShapesGroup)::DownCast( theGroupsList.Value( anIndex ) );
 179     if( aGeomGroup.IsNull() )
 180       continue;
 181
 182     TopTools_SequenceOfShape aGroupShapes;
 183     aGeomGroup->GetShapes( aGroupShapes );
 184     for( int i = 1, aNbShapes = aGroupShapes.Length(); i <= aNbShapes; i++ ) {
 185       const TopoDS_Shape& aGroupShape = aGroupShapes.Value( i );
 186       if ( aGroupShape.IsNull() )
 187         continue;
 188
 189       if ( aGroupShape.ShapeType() == TopAbs_COMPOUND ) {
 190         TopExp_Explorer anExp( aGroupShape, TopAbs_EDGE );
 191         for ( ; anExp.More(); anExp.Next() ) {
 192           const TopoDS_Edge& anEdge = TopoDS::Edge( anExp.Current() );
 193           if ( !anEdge.IsNull() ) {
 194             SplitData aSplitData( SplitData::Data_Edge, anEdge, aGeomGroup->GetName() );
 195             anInputGroupList.append( aSplitData );
 196           }
 197         }
 198       } else {
 199         SplitData aSplitData( SplitData::Data_Edge, aGroupShape, aGeomGroup->GetName() );
 200         anInputGroupList.append( aSplitData );
 201 #ifdef DEB_SPLIT_TO_ZONES
 202         QString aStr = aSplitData.ObjectNames.join(" ");
 203         cout << "Group# = " <<anIndex <<" Nmb in grp = " << i<< " Names = "<<aStr.toStdString() << " size = " <<aSplitData.ObjectNames.size() <<endl;
 204 #endif
 205       }
 206     }
 207   }
 208   // If only one shape is given we don't split it
 209   // algorithm just returns the unpacked input data
 210   bool limplus1Object(false);
 211   if(theObjectList.Size() == 1 )  {
 212     if(thePolyline.IsNull()) {
 213       anOutputSplitDataList.append(anInputSplitDataList);
 214       if(!theGroupsList.IsEmpty() )
 215         anOutputSplitDataList.append(anInputGroupList);
 216       return anOutputSplitDataList;
 217     } else
 218         limplus1Object = true;// size =1 && hasLimits
 219   }
 220   HYDROData_DataMapOfShapeListOfString aDM3;
 221   if(!anInputGroupList.isEmpty()) {// Old edge ==> List_Of_Names
 222     QStringList aListOfNames;
 223     for (int i=0;i < anInputGroupList.size() ;i++) {
 224       const TopoDS_Shape& aSh = anInputGroupList.at(i).Shape;
 225       aDM3.Add(aSh, anInputGroupList.at(i).ObjectNames);
 226     }
 227   }
 228
 229   // Step 1. Prepare Partition structures.
 230   TopoDS_Shape aResult;
 231   BOPCol_ListOfShape aLS;
 232   QStringList aListOfNames;
 233   TopoDS_Compound aCmp;
 234   BRep_Builder aBB;
 235   aBB.MakeCompound(aCmp);
 236   for (int i=0;i < anInputSplitDataList.size() ;i++) {
 237     const TopoDS_Shape& aSh = anInputSplitDataList.at(i).Shape;
 238     aDM3.Add(aSh, anInputSplitDataList.at(i).ObjectNames);
 239     aLS.Append(aSh);
 240     aBB.Add(aCmp,aSh);
 241 #ifdef DEB_SPLIT_TO_ZONES
 242     //TCollection_AsciiString aName = aNam + i + ".brep";
 243     //BRepTools::Write(aSh, aName.ToCString());
 244 #endif
 245   }
 246 #ifdef DEB_SPLIT_TO_ZONES_CHECK_PARTITION
 247   TCollection_AsciiString aNameBefore = fileNameBefore + ".brep";
 248   BRepTools::Write(aCmp, aNameBefore.ToCString());
 249 #endif
 250
 251   HYDROData_DataMapOfShapeListOfShape aDM1;
 252   if(anInputSplitDataList.size() > 1) {
 253     HYDROData_Transform splitTool;
 254     const Standard_Integer anErr = SplitFaces(aCmp, splitTool);
 255     if(anErr)
 256       return anOutputSplitDataList;
 257     aResult = splitTool.Shape();
 258     if (aResult.IsNull())
 259       return anOutputSplitDataList;
 260     BRepCheck_Analyzer aCheck (aResult);
 261     if(!aCheck.IsValid()) {
 262   #ifdef DEB_SPLIT_TO_ZONES
 263       cout << "result is not valid" <<endl;
 264       BRepTools::Write(aResult, "SplitFacesNV.brep");
 265   #endif
 266       return anOutputSplitDataList;
 267     }
 268   #ifdef DEB_SPLIT_TO_ZONES
 269     BRepTools::Write(aResult, "SplitFacesV.brep");
 270   #endif
 271
 272     // Step 3. Collect history
 273     //HYDROData_DataMapOfShapeListOfShape aDM1;
 274     BOPCol_ListIteratorOfListOfShape anIt(aLS);
 275   #ifdef DEB_SPLIT_TO_ZONES
 276     TCollection_AsciiString aNamM ("EdgM_");
 277     TCollection_AsciiString aNamG ("EdgG_");
 278   #endif
 279     for (int i =1;anIt.More();anIt.Next(),i++) {
 280       Standard_Boolean foundF(Standard_False);
 281       const TopTools_ListOfShape& aListOfNew = splitTool.Modified(anIt.Value());
 282       if(!aListOfNew.IsEmpty())
 283         foundF = Standard_True;
 284
 285       TopTools_ListOfShape aList;
 286
 287       TopTools_ListIteratorOfListOfShape it(aListOfNew);
 288       for(;it.More();it.Next())
 289         aList.Append(it.Value());
 290           /* *********************************************************************
 291       // Bug in History: partition should give only modified entities! => temporary solution is used
 292       //const TopTools_ListOfShape& aListOfGen = splitTool.Generated(anIt.Value());
 293       //if(!aListOfGen.IsEmpty())
 294         //foundF = Standard_True;
 295       //it.Initialize(aListOfGen);
 296       //for(;it.More();it.Next())
 297        // aList.Append(it.Value());
 298            ********************************************************************* */
 299       if(!foundF) // face is not modified
 300         aList.Append (anIt.Value());
 301       aDM1.Add(anIt.Value(), aList);
 302   #ifdef DEB_SPLIT_TO_ZONES
 303       TCollection_AsciiString aName;
 304   #endif
 305       if(!anInputGroupList.isEmpty() ) { // 1
 306         TopExp_Explorer exp (anIt.Value(), TopAbs_EDGE);
 307         for (int j =1;exp.More();exp.Next(),j++) {
 308           aList.Clear();
 309           Standard_Boolean foundE(Standard_False);
 310           const TopTools_ListOfShape& aListM = splitTool.Modified(exp.Current());
 311           //cout << "NB_EDGE_M = " << aListM.Extent() <<endl;
 312           if(aListM.Extent()) foundE = Standard_True;
 313           it.Initialize(aListM);
 314           for(int k=1;it.More();it.Next(),k++) {
 315             aList.Append(it.Value());
 316   #ifdef DEB_SPLIT_TO_ZONES
 317             //aName = aNamM + i + j +k +".brep";
 318             //BRepTools::Write(it.Value(),aName.ToCString());
 319   #endif
 320           }
 321           /* *********************************************************************
 322                   //const TopTools_ListOfShape& aListG = splitTool.Generated(exp.Current());
 323           //if(aListG.Extent()) foundE = Standard_True;
 324           //it.Initialize(aListG);
 325           //for(int k=1;it.More();it.Next(),k++)
 326             //aList.Append(it.Value());
 327           //cout << "NB_EDGE = " << aList.Extent() <<endl;
 328                   ************************************************************************** */
 329           if(!foundE) {
 330             aList.Append (exp.Current());
 331   #ifdef DEB_SPLIT_TO_ZONES
 332             aName = aNamG + i + j +".brep";
 333             BRepTools::Write(exp.Current(),aName.ToCString());
 334             cout << aName.ToCString()<< " = " << exp.Current().TShape() <<endl;
 335   #endif
 336           }
 337           aDM1.Add(exp.Current(), aList);
 338         }
 339       }
 340     }
 341   } else {
 342       aResult = anInputSplitDataList.at(0).Shape; // get single input shape
 343   }
 344
 345   // aDM2: NewShape ==> ListOfOldShapes
 346   HYDROData_DataMapOfShapeListOfShape aDM2;
 347   // make limiting face
 348   HYDROData_DataMapOfShapeListOfShape aDM4;
 349   Standard_Boolean hasLimits(Standard_False);
 350   QString aBndName;
 351   HYDROData_MapOfShape aBndView;
 352   if (! thePolyline.IsNull()) {
 353     const TopoDS_Wire aBndWire = TopoDS::Wire(thePolyline->GetShape());
 354     if(!aBndWire.IsNull()) {
 355       TopoDS_Face limFace;
 356       if(buildLimFace(aBndWire, limFace)) {
 357         TopoDS_Shape aComResult;
 358         BRepAlgoAPI_Common mkCom(aResult, limFace);
 359         if(mkCom.IsDone()) {
 360           aComResult = mkCom.Shape();
 361           BRepCheck_Analyzer aCheck (aComResult);
 362           if(aCheck.IsValid()) {
 363 #ifdef DEB_SPLIT_TO_ZONES
 364             BRepTools::Write(aComResult,"CommonV.brep");
 365             BRepTools::Write(limFace,"limFace.brep");
 366 #endif
 367             aBndName = thePolyline->GetName();
 368             hasLimits = Standard_True; // DM2 should be filled here
 369             TopExp_Explorer exp (limFace, TopAbs_EDGE);
 370             for (int i =1;exp.More();exp.Next(),i++) {
 371               const TopoDS_Shape& anEdge = exp.Current();
 372               if(anEdge.IsNull()) continue;
 373               aBndView.Add(anEdge);
 374               QStringList aListOfNames;
 375               aListOfNames.append(aBndName);
 376               aDM3.Add(anEdge, aListOfNames);
 377               TopTools_ListOfShape aList;
 378               aList.Append(anEdge);
 379               aDM1.Add(anEdge,aList);
 380             }
 381             HYDROData_MapOfShape aView;
 382             exp.Init (aResult, TopAbs_FACE);
 383             for (int i =1;exp.More();exp.Next(),i++) {
 384               const TopoDS_Shape& aFace = exp.Current();
 385               if(!aFace.IsNull()) {
 386                 const TopTools_ListOfShape& aListOfNew = mkCom.Modified(aFace);
 387                 //cout << "Modified: " << aListOfNew.Extent() <<endl;
 388                 if(!aListOfNew.IsEmpty()) {
 389                   aDM4.Add(aFace, aListOfNew);
 390 #ifdef DEB_SPLIT_TO_ZONES
 391                   //TCollection_AsciiString aName = aNam + i + ".brep";
 392                   //BRepTools::Write(aListOfNew.Last(), aName.ToCString());
 393 #endif
 394                 }
 395                 else {
 396                     if(!mkCom.IsDeleted(aFace)) {
 397                       const TopTools_ListOfShape& aListOfGen = mkCom.Generated(aFace);
 398                       if(!aListOfGen.IsEmpty()) {
 399                         /* aDM4.Bind(aFace, aListOfGen); ???   */
 400 #ifdef DEB_SPLIT_TO_ZONES
 401                         //TCollection_AsciiString aName = aNam + i + "g.brep";
 402                         //BRepTools::Write(aListOfGen.Last(), aName.ToCString());
 403 #endif
 404                       }
 405                        else {
 406                         TopTools_ListOfShape aList;
 407                         aList.Append(aFace);
 408                         aDM4.Add(aFace, aList); //the same face - not modified
 409                        }
 410                     }
 411              }
 412              TopExp_Explorer exp2 (aFace, TopAbs_EDGE);
 413              for (int j =1;exp2.More();exp2.Next(),j++) {
 414                const TopoDS_Shape& anEdge = exp2.Current();
 415                if(!anEdge.IsNull()) {
 416                  if(aView.Contains(anEdge)) continue;
 417                  aView.Add(anEdge);
 418                  const TopTools_ListOfShape& aListOfNewEd = mkCom.Modified(anEdge);
 419                  if(!aListOfNewEd.IsEmpty())
 420                    aDM4.Add(anEdge, aListOfNewEd);
 421                  else {
 422                    if(!mkCom.IsDeleted(anEdge)) {
 423                      const TopTools_ListOfShape& aListOfGenEd = mkCom.Generated(anEdge);
 424                      if(!aListOfGenEd.IsEmpty()) {
 425                        /* aDM4.Bind(anEdge, aListOfGenEd); ???*/
 426                      } else {
 427                        TopTools_ListOfShape aList;
 428                        aList.Append(anEdge);
 429                        aDM4.Add(anEdge, aList);//the same edge - not modified
 430                      }
 431                    }
 432                  }
 433                }
 434              }
 435             }
 436             } //end DM4 filling (phase 1)
 437             //cout << "DM4 Ext = " <<aDM4.Extent() <<endl;
 438          // phase 2 (from tool)
 439 #ifdef DEB_SPLIT_TO_ZONES
 440          TCollection_AsciiString aNam("BndEd_");
 441 #endif
 442          TopExp_Explorer expt (limFace, TopAbs_EDGE);
 443          for(int i =1;expt.More();expt.Next(),i++) {
 444            const TopoDS_Shape& anEdge = expt.Current();
 445            if(!anEdge.IsNull()) {
 446              const TopTools_ListOfShape& aListOfNewEd = mkCom.Modified(anEdge);
 447 #ifdef DEB_SPLIT_TO_ZONES
 448              TopTools_ListIteratorOfListOfShape itl(aListOfNewEd);
 449              for(int j=1;itl.More();itl.Next(),j++) {
 450                TCollection_AsciiString aName = aNam + i + "_" + j + ".brep";
 451                BRepTools::Write(itl.Value(), aName.ToCString());
 452                cout <<aName.ToCString()<<" = "<< itl.Value().TShape() <<endl;
 453              }
 454 #endif
 455              if(!aListOfNewEd.IsEmpty())
 456                aDM4.Add(anEdge, aListOfNewEd);
 457              else {
 458                if(!mkCom.IsDeleted(anEdge)) {
 459                  const TopTools_ListOfShape& aListOfGenEd = mkCom.Generated(anEdge);
 460                  if(!aListOfGenEd.IsEmpty()) {
 461                    /* aDM4.Bind(anEdge, aListOfGenEd); ??? */
 462                  } else {
 463                    TopTools_ListOfShape aList;
 464                    aList.Append(anEdge);
 465                    aDM4.Add(anEdge, aList);//the same edge - not modified
 466                  }
 467                }
 468              }
 469            }
 470          }
 471          //cout << "DM4 Ext = " <<aDM4.Extent() <<endl;
 472          if(limplus1Object) {
 473            // fill DM1 (old - new) and DM2 (new - old)
 474            HYDROData_DataMapIteratorOfDataMapOfShapeListOfShape mIt(aDM4);
 475            for(;mIt.More();mIt.Next()) {
 476              const TopoDS_Shape& aKey = mIt.Key();//old
 477              TopTools_ListOfShape aList;
 478              aList.Append(aKey);
 479              const TopTools_ListOfShape& aListOfNew = mIt.Value();
 480              aDM1.Add(aKey, aListOfNew);
 481              TopTools_ListIteratorOfListOfShape it(aListOfNew);
 482              for(;it.More();it.Next()) {
 483                if(!aDM2.Contains(it.Value()))
 484                  aDM2.Add(it.Value(), aList);
 485                else {
 486                  TopTools_ListOfShape& aList = aDM2.ChangeFromKey(it.Value());
 487                  aList.Prepend(aKey);
 488                }
 489              }
 490            }
 491          } else {
 492            HYDROData_DataMapIteratorOfDataMapOfShapeListOfShape aMIt(aDM1);
 493            //DM4 contains Old - New after common op. DM1: old - new after Split op.
 494            for(;aMIt.More();aMIt.Next()) {
 495              const TopoDS_Shape& aKey = aMIt.Key();
 496              TopTools_ListOfShape aList;
 497              aList.Append(aKey);
 498              const TopTools_ListOfShape& aListOfNew = aMIt.Value();
 499              TopTools_ListIteratorOfListOfShape it(aListOfNew);
 500              for(;it.More();it.Next()) {// iterate new: b1, b2, b3...
 501                if(!aDM4.Contains(it.Value())) // bi - is deleted
 502                  continue; // go to the next bi
 503                else {
 504                  const TopTools_ListOfShape& aListOfNew4 = aDM4.FindFromKey(it.Value());
 505                   TopTools_ListIteratorOfListOfShape it4(aListOfNew4); // {c1, c2, c3,...}
 506                   for(;it4.More();it4.Next()) {
 507                     if(!aDM2.Contains(it4.Value()))
 508                       aDM2.Add(it4.Value(), aList);
 509                     else {
 510                       TopTools_ListOfShape& aList = aDM2.ChangeFromKey(it4.Value());
 511                       aList.Prepend(aKey);
 512                     }
 513                   }
 514                }
 515              }
 516            }
 517          }
 518          } else {
 519              hasLimits = Standard_False;
 520 #ifdef DEB_SPLIT_TO_ZONES
 521              BRepTools::Write(aComResult,"CommonNV.brep");
 522 #endif
 523           }
 524        }
 525       }
 526     }
 527   }// end limits processing
 528   if(!hasLimits) {
 529     HYDROData_DataMapIteratorOfDataMapOfShapeListOfShape aMIt(aDM1);
 530     for(;aMIt.More();aMIt.Next()) {
 531       const TopoDS_Shape& aKey = aMIt.Key();
 532       TopTools_ListOfShape aList;
 533       aList.Append(aKey);
 534       const TopTools_ListOfShape& aListOfNew = aMIt.Value();
 535       TopTools_ListIteratorOfListOfShape it(aListOfNew);
 536       for(;it.More();it.Next()) {
 537         if(!aDM2.Contains(it.Value()))
 538           aDM2.Add(it.Value(), aList);
 539         else {
 540           TopTools_ListOfShape& aList = aDM2.ChangeFromKey(it.Value());
 541           aList.Prepend(aKey);
 542         }
 543       }
 544     }
 545   }
 546   // Step 4. Fill output structure.
 547 #ifdef DEB_SPLIT_TO_ZONES
 548   TCollection_AsciiString aNam4 ("SC_");
 549 #endif
 550   HYDROData_DataMapIteratorOfDataMapOfShapeListOfShape aMIt(aDM2);
 551   for(int i =1;aMIt.More();aMIt.Next(),i++) {
 552     SplitData aDestSplitData;
 553     const TopoDS_Shape& aKey = aMIt.Key(); //new
 554     aDestSplitData.Shape = aKey;
 555     if(aKey.ShapeType() == TopAbs_FACE)
 556       aDestSplitData.Type = SplitData::Data_Zone;
 557     else {
 558       aDestSplitData.Type = SplitData::Data_Edge;
 559 #ifdef DEB_SPLIT_TO_ZONES
 560       TCollection_AsciiString aName = aNam4 + i + ".brep";
 561       BRepTools::Write(aKey,aName.ToCString());
 562 #endif
 563     }
 564
 565     QStringList aListOfNames; // names processing
 566     const TopTools_ListOfShape& aListOfOld = aMIt.Value();
 567     TopTools_ListIteratorOfListOfShape it(aListOfOld);
 568     for(int j =1;it.More();it.Next(),j++) {
 569       const TopoDS_Shape& aSh = it.Value(); //old
 570       if(aDM3.Contains(aSh)) {
 571         const QStringList& ObjectNames = aDM3.FindFromKey(aSh);
 572         aListOfNames.append(ObjectNames);
 573 #ifdef DEB_SPLIT_TO_ZONES
 574         TCollection_AsciiString aName = aNam4 + i +"_" + j + ".brep";
 575         BRepTools::Write(aSh ,aName.ToCString());
 576 #endif
 577       }    else {
 578 #ifdef DEB_SPLIT_TO_ZONES
 579         TCollection_AsciiString aName = aNam4 +"__" + i +"_" + j + ".brep";
 580         BRepTools::Write(aSh ,aName.ToCString());
 581         cout <<aName.ToCString()<<" = "<< aSh.TShape() <<endl;
 582 #endif
 583           if(aBndView.Contains(aSh) && hasLimits) {
 584             aListOfNames.append(aBndName);
 585 #ifdef DEB_SPLIT_TO_ZONES
 586             cout << " BndName = "<<aBndName.toStdString() <<endl;
 587 #endif
 588           }
 589       }
 590     }
 591
 592     aDestSplitData.ObjectNames = aListOfNames;
 593     anOutputSplitDataList.append(aDestSplitData);
 594 #ifdef DEB_SPLIT_TO_ZONES
 595     QString aStr = aDestSplitData.ObjectNames.join(" ");
 596     cout << "New# = " << i<< " Names = "<<aStr.toStdString() << " size = " <<aDestSplitData.ObjectNames.size() <<endl;
 597 #endif
 598   }
 599
 600   return anOutputSplitDataList;
 601 }
 602