// Copyright (C) 2007-2012  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 library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
// License as published by the Free Software Foundation; either
// version 2.1 of the License.
//
// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
//
// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
//

// File:        GEOMAlgo_BuilderTools.cxx
// Author:      Peter KURNEV

#include <GEOMAlgo_BuilderTools.hxx>

#include <Basics_OCCTVersion.hxx>

#include <TColStd_Array1OfReal.hxx>

#include <gp_Pnt2d.hxx>
#include <gp_Vec.hxx>
#include <gp_Dir.hxx>
#include <gp_Pnt.hxx>

#include <TColgp_Array1OfPnt.hxx>

#include <Poly_Triangulation.hxx>
#include <Poly_Array1OfTriangle.hxx>
#include <Poly_Triangle.hxx>

#include <Geom2d_Curve.hxx>
#include <Geom2dInt_Geom2dCurveTool.hxx>

#include <TopLoc_Location.hxx>
#include <TopAbs_Orientation.hxx>

#include <TopoDS_Face.hxx>
#include <TopoDS_Iterator.hxx>
#include <TopoDS_Wire.hxx>
#include <TopoDS.hxx>
#include <TopoDS_Edge.hxx>
#include <TopExp_Explorer.hxx>

#if OCC_VERSION_LARGE > 0x06050100 // for OCC-6.5.2 and higher version
#include <TopExp.hxx>
#include <TopTools_IndexedDataMapOfShapeListOfShape.hxx>
#endif

#include <BRep_Tool.hxx>
#include <BRepBndLib.hxx>
#include <BRepMesh_FastDiscret.hxx>
#include <Bnd_Box.hxx>
#include <BRepAdaptor_Curve2d.hxx>

static
  Standard_Integer ComputeProps(const TopoDS_Face& aF,
                                Standard_Real& aA,
                                Standard_Real& aV);
static
  void BuildTriangulation(const TopoDS_Face& aF);

//=======================================================================
//function : IsHole
//purpose  :
//=======================================================================
  Standard_Boolean GEOMAlgo_BuilderTools::IsHole(const TopoDS_Shape& aW,
                                                 const TopoDS_Shape& aFace)
{
  Standard_Boolean bIsHole;
  Standard_Integer i, aNbS;
  Standard_Real aT1, aT2, aS;
  Standard_Real aU1, aU2, aU, dU;
  Standard_Real aX1, aY1, aX0, aY0;
  TopAbs_Orientation aOr;

  gp_Pnt2d aP2D0, aP2D1;
  Handle(Geom2d_Curve) aC2D;
  TopoDS_Face aF, aFF;
  TopoDS_Iterator aItW;
  //
  bIsHole=Standard_False;
  //
  aF=TopoDS::Face(aFace);
  aFF=aF;
  aFF.Orientation(TopAbs_FORWARD);
  //
  aS=0.;
  aItW.Initialize(aW);
  for (; aItW.More(); aItW.Next()) {
    const TopoDS_Edge& aE=TopoDS::Edge(aItW.Value());
    aOr=aE.Orientation();
    if (!(aOr==TopAbs_FORWARD ||
          aOr==TopAbs_REVERSED)) {
      continue;
    }
    //
    aC2D=BRep_Tool::CurveOnSurface(aE, aFF, aT1, aT2);
    if (aC2D.IsNull()) {
      break; //xx
    }
    //
    BRepAdaptor_Curve2d aBAC2D(aE, aFF);
    aNbS=Geom2dInt_Geom2dCurveTool::NbSamples(aBAC2D);
    if (aNbS>2) {
      aNbS*=4;
    }
    //
    dU=(aT2-aT1)/(Standard_Real)(aNbS-1);
    aU =aT1;
    aU1=aT1;
    aU2=aT2;
    if (aOr==TopAbs_REVERSED) {
      aU =aT2;
      aU1=aT2;
      aU2=aT1;
      dU=-dU;
    }
    //
    aC2D->D0(aU, aP2D0);
    for(i=2; i<=aNbS; i++) {
      aU=aU1+(i-1)*dU;
      aC2D->D0(aU, aP2D1);
      aP2D0.Coord(aX0, aY0);
      aP2D1.Coord(aX1, aY1);
      //
      aS=aS+(aY0+aY1)*(aX1-aX0);
      //
      aP2D0=aP2D1;
    }
  }//for (; aItW.More(); aItW.Next()) {
  bIsHole=(aS>0.);
  return bIsHole;
}
//=======================================================================
//function : IsHole
//purpose  :
//=======================================================================
  Standard_Boolean GEOMAlgo_BuilderTools::IsHole(const TopoDS_Shape& aShell)
{
  Standard_Boolean bIsHole;
  Standard_Integer iRet;
  Standard_Real aAi, aA, aV, aVi;
  TopExp_Explorer aExp;
  //
  aA=0.;
  aV=0.;
  aExp.Init(aShell, TopAbs_FACE);
  for (; aExp.More(); aExp.Next()) {
    const TopoDS_Face& aF=TopoDS::Face(aExp.Current());
    iRet=ComputeProps(aF, aAi, aVi);
    if (!iRet) {
      aA+=aAi;
      aV+=aVi;
    }
  }
  //
  bIsHole=aV<0.;
  return bIsHole;
}
//=======================================================================
//function : ComputeProps
//purpose  :
//=======================================================================
Standard_Integer ComputeProps(const TopoDS_Face& aF,
                              Standard_Real& aA,
                              Standard_Real& aV)
{
  Standard_Integer j, i, i1, i2, aNbNodes, aNbTrigs, n[3];
  Standard_Real aAi, aVi;
  gp_Pnt aP[3], aGC, aGC1;
  TopLoc_Location aLoc;
  TopAbs_Orientation aOr;
  Handle(Poly_Triangulation) aTri;
  //
  aA=0.;
  aV=0.;
  //
  aTri=BRep_Tool::Triangulation(aF, aLoc);
  if(aTri.IsNull()) {
    BuildTriangulation(aF);
    aTri=BRep_Tool::Triangulation(aF, aLoc);
    if(aTri.IsNull()) {
      return 1;// a face is without triangulation
    }
  }
  //
  aNbNodes=aTri->NbNodes();
  aNbTrigs=aTri->NbTriangles();
  if (!aNbTrigs){
    return 2;//no triangles
  }
  //
  aOr=aF.Orientation();
  //
  const TColgp_Array1OfPnt& aNodes=aTri->Nodes();
  const Poly_Array1OfTriangle& aTriangles=aTri->Triangles();
  //
  i1=aTriangles.Lower();
  i2=aTriangles.Upper();
  //
  for (i=i1; i<=i2; ++i){
    const Poly_Triangle& aTriangle=aTriangles.Value(i);
    aTriangle.Get(n[0], n[1], n[2]);
    aGC.SetCoord(0.,0.,0.);
    for (j=0; j<3; ++j) {
      aP[j]=aNodes.Value(n[j]);
      aGC.ChangeCoord()+=aP[j].XYZ();
    }
    aGC.ChangeCoord()*=0.333333333333;
    //
    // Normal
    gp_Vec aV01(aP[0], aP[1]);
    gp_Vec aV12(aP[1], aP[2]);
    gp_Vec aVN=aV01^aV12;
    aAi=aVN.Magnitude();
    aA=aA+aAi;
    //
    if (aAi>0.0000001) {
      Standard_Real aSx, aZx;
      gp_Dir aDN(aVN);
      if (aOr==TopAbs_REVERSED) {
        aDN.Reverse();
      }
      //
      aSx=aAi*aDN.Z();
      aZx=aGC.Z();
      aVi=aZx*aSx;
      aV=aV+aVi;
    }
  }
  return 0;
}
//=======================================================================
//function : BuildTriangulation
//purpose  :
//=======================================================================
void BuildTriangulation(const TopoDS_Face& aF)
{
  Standard_Boolean bWithShare;
  Standard_Real aDiscret, aXmin, aYmin, aZmin, aXmax, aYmax, aZmax;
  Standard_Real dX, dY, dZ, dMax, aCoeff, aAngle;
  Bnd_Box aBox;
  //
  bWithShare=Standard_False;
  aAngle=0.5;
  //
  BRepBndLib::Add(aF, aBox);
  //
  // aDiscret
  aBox.Get(aXmin, aYmin, aZmin, aXmax, aYmax, aZmax);
  dX=aXmax-aXmin;
  dY=aYmax-aYmin;
  dZ=aZmax-aZmin;
  dMax=dX;
  if (dY>dMax) {
    dMax=dY;
  }
  if (dZ>dMax) {
    dMax=dZ;
  }
  //
  aCoeff=0.1;
  aDiscret=aCoeff*dMax;
  //
  BRepMesh_FastDiscret aMesher(aDiscret,
                               aAngle,
                               aBox,
                               bWithShare,
                               Standard_True,
                               Standard_False,
                               Standard_True);

#if OCC_VERSION_LARGE > 0x06050100 // for OCC-6.5.2 and higher version
  TopTools_IndexedDataMapOfShapeListOfShape anAncestors;
  TopExp::MapShapesAndAncestors(aF, TopAbs_EDGE, TopAbs_FACE, anAncestors);
  aMesher.Add(aF, anAncestors);
#else
  aMesher.Add(aF);
#endif
}
