// File:	TopOpeBRepDS_ProcessInterferencesTool.cxx
// Created:	Fri Feb 14 11:18:00 1997
// Author:	Jean Yves LEBEY
//		<jyl@bistrox.paris1.matra-dtv.fr>

#include <TopOpeBRepDS_InterferenceTool.hxx>
#include <TopOpeBRepDS_ProcessInterferencesTool.hxx>
#include <TopOpeBRepDS_EdgeVertexInterference.hxx>
#include <TopOpeBRepDS_InterferenceIterator.hxx>
#include <TopOpeBRepDS_TKI.hxx>
#include <TopoDS.hxx>
#include <BRepLProp_SLProps.hxx>
#include <BRepAdaptor_Surface.hxx>
//#include <BRepAdaptor_Curve2d.hxx>
#include <BRep_Tool.hxx>
#include <gp_Pnt2d.hxx>
#include <gp_Vec.hxx>
#include <gp_Dir.hxx>
#include <Precision.hxx>
#include <TColStd_MapOfInteger.hxx>
#include <Standard_ProgramError.hxx>
#include <TopOpeBRepDS_define.hxx>
#include <TopOpeBRepTool_TOOL.hxx>
#include <TopOpeBRepTool_EXPORT.hxx>
#include <TopOpeBRepDS_EXPORT.hxx>

#ifdef DEB
Standard_IMPORT Standard_Boolean TopOpeBRepDS_GettracePEI();
Standard_IMPORT Standard_Boolean TopOpeBRepDS_GettracePI();
Standard_IMPORT Standard_Boolean TopOpeBRepDS_GettraceSPSX(const Standard_Integer);
static Standard_Boolean TRCE(const Standard_Integer EIX) {
  Standard_Boolean b1 = TopOpeBRepDS_GettracePEI();
  Standard_Boolean b2 = TopOpeBRepDS_GettracePI();
  Standard_Boolean b3 = TopOpeBRepDS_GettraceSPSX(EIX);
  return (b1 || b2 || b3);
}
Standard_EXPORT void debredudoub(const Standard_Integer /*i*/) {}
Standard_EXPORT void debredudoub1(const Standard_Integer /*i*/) {}
Standard_EXPORT void debredudoub2(const Standard_Integer /*i*/) {} 
Standard_EXPORT void debredunk(const Standard_Integer /*i*/) {}
#endif

//-----------------------------------------------------------------------
Standard_EXPORT Handle(TopOpeBRepDS_Interference) MakeCPVInterference
(const TopOpeBRepDS_Transition& T, // transition
 const Standard_Integer SI, // curve/edge index
 const Standard_Integer GI, // point/vertex index
 const Standard_Real P, // parameter of G on S
 const TopOpeBRepDS_Kind GK) // POINT/VERTEX
//-----------------------------------------------------------------------
{
  Handle(TopOpeBRepDS_Interference) I;
  TopOpeBRepDS_Kind SK = TopOpeBRepDS_CURVE;
  I = TopOpeBRepDS_InterferenceTool::MakeCurveInterference(T,SK,SI,GK,GI,P);
  return I;
}

//-----------------------------------------------------------------------
Standard_EXPORT Handle(TopOpeBRepDS_Interference) MakeEPVInterference
(const TopOpeBRepDS_Transition& T, // transition
 const Standard_Integer SI, // curve/edge index
 const Standard_Integer GI, // point/vertex index
 const Standard_Real P, // parameter of G on S
 const TopOpeBRepDS_Kind GK,
 const Standard_Boolean B) // G is a vertex (or not) of the interference master
//-----------------------------------------------------------------------
{
  Handle(TopOpeBRepDS_Interference) I;
  TopOpeBRepDS_Kind SK = TopOpeBRepDS_EDGE;
  if      ( GK == TopOpeBRepDS_POINT ) {
    I = TopOpeBRepDS_InterferenceTool::MakeEdgeInterference(T,SK,SI,GK,GI,P);
  }
  else if ( GK == TopOpeBRepDS_VERTEX ) {
    I = TopOpeBRepDS_InterferenceTool::MakeEdgeVertexInterference(T,SI,GI,B,TopOpeBRepDS_UNSHGEOMETRY,P);
    I->GeometryType(GK);
  }
  return I;
}

//-----------------------------------------------------------------------
Standard_EXPORT Handle(TopOpeBRepDS_Interference) MakeEPVInterference
(const TopOpeBRepDS_Transition& T, // transition
 const Standard_Integer S, // curve/edge index
 const Standard_Integer G, // point/vertex index
 const Standard_Real P, // parameter of G on S
 const TopOpeBRepDS_Kind GK, // POINT/VERTEX
 const TopOpeBRepDS_Kind SK,
 const Standard_Boolean B) // G is a vertex (or not) of the interference master
//-----------------------------------------------------------------------
{
  Handle(TopOpeBRepDS_Interference) I = ::MakeEPVInterference(T,S,G,P,GK,B);
  I->SupportType(SK);
  return I;
}

//------------------------------------------------------
Standard_EXPORT Standard_Boolean FUN_hasStateShape
//------------------------------------------------------
// FUN_hasStateShape : True if transition T is (state,shape) before or after
(const TopOpeBRepDS_Transition& T,
 const TopAbs_State state,
 const TopAbs_ShapeEnum shape)
{
  TopAbs_State     staB = T.Before(), staA = T.After();
  TopAbs_ShapeEnum shaB = T.ShapeBefore(), shaA = T.ShapeAfter();
  Standard_Boolean B = ( staB == state ) && ( shaB == shape);
  Standard_Boolean A = ( staA == state ) && ( shaA == shape);
  Standard_Boolean result = B || A;
  return result;
}

//------------------------------------------------------
Standard_EXPORT Standard_Integer FUN_selectTRASHAinterference
//------------------------------------------------------
// selection des interf (sha2) de L1 dans L2. retourne leur nombre.
(TopOpeBRepDS_ListOfInterference& L1,const TopAbs_ShapeEnum sha2, TopOpeBRepDS_ListOfInterference& L2)
{
  TopOpeBRepDS_ListIteratorOfListOfInterference it1(L1);
  while (it1.More()) {
    Handle(TopOpeBRepDS_Interference) I1 = it1.Value();
    const TopOpeBRepDS_Transition& T1 = I1->Transition();

    TopAbs_ShapeEnum shab = T1.ShapeBefore(), shaa = T1.ShapeAfter();
    Standard_Boolean sel=((shab==sha2)||(shaa==sha2));
    if (sel) {
      L2.Append(I1);
      L1.Remove(it1);
    }
    else it1.Next();
  }
  Standard_Integer n2 = L2.Extent();
  return n2;
}

//------------------------------------------------------
Standard_EXPORT Standard_Integer FUN_selectITRASHAinterference
//------------------------------------------------------
// selection des interf (sha2) de L1 dans L2. retourne leur nombre.
(TopOpeBRepDS_ListOfInterference& L1,const Standard_Integer Index, TopOpeBRepDS_ListOfInterference& L2)
{
  if (Index == 0) return 0;
  TopOpeBRepDS_ListIteratorOfListOfInterference it1(L1);
  while (it1.More()) {
    Handle(TopOpeBRepDS_Interference) I1 = it1.Value();
    const TopOpeBRepDS_Transition& T1 = I1->Transition();
    Standard_Integer Ind = T1.Index();
    if (Ind == Index) {
      L2.Append(I1);
      L1.Remove(it1);
    }
    else it1.Next();
  }
  Standard_Integer n2 = L2.Extent();
  return n2;
}

//------------------------------------------------------
Standard_EXPORT Standard_Integer FUN_selectTRAUNKinterference
//------------------------------------------------------
// selection des interf de transition (UNK,UNK) de L1 dans L2. retourne leur nombre.
(TopOpeBRepDS_ListOfInterference& L1,TopOpeBRepDS_ListOfInterference& L2)
{
  TopOpeBRepDS_ListIteratorOfListOfInterference it1(L1);
  while (it1.More()) {
    Handle(TopOpeBRepDS_Interference) I1 = it1.Value();
    const TopOpeBRepDS_Transition& T1 = I1->Transition();
    Standard_Boolean sel = T1.IsUnknown();
    if (sel) {
      L2.Append(I1);
      L1.Remove(it1);
    }
    else it1.Next();
  }
  Standard_Integer n2 = L2.Extent();
  return n2;
}

//------------------------------------------------------
Standard_EXPORT Standard_Integer FUN_selectTRAORIinterference
//------------------------------------------------------
// selection des interf d'orientation donnee O de L1 dans L2. retourne leur nombre.
(TopOpeBRepDS_ListOfInterference& L1, const TopAbs_Orientation O, TopOpeBRepDS_ListOfInterference& L2)
{
  TopOpeBRepDS_ListIteratorOfListOfInterference it1(L1);
  while (it1.More()) {
    Handle(TopOpeBRepDS_Interference) I1 = it1.Value();
    const TopOpeBRepDS_Transition& T1 = I1->Transition();
    TopAbs_Orientation ori = T1.Orientation(TopAbs_IN);
    Standard_Boolean sel = (ori ==O);
    if (sel) {
      L2.Append(I1);
      L1.Remove(it1);
    }
    else it1.Next();
  }
  Standard_Integer n2 = L2.Extent();
  return n2;
}

//------------------------------------------------------
Standard_EXPORT Standard_Integer FUN_selectGKinterference
//------------------------------------------------------
(TopOpeBRepDS_ListOfInterference& L1,const TopOpeBRepDS_Kind GK,TopOpeBRepDS_ListOfInterference& L2)
{
  TopOpeBRepDS_ListIteratorOfListOfInterference it1(L1);
  while (it1.More()) {
    TopOpeBRepDS_Kind gki = it1.Value()->GeometryType();
    if (gki == GK) {
      L2.Append(it1.Value());
      L1.Remove(it1);
    }
    else it1.Next();
  }
  Standard_Integer n2 = L2.Extent();
  return n2;
}

//------------------------------------------------------
Standard_EXPORT Standard_Integer FUN_selectSKinterference
//------------------------------------------------------
// selection des interf de type de support SK de L1 dans L2. retourne leur nombre.
(TopOpeBRepDS_ListOfInterference& L1,const TopOpeBRepDS_Kind SK,TopOpeBRepDS_ListOfInterference& L2)
{
  TopOpeBRepDS_ListIteratorOfListOfInterference it1(L1);
  while (it1.More()) {
    TopOpeBRepDS_Kind ski = it1.Value()->SupportType();
    if (ski == SK) {
      L2.Append(it1.Value());
      L1.Remove(it1);
    }
    else it1.Next();
  }
  Standard_Integer n2 = L2.Extent();
  return n2;
}

//------------------------------------------------------
Standard_EXPORT Standard_Integer FUN_selectGIinterference
//------------------------------------------------------
// selection des interf d' index de geometrie GI de L1 dans L2. retourne leur nombre.
(TopOpeBRepDS_ListOfInterference& L1,const Standard_Integer GI,TopOpeBRepDS_ListOfInterference& L2)
{
  if (GI == 0) return 0;
  TopOpeBRepDS_ListIteratorOfListOfInterference it1(L1);
  while (it1.More()) {
    Standard_Integer gi = it1.Value()->Geometry();
    if (gi == GI) {
      L2.Append(it1.Value());
      L1.Remove(it1);
    }
    else it1.Next();
  }
  Standard_Integer n2 = L2.Extent();
  return n2;
}

//------------------------------------------------------
Standard_EXPORT Standard_Integer FUN_selectSIinterference
//------------------------------------------------------
// selection des interf d' index de support SI de L1 dans L2. retourne leur nombre.
(TopOpeBRepDS_ListOfInterference& L1,const Standard_Integer SI,TopOpeBRepDS_ListOfInterference& L2)
{
  if (SI == 0) return 0;
  TopOpeBRepDS_ListIteratorOfListOfInterference it1(L1);
  while (it1.More()) {
    Standard_Integer si = it1.Value()->Support();
    if (si == SI) {
      L2.Append(it1.Value());
      L1.Remove(it1);
    }
    else it1.Next();
  }
  Standard_Integer n2 = L2.Extent();
  return n2;
}

//------------------------------------------------------
Standard_EXPORT Standard_Boolean FUN_interfhassupport
//------------------------------------------------------
// FUN_interfhassupport : True si shape support de l 'interf I est IsSame(S).
(const TopOpeBRepDS_DataStructure& DS,const Handle(TopOpeBRepDS_Interference)& I,const TopoDS_Shape& S)
{
  Standard_Boolean h = Standard_True;
  const Standard_Integer index = I->Support();
  const TopoDS_Shape& SofI = DS.Shape(index);
  if ( SofI.IsSame(S) ) h =  Standard_True;
  else                  h =  Standard_False;
  return h;
}

//------------------------------------------------------
Standard_EXPORT Standard_Boolean FUN_transitionEQUAL
//------------------------------------------------------
(const TopOpeBRepDS_Transition& T1,const TopOpeBRepDS_Transition& T2)
{
  Standard_Boolean id1 = FUN_transitionSTATEEQUAL(T1,T2);
  Standard_Boolean id2 = FUN_transitionSHAPEEQUAL(T1,T2);
  Standard_Boolean id3 = FUN_transitionINDEXEQUAL(T1,T2);
  Standard_Boolean id = id1 && id2 && id3;
  return id;
}

//------------------------------------------------------
Standard_EXPORT Standard_Boolean FUN_transitionSTATEEQUAL
//------------------------------------------------------
(const TopOpeBRepDS_Transition& T1,const TopOpeBRepDS_Transition& T2)
{
  Standard_Boolean id = T1.Before()==T2.Before() && T1.After()==T2.After();
  return id;
}

//------------------------------------------------------
Standard_EXPORT Standard_Boolean FUN_transitionSHAPEEQUAL
//------------------------------------------------------
(const TopOpeBRepDS_Transition& T1,const TopOpeBRepDS_Transition& T2)
{
  Standard_Boolean id = T1.ShapeBefore()==T2.ShapeBefore() && T1.ShapeAfter()==T2.ShapeAfter();
  return id;
}

//------------------------------------------------------
Standard_EXPORT Standard_Boolean FUN_transitionINDEXEQUAL
//------------------------------------------------------
(const TopOpeBRepDS_Transition& T1,const TopOpeBRepDS_Transition& T2)
{
  Standard_Boolean id = T1.IndexBefore()==T2.IndexBefore() && T1.IndexAfter()==T2.IndexAfter();
  return id;
}

//------------------------------------------------------
Standard_EXPORT void FUN_reducedoublons
//------------------------------------------------------
(TopOpeBRepDS_ListOfInterference& LI,const TopOpeBRepDS_DataStructure& BDS,const Standard_Integer SIX)    
{
#ifdef DEB
  Standard_Boolean TRC=TRCE(SIX); if (TRC) debredudoub(SIX);
  Standard_Boolean modif = Standard_False;
#endif
  const TopoDS_Shape& E = BDS.Shape(SIX);
  TopOpeBRepDS_ListIteratorOfListOfInterference it1;
  
  // process interferences of LI with VERTEX geometry
  
  it1.Initialize(LI);
  while (it1.More() ) {
    Handle(TopOpeBRepDS_Interference)& I1 = it1.Value();
    const TopOpeBRepDS_Transition& T1 = I1->Transition();
    TopOpeBRepDS_Kind GT1,ST1; Standard_Integer G1,S1; FDS_data(I1,GT1,G1,ST1,S1);
    TopAbs_ShapeEnum tsb1,tsa1; Standard_Integer isb1,isa1;
    FDS_Tdata(I1,tsb1,isb1,tsa1,isa1);
    
#ifdef DEB
    if (TRC) debredudoub1(SIX);
#endif
    
    TopOpeBRepDS_ListIteratorOfListOfInterference it2(it1); it2.Next();
    while ( it2.More() ) {

      const Handle(TopOpeBRepDS_Interference)& I2 = it2.Value();
      const TopOpeBRepDS_Transition& T2 = I2->Transition();
      TopOpeBRepDS_Kind GT2,ST2; Standard_Integer G2,S2; FDS_data(I2,GT2,G2,ST2,S2);
      TopAbs_ShapeEnum tsb2,tsa2; Standard_Integer isb2,isa2;
      FDS_Tdata(I2,tsb2,isb2,tsa2,isa2);
      
#ifdef DEB
      if (TRC) debredudoub2(SIX);
#endif
      
      Standard_Boolean idGS = (GT2 == GT1 && G2 == G1 && ST2 == ST1 && S2 == S1);
      if (idGS) {

	Standard_Boolean id1 = FUN_transitionSTATEEQUAL(T1,T2);
	Standard_Boolean id2 = FUN_transitionSHAPEEQUAL(T1,T2);
	Standard_Boolean id3 = FUN_transitionINDEXEQUAL(T1,T2);
	Standard_Boolean idT = id1 && id2 && id3;
	
	if ( idT ) { 
	  const Handle(TopOpeBRepDS_EdgeVertexInterference) EVI1 = Handle(TopOpeBRepDS_EdgeVertexInterference)::DownCast(I1);
	  const Handle(TopOpeBRepDS_EdgeVertexInterference) EVI2 = Handle(TopOpeBRepDS_EdgeVertexInterference)::DownCast(I2);
	  Standard_Boolean evi = (!EVI1.IsNull()) && (!EVI2.IsNull());
	  if (evi) {
	    //xpu130898 : if vG is on E's closing vertex -> filter
	    const TopoDS_Vertex& vG = TopoDS::Vertex(BDS.Shape(G1));
	    TopoDS_Vertex OOv;  Standard_Boolean hasOO = FUN_ds_getoov(vG,BDS,OOv);
	    TopoDS_Vertex vclo; Standard_Boolean Eclosed = TopOpeBRepTool_TOOL::ClosedE(TopoDS::Edge(E),vclo);
	    Standard_Boolean onvclo = Standard_False;
	    if (Eclosed) {
	      onvclo = vG.IsSame(vG);
	      if (hasOO && !onvclo) onvclo = vG.IsSame(OOv);
	    }
	    if (onvclo) idT = Standard_True;
	    else {
	      //xpu100697 : if interf are EVI compare parameters
	      Standard_Real tolE = FUN_tool_maxtol(E); Standard_Real t = Precision::Parametric(tolE);
	      Standard_Real t1 = EVI1->Parameter(); Standard_Real t2 = EVI2->Parameter();
	      Standard_Real dd = t1-t2;
	      Standard_Boolean samepar = (Abs(dd) <= t);
	      idT = samepar;
	    }
	  } // evi
	}	
	if (idT) {
	  // les 2 interferences I1 et I2 sont identiques : on supprime I2
#ifdef DEB
	  if(TRC){cout<<"shape "<<SIX<<" : doublon ";I1->Dump(cout);cout<<endl;}
	  modif = Standard_True;
#endif
	  LI.Remove(it2);
	}
	else it2.Next();
      }
      else it2.Next();
    }
    it1.Next();
  } // it1.More()
  
#ifdef DEB
  if(TRC && !modif){ FDS_dumpLI(LI,"sans doublon : "); }
#endif  
} // reducedoublons

//------------------------------------------------------
Standard_EXPORT void FUN_unkeepUNKNOWN
//------------------------------------------------------
  (TopOpeBRepDS_ListOfInterference& LI,
   TopOpeBRepDS_DataStructure& /*BDS*/,
#ifdef DEB
   const Standard_Integer SIX)    
#else
   const Standard_Integer)    
#endif
{
#ifdef DEB
  Standard_Boolean TRC=TRCE(SIX); if (TRC) debredunk(SIX);
  Standard_Boolean modif = Standard_False;
#endif
//                 BDS.Shape(SIX);
  TopOpeBRepDS_ListIteratorOfListOfInterference it1;
  
  // process interferences of LI with UNKNOWN transition
  
  it1.Initialize(LI);
  while (it1.More() ) {
    Handle(TopOpeBRepDS_Interference)& I1 = it1.Value();
    const TopOpeBRepDS_Transition& T1 = I1->Transition();
    Standard_Boolean isunk = T1.IsUnknown();

    if (isunk) {
#ifdef DEB
      if(TRC) {
	if (isunk) debredunk(SIX);
	cout<<"shape "<<SIX<<" : UNKNOWN transition ";I1->Dump(cout);cout<<endl;
      }
#endif
      LI.Remove(it1);
    }
    else it1.Next();
  } // it1.More()
  
#ifdef DEB
  if(TRC && !modif){ FDS_dumpLI(LI,"sans UNKNOWN : "); }
#endif  
} // unkeepUNKNOWN

// -----------------------------------------------------------
static Standard_Integer FUN_select3dI(const Standard_Integer SIX, TopOpeBRepDS_DataStructure& BDS,
			 TopOpeBRepDS_ListOfInterference& lFE, TopOpeBRepDS_ListOfInterference& lFEresi, TopOpeBRepDS_ListOfInterference& l3dFE)
//------------------------------------------------------
{
  l3dFE.Clear();
  lFEresi.Clear();
  Standard_Integer n3d = 0;
#ifdef DEB
  Standard_Boolean TRC=TRCE(SIX); 
#endif

  Standard_Integer nFE = lFE.Extent();
  if (nFE <= 1) return n3d;

  const TopoDS_Edge& E = TopoDS::Edge(BDS.Shape(SIX));
  Standard_Integer rankE  = BDS.AncestorRank(E);
  TopoDS_Shape OOv; Standard_Integer OOG = 0, Gsta = 0;

  // attached to a given edge :
  // <lFE> = {I=(T(face),G=POINT,S=EDGE)}/ {I=(T(face),G=VERTEX,S=EDGE)}
  //             = set of interferences of same geometry kind.
  //    -> <l3dFE> + <lFE>
  // <l3dFE> = sets of interferences with same geometry,
  //           describing 3d complex transitions.

  TopOpeBRepDS_ListIteratorOfListOfInterference it1(lFE);
  while (it1.More()) {

    Standard_Boolean complex3d=Standard_False;
    const Handle(TopOpeBRepDS_Interference)& I1 = it1.Value();
    TopOpeBRepDS_Kind GT1,ST1; Standard_Integer G1,S1; FDS_data(I1,GT1,G1,ST1,S1);
    TopAbs_ShapeEnum SB1,SA1; Standard_Integer IB1,IA1; FDS_Tdata(I1,SB1,IB1,SA1,IA1);
    Standard_Boolean vertex1 = (GT1 == TopOpeBRepDS_VERTEX);
    
    TopOpeBRepDS_ListIteratorOfListOfInterference it2(it1);
    if (it2.More()) it2.Next();
    else break;

    // <Gsta>, <OOv>
    if (vertex1) {
      TopoDS_Vertex vG1 = TopoDS::Vertex(BDS.Shape(G1)); Standard_Integer rankvG1 = BDS.AncestorRank(vG1);
      Standard_Boolean G1hsd = FUN_ds_getVsdm(BDS,G1,OOG);	
      if (rankvG1 != rankE) // vG1 not on E
	{OOv = vG1; Gsta = G1hsd? 3: 2;}
      else // vG on E
	{if (G1hsd) OOv = BDS.Shape(OOG); Gsta = G1hsd? 3: 1;}
    }

    while (it2.More()){
      const Handle(TopOpeBRepDS_Interference)& I2 = it2.Value();
      TopOpeBRepDS_Kind GT2,ST2; Standard_Integer G2,S2; FDS_data(I2,GT2,G2,ST2,S2);
      TopAbs_ShapeEnum SB2,SA2; Standard_Integer IB2,IA2; FDS_Tdata(I2,SB2,IB2,SA2,IA2);
	
      Standard_Boolean sameG = (GT2 == GT1);
      if (!sameG) break;
      sameG = (G2 == G1) || (G2 == OOG);

      const TopoDS_Face& F1 = TopoDS::Face(BDS.Shape(IB1));
      const TopoDS_Face& F2 = TopoDS::Face(BDS.Shape(IB2));

      // same domain faces, -> 2d interferences
      Standard_Boolean sdmFT = (IB1==IB2) || FUN_ds_sdm(BDS,F1,F2);
//      if (sdmFT) {it2.Next(); continue;}      
      if (sdmFT) {lFEresi.Append(I2); lFE.Remove(it2); continue;}  

      // a. if (S2 == S1) ->ok
      // b. else :
      //     if (G is POINT on Fi i = 1,2) ->continue
      //     else : look for E1(edge S1) on F1(IB1) shared by F2(IB2)
      //            with bound G.
      Standard_Boolean sameS = (ST2 == ST1) && (S2 == S1);

      Standard_Boolean sameSorsharedEbyTRASHA = sameS;
      Standard_Boolean hasOOv = (Gsta > 1);
      if (!sameSorsharedEbyTRASHA && hasOOv) {
	TopoDS_Shape Eshared; Standard_Boolean foundsh = FUN_tool_Eshared(OOv,F1,F2,Eshared);
	if (!foundsh) {it2.Next(); continue;}

	sameSorsharedEbyTRASHA = Standard_True;
	if (!BDS.HasShape(Eshared)) 
	  {Standard_Integer OOrank = BDS.AncestorRank(OOv); BDS.AddShape(Eshared,OOrank);}
	S1 = BDS.Shape(Eshared); S2 = S1;
#ifdef DEB
	if (TRC) {cout<<"    finding out E shared by F"<<IB1<<" and F"<<IB2;
		  cout<<" new support is E"<<S1<<endl;}
#endif
      } 	

      if (sameSorsharedEbyTRASHA) { // xpu : 09-03-98
	Standard_Boolean sdm = FUN_ds_sdm(BDS,BDS.Shape(SIX),BDS.Shape(S1));
	if (sdm) {
#ifdef DEB
	  if (TRC) cout<<"  NO I3d : e"<<SIX<<" same domain with e"<<S1<<endl;
#endif
	  it2.Next(); continue;
	}
      } // xpu : 09-03-98

      if (sameSorsharedEbyTRASHA) {l3dFE.Append(I2); lFE.Remove(it2); complex3d=Standard_True;}
      else                        it2.Next();
    } // it2(it1);

    if (complex3d) {l3dFE.Append(I1); lFE.Remove(it1);}
    else           it1.Next();
  } // it1(lFE);

  // xpu170898 : cto009H1 (I1(T1(F),G,S1),I1'(T2(F),G,S2), I2)
  //             (I1,I2) -> to reduce => I1' to delete
  n3d = l3dFE.Extent();
  Standard_Integer nresi = lFEresi.Extent();
  if (nresi == 0) return n3d;

  if      (n3d != 0) {
    TopOpeBRepDS_ListIteratorOfListOfInterference it(lFEresi);
    while (it.More()){
      const Handle(TopOpeBRepDS_Interference)& I = it.Value();
      Standard_Integer TRA = I->Transition().Index();
      TopOpeBRepDS_ListOfInterference lcopy; FDS_assign(l3dFE,lcopy);
      TopOpeBRepDS_ListOfInterference lfound; Standard_Integer nfound = FUN_selectITRASHAinterference(lcopy,TRA,lfound);
      if (nfound == 0) lFE.Remove(it);
      else             it.Next();
    }
  }
  else lFE.Append(lFEresi);
  return n3d;
}

static Standard_Boolean FUN_find3dISEsameISF(const Handle(TopOpeBRepDS_Interference)& I1, TopOpeBRepDS_ListIteratorOfListOfInterference& it2)
// I1 = (T0,G0,S)
// looking among interferences of <it2> for I2/
// I2 = (T0,G0,S')     
{
  TopAbs_Orientation O1 = I1->Transition().Orientation(TopAbs_IN);
  TopOpeBRepDS_Kind GT1,ST1; Standard_Integer G1,S1; FDS_data(I1,GT1,G1,ST1,S1);
  TopAbs_ShapeEnum SB1,SA1; Standard_Integer IB1,IA1; FDS_Tdata(I1,SB1,IB1,SA1,IA1);
  
  while ( it2.More()){
    Handle(TopOpeBRepDS_Interference)& I2 = it2.Value();
    TopAbs_Orientation O2 = I2->Transition().Orientation(TopAbs_IN);
    TopOpeBRepDS_Kind GT2,ST2; Standard_Integer G2,S2; FDS_data(I2,GT2,G2,ST2,S2);
    TopAbs_ShapeEnum SB2,SA2; Standard_Integer IB2,IA2; FDS_Tdata(I2,SB2,IB2,SA2,IA2);
    Standard_Boolean sameT = (SB1 == SB2) && (IB1 == IB2) && (O1 == O2);
    Standard_Boolean sameG = (GT1 == GT2) && (G1 == G2); 
    if (sameT && sameG) return Standard_True;
    else                it2.Next();
  } // it2
  return Standard_False;
}

// ------------------------------------------------------------
static Standard_Integer FUN_select3dISEsameISF
(TopOpeBRepDS_ListOfInterference& lFE, TopOpeBRepDS_ListOfInterference& l3dFE, TopOpeBRepDS_ListOfInterference& l3dFEresi,
 TopOpeBRepDS_ListOfInterference& lF, TopOpeBRepDS_ListOfInterference& l3dF)
//-------------------------------------------------------------
{
  //    I3d in <l3dFE>, I3d = (T(face),G=POINT/VERTEX,S=EDGE)
  // <->I3d' in <l3dF>,  I3d'= (T(face),G=POINT/VERTEX,S=FACE)

  // <l3dFEold> -> {I3d} + <l3dFEnew>
  // <lFEold>   -> <lFE> + (<l3dFEold> - <l3dFEnew>)
  // <lFold>    -> {I3d'}+ <lFnew>;    {I3d'}=<l3dF>

  TopOpeBRepDS_ListIteratorOfListOfInterference it1(l3dFE);
  while (it1.More()) {
    const Handle(TopOpeBRepDS_Interference)& I1 = it1.Value();
//    TopAbs_Orientation O1 = I1->Transition().Orientation(TopAbs_IN);
//    TopOpeBRepDS_Kind GT1,ST1; Standard_Integer G1,S1; FDS_data(I1,GT1,G1,ST1,S1);
//    TopAbs_ShapeEnum SB1,SA1; Standard_Integer IB1,IA1; FDS_Tdata(I1,SB1,IB1,SA1,IA1);
    
//    Standard_Boolean found = Standard_False;
//    TopOpeBRepDS_ListIteratorOfListOfInterference it2(lF);
//    while ( it2.More()){
//      Handle(TopOpeBRepDS_Interference)& I2 = it2.Value();
//      TopAbs_Orientation O2 = I2->Transition().Orientation(TopAbs_IN);
//      TopOpeBRepDS_Kind GT2,ST2; Standard_Integer G2,S2; FDS_data(I2,GT2,G2,ST2,S2);
//      TopAbs_ShapeEnum SB2,SA2; Standard_Integer IB2,IA2; FDS_Tdata(I2,SB2,IB2,SA2,IA2);
//      Standard_Boolean sameT = (SB1 == SB2) && (IB1 == IB2) && (O1 == O2);
//      Standard_Boolean sameG = (GT1 == GT2) && (G1 == G2); 
//      if (sameT && sameG) {found = Standard_True; l3dF.Append(I2); lF.Remove(it2);break;}
//      else                it2.Next();
//    } // it2
//    if (found) it1.Next();
//    else       {lFE.Append(I1); l3dFE.Remove(it1);}

    // {IFE = (T0,G0,Sedge); IFF = (T0,G0,Sface)}
    // => l3dFE += IFE, l3dF += iFF
    TopOpeBRepDS_ListIteratorOfListOfInterference it2(lF);
    Standard_Boolean found = FUN_find3dISEsameISF(I1,it2);

    if (found) {
      l3dF.Append(it2.Value()); 
      lF.Remove(it2); 
      it1.Next();
    }
    else       {
      // xpu : 16-02-98 :      
      // IFE1 = (T0,G0,Sedge1), IFE2 = (T0,G0,Sedge2)
      // IFF = (T0,G0,Sface)
      // => l3dFE += IFE1, l3dF += iFF, 
      //    l3dFEresi += IFE2
      TopOpeBRepDS_ListIteratorOfListOfInterference it3(l3dF);
      found = FUN_find3dISEsameISF(I1,it3);
      
      if (found) l3dFEresi.Append(I1);
      else       lFE.Append(I1);
      l3dFE.Remove(it1);

    }
  } // it1
  return l3dF.Extent();
}

//------------------------------------------------------
Standard_EXPORT Standard_Integer FUN_select2dI
(const Standard_Integer SIX, TopOpeBRepDS_DataStructure& BDS,const TopAbs_ShapeEnum TRASHAk,TopOpeBRepDS_ListOfInterference& lI,TopOpeBRepDS_ListOfInterference& l2dI)
//------------------------------------------------------
// <lI> -> <l2dI> + <lI>
// TRASHAk = TopAbs_FACE : <l2dI> = {I2d = (T on same domain faces, G,S)} 
// TRASHAk = TopAbs_EDGE : <l2dI> = {I2d = (T on same edge, G,S)} 
{
  l2dI.Clear();
  Standard_Integer n2d = 0;
  Standard_Integer nFE = lI.Extent();
  if (nFE <=1) return n2d;
  
  Standard_Boolean TonFace = (TRASHAk == TopAbs_FACE);

  // xpu201098 : check FEI is not 3dFEI (cto904F6,SIX10,G6)
  TColStd_MapOfInteger mapftra;
  TopOpeBRepDS_ListOfInterference lIE; FDS_copy(BDS.ShapeInterferences(SIX),lIE);
  TopOpeBRepDS_ListOfInterference l3dF;
#ifdef DEB
//  Standard_Integer n3d =
#endif
            FUN_selectSKinterference(lIE,TopOpeBRepDS_FACE,l3dF);
  for (TopOpeBRepDS_ListIteratorOfListOfInterference itt(l3dF); itt.More(); itt.Next()) mapftra.Add(itt.Value()->Support());
  TopOpeBRepDS_ListOfInterference lII; TopOpeBRepDS_ListIteratorOfListOfInterference it1(lI);
  while (it1.More()) {
    const Handle(TopOpeBRepDS_Interference)& I = it1.Value();
    Standard_Integer ITRASHA = I->Transition().Index();
    if (mapftra.Contains(ITRASHA)) {it1.Next();}
    else                           {lII.Append(I);lI.Remove(it1);}
  }

  it1.Initialize(lII);
  while (it1.More()) {
    const Handle(TopOpeBRepDS_Interference)& I1 = it1.Value();
    TopOpeBRepDS_Kind GT1,ST1;  Standard_Integer G1,S1;   FDS_data(I1,GT1,G1,ST1,S1);
    TopAbs_ShapeEnum SB1,SA1; Standard_Integer IB1,IA1; FDS_Tdata(I1,SB1,IB1,SA1,IA1);

    if (SB1 != TRASHAk) {it1.Next(); continue;}
    
    Standard_Boolean complex2d = Standard_False;
    TopOpeBRepDS_ListIteratorOfListOfInterference it2(it1);
    if (it2.More()) it2.Next();
    else break;
    while (it2.More()){
      const Handle(TopOpeBRepDS_Interference)& I2 = it2.Value();
      TopOpeBRepDS_Kind GT2,ST2;  Standard_Integer G2,S2;   FDS_data(I2,GT2,G2,ST2,S2);
      TopAbs_ShapeEnum SB2,SA2; Standard_Integer IB2,IA2; FDS_Tdata(I2,SB2,IB2,SA2,IA2);

      Standard_Boolean cond = (SB1 == SB2) && (IB1 == IB2) && (IA1 == IA2);
      if (!cond) {it2.Next(); continue;}
      
      Standard_Boolean sameG = (GT2 == GT1) && (G2 == G1);
      Standard_Boolean sameST = (ST2 == ST1);
      cond = sameG && sameST;
      if (!cond) {it2.Next(); continue;}

      complex2d = (IB1==IB2);
      if (TonFace) {
	const TopoDS_Face& F1 = TopoDS::Face(BDS.Shape(IB1));
	const TopoDS_Face& F2 = TopoDS::Face(BDS.Shape(IB2));
	complex2d = complex2d || FUN_ds_sdm(BDS,F1,F2);
      }
      if (complex2d) {l2dI.Append(I2); lII.Remove(it2);}
      else           it2.Next(); 
    } // it2
    if (complex2d) {l2dI.Append(I1); lII.Remove(it1);}
    else           it1.Next();
  } // it1
  lI.Append(lII);
  n2d = l2dI.Extent();
  return n2d;
}

//------------------------------------------------------
Standard_EXPORT Standard_Integer FUN_selectpure2dI
(const TopOpeBRepDS_ListOfInterference& lF, TopOpeBRepDS_ListOfInterference& lFE, TopOpeBRepDS_ListOfInterference& l2dFE)
//------------------------------------------------------
// <lFE> -> <lFE> + <l2dFE>
// <l2dFE> = {I2dFE = (T(face),G,S=edge) / we cannot find (T(face),G,face)}
{
  l2dFE.Clear();
  TopOpeBRepDS_ListIteratorOfListOfInterference itFE(lFE);
  while (itFE.More()) {
    const Handle(TopOpeBRepDS_Interference)& IFE = itFE.Value();
    const TopOpeBRepDS_Transition& TFE = IFE->Transition();
    Standard_Integer IB = TFE.IndexBefore(), IA = TFE.IndexAfter();
    if (IB != IA) {itFE.Next(); continue;}

    Standard_Boolean foundIF = Standard_False;
    for (TopOpeBRepDS_ListIteratorOfListOfInterference itF(lF); itF.More(); itF.Next()){
      const Handle(TopOpeBRepDS_Interference)& IF = itF.Value();
      Standard_Integer S = IF->Support();      
      if (S == IB) {foundIF = Standard_True; break;}
    }
    if (foundIF) {itFE.Next(); continue;}

    l2dFE.Append(IFE); lFE.Remove(itFE);
  }
  Standard_Integer n2dFE = l2dFE.Extent();
  return n2dFE;
} // selectpure2dI

//------------------------------------------------------
Standard_EXPORT Standard_Integer FUN_select1dI(const Standard_Integer SIX, TopOpeBRepDS_DataStructure& BDS,TopOpeBRepDS_ListOfInterference& LI,TopOpeBRepDS_ListOfInterference& l1dI)
//------------------------------------------------------
// attached to edge E : <lI>={I=(T(ES),G,ES)} -> <l1dI> + <lI>
// l1dI contains {(I1,I2) / I1=(T1(ES1),VG,ES1), I2=(T2(ES2),VG,ES2), 
//                          E sdm ES1, ES2}     
{
  l1dI.Clear();
  Standard_Integer n1d = 0;
  Standard_Integer nFE = LI.Extent();
  if (nFE <=1) return n1d;

  TopOpeBRepDS_ListOfInterference newLI;
  const TopoDS_Shape& EIX = BDS.Shape(SIX);
  TopOpeBRepDS_TKI tki; tki.FillOnGeometry(LI);
  for (tki.Init(); tki.More(); tki.Next()) {
    TopOpeBRepDS_Kind K; Standard_Integer G; tki.Value(K,G);
    TopOpeBRepDS_ListOfInterference& loi = tki.ChangeValue(K,G); TopOpeBRepDS_ListOfInterference Rloi;
    Standard_Integer nloi = loi.Extent();

    Standard_Boolean ok =   (K == TopOpeBRepDS_VERTEX);
    ok = ok && (nloi > 1);
    if (!ok) {newLI.Append(loi); continue;}

    TopOpeBRepDS_ListIteratorOfListOfInterference it1(loi);
    while (it1.More()) {
      const Handle(TopOpeBRepDS_Interference)& I1 = it1.Value();
      TopOpeBRepDS_Kind GT1,ST1;  Standard_Integer G1,S1;   FDS_data(I1,GT1,G1,ST1,S1);
      TopAbs_ShapeEnum SB1,SA1; Standard_Integer IB1,IA1; FDS_Tdata(I1,SB1,IB1,SA1,IA1);

      Standard_Boolean cond1 = (SB1 == TopAbs_EDGE) && (IB1 == IA1);      
      if (!cond1) {newLI.Append(I1); it1.Next(); continue;}
      cond1 = FUN_ds_sdm(BDS,EIX,BDS.Shape(S1));    
      if (!cond1) {newLI.Append(I1); it1.Next(); continue;}
        
      Standard_Boolean complex1d = Standard_False;
      TopOpeBRepDS_ListIteratorOfListOfInterference it2(it1);
      if (it2.More()) it2.Next();
      else break;
      while (it2.More()){
	const Handle(TopOpeBRepDS_Interference)& I2 = it2.Value();
	TopOpeBRepDS_Kind GT2,ST2;  Standard_Integer G2,S2;   FDS_data(I2,GT2,G2,ST2,S2);
	TopAbs_ShapeEnum SB2,SA2; Standard_Integer IB2,IA2; FDS_Tdata(I2,SB2,IB2,SA2,IA2);

	Standard_Boolean cond2 = (SB1 == SB2) && (IB2 == IA2);
	if (!cond2) {newLI.Append(I2); it2.Next(); continue;}
	complex1d = FUN_ds_sdm(BDS,EIX,BDS.Shape(S2)); 
	if (!complex1d) {newLI.Append(I2); it2.Next(); continue;}
//                 LI.Extent(); //deb
	l1dI.Append(I2); it2.Next();   
      } // it2
      if (complex1d) {l1dI.Append(I1); it1.Next();}
      else           it1.Next();
    } // it1
  } // tki
  LI.Clear(); LI.Append(newLI);
  n1d = l1dI.Extent();
  return n1d;
} // select1dI

//------------------------------------------------------
Standard_EXPORT void FUN_select3dinterference
(const Standard_Integer SIX, TopOpeBRepDS_DataStructure& BDS,
 TopOpeBRepDS_ListOfInterference& lF, TopOpeBRepDS_ListOfInterference& l3dF, 
 TopOpeBRepDS_ListOfInterference& lFE, TopOpeBRepDS_ListOfInterference& lFEresi, TopOpeBRepDS_ListOfInterference& l3dFE, TopOpeBRepDS_ListOfInterference& l3dFEresi,
 TopOpeBRepDS_ListOfInterference& l2dFE)  
//------------------------------------------------------
{
  const TopoDS_Edge& E = TopoDS::Edge(BDS.Shape(SIX));
  Standard_Boolean isdg = BRep_Tool::Degenerated(E);
  if (isdg) return;

  l3dF.Clear(); l3dFE.Clear();
  // attached to an edge the list of interferences :
  // IN  : <lF>  = {I=(T(face),G=POINT/VERTEX,S=FACE)}
  //       <lFE> = {I=(T(face),G=POINT/VERTEX,S=EDGE)}

  // OUT : <lFE> -> <l3dFE> = {I describing 3d interferences on S=EDGE}
  //                   is the list of 3d interferences to reduce
  //                +<lFE> +lFEresi
  //       <lF>  -> <l3dF> = {I describing 3d interferences on S=FACE}
  //                +<lF>
  //     I  = (T(face),G=POINT/VERTEX,S=EDGE) in <l3dFE>
  // <-> I' = (T(face),G=POINT/VERTEX,S=FACE) in <l3dF>

#ifdef DEB
//  Standard_Integer n3dFE =
#endif
              ::FUN_select3dI(SIX,BDS,lFE,lFEresi,l3dFE);
#ifdef DEB
//  Standard_Integer n3dFF =
#endif
              ::FUN_select3dISEsameISF(lFE,l3dFE,l3dFEresi,lF,l3dF);
#ifdef DEB
//  Standard_Integer n2dFE =
#endif
              FUN_select2dI(SIX,BDS,TopAbs_FACE,lFE,l2dFE); 
}
