mirror of
https://git.dev.opencascade.org/repos/occt.git
synced 2025-05-16 10:54:53 +03:00
01a6e62bc2
BRepMesh_Classifier: Two-pass approach for intersection check with possibility to run it in parallel mode. First pass - bounding boxes of segments are checked for overlapping; Second pass - intersection point is calculated in case if overlapping is detected. Make NCollection_UBTree::ChangeLastNode() exported due to compilation error on Linux platform. Reason: method does not depend on template parameters, so it should be available. Revert previous change and try to use another trick for Linux Fix compilation warning on MacOS: remove redundant constant Fix regressions: do not consider insignificant loops in case of self intersections on the same wire. More sugar solution for compilation errors on NCollection_EBTree on Linux Test cases for issue CR24968
142 lines
4.1 KiB
C++
142 lines
4.1 KiB
C++
// Created on: 1997-06-26
|
|
// Created by: Laurent PAINNOT
|
|
// Copyright (c) 1997-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 <BRepMesh_Classifier.hxx>
|
|
|
|
#include <Precision.hxx>
|
|
#include <gp_Pnt2d.hxx>
|
|
#include <CSLib_Class2d.hxx>
|
|
#include <TColgp_Array1OfPnt2d.hxx>
|
|
|
|
//=======================================================================
|
|
//function : Constructor
|
|
//purpose :
|
|
//=======================================================================
|
|
BRepMesh_Classifier::BRepMesh_Classifier()
|
|
{
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : Perform
|
|
//purpose :
|
|
//=======================================================================
|
|
TopAbs_State BRepMesh_Classifier::Perform(const gp_Pnt2d& thePoint) const
|
|
{
|
|
Standard_Boolean isOut = Standard_False;
|
|
Standard_Integer aNb = myTabClass.Length();
|
|
|
|
for (Standard_Integer i = 1; i <= aNb; i++)
|
|
{
|
|
Standard_Integer aCur = ((CSLib_Class2d*)myTabClass(i))->SiDans(thePoint);
|
|
if (aCur == 0)
|
|
{
|
|
// Point is ON, but mark it as OUT
|
|
isOut = Standard_True;
|
|
}
|
|
else
|
|
isOut = myTabOrient(i) ? (aCur == -1) : (aCur == 1);
|
|
|
|
if (isOut)
|
|
return TopAbs_OUT;
|
|
}
|
|
|
|
return TopAbs_IN;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : RegisterWire
|
|
//purpose :
|
|
//=======================================================================
|
|
void BRepMesh_Classifier::RegisterWire(
|
|
const NCollection_Sequence<gp_Pnt2d>& theWire,
|
|
const Standard_Real theTolUV,
|
|
const Standard_Real theUmin,
|
|
const Standard_Real theUmax,
|
|
const Standard_Real theVmin,
|
|
const Standard_Real theVmax)
|
|
{
|
|
const Standard_Integer aNbPnts = theWire.Length();
|
|
if (aNbPnts < 2)
|
|
return;
|
|
|
|
// Accumulate angle
|
|
TColgp_Array1OfPnt2d aPClass(1, aNbPnts);
|
|
Standard_Real anAngle = 0.0;
|
|
gp_Pnt2d p1 = theWire(1), p2 = theWire(2), p3;
|
|
aPClass(1) = p1;
|
|
aPClass(2) = p2;
|
|
|
|
const Standard_Real aAngTol = Precision::Angular();
|
|
const Standard_Real aSqConfusion =
|
|
Precision::PConfusion() * Precision::PConfusion();
|
|
|
|
for (Standard_Integer i = 1; i <= aNbPnts; i++)
|
|
{
|
|
Standard_Integer ii = i + 2;
|
|
if (ii > aNbPnts)
|
|
{
|
|
p3 = aPClass(ii - aNbPnts);
|
|
}
|
|
else
|
|
{
|
|
p3 = theWire.Value(ii);
|
|
aPClass(ii) = p3;
|
|
}
|
|
|
|
gp_Vec2d A(p1,p2), B(p2,p3);
|
|
if (A.SquareMagnitude() > aSqConfusion &&
|
|
B.SquareMagnitude() > aSqConfusion)
|
|
{
|
|
const Standard_Real aCurAngle = A.Angle(B);
|
|
const Standard_Real aCurAngleAbs = Abs(aCurAngle);
|
|
// Check if vectors are opposite
|
|
if (aCurAngleAbs > aAngTol && (M_PI - aCurAngleAbs) > aAngTol)
|
|
{
|
|
anAngle += aCurAngle;
|
|
p1 = p2;
|
|
}
|
|
}
|
|
p2 = p3;
|
|
}
|
|
// Check for zero angle - treat self intersecting wire as outer
|
|
if (Abs(anAngle) < aAngTol)
|
|
anAngle = 0.0;
|
|
|
|
myTabClass.Append( (void *)new CSLib_Class2d(aPClass,
|
|
theTolUV, theTolUV, theUmin, theVmin, theUmax, theVmax) );
|
|
myTabOrient.Append( !(anAngle < 0.0) );
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : Destroy
|
|
//purpose :
|
|
//=======================================================================
|
|
void BRepMesh_Classifier::Destroy()
|
|
{
|
|
Standard_Integer aNb = myTabClass.Length();
|
|
for (Standard_Integer i = 1; i <= aNb; i++)
|
|
{
|
|
if (myTabClass(i))
|
|
{
|
|
delete ((CSLib_Class2d*)myTabClass(i));
|
|
myTabClass(i) = NULL;
|
|
}
|
|
}
|
|
|
|
myTabClass.Clear();
|
|
myTabOrient.Clear();
|
|
}
|