OCCT

// Created on: 1994-03-30
// Created by: Laurent BUCHARD
// Copyright (c) 1994-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.

#define MARCHEPASSIUNESEULEFACE 0

#ifdef OCCT_DEBUG
  #define LBRCOMPT 0
#else
  #define LBRCOMPT 0
#endif

#if LBRCOMPT
  #include <stdio.h>

class StatistiquesBRepClass3d
{
public:
  long unsigned NbConstrVide;
  long unsigned NbLoad;
  long unsigned NbConstrShape;
  long unsigned NbConstrShapePnt;
  long unsigned NbPerform;
  long unsigned NbPerformRejection;
  long unsigned NbPerformInfinitePoint;
  long unsigned NbDestroy;

public:
  StatistiquesBRepClass3d()
  {
    NbConstrVide = NbLoad = NbConstrShape = NbConstrShapePnt = NbPerform = NbPerformInfinitePoint =
      NbDestroy                                                          = 0;
  }

  ~StatistiquesBRepClass3d()
  {
    printf("\n--- Statistiques BRepClass3d:\n");

    printf("\nConstructeurVide    : %10lu", NbConstrVide);
    printf("\nConstructeur(Shape) : %10lu", NbConstrShape);
    printf("\nLoad(Shape)         : %10lu", NbLoad);
    printf("\nPerform(pnt3d)      : %10lu", NbPerform);
    printf("\nPerform(pnt3d) REJ  : %10lu", NbPerformRejection);
    printf("\nPerformInfinitePoint: %10lu", NbPerformInfinitePoint);
    printf("\nDestroy             : %10lu", NbDestroy);
  }
};

static StatistiquesBRepClass3d STAT;
#endif

#include <BRepClass3d_SolidClassifier.hxx>
#include <gp_Pnt.hxx>
#include <TopoDS_Shape.hxx>

BRepClass3d_SolidClassifier::BRepClass3d_SolidClassifier()
{
  aSolidLoaded = isaholeinspace = false;
#if LBRCOMPT
  STAT.NbConstrVide++;
#endif
}

void BRepClass3d_SolidClassifier::Load(const TopoDS_Shape& S)
{

#if LBRCOMPT
  STAT.NbLoad++;
#endif

  if (aSolidLoaded)
  {
    explorer.Destroy();
  }
  explorer.InitShape(S);
  aSolidLoaded = true;

#if MARCHEPASSIUNESEULEFACE
  PerformInfinitePoint(1e-7);
  if (State() == TopAbs_OUT)
  {
    isaholeinspace = false;
  }
  else
  {
    isaholeinspace = true;
  }
#endif
}

BRepClass3d_SolidClassifier::BRepClass3d_SolidClassifier(const TopoDS_Shape& S)
    : aSolidLoaded(true),
      explorer(S),
      isaholeinspace(false)
{
#if LBRCOMPT
  STAT.NbConstrShape++;
#endif
#if MARCHEPASSIUNESEULEFACE
  PerformInfinitePoint(1e-7);
  if (State() == TopAbs_OUT)
  {
    isaholeinspace = false;
  }
  else
  {
    isaholeinspace = true;
  }
#endif
}

BRepClass3d_SolidClassifier::BRepClass3d_SolidClassifier(const TopoDS_Shape& S,
                                                         const gp_Pnt&       P,
                                                         const double        Tol)
    : explorer(S)
{
#if LBRCOMPT
  STAT.NbConstrShapePnt++;
#endif
  aSolidLoaded = true;
#if MARCHEPASSIUNESEULEFACE
  PerformInfinitePoint(1e-7);
  if (State() == TopAbs_OUT)
  {
    isaholeinspace = false;
  }
  else
  {
    isaholeinspace = true;
  }

  if (isaholeinspace == false)
  {
    if (explorer.Box().IsOut(P))
    {
      ForceOut();
    }
    else
    {
      Perform(P, Tol);
    }
  }
  else
  {
    if (explorer.Box().IsOut(P))
    {
      ForceIn();
    }
    else
    {
      Perform(P, Tol);
    }
  }
#else
  Perform(P, Tol);
#endif
}

void BRepClass3d_SolidClassifier::Perform(const gp_Pnt& P, const double Tol)
{
#if LBRCOMPT
  STAT.NbPerform++;
#endif
#if MARCHEPASSIUNESEULEFACE
  if (aSolidLoaded)
  {
    if (isaholeinspace == false)
    {
      if (explorer.Box().IsOut(P))
      {
        ForceOut();
  #if LBRCOMPT
        STAT.NbPerformRejection++;
  #endif
      }
      else
      {
        BRepClass3d_SClassifier::Perform(explorer, P, Tol);
      }
    }
    else
    {
      if (explorer.Box().IsOut(P))
      {
        ForceIn();
  #if LBRCOMPT
        STAT.NbPerformRejection++;
  #endif
      }
      else
      {
        BRepClass3d_SClassifier::Perform(explorer, P, Tol);
      }
    }
  }
#else
  BRepClass3d_SClassifier::Perform(explorer, P, Tol);
#endif
}

void BRepClass3d_SolidClassifier::PerformInfinitePoint(const double Tol)
{
#if LBRCOMPT
  STAT.NbPerformInfinitePoint++;
#endif

  if (aSolidLoaded)
  {
    BRepClass3d_SClassifier::PerformInfinitePoint(explorer, Tol);
    isaholeinspace = State() != TopAbs_OUT;
  }
}

void BRepClass3d_SolidClassifier::Destroy()
{
#if LBRCOMPT
  STAT.NbDestroy++;
#endif

  if (aSolidLoaded)
  {
    explorer.Destroy();
    aSolidLoaded = false;
  }
}