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

// Created by: Peter KURNEV
// Copyright (c) 2010-2014 OPEN CASCADE SAS
// Copyright (c) 2007-2010 CEA/DEN, EDF R&D, OPEN CASCADE
// Copyright (c) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN, CEDRAT,
//                         EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
//
// 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 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_PaveFiller.ixx>

#include <NCollection_IncAllocator.hxx>

#include <Bnd_Box.hxx>

#include <TopoDS_Vertex.hxx>
#include <TopoDS_Edge.hxx>
#include <TopoDS_Face.hxx>
#include <BRep_Tool.hxx>

#include <IntTools_EdgeFace.hxx>
#include <IntTools_Range.hxx>
#include <IntTools_SequenceOfCommonPrts.hxx>
#include <IntTools_CommonPrt.hxx>
#include <BOPTools_AlgoTools.hxx>

#include <BOPCol_MapOfInteger.hxx>

#include <BOPInt_Context.hxx>

#include <BOPDS_Interf.hxx>
#include <BOPDS_Iterator.hxx>
#include <BOPDS_PaveBlock.hxx>
#include <BOPDS_MapOfPaveBlock.hxx>
#include <BOPDS_DataMapOfPaveBlockListOfInteger.hxx>
#include <BOPDS_CommonBlock.hxx>
#include <BOPDS_Pave.hxx>

#include <BOPTools_AlgoTools.hxx>
#include <BOPDS_CoupleOfPaveBlocks.hxx>
#include <BRepBndLib.hxx>
#include <BOPAlgo_Tools.hxx>
#include <BOPInt_Tools.hxx>
#include <BRepAdaptor_Curve.hxx>
#include <BRep_Builder.hxx>
#include <GeomAPI_ProjectPointOnSurf.hxx>


//=======================================================================
//function : PerformEF
//purpose  : 
//=======================================================================
  void BOPAlgo_PaveFiller::PerformEF()
{
  Standard_Integer iSize;
  //
  myErrorStatus=0;
  //
  myIterator->Initialize(TopAbs_EDGE, TopAbs_FACE);
  iSize=myIterator->ExpectedLength();
  if (!iSize) {
    return; 
  }
  //----------------------------------------------------------------------
  Standard_Boolean bJustAdd, bV[2];
  Standard_Integer nE, nF, aDiscretize, i, aNbCPrts, iX, nV[2];
  Standard_Real aTolE, aTolF, aTS1, aTS2, aT1, aT2, aDeflection;
  Handle(NCollection_IncAllocator) aAllocator;
  TopAbs_ShapeEnum aType;
  BOPDS_ListIteratorOfListOfPaveBlock aIt;
  //-----------------------------------------------------scope f
  //
  BRep_Builder aBB;
  //
  aAllocator=new NCollection_IncAllocator();
  
  BOPCol_MapOfInteger aMIEFC(100, aAllocator);
  BOPDS_IndexedDataMapOfShapeCoupleOfPaveBlocks aMVCPB(100, aAllocator);
  BOPDS_IndexedDataMapOfPaveBlockListOfInteger aMPBLI(100, aAllocator);
  //
  aDiscretize=35;
  aDeflection=0.01;
  //
  BOPDS_VectorOfInterfEF& aEFs=myDS->InterfEF();
  aEFs.SetStartSize(iSize);
  aEFs.SetIncrement(iSize);
  aEFs.Init();
  //
  for (; myIterator->More(); myIterator->Next()) {
    myIterator->Value(nE, nF, bJustAdd);
    if(bJustAdd) {
      continue;
    }
    //
    const BOPDS_ShapeInfo& aSIE=myDS->ShapeInfo(nE);
    if (aSIE.HasFlag()){//degenerated 
      continue;
    }
    //
    const TopoDS_Edge& aE=(*(TopoDS_Edge *)(&aSIE.Shape()));
    const TopoDS_Face& aF=(*(TopoDS_Face *)(&myDS->Shape(nF)));
    const Bnd_Box& aBBF=myDS->ShapeInfo(nF).Box(); 
    //
    BOPDS_FaceInfo& aFI=myDS->ChangeFaceInfo(nF);
    const BOPDS_IndexedMapOfPaveBlock& aMPBF=aFI.PaveBlocksOn();
    const BOPCol_MapOfInteger& aMIFOn=aFI.VerticesOn();
    const BOPCol_MapOfInteger& aMIFIn=aFI.VerticesIn();
    //
    aTolE=BRep_Tool::Tolerance(aE);
    aTolF=BRep_Tool::Tolerance(aF);
    //
    BOPDS_ListOfPaveBlock& aLPB=myDS->ChangePaveBlocks(nE);
    aIt.Initialize(aLPB);
    for (; aIt.More(); aIt.Next()) {
      Handle(BOPDS_PaveBlock)& aPB=aIt.ChangeValue();
      //
      const Handle(BOPDS_PaveBlock) aPBR=myDS->RealPaveBlock(aPB);
      if (aMPBF.Contains(aPBR)) {
        continue;
      }
      //
      if (!aPB->HasShrunkData()) {
        FillShrunkData(aPB);
        if (myWarningStatus) {
          continue;
        }
      }
      //
      Bnd_Box aBBE;
      aPB->ShrunkData(aTS1, aTS2, aBBE);
      //
      if (aBBF.IsOut (aBBE)) {
        continue;
      }
      //
      // -----------f
      IntTools_EdgeFace aEdgeFace;
      //
      aEdgeFace.SetEdge (aE);
      aEdgeFace.SetFace (aF);
      aEdgeFace.SetTolE (aTolE);
      aEdgeFace.SetTolF (aTolF);
      aEdgeFace.SetDiscretize (aDiscretize);
      aEdgeFace.SetDeflection (aDeflection);
      aEdgeFace.SetContext(myContext);
      //
      IntTools_Range aSR(aTS1, aTS2);
      IntTools_Range anewSR=aSR;
      BOPTools_AlgoTools::CorrectRange(aE, aF, aSR, anewSR);
      //
      aPB->Range(aT1, aT2);
      IntTools_Range aPBRange(aT1, aT2);
      aSR = aPBRange;
      BOPTools_AlgoTools::CorrectRange(aE, aF, aSR, aPBRange);
      //
      aEdgeFace.SetRange (aPBRange);
      //
      aEdgeFace.Perform();
      if (!aEdgeFace.IsDone()) {
        continue;
      }
      //
      aPB->Indices(nV[0], nV[1]);
      //
      const IntTools_SequenceOfCommonPrts& aCPrts=aEdgeFace.CommonParts();
      aNbCPrts=aCPrts.Length();
      for (i=1; i<=aNbCPrts; ++i) {
        const IntTools_CommonPrt& aCPart=aCPrts(i);
        aType=aCPart.Type();
        switch (aType) {
        case TopAbs_VERTEX:  {
          Standard_Boolean bIsOnPave[2];
          Standard_Integer j;
          Standard_Real aT, aTolToDecide; 
          TopoDS_Vertex aVnew;
          
          BOPInt_Tools::VertexParameter(aCPart, aT);
          BOPTools_AlgoTools::MakeNewVertex(aE, aT, aF, aVnew);
          //
          const IntTools_Range& aR=aCPart.Range1();
          aTolToDecide=5.e-8;
          //
          IntTools_Range aR1(aT1, anewSR.First()), aR2(anewSR.Last(), aT2);
          //
          bIsOnPave[0]=BOPInt_Tools::IsInRange(aR1, aR, aTolToDecide); 
          bIsOnPave[1]=BOPInt_Tools::IsInRange(aR2, aR, aTolToDecide); 
          //
          if (bIsOnPave[0] && bIsOnPave[1]) {
            bV[0]=CheckFacePaves(nV[0], aMIFOn, aMIFIn);
            bV[1]=CheckFacePaves(nV[1], aMIFOn, aMIFIn);
            if (bV[0] && bV[1]) {
              iX=aEFs.Append()-1;
              IntTools_CommonPrt aCP = aCPart;
              aCP.SetType(TopAbs_EDGE);
              BOPDS_InterfEF& aEF=aEFs(iX);
              aEF.SetIndices(nE, nF);
              aEF.SetCommonPart(aCP);
              myDS->AddInterf(nE, nF);
              // 3          
              BOPAlgo_Tools::FillMap(aPB, nF, aMPBLI, aAllocator);
              break;
            }
          }
          for (j=0; j<2; ++j) {
            if (bIsOnPave[j]) {
              bV[j]=CheckFacePaves(nV[j], aMIFOn, aMIFIn);
              if (bV[j]) {
                const TopoDS_Vertex& aV = (*(TopoDS_Vertex *)(&myDS->Shape(nV[j])));
                BOPTools_AlgoTools::UpdateVertex(aE, aT, aV);
                BOPDS_ShapeInfo& aSIDS=myDS->ChangeShapeInfo(nV[j]);
                Bnd_Box& aBoxDS=aSIDS.ChangeBox();
                BRepBndLib::Add(aV, aBoxDS);
              }
              else {
                bIsOnPave[j] = ForceInterfVF(nV[j], nF);
              }
            }
          }
          //
          if (!bIsOnPave[0] && !bIsOnPave[1]) {
            if (CheckFacePaves(aVnew, aMIFOn)) {
              continue;
            }
            //
            const gp_Pnt& aPnew = BRep_Tool::Pnt(aVnew);
            if (!myContext->IsValidPointForFace(aPnew, aF, aTolE+aTolF)) {
              continue;
            }
            //
            aBB.UpdateVertex(aVnew, aTolE);
            //
            aMIEFC.Add(nF);
            // 1
            iX=aEFs.Append()-1;
            BOPDS_InterfEF& aEF=aEFs(iX);
            aEF.SetIndices(nE, nF);
            aEF.SetCommonPart(aCPart);
            // 2
            myDS->AddInterf(nE, nF);
            // 3
            BOPDS_CoupleOfPaveBlocks aCPB;
            //
            aCPB.SetPaveBlocks(aPB, aPB);
            aCPB.SetIndexInterf(iX);
            aMVCPB.Add(aVnew, aCPB);
          }
        }
          break;
        case TopAbs_EDGE:  {
          aMIEFC.Add(nF);
          //
          // 1
          iX=aEFs.Append()-1;
          BOPDS_InterfEF& aEF=aEFs(iX);
          aEF.SetIndices(nE, nF);
          //
          bV[0]=CheckFacePaves(nV[0], aMIFOn, aMIFIn);
          bV[1]=CheckFacePaves(nV[1], aMIFOn, aMIFIn);
          if (!bV[0] || !bV[1]) {
            myDS->AddInterf(nE, nF);
            break;
          }
          //update tolerance of edge if needed
          if (aTolE < aTolF) {
            myDS->UpdateEdgeTolerance(nE, aTolF);
            aTolE = aTolF;
          }
          aEF.SetCommonPart(aCPart);
          // 2
          myDS->AddInterf(nE, nF);
          // 3
          BOPAlgo_Tools::FillMap(aPB, nF, aMPBLI, aAllocator);
          
        }
          break; 
        default:
          break; 
        }//switch (aType) {
      }//for (i=1; i<=aNbCPrts; ++i) {
      // -----------t
    }//for (; aIt.More(); aIt.Next()) {
  }//for (; myIterator->More(); myIterator->Next()) {
  // 
  //=========================================
  // post treatment
  //=========================================
  BOPAlgo_Tools::PerformCommonBlocks(aMPBLI, aAllocator, myDS);
  PerformVerticesEF(aMVCPB, aAllocator);
  //
  // Update FaceInfoIn for all faces having EF common parts
  BOPCol_MapIteratorOfMapOfInteger aItMI;
  aItMI.Initialize(aMIEFC);
  for (; aItMI.More(); aItMI.Next()) {
    nF=aItMI.Value();
    myDS->UpdateFaceInfoIn(nF);
  }
  // Refine FaceInfoOn to remove all formal pave blocks 
  // made during EF processing 
  //myDS->RefineFaceInfoOn();
  //-----------------------------------------------------scope t
  aMIEFC.Clear();
  aMVCPB.Clear();
  aMPBLI.Clear();
  aAllocator.Nullify();
  //
  
}
//=======================================================================
//function : PerformVertices1
//purpose  : 
//=======================================================================
  Standard_Integer BOPAlgo_PaveFiller::PerformVerticesEF
    (BOPDS_IndexedDataMapOfShapeCoupleOfPaveBlocks& theMVCPB,
     Handle(NCollection_BaseAllocator)& theAllocator)
{
  Standard_Integer aNbV, iRet;
  //
  iRet=0;
  aNbV=theMVCPB.Extent();
  if (!aNbV) {
    return iRet;
  }
  //
  Standard_Integer nVx, nVSD, iV, iErr, nE, iFlag, iX, i, aNbPBLI;
  Standard_Real aT;
  TopoDS_Shape aV;
  BOPCol_ListIteratorOfListOfShape aItLS;
  BOPCol_ListIteratorOfListOfInteger aItLI;
  BOPDS_PDS aPDS;
  BOPDS_ShapeInfo aSI;
  BOPDS_Pave aPave;
  //
  BOPCol_ListOfShape aLS(theAllocator);
  BOPCol_DataMapOfShapeInteger aMVI(100, theAllocator);
  BOPDS_IndexedDataMapOfPaveBlockListOfInteger aMPBLI(100, theAllocator);
  BOPAlgo_PaveFiller aPF(theAllocator); 
  //
  aSI.SetShapeType(TopAbs_VERTEX);
  BOPDS_VectorOfInterfEF& aEFs=myDS->InterfEF();
  //
  // 1 prepare arguments
  for (i=1; i<=aNbV; ++i) {
    const TopoDS_Shape& aS=theMVCPB.FindKey(i);
    aLS.Append(aS);
  }
  //
  // 2 Fuse vertices
  aPF.SetArguments(aLS);
  aPF.Perform();
  iErr=aPF.ErrorStatus();
  if (iErr) {
    iRet=1;
    return iRet;
  }
  aPDS=aPF.PDS();
  //
  // 3 Add new vertices to theDS; 
  // 4 Map PaveBlock/ListOfVertices to add to this PaveBlock ->aMPBLI
  aItLS.Initialize(aLS);
  for (; aItLS.More(); aItLS.Next()) {
    const TopoDS_Shape& aVx=aItLS.Value();
    nVx=aPDS->Index(aVx);
    //
    if (aPDS->HasShapeSD(nVx, nVSD)) {
      aV=aPDS->Shape(nVSD);
    }
    else {
      aV=aVx;
    }
    // index of new vertex in theDS -> iV
    if (!aMVI.IsBound(aV)) {
      aSI.SetShape(aV);
      iV=myDS->Append(aSI);
      //
      BOPDS_ShapeInfo& aSIDS=myDS->ChangeShapeInfo(iV);
      Bnd_Box& aBox=aSIDS.ChangeBox();
      BRepBndLib::Add(aV, aBox);
      //
      aMVI.Bind(aV, iV);
    }
    else {
      iV=aMVI.Find(aV);
    }
    //
    BOPDS_CoupleOfPaveBlocks &aCPB=theMVCPB.ChangeFromKey(aVx);
    aCPB.SetIndex(iV);
    // update EF interference
    iX=aCPB.IndexInterf();
    BOPDS_InterfEF& aEF=aEFs(iX);
    aEF.SetIndexNew(iV);
    // map aMPBLI
    const Handle(BOPDS_PaveBlock)& aPB=aCPB.PaveBlock1();
    if (aMPBLI.Contains(aPB)) {
      BOPCol_ListOfInteger& aLI=aMPBLI.ChangeFromKey(aPB);
      aLI.Append(iV);
    }
    else {
      BOPCol_ListOfInteger aLI(theAllocator);
      aLI.Append(iV);
      aMPBLI.Add(aPB, aLI);
    }
  }
  //
  // 5 
  // 5.1  Compute Extra Paves and 
  // 5.2. Add Extra Paves to the PaveBlocks
  aNbPBLI=aMPBLI.Extent();
  for (i=1; i<=aNbPBLI; ++i) {
    Handle(BOPDS_PaveBlock) aPB=aMPBLI.FindKey(i);
    const BOPCol_ListOfInteger& aLI=aMPBLI.FindFromIndex(i);
    nE=aPB->OriginalEdge();
    const TopoDS_Edge& aE=(*(TopoDS_Edge *)(&myDS->Shape(nE)));
    // 
    aItLI.Initialize(aLI);
    for (; aItLI.More(); aItLI.Next()) {
      nVx=aItLI.Value();
      const TopoDS_Vertex& aVx=(*(TopoDS_Vertex *)(&myDS->Shape(nVx)));
      //
      iFlag=myContext->ComputeVE (aVx, aE, aT);
      if (!iFlag) {
        aPave.SetIndex(nVx);
        aPave.SetParameter(aT);
        aPB->AppendExtPave(aPave);
      }
    }
  }
  // 6  Split PaveBlocksa
  for (i=1; i<=aNbPBLI; ++i) {
    Handle(BOPDS_PaveBlock) aPB=aMPBLI.FindKey(i);
    nE=aPB->OriginalEdge();
    // 3
    if (!myDS->IsCommonBlock(aPB)) {
      myDS->UpdatePaveBlock(aPB);
    }
    else {
      const Handle(BOPDS_CommonBlock)& aCB=myDS->CommonBlock(aPB);
      myDS->UpdateCommonBlock(aCB);
    }    
  }//for (; aItMPBLI.More(); aItMPBLI.Next()) {
  // 
  return iRet;
}
//=======================================================================
// function: CheckFacePaves
// purpose: 
//=======================================================================
  Standard_Boolean BOPAlgo_PaveFiller::CheckFacePaves (const Standard_Integer nVx,
                                                       const BOPCol_MapOfInteger& aMIFOn,
                                                       const BOPCol_MapOfInteger& aMIFIn)
{
  Standard_Boolean bRet;
  Standard_Integer nV;
  BOPCol_MapIteratorOfMapOfInteger aIt;
  //
  bRet=Standard_False;
  //
  aIt.Initialize(aMIFOn);
  for (; aIt.More(); aIt.Next()) {
    nV=aIt.Value();
    if (nV==nVx) {
      bRet=!bRet;
      return bRet;
    }
  }
  aIt.Initialize(aMIFIn);
  for (; aIt.More(); aIt.Next()) {
    nV=aIt.Value();
    if (nV==nVx) {
      bRet=!bRet;
      return bRet;
    }
  }
  //
  return bRet;
}
//=======================================================================
// function: CheckFacePaves
// purpose: 
//=======================================================================
  Standard_Boolean BOPAlgo_PaveFiller::CheckFacePaves (const TopoDS_Vertex& aVnew,
                                                       const BOPCol_MapOfInteger& aMIF)
{
  Standard_Boolean bRet;
  Standard_Integer nV, iFlag;
  BOPCol_MapIteratorOfMapOfInteger aIt;
  //
  bRet=Standard_True;
  //
  aIt.Initialize(aMIF);
  for (; aIt.More(); aIt.Next()) {
    nV=aIt.Value();
    const TopoDS_Vertex& aV=(*(TopoDS_Vertex *)(&myDS->Shape(nV)));
    iFlag=BOPTools_AlgoTools::ComputeVV(aVnew, aV);
    if (!iFlag) {
      return bRet;
    }
  }
  //
  return !bRet;
}
//=======================================================================
//function : ForceInterfVF
//purpose  : 
//=======================================================================
Standard_Boolean BOPAlgo_PaveFiller::ForceInterfVF(const Standard_Integer nV, 
                                                   const Standard_Integer nF)
{
  Standard_Boolean bRet;
  //
  bRet = Standard_False;
  const TopoDS_Vertex& aV = *(TopoDS_Vertex*)&myDS->Shape(nV);
  const TopoDS_Face&   aF = *(TopoDS_Face*)  &myDS->Shape(nF);
  //
  GeomAPI_ProjectPointOnSurf& aProj = myContext->ProjPS(aF);
  const gp_Pnt& aP = BRep_Tool::Pnt(aV);
  aProj.Perform(aP);
  if (!aProj.IsDone()) {
    return bRet;
  }
  Standard_Real aDist, U, V;
  //
  aDist=aProj.LowerDistance();
  aProj.LowerDistanceParameters(U, V);
  //
  gp_Pnt2d aP2d(U, V);
  bRet = myContext->IsPointInFace (aF, aP2d);
  if (bRet) {
    Standard_Integer i;
    BRep_Builder aBB;
    //
    BOPDS_VectorOfInterfVF& aVFs=myDS->InterfVF();
    i=aVFs.Append()-1;
    BOPDS_InterfVF& aVF=aVFs(i);
    aVF.SetIndices(nV, nF);
    aVF.SetUV(U, V);
    //
    myDS->AddInterf(nV, nF);
    //
    aBB.UpdateVertex(aV, aDist);
    BOPDS_ShapeInfo& aSIDS=myDS->ChangeShapeInfo(nV);
    Bnd_Box& aBoxDS=aSIDS.ChangeBox();
    BRepBndLib::Add(aV, aBoxDS);
    //
    BOPDS_FaceInfo& aFI=myDS->ChangeFaceInfo(nF);
    BOPCol_MapOfInteger& aMVIn=aFI.ChangeVerticesIn();
    aMVIn.Add(nV);
  }
  //
  return bRet;
}
Commit	Line	Data
4e57c75e	1	// Created by: Peter KURNEV
973c2be1	2	// Copyright (c) 2010-2014 OPEN CASCADE SAS
4e57c75e	3	// Copyright (c) 2007-2010 CEA/DEN, EDF R&D, OPEN CASCADE
	4	// Copyright (c) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN, CEDRAT,
	5	// EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
	6	//
973c2be1	7	// This file is part of Open CASCADE Technology software library.
4e57c75e	8	//
973c2be1	9	// This library is free software; you can redistribute it and / or modify it
	10	// under the terms of the GNU Lesser General Public version 2.1 as published
	11	// by the Free Software Foundation, with special exception defined in the file
	12	// OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
	13	// distribution for complete text of the license and disclaimer of any warranty.
4e57c75e	14	//
973c2be1	15	// Alternatively, this file may be used under the terms of Open CASCADE
973c2be1	16	// commercial license or contractual agreement.
4e57c75e	17
	18	#include <BOPAlgo_PaveFiller.ixx>
	19
	20	#include <NCollection_IncAllocator.hxx>
	21
	22	#include <Bnd_Box.hxx>
	23
	24	#include <TopoDS_Vertex.hxx>
	25	#include <TopoDS_Edge.hxx>
	26	#include <TopoDS_Face.hxx>
	27	#include <BRep_Tool.hxx>
	28
	29	#include <IntTools_EdgeFace.hxx>
	30	#include <IntTools_Range.hxx>
	31	#include <IntTools_SequenceOfCommonPrts.hxx>
	32	#include <IntTools_CommonPrt.hxx>
	33	#include <BOPTools_AlgoTools.hxx>
	34
	35	#include <BOPCol_MapOfInteger.hxx>
	36
	37	#include <BOPInt_Context.hxx>
	38
	39	#include <BOPDS_Interf.hxx>
	40	#include <BOPDS_Iterator.hxx>
	41	#include <BOPDS_PaveBlock.hxx>
	42	#include <BOPDS_MapOfPaveBlock.hxx>
	43	#include <BOPDS_DataMapOfPaveBlockListOfInteger.hxx>
	44	#include <BOPDS_CommonBlock.hxx>
	45	#include <BOPDS_Pave.hxx>
	46
	47	#include <BOPTools_AlgoTools.hxx>
	48	#include <BOPDS_CoupleOfPaveBlocks.hxx>
	49	#include <BRepBndLib.hxx>
	50	#include <BOPAlgo_Tools.hxx>
	51	#include <BOPInt_Tools.hxx>
	52	#include <BRepAdaptor_Curve.hxx>
	53	#include <BRep_Builder.hxx>
b4109929	54	#include <GeomAPI_ProjectPointOnSurf.hxx>
4e57c75e	55
	56
	57	//=======================================================================
	58	//function : PerformEF
	59	//purpose :
	60	//=======================================================================
	61	void BOPAlgo_PaveFiller::PerformEF()
	62	{
	63	Standard_Integer iSize;
	64	//
	65	myErrorStatus=0;
	66	//
	67	myIterator->Initialize(TopAbs_EDGE, TopAbs_FACE);
	68	iSize=myIterator->ExpectedLength();
	69	if (!iSize) {
	70	return;
	71	}
	72	//----------------------------------------------------------------------
b4109929	73	Standard_Boolean bJustAdd, bV[2];
	74	Standard_Integer nE, nF, aDiscretize, i, aNbCPrts, iX, nV[2];
	75	Standard_Real aTolE, aTolF, aTS1, aTS2, aT1, aT2, aDeflection;
4e57c75e	76	Handle(NCollection_IncAllocator) aAllocator;
	77	TopAbs_ShapeEnum aType;
	78	BOPDS_ListIteratorOfListOfPaveBlock aIt;
	79	//-----------------------------------------------------scope f
	80	//
	81	BRep_Builder aBB;
	82	//
	83	aAllocator=new NCollection_IncAllocator();
	84
	85	BOPCol_MapOfInteger aMIEFC(100, aAllocator);
	86	BOPDS_IndexedDataMapOfShapeCoupleOfPaveBlocks aMVCPB(100, aAllocator);
	87	BOPDS_IndexedDataMapOfPaveBlockListOfInteger aMPBLI(100, aAllocator);
	88	//
	89	aDiscretize=35;
	90	aDeflection=0.01;
	91	//
	92	BOPDS_VectorOfInterfEF& aEFs=myDS->InterfEF();
	93	aEFs.SetStartSize(iSize);
	94	aEFs.SetIncrement(iSize);
	95	aEFs.Init();
	96	//
	97	for (; myIterator->More(); myIterator->Next()) {
	98	myIterator->Value(nE, nF, bJustAdd);
	99	if(bJustAdd) {
	100	continue;
	101	}
	102	//
	103	const BOPDS_ShapeInfo& aSIE=myDS->ShapeInfo(nE);
	104	if (aSIE.HasFlag()){//degenerated
	105	continue;
	106	}
	107	//
	108	const TopoDS_Edge& aE=((TopoDS_Edge )(&aSIE.Shape()));
	109	const TopoDS_Face& aF=((TopoDS_Face )(&myDS->Shape(nF)));
	110	const Bnd_Box& aBBF=myDS->ShapeInfo(nF).Box();
	111	//
	112	BOPDS_FaceInfo& aFI=myDS->ChangeFaceInfo(nF);
	113	const BOPDS_IndexedMapOfPaveBlock& aMPBF=aFI.PaveBlocksOn();
	114	const BOPCol_MapOfInteger& aMIFOn=aFI.VerticesOn();
	115	const BOPCol_MapOfInteger& aMIFIn=aFI.VerticesIn();
	116	//
	117	aTolE=BRep_Tool::Tolerance(aE);
	118	aTolF=BRep_Tool::Tolerance(aF);
	119	//
	120	BOPDS_ListOfPaveBlock& aLPB=myDS->ChangePaveBlocks(nE);
	121	aIt.Initialize(aLPB);
	122	for (; aIt.More(); aIt.Next()) {
	123	Handle(BOPDS_PaveBlock)& aPB=aIt.ChangeValue();
	124	//
5a77460e	125	const Handle(BOPDS_PaveBlock) aPBR=myDS->RealPaveBlock(aPB);
4e57c75e	126	if (aMPBF.Contains(aPBR)) {
	127	continue;
	128	}
	129	//
	130	if (!aPB->HasShrunkData()) {
	131	FillShrunkData(aPB);
	132	if (myWarningStatus) {
	133	continue;
	134	}
	135	}
	136	//
	137	Bnd_Box aBBE;
	138	aPB->ShrunkData(aTS1, aTS2, aBBE);
	139	//
	140	if (aBBF.IsOut (aBBE)) {
	141	continue;
	142	}
	143	//
	144	// -----------f
	145	IntTools_EdgeFace aEdgeFace;
	146	//
	147	aEdgeFace.SetEdge (aE);
	148	aEdgeFace.SetFace (aF);
	149	aEdgeFace.SetTolE (aTolE);
	150	aEdgeFace.SetTolF (aTolF);
	151	aEdgeFace.SetDiscretize (aDiscretize);
	152	aEdgeFace.SetDeflection (aDeflection);
	153	aEdgeFace.SetContext(myContext);
	154	//
	155	IntTools_Range aSR(aTS1, aTS2);
	156	IntTools_Range anewSR=aSR;
4e57c75e	157	BOPTools_AlgoTools::CorrectRange(aE, aF, aSR, anewSR);
b4109929	158	//
	159	aPB->Range(aT1, aT2);
	160	IntTools_Range aPBRange(aT1, aT2);
	161	aSR = aPBRange;
	162	BOPTools_AlgoTools::CorrectRange(aE, aF, aSR, aPBRange);
	163	//
	164	aEdgeFace.SetRange (aPBRange);
4e57c75e	165	//
	166	aEdgeFace.Perform();
	167	if (!aEdgeFace.IsDone()) {
	168	continue;
	169	}
	170	//
b4109929	171	aPB->Indices(nV[0], nV[1]);
b4109929	172	//
4e57c75e	173	const IntTools_SequenceOfCommonPrts& aCPrts=aEdgeFace.CommonParts();
	174	aNbCPrts=aCPrts.Length();
	175	for (i=1; i<=aNbCPrts; ++i) {
	176	const IntTools_CommonPrt& aCPart=aCPrts(i);
	177	aType=aCPart.Type();
	178	switch (aType) {
	179	case TopAbs_VERTEX: {
b4109929	180	Standard_Boolean bIsOnPave[2];
b4109929	181	Standard_Integer j;
4e57c75e	182	Standard_Real aT, aTolToDecide;
	183	TopoDS_Vertex aVnew;
	184
	185	BOPInt_Tools::VertexParameter(aCPart, aT);
	186	BOPTools_AlgoTools::MakeNewVertex(aE, aT, aF, aVnew);
	187	//
	188	const IntTools_Range& aR=aCPart.Range1();
	189	aTolToDecide=5.e-8;
4e57c75e	190	//
b4109929	191	IntTools_Range aR1(aT1, anewSR.First()), aR2(anewSR.Last(), aT2);
4e57c75e	192	//
b4109929	193	bIsOnPave[0]=BOPInt_Tools::IsInRange(aR1, aR, aTolToDecide);
	194	bIsOnPave[1]=BOPInt_Tools::IsInRange(aR2, aR, aTolToDecide);
	195	//
	196	if (bIsOnPave[0] && bIsOnPave[1]) {
	197	bV[0]=CheckFacePaves(nV[0], aMIFOn, aMIFIn);
	198	bV[1]=CheckFacePaves(nV[1], aMIFOn, aMIFIn);
	199	if (bV[0] && bV[1]) {
4e57c75e	200	iX=aEFs.Append()-1;
	201	IntTools_CommonPrt aCP = aCPart;
	202	aCP.SetType(TopAbs_EDGE);
	203	BOPDS_InterfEF& aEF=aEFs(iX);
	204	aEF.SetIndices(nE, nF);
	205	aEF.SetCommonPart(aCP);
	206	myDS->AddInterf(nE, nF);
	207	// 3
	208	BOPAlgo_Tools::FillMap(aPB, nF, aMPBLI, aAllocator);
	209	break;
	210	}
	211	}
b4109929	212	for (j=0; j<2; ++j) {
	213	if (bIsOnPave[j]) {
	214	bV[j]=CheckFacePaves(nV[j], aMIFOn, aMIFIn);
	215	if (bV[j]) {
	216	const TopoDS_Vertex& aV = ((TopoDS_Vertex )(&myDS->Shape(nV[j])));
	217	BOPTools_AlgoTools::UpdateVertex(aE, aT, aV);
	218	BOPDS_ShapeInfo& aSIDS=myDS->ChangeShapeInfo(nV[j]);
	219	Bnd_Box& aBoxDS=aSIDS.ChangeBox();
	220	BRepBndLib::Add(aV, aBoxDS);
	221	}
	222	else {
	223	bIsOnPave[j] = ForceInterfVF(nV[j], nF);
	224	}
4e57c75e	225	}
4e57c75e	226	}
4e57c75e	227	//
b4109929	228	if (!bIsOnPave[0] && !bIsOnPave[1]) {
b4109929	229	if (CheckFacePaves(aVnew, aMIFOn)) {
4e57c75e	230	continue;
4e57c75e	231	}
b4109929	232	//
b4109929	233	const gp_Pnt& aPnew = BRep_Tool::Pnt(aVnew);
7eed5d29	234	if (!myContext->IsValidPointForFace(aPnew, aF, aTolE+aTolF)) {
4e57c75e	235	continue;
4e57c75e	236	}
b4109929	237	//
4e57c75e	238	aBB.UpdateVertex(aVnew, aTolE);
	239	//
	240	aMIEFC.Add(nF);
	241	// 1
	242	iX=aEFs.Append()-1;
	243	BOPDS_InterfEF& aEF=aEFs(iX);
	244	aEF.SetIndices(nE, nF);
	245	aEF.SetCommonPart(aCPart);
	246	// 2
	247	myDS->AddInterf(nE, nF);
	248	// 3
	249	BOPDS_CoupleOfPaveBlocks aCPB;
	250	//
	251	aCPB.SetPaveBlocks(aPB, aPB);
	252	aCPB.SetIndexInterf(iX);
	253	aMVCPB.Add(aVnew, aCPB);
	254	}
	255	}
	256	break;
	257	case TopAbs_EDGE: {
	258	aMIEFC.Add(nF);
	259	//
	260	// 1
	261	iX=aEFs.Append()-1;
	262	BOPDS_InterfEF& aEF=aEFs(iX);
	263	aEF.SetIndices(nE, nF);
	264	//
b4109929	265	bV[0]=CheckFacePaves(nV[0], aMIFOn, aMIFIn);
	266	bV[1]=CheckFacePaves(nV[1], aMIFOn, aMIFIn);
	267	if (!bV[0] \|\| !bV[1]) {
4e57c75e	268	myDS->AddInterf(nE, nF);
	269	break;
	270	}
	271	//update tolerance of edge if needed
	272	if (aTolE < aTolF) {
	273	myDS->UpdateEdgeTolerance(nE, aTolF);
	274	aTolE = aTolF;
	275	}
	276	aEF.SetCommonPart(aCPart);
	277	// 2
	278	myDS->AddInterf(nE, nF);
	279	// 3
	280	BOPAlgo_Tools::FillMap(aPB, nF, aMPBLI, aAllocator);
	281
	282	}
	283	break;
	284	default:
	285	break;
	286	}//switch (aType) {
	287	}//for (i=1; i<=aNbCPrts; ++i) {
	288	// -----------t
	289	}//for (; aIt.More(); aIt.Next()) {
	290	}//for (; myIterator->More(); myIterator->Next()) {
	291	//
	292	//=========================================
	293	// post treatment
	294	//=========================================
5a77460e	295	BOPAlgo_Tools::PerformCommonBlocks(aMPBLI, aAllocator, myDS);
4e57c75e	296	PerformVerticesEF(aMVCPB, aAllocator);
	297	//
	298	// Update FaceInfoIn for all faces having EF common parts
	299	BOPCol_MapIteratorOfMapOfInteger aItMI;
	300	aItMI.Initialize(aMIEFC);
	301	for (; aItMI.More(); aItMI.Next()) {
	302	nF=aItMI.Value();
	303	myDS->UpdateFaceInfoIn(nF);
	304	}
	305	// Refine FaceInfoOn to remove all formal pave blocks
	306	// made during EF processing
	307	//myDS->RefineFaceInfoOn();
	308	//-----------------------------------------------------scope t
	309	aMIEFC.Clear();
	310	aMVCPB.Clear();
	311	aMPBLI.Clear();
	312	aAllocator.Nullify();
	313	//
	314
	315	}
	316	//=======================================================================
	317	//function : PerformVertices1
	318	//purpose :
	319	//=======================================================================
	320	Standard_Integer BOPAlgo_PaveFiller::PerformVerticesEF
	321	(BOPDS_IndexedDataMapOfShapeCoupleOfPaveBlocks& theMVCPB,
	322	Handle(NCollection_BaseAllocator)& theAllocator)
	323	{
	324	Standard_Integer aNbV, iRet;
	325	//
	326	iRet=0;
	327	aNbV=theMVCPB.Extent();
	328	if (!aNbV) {
	329	return iRet;
	330	}
	331	//
	332	Standard_Integer nVx, nVSD, iV, iErr, nE, iFlag, iX, i, aNbPBLI;
	333	Standard_Real aT;
	334	TopoDS_Shape aV;
	335	BOPCol_ListIteratorOfListOfShape aItLS;
	336	BOPCol_ListIteratorOfListOfInteger aItLI;
	337	BOPDS_PDS aPDS;
	338	BOPDS_ShapeInfo aSI;
	339	BOPDS_Pave aPave;
	340	//
	341	BOPCol_ListOfShape aLS(theAllocator);
	342	BOPCol_DataMapOfShapeInteger aMVI(100, theAllocator);
	343	BOPDS_IndexedDataMapOfPaveBlockListOfInteger aMPBLI(100, theAllocator);
	344	BOPAlgo_PaveFiller aPF(theAllocator);
	345	//
	346	aSI.SetShapeType(TopAbs_VERTEX);
	347	BOPDS_VectorOfInterfEF& aEFs=myDS->InterfEF();
	348	//
	349	// 1 prepare arguments
	350	for (i=1; i<=aNbV; ++i) {
	351	const TopoDS_Shape& aS=theMVCPB.FindKey(i);
	352	aLS.Append(aS);
	353	}
	354	//
	355	// 2 Fuse vertices
	356	aPF.SetArguments(aLS);
	357	aPF.Perform();
	358	iErr=aPF.ErrorStatus();
	359	if (iErr) {
360	iRet=1;
361	return iRet;
362	}
363	aPDS=aPF.PDS();
364	//
365	// 3 Add new vertices to theDS;
366	// 4 Map PaveBlock/ListOfVertices to add to this PaveBlock ->aMPBLI
367	aItLS.Initialize(aLS);
368	for (; aItLS.More(); aItLS.Next()) {
369	const TopoDS_Shape& aVx=aItLS.Value();
370	nVx=aPDS->Index(aVx);
371	//
372	if (aPDS->HasShapeSD(nVx, nVSD)) {
373	aV=aPDS->Shape(nVSD);
374	}
375	else {
376	aV=aVx;
377	}
378	// index of new vertex in theDS -> iV
379	if (!aMVI.IsBound(aV)) {
380	aSI.SetShape(aV);
381	iV=myDS->Append(aSI);
382	//
383	BOPDS_ShapeInfo& aSIDS=myDS->ChangeShapeInfo(iV);
384	Bnd_Box& aBox=aSIDS.ChangeBox();
385	BRepBndLib::Add(aV, aBox);
386	//
387	aMVI.Bind(aV, iV);
388	}
389	else {
390	iV=aMVI.Find(aV);
391	}
392	//
393	BOPDS_CoupleOfPaveBlocks &aCPB=theMVCPB.ChangeFromKey(aVx);
394	aCPB.SetIndex(iV);
395	// update EF interference
396	iX=aCPB.IndexInterf();
397	BOPDS_InterfEF& aEF=aEFs(iX);
398	aEF.SetIndexNew(iV);
399	// map aMPBLI
400	const Handle(BOPDS_PaveBlock)& aPB=aCPB.PaveBlock1();
401	if (aMPBLI.Contains(aPB)) {
402	BOPCol_ListOfInteger& aLI=aMPBLI.ChangeFromKey(aPB);
403	aLI.Append(iV);
404	}
405	else {
406	BOPCol_ListOfInteger aLI(theAllocator);
407	aLI.Append(iV);
408	aMPBLI.Add(aPB, aLI);
409	}
410	}
411	//
412	// 5
413	// 5.1 Compute Extra Paves and
414	// 5.2. Add Extra Paves to the PaveBlocks
415	aNbPBLI=aMPBLI.Extent();
416	for (i=1; i<=aNbPBLI; ++i) {
417	Handle(BOPDS_PaveBlock) aPB=aMPBLI.FindKey(i);
418	const BOPCol_ListOfInteger& aLI=aMPBLI.FindFromIndex(i);
419	nE=aPB->OriginalEdge();
420	const TopoDS_Edge& aE=((TopoDS_Edge )(&myDS->Shape(nE)));
421	//
422	aItLI.Initialize(aLI);
423	for (; aItLI.More(); aItLI.Next()) {
424	nVx=aItLI.Value();
425	const TopoDS_Vertex& aVx=((TopoDS_Vertex )(&myDS->Shape(nVx)));
426	//
427	iFlag=myContext->ComputeVE (aVx, aE, aT);
428	if (!iFlag) {
429	aPave.SetIndex(nVx);
430	aPave.SetParameter(aT);
431	aPB->AppendExtPave(aPave);
432	}
433	}
434	}
435	// 6 Split PaveBlocksa
436	for (i=1; i<=aNbPBLI; ++i) {
437	Handle(BOPDS_PaveBlock) aPB=aMPBLI.FindKey(i);
438	nE=aPB->OriginalEdge();
439	// 3
5a77460e	440	if (!myDS->IsCommonBlock(aPB)) {
4e57c75e	441	myDS->UpdatePaveBlock(aPB);
	442	}
	443	else {
5a77460e	444	const Handle(BOPDS_CommonBlock)& aCB=myDS->CommonBlock(aPB);
4e57c75e	445	myDS->UpdateCommonBlock(aCB);
	446	}
	447	}//for (; aItMPBLI.More(); aItMPBLI.Next()) {
	448	//
	449	return iRet;
	450	}
	451	//=======================================================================
	452	// function: CheckFacePaves
	453	// purpose:
	454	//=======================================================================
	455	Standard_Boolean BOPAlgo_PaveFiller::CheckFacePaves (const Standard_Integer nVx,
	456	const BOPCol_MapOfInteger& aMIFOn,
	457	const BOPCol_MapOfInteger& aMIFIn)
	458	{
	459	Standard_Boolean bRet;
	460	Standard_Integer nV;
	461	BOPCol_MapIteratorOfMapOfInteger aIt;
	462	//
	463	bRet=Standard_False;
	464	//
	465	aIt.Initialize(aMIFOn);
	466	for (; aIt.More(); aIt.Next()) {
	467	nV=aIt.Value();
	468	if (nV==nVx) {
	469	bRet=!bRet;
	470	return bRet;
	471	}
	472	}
	473	aIt.Initialize(aMIFIn);
	474	for (; aIt.More(); aIt.Next()) {
	475	nV=aIt.Value();
	476	if (nV==nVx) {
	477	bRet=!bRet;
	478	return bRet;
	479	}
	480	}
	481	//
	482	return bRet;
	483	}
	484	//=======================================================================
	485	// function: CheckFacePaves
	486	// purpose:
	487	//=======================================================================
	488	Standard_Boolean BOPAlgo_PaveFiller::CheckFacePaves (const TopoDS_Vertex& aVnew,
	489	const BOPCol_MapOfInteger& aMIF)
	490	{
	491	Standard_Boolean bRet;
	492	Standard_Integer nV, iFlag;
	493	BOPCol_MapIteratorOfMapOfInteger aIt;
	494	//
	495	bRet=Standard_True;
	496	//
	497	aIt.Initialize(aMIF);
	498	for (; aIt.More(); aIt.Next()) {
	499	nV=aIt.Value();
	500	const TopoDS_Vertex& aV=((TopoDS_Vertex )(&myDS->Shape(nV)));
	501	iFlag=BOPTools_AlgoTools::ComputeVV(aVnew, aV);
	502	if (!iFlag) {
	503	return bRet;
	504	}
	505	}
	506	//
	507	return !bRet;
	508	}
b4109929	509	//=======================================================================
	510	//function : ForceInterfVF
	511	//purpose :
	512	//=======================================================================
	513	Standard_Boolean BOPAlgo_PaveFiller::ForceInterfVF(const Standard_Integer nV,
	514	const Standard_Integer nF)
	515	{
	516	Standard_Boolean bRet;
	517	//
	518	bRet = Standard_False;
	519	const TopoDS_Vertex& aV = (TopoDS_Vertex)&myDS->Shape(nV);
	520	const TopoDS_Face& aF = (TopoDS_Face) &myDS->Shape(nF);
	521	//
	522	GeomAPI_ProjectPointOnSurf& aProj = myContext->ProjPS(aF);
	523	const gp_Pnt& aP = BRep_Tool::Pnt(aV);
	524	aProj.Perform(aP);
	525	if (!aProj.IsDone()) {
	526	return bRet;
	527	}
	528	Standard_Real aDist, U, V;
	529	//
	530	aDist=aProj.LowerDistance();
	531	aProj.LowerDistanceParameters(U, V);
	532	//
	533	gp_Pnt2d aP2d(U, V);
	534	bRet = myContext->IsPointInFace (aF, aP2d);
	535	if (bRet) {
	536	Standard_Integer i;
	537	BRep_Builder aBB;
	538	//
	539	BOPDS_VectorOfInterfVF& aVFs=myDS->InterfVF();
	540	i=aVFs.Append()-1;
	541	BOPDS_InterfVF& aVF=aVFs(i);
	542	aVF.SetIndices(nV, nF);
	543	aVF.SetUV(U, V);
	544	//
	545	myDS->AddInterf(nV, nF);
	546	//
	547	aBB.UpdateVertex(aV, aDist);
	548	BOPDS_ShapeInfo& aSIDS=myDS->ChangeShapeInfo(nV);
	549	Bnd_Box& aBoxDS=aSIDS.ChangeBox();
	550	BRepBndLib::Add(aV, aBoxDS);
	551	//
	552	BOPDS_FaceInfo& aFI=myDS->ChangeFaceInfo(nF);
	553	BOPCol_MapOfInteger& aMVIn=aFI.ChangeVerticesIn();
	554	aMVIn.Add(nV);
	555	}
	556	//
	557	return bRet;
	558	}
	559