0025412: BRepMesh_FastDiscretFace::tryToInsertAnalyticVertex() - fix invalid comparison

[occt.git] / src / BOPAlgo / BOPAlgo_BOP_1.cxx

// Created by: Peter KURNEV
// Copyright (c) 1999-2014 OPEN CASCADE SAS
//
// This file is part of Open CASCADE Technology software library.
//
// This library is free software; you can redistribute it and/or modify it under
// the terms of the GNU Lesser General Public License version 2.1 as published
// by the Free Software Foundation, with special exception defined in the file
// OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
// distribution for complete text of the license and disclaimer of any warranty.
//
// Alternatively, this file may be used under the terms of Open CASCADE
// commercial license or contractual agreement.

#include <BOPAlgo_BOP.ixx>

#include <BOPDS_DS.hxx>
#include <BOPDS_PaveBlock.hxx>
#include <BOPDS_CommonBlock.hxx>
#include <BOPDS_FaceInfo.hxx>
#include <BOPDS_ListOfPaveBlock.hxx>
#include <BOPDS_MapOfPaveBlock.hxx>
#include <BOPDS_VectorOfListOfPaveBlock.hxx>
#include <BOPDS_VectorOfFaceInfo.hxx>

#include <BOPTools.hxx>
#include <BOPTools_AlgoTools.hxx>

#include <BOPCol_MapOfShape.hxx>

#include <BRep_Builder.hxx>
#include <TopExp_Explorer.hxx>


//=======================================================================
//function : BuildSection
//purpose  : 
//=======================================================================
void BOPAlgo_BOP::BuildSection()
{

  Standard_Integer i, j, k, nE, nF1, nF2, aNbPB, aNbFF;
  Standard_Boolean bFlag;
  TopoDS_Shape aRC;
  BRep_Builder aBB;
  BOPCol_MapOfShape aME;
  BOPCol_IndexedMapOfShape aME1, aME2;
  //
  myErrorStatus=0;
  //
  BOPTools_AlgoTools::MakeContainer(TopAbs_COMPOUND, aRC);  
  BOPDS_VectorOfInterfFF& aFFs=myDS->InterfFF();
  aNbFF=aFFs.Extent();
  //
  for (i=0; i<aNbFF; ++i) {
    BOPDS_InterfFF& aFF=aFFs(i);
    aFF.Indices(nF1, nF2);
    const BOPDS_FaceInfo& aFI1=myDS->FaceInfo(nF1);
    const BOPDS_FaceInfo& aFI2=myDS->FaceInfo(nF2);
    //
    const BOPDS_IndexedMapOfPaveBlock& aMPBIn1=aFI1.PaveBlocksIn();
    const BOPDS_IndexedMapOfPaveBlock& aMPBOn1=aFI1.PaveBlocksOn();
    const BOPDS_IndexedMapOfPaveBlock& aMPBSc1=aFI1.PaveBlocksSc();
    //
    const BOPDS_IndexedMapOfPaveBlock& aMPBIn2=aFI2.PaveBlocksIn();
    const BOPDS_IndexedMapOfPaveBlock& aMPBOn2=aFI2.PaveBlocksOn();
    //
    //1. Section edges
    aNbPB = aMPBSc1.Extent();
    for (j=1; j<=aNbPB; ++j) {
      const Handle(BOPDS_PaveBlock)& aPB=aMPBSc1(j);
      nE=aPB->Edge();
      const TopoDS_Shape& aE=myDS->Shape(nE);
      if (aME.Add(aE)) {
        aBB.Add(aRC, aE);
      }
    }
    //2. Common edges
    BOPDS_IndexedMapOfPaveBlock aMPB[4] = {aMPBOn2, aMPBIn1, aMPBIn2, aMPBOn1};
    for (k = 0; k < 3; ++k) {
      aNbPB = aMPB[k].Extent();
      for (j=1; j<=aNbPB; ++j) {
        const Handle(BOPDS_PaveBlock)& aPB=aMPB[k](j);
        bFlag = (k==0) ? aMPB[3].Contains(aPB) :
          (aMPB[k-1].Contains(aPB) || aMPB[k+1].Contains(aPB));
        if (bFlag) {
          nE=aPB->Edge();
          const TopoDS_Shape& aE=myDS->Shape(nE);
          if (aME.Add(aE)) {
            aBB.Add(aRC, aE);
          }
        }
      }
    }
    //3. Shared edges
    aME1.Clear();
    aME2.Clear();
    //
    const TopoDS_Face& aF1=(*(TopoDS_Face *)(&myDS->Shape(nF1)));
    const TopoDS_Face& aF2=(*(TopoDS_Face *)(&myDS->Shape(nF2)));
    //
    BOPTools::MapShapes(aF1, TopAbs_EDGE, aME1);
    BOPTools::MapShapes(aF2, TopAbs_EDGE, aME2);
    //
    aNbPB = aME1.Extent();
    for (j=1; j<=aNbPB; ++j) {
      const TopoDS_Shape& aE = aME1(j);
      if (aME2.Contains(aE)) {
        if (aME.Add(aE)) {
          aBB.Add(aRC, aE);
        }
      }
    }
  }
  //
  myShape=aRC;
}
//=======================================================================
//function : Generated
//purpose  : 
//=======================================================================
const TopTools_ListOfShape& BOPAlgo_BOP::Generated
  (const TopoDS_Shape& theS)
{
  myHistShapes.Clear();
  if (theS.IsNull() || (myOperation != BOPAlgo_SECTION)) {
    return myHistShapes;
  }
  //
  TopAbs_ShapeEnum aType = theS.ShapeType();
  if (aType != TopAbs_FACE) {
    return myHistShapes;
  }
  //
  Standard_Integer nS = myDS->Index(theS);
  if (nS < 0) {
    return myHistShapes;
  }
  //
  if (!myDS->HasFaceInfo(nS)) {
    return myHistShapes;
  }
  //
  //collect section edges of the face theS
  Standard_Integer i, aNb, nSp;
  //
  const BOPDS_FaceInfo& aFI = myDS->FaceInfo(nS);
  const BOPDS_IndexedMapOfPaveBlock& aMPBSc = aFI.PaveBlocksSc();
  aNb = aMPBSc.Extent();
  for (i = 1; i <= aNb; ++i) {
    const Handle(BOPDS_PaveBlock)& aPB = aMPBSc(i);
    nSp = aPB->Edge();
    const TopoDS_Shape& aSp = myDS->Shape(nSp);
    myHistShapes.Append(aSp);
  }
  //
  return myHistShapes;
}