OCCT

// Created on: 1993-12-15
// Created by: Christophe MARION
// Copyright (c) 1993-1999 Matra Datavision
// 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.

// modified by MPS (june 96) : on utilise BRepClass_FaceClassifier seulement
//  si IsDone de Extrema est vrai

#include <BRepExtrema_ExtPF.hxx>

#include <BRep_Tool.hxx>
#include <BRepTools.hxx>
#include <BRepClass_FaceClassifier.hxx>
#include <gp_Pnt2d.hxx>
#include <BRepAdaptor_Surface.hxx>
#include <Precision.hxx>

//=================================================================================================

BRepExtrema_ExtPF::BRepExtrema_ExtPF(const TopoDS_Vertex&  TheVertex,
                                     const TopoDS_Face&    TheFace,
                                     const Extrema_ExtFlag TheFlag,
                                     const Extrema_ExtAlgo TheAlgo)
{
  Initialize(TheFace, TheFlag, TheAlgo);
  Perform(TheVertex, TheFace);
}

//=================================================================================================

void BRepExtrema_ExtPF::Initialize(const TopoDS_Face&    TheFace,
                                   const Extrema_ExtFlag TheFlag,
                                   const Extrema_ExtAlgo TheAlgo)
{
  // cette surface doit etre en champ. Extrema ne fait
  // pas de copie et prend seulement un pointeur dessus.
  mySurf.Initialize(TheFace, false);

  if (mySurf.GetType() == GeomAbs_OtherSurface)
  {
    return; // protect against non-geometric type (e.g. triangulation)
  }

  double Tol = std::min(BRep_Tool::Tolerance(TheFace), Precision::Confusion());
  double aTolU, aTolV;
  aTolU = std::max(mySurf.UResolution(Tol), Precision::PConfusion());
  aTolV = std::max(mySurf.VResolution(Tol), Precision::PConfusion());
  double U1, U2, V1, V2;
  BRepTools::UVBounds(TheFace, U1, U2, V1, V2);
  myExtPS.SetFlag(TheFlag);
  myExtPS.SetAlgo(TheAlgo);
  myExtPS.Initialize(mySurf, U1, U2, V1, V2, aTolU, aTolV);
}

//=================================================================================================

void BRepExtrema_ExtPF::Perform(const TopoDS_Vertex& TheVertex, const TopoDS_Face& TheFace)
{
  mySqDist.Clear();
  myPoints.Clear();

  const gp_Pnt P = BRep_Tool::Pnt(TheVertex);
  if (mySurf.GetType() == GeomAbs_OtherSurface)
  {
    return; // protect against non-geometric type (e.g. triangulation)
  }

  myExtPS.Perform(P);

  // Exploration of points and classification
  if (myExtPS.IsDone())
  {
    BRepClass_FaceClassifier classifier;
    double                   U1, U2;
    const double             Tol = BRep_Tool::Tolerance(TheFace);
    for (int i = 1; i <= myExtPS.NbExt(); i++)
    {
      myExtPS.Point(i).Parameter(U1, U2);
      const gp_Pnt2d Puv(U1, U2);
      classifier.Perform(TheFace, Puv, Tol);
      const TopAbs_State state = classifier.State();
      if (state == TopAbs_ON || state == TopAbs_IN)
      {
        mySqDist.Append(myExtPS.SquareDistance(i));
        myPoints.Append(myExtPS.Point(i));
      }
    }
  }
}