OCCT

// Created on: 1993-07-08
// Created by: Jean Yves LEBEY
// 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.

#include <BRep_Tool.hxx>
#include <BRepBndLib.hxx>
#include <gp_Pnt.hxx>
#include <Precision.hxx>
#include <Standard_ProgramError.hxx>
#include <Standard_Type.hxx>
#include <TopExp_Explorer.hxx>
#include <TopoDS.hxx>
#include <TopoDS_Shape.hxx>
#include <TopOpeBRepTool_box.hxx>
#include <TopAbs_ShapeEnum.hxx>
#include <TopAbs_Orientation.hxx>
#include <TopAbs_State.hxx>
#include <TopTools_ShapeMapHasher.hxx>
#include <NCollection_Map.hxx>
#include <NCollection_List.hxx>
#include <NCollection_IndexedMap.hxx>
#include <NCollection_DataMap.hxx>
#include <Standard_Integer.hxx>
#include <NCollection_IndexedDataMap.hxx>
#include <TopoDS_Face.hxx>
#include <TopoDS_Edge.hxx>
#include <TopoDS_Vertex.hxx>
#include <TCollection_AsciiString.hxx>
#include <TopOpeBRepTool_HBoxTool.hxx>

IMPLEMENT_STANDARD_RTTIEXT(TopOpeBRepTool_HBoxTool, Standard_Transient)

#ifdef OCCT_DEBUG
  #define TBOX TopOpeBRepTool_GettraceBOX()
#endif

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

TopOpeBRepTool_HBoxTool::TopOpeBRepTool_HBoxTool() = default;

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

void TopOpeBRepTool_HBoxTool::Clear()
{
  myIMS.Clear();
}

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

void TopOpeBRepTool_HBoxTool::AddBoxes(const TopoDS_Shape&    S,
                                       const TopAbs_ShapeEnum TS,
                                       const TopAbs_ShapeEnum TA)
{
  TopExp_Explorer ex;
  for (ex.Init(S, TS, TA); ex.More(); ex.Next())
  {
    //  for (TopExp_Explorer ex(S,TS,TA);ex.More();ex.Next()) {
    const TopoDS_Shape& ss = ex.Current();
    AddBox(ss);
  }
}

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

void TopOpeBRepTool_HBoxTool::AddBox(const TopoDS_Shape& S)
{
#ifdef OCCT_DEBUG
  TopAbs_ShapeEnum t =
#endif
    S.ShapeType();
  bool hasb = HasBox(S);
  if (hasb)
  {
    return;
  }

  Bnd_Box B;
  ComputeBox(S, B);
  myIMS.Add(S, B);
#ifdef OCCT_DEBUG
  if (TBOX)
  {
    DumpB(Box(S));
    std::cout << "; # HBT::AddBox " << Index(S) << " : ";
    TopAbs::Print(t, std::cout);
    std::cout << "\n";
    std::cout << std::endl;
  }
#endif
}

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

void TopOpeBRepTool_HBoxTool::ComputeBox(const TopoDS_Shape& S, Bnd_Box& B)
{
  TopAbs_ShapeEnum t = S.ShapeType();
  if (t == TopAbs_FACE)
  {
    BRepBndLib::Add(S, B);
  }
  else if (t == TopAbs_EDGE)
  {
    BRepBndLib::Add(S, B);
    // modified by NIZHNY-MZV  Tue Sep 21 14:04:33 1999
  }
  else if (t == TopAbs_SOLID)
  {
    BRepBndLib::Add(S, B);
  }
  else if (t == TopAbs_COMPOUND)
  {
    BRepBndLib::Add(S, B);
    // modified by NIZHNY-MZV  Wed Apr  5 10:05:53 2000
  }
  else if (t == TopAbs_VERTEX)
  {
    BRepBndLib::Add(S, B);
  }
  else
  {
    throw Standard_ProgramError("HBT::ComputeBox : invalid type");
  }
}

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

void TopOpeBRepTool_HBoxTool::ComputeBoxOnVertices(const TopoDS_Shape& S, Bnd_Box& B)
{
  TopExp_Explorer ex(S, TopAbs_VERTEX);
  if (!ex.More())
  {
    B.Update(-1.e5, -1.e5, -1.e5, 1.e5, 1.e5, 1.e5);
    return;
  }
  double tol = Precision::Confusion();
  for (; ex.More(); ex.Next())
  {
    double x, y, z;
    BRep_Tool::Pnt(TopoDS::Vertex(ex.Current())).Coord(x, y, z);
    B.Update(x, y, z);
    tol = std::max(tol, BRep_Tool::Tolerance(TopoDS::Vertex(ex.Current())));
  }
  B.Enlarge(tol);
}

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

const Bnd_Box& TopOpeBRepTool_HBoxTool::Box(const TopoDS_Shape& S)
{
  bool hb = HasBox(S);
  if (!hb)
  {
    throw Standard_ProgramError("HBT::Box1");
  }

  const Bnd_Box& B = myIMS.FindFromKey(S);
  return B;
}

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

const Bnd_Box& TopOpeBRepTool_HBoxTool::Box(const int I) const
{
  int iu = Extent();
  int hb = (I >= 1 && I <= iu);
  if (!hb)
  {
    throw Standard_ProgramError("HBT::Box2");
  }
  const Bnd_Box& B = myIMS.FindFromIndex(I);
  return B;
}

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

bool TopOpeBRepTool_HBoxTool::HasBox(const TopoDS_Shape& S) const
{
  bool hb = myIMS.Contains(S);
  return hb;
}

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

const TopoDS_Shape& TopOpeBRepTool_HBoxTool::Shape(const int I) const
{
  int iu = Extent();
  int hs = (I >= 1 && I <= iu);
  if (!hs)
  {
    throw Standard_ProgramError("HBT::Box4");
  }
  const TopoDS_Shape& S = myIMS.FindKey(I);
  return S;
}

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

int TopOpeBRepTool_HBoxTool::Index(const TopoDS_Shape& S) const
{
  int i = myIMS.FindIndex(S);
  return i;
}

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

int TopOpeBRepTool_HBoxTool::Extent() const
{
  int n = myIMS.Extent();
  return n;
}

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

NCollection_IndexedDataMap<TopoDS_Shape, Bnd_Box>& TopOpeBRepTool_HBoxTool::ChangeIMS()
{
  return myIMS;
}

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

const NCollection_IndexedDataMap<TopoDS_Shape, Bnd_Box>& TopOpeBRepTool_HBoxTool::IMS() const
{
  return myIMS;
}

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

void TopOpeBRepTool_HBoxTool::DumpB
#ifdef OCCT_DEBUG
  (const Bnd_Box& B)
{
  if (B.IsVoid())
    std::cout << "# IsVoid";
  else if (B.IsWhole())
    std::cout << "# IsWhole";
  else
  {
    double x, y, z, X, Y, Z;
    B.Get(x, y, z, X, Y, Z);
    std::cout << "bounding " << x << " " << y << " " << z << " " << X << " " << Y << " " << Z;
    std::cout.flush();
  }
}
#else
  (const Bnd_Box&)
{
}
#endif