src/BRepMesh/BRepMesh_Classifier.cxx

   1 // Created on: 1997-06-26
   2 // Created by: Laurent PAINNOT
   3 // Copyright (c) 1997-1999 Matra Datavision
   4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
   5 //
   6 // This file is part of Open CASCADE Technology software library.
   7 //
   8 // This library is free software; you can redistribute it and/or modify it under
   9 // the terms of the GNU Lesser General Public License version 2.1 as published
  10 // by the Free Software Foundation, with special exception defined in the file
  11 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
  12 // distribution for complete text of the license and disclaimer of any warranty.
  13 //
  14 // Alternatively, this file may be used under the terms of Open CASCADE
  15 // commercial license or contractual agreement.
  16
  17 #include <BRepMesh_Classifier.hxx>
  18
  19 #include <Precision.hxx>
  20 #include <gp_Pnt2d.hxx>
  21 #include <CSLib_Class2d.hxx>
  22 #include <TColgp_Array1OfPnt2d.hxx>
  23
  24 //=======================================================================
  25 //function : Constructor
  26 //purpose  :
  27 //=======================================================================
  28 BRepMesh_Classifier::BRepMesh_Classifier()
  29 {
  30 }
  31
  32 //=======================================================================
  33 //function : Perform
  34 //purpose  :
  35 //=======================================================================
  36 TopAbs_State BRepMesh_Classifier::Perform(const gp_Pnt2d& thePoint) const
  37 {
  38   Standard_Boolean isOut = Standard_False;
  39   Standard_Integer aNb   = myTabClass.Length();
  40
  41   for (Standard_Integer i = 1; i <= aNb; i++)
  42   {
  43     Standard_Integer aCur = ((CSLib_Class2d*)myTabClass(i))->SiDans(thePoint);
  44     if (aCur == 0)
  45     {
  46       // Point is ON, but mark it as OUT
  47       isOut = Standard_True;
  48     }
  49     else
  50       isOut = myTabOrient(i) ? (aCur == -1) : (aCur == 1);
  51
  52     if (isOut)
  53       return TopAbs_OUT;
  54   }
  55
  56   return TopAbs_IN;
  57 }
  58
  59 //=======================================================================
  60 //function : RegisterWire
  61 //purpose  :
  62 //=======================================================================
  63 void BRepMesh_Classifier::RegisterWire(
  64   const NCollection_Sequence<gp_Pnt2d>& theWire,
  65   const Standard_Real                   theTolUV,
  66   const Standard_Real                   theUmin,
  67   const Standard_Real                   theUmax,
  68   const Standard_Real                   theVmin,
  69   const Standard_Real                   theVmax)
  70 {
  71   const Standard_Integer aNbPnts = theWire.Length();
  72   if (aNbPnts < 2)
  73     return;
  74
  75   // Accumulate angle
  76   TColgp_Array1OfPnt2d aPClass(1, aNbPnts);
  77   Standard_Real anAngle = 0.0;
  78   gp_Pnt2d p1 = theWire(1), p2 = theWire(2), p3;
  79   aPClass(1) = p1;
  80   aPClass(2) = p2;
  81
  82   const Standard_Real aAngTol = Precision::Angular();
  83   const Standard_Real aSqConfusion =
  84     Precision::PConfusion() * Precision::PConfusion();
  85
  86   for (Standard_Integer i = 1; i <= aNbPnts; i++)
  87   {
  88     Standard_Integer ii = i + 2;
  89     if (ii > aNbPnts)
  90     {
  91       p3 = aPClass(ii - aNbPnts);
  92     }
  93     else
  94     {
  95       p3 = theWire.Value(ii);
  96       aPClass(ii) = p3;
  97     }
  98
  99     gp_Vec2d A(p1,p2), B(p2,p3);
 100     if (A.SquareMagnitude() > aSqConfusion &&
 101         B.SquareMagnitude() > aSqConfusion)
 102     {
 103       const Standard_Real aCurAngle    = A.Angle(B);
 104       const Standard_Real aCurAngleAbs = Abs(aCurAngle);
 105       // Check if vectors are opposite
 106       if (aCurAngleAbs > aAngTol && (M_PI - aCurAngleAbs) > aAngTol)
 107       {
 108         anAngle += aCurAngle;
 109         p1 = p2;
 110       }
 111     }
 112     p2 = p3;
 113   }
 114   // Check for zero angle - treat self intersecting wire as outer
 115   if (Abs(anAngle) < aAngTol)
 116     anAngle = 0.0;
 117
 118   myTabClass.Append( (void *)new CSLib_Class2d(aPClass,
 119     theTolUV, theTolUV, theUmin, theVmin, theUmax, theVmax) );
 120   myTabOrient.Append( !(anAngle < 0.0) );
 121 }
 122
 123 //=======================================================================
 124 //function : Destroy
 125 //purpose  :
 126 //=======================================================================
 127 void BRepMesh_Classifier::Destroy()
 128 {
 129   Standard_Integer aNb = myTabClass.Length();
 130   for (Standard_Integer i = 1; i <= aNb; i++)
 131   {
 132     if (myTabClass(i))
 133     {
 134       delete ((CSLib_Class2d*)myTabClass(i));
 135       myTabClass(i) = NULL;
 136     }
 137   }
 138
 139   myTabClass.Clear();
 140   myTabOrient.Clear();
 141 }