Salome HOME
982e21921f602a39a10d68c1fa57ae6af7f5b6de
[modules/geom.git] / src / GEOMAlgo / GEOMAlgo_VertexSolid.cxx
1 // Copyright (C) 2007-2014  CEA/DEN, EDF R&D, OPEN CASCADE
2 //
3 // Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
5 //
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License, or (at your option) any later version.
10 //
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14 // Lesser General Public License for more details.
15 //
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
19 //
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
21 //
22
23 // File:        GEOMAlgo_VertexSolid.cxx
24 // Created:     Wed Jan 12 16:36:40 2005
25 // Author:      Peter KURNEV
26 //              <pkv@irinox>
27 //
28 #include <GEOMAlgo_VertexSolid.hxx>
29
30 #include <Basics_OCCTVersion.hxx>
31
32 #include <gp_Pnt.hxx>
33
34 #include <TopAbs_ShapeEnum.hxx>
35 #include <TopAbs_State.hxx>
36
37 #include <TopTools_ListIteratorOfListOfShape.hxx>
38 #include <TopTools_ListOfShape.hxx>
39 #include <TopTools_IndexedMapOfShape.hxx>
40
41 #include <TopoDS.hxx>
42 #include <TopoDS_Solid.hxx>
43 #include <TopoDS_Vertex.hxx>
44
45 #include <TopExp.hxx>
46
47 #include <BRep_Tool.hxx>
48 #include <BRepClass3d_SolidClassifier.hxx>
49 //
50 #include <BOPCol_ListOfShape.hxx>
51 #if OCC_VERSION_LARGE > 0x06070100
52 #include <IntTools_Context.hxx>
53 #else
54 #include <BOPInt_Context.hxx>
55 #endif
56 //
57 #include <BOPDS_DS.hxx>
58 #include <BOPDS_IndexRange.hxx>
59 #include <BOPDS_VectorOfInterfVV.hxx>
60 #include <BOPDS_VectorOfInterfVE.hxx>
61 #include <BOPDS_VectorOfInterfVF.hxx>
62 #include <BOPDS_Interf.hxx>
63
64 //=======================================================================
65 //function : GEOMAlgo_VertexSolid
66 //purpose  :
67 //=======================================================================
68 GEOMAlgo_VertexSolid::GEOMAlgo_VertexSolid()
69 :
70   GEOMAlgo_ShapeSolid()
71 {
72 }
73 //=======================================================================
74 //function : ~
75 //purpose  :
76 //=======================================================================
77 GEOMAlgo_VertexSolid::~GEOMAlgo_VertexSolid()
78 {
79 }
80 //=======================================================================
81 // function: Perform
82 // purpose:
83 //=======================================================================
84 void GEOMAlgo_VertexSolid::Perform()
85 {
86   myErrorStatus=0;
87   //
88   try {
89     if (myDSFiller==NULL) {
90       myErrorStatus=10;
91       return;
92     }
93     if(myDSFiller->ErrorStatus()) {
94       myErrorStatus=11;
95       return;
96     }
97     //
98     Standard_Integer aNbF, aNbArgs;
99     TopTools_IndexedMapOfShape aM;
100     //
101     const BOPDS_DS& aDS=myDSFiller->DS();
102     const BOPCol_ListOfShape& aLS=aDS.Arguments();
103     aNbArgs=aLS.Extent();
104     if (aNbArgs!=2) {
105       myErrorStatus=14;
106       return;
107     }
108     
109     const TopoDS_Shape& aObj=aLS.First();
110     //
111     TopExp::MapShapes(aObj, TopAbs_FACE, aM);
112     aNbF=aM.Extent();
113     myRank=(aNbF) ? 1 : 0;
114     //
115     BuildResult();
116   }
117   //
118   catch (Standard_Failure) {
119     myErrorStatus = 12;
120   }
121 }
122 //=======================================================================
123 // function: BuildResult
124 // purpose:
125 //=======================================================================
126 void GEOMAlgo_VertexSolid::BuildResult()
127 {
128   Standard_Integer i, iBeg, iEnd, aNbVV, aNbVE, aNbVF, j, iFound, aNbRanges;
129   Standard_Real aTol;
130   TopAbs_State aSt;
131   TopAbs_ShapeEnum aType;
132   gp_Pnt aP3D;
133   //
134   myLSIN.Clear();
135   myLSOUT.Clear();
136   myLSON.Clear();
137   //
138   const BOPDS_DS& aDS=myDSFiller->DS();
139   BOPDS_DS* pDS=(BOPDS_DS*)&aDS;
140   //
141   BOPDS_VectorOfInterfVV& aVVs=pDS->InterfVV();
142   BOPDS_VectorOfInterfVE& aVEs=pDS->InterfVE();
143   BOPDS_VectorOfInterfVF& aVFs=pDS->InterfVF();
144   //
145   const BOPCol_ListOfShape& aLS=aDS.Arguments();
146   const TopoDS_Shape& aObj=aLS.First();
147   //
148   const TopoDS_Shape& aTool=aLS.Last();
149   const TopoDS_Solid& aSolid=(myRank==0) ? TopoDS::Solid(aTool) : TopoDS::Solid(aObj);
150   //
151 #if OCC_VERSION_LARGE > 0x06070100
152   Handle(IntTools_Context) aCtx=myDSFiller->Context();
153 #else
154   Handle(BOPInt_Context) aCtx=myDSFiller->Context();
155 #endif
156   BRepClass3d_SolidClassifier& aSC=aCtx->SolidClassifier(aSolid);
157   //
158   aNbRanges=aDS.NbRanges();
159   const BOPDS_IndexRange& aRange=aDS.Range(myRank);
160   aRange.Indices(iBeg, iEnd);
161   //
162   for (i=iBeg; i<=iEnd; ++i) {
163     const TopoDS_Shape& aS=aDS.Shape(i);
164     aType=aS.ShapeType();
165     if (aType!=TopAbs_VERTEX) {
166       continue; 
167     }
168     //
169     const TopoDS_Vertex& aV=TopoDS::Vertex(aS);
170     //
171     iFound=0;
172     //
173     // 1
174     aNbVV=aVVs.Extent();
175     for (j=0; j<aNbVV; ++j) {
176       BOPDS_InterfVV& aVV=aVVs(j);
177       if (aVV.Contains(i)) {
178         myLSON.Append(aV);
179         iFound=1;
180         break;
181       }
182     }
183     if (iFound) {
184       continue; 
185     }
186     // 2
187     aNbVE=aVEs.Extent();
188     for (j=0; j<aNbVE; ++j) {
189       BOPDS_InterfVE& aVE=aVEs(j);
190       if (aVE.Contains(i)) {
191         myLSON.Append(aV);
192         iFound=1;
193         break;
194       }
195     }
196     if (iFound) {
197       continue; 
198     }
199     // 3
200     aNbVF=aVFs.Extent();
201     for (j=0; j<aNbVF; ++j) {
202       BOPDS_InterfVF& aVF=aVFs(j);
203       if (aVF.Contains(i)) {
204         myLSON.Append(aV);
205         iFound=1;
206         break;
207       }
208     }
209     if (iFound) {
210       continue; 
211     }
212     //
213     // 4
214     aP3D=BRep_Tool::Pnt(aV);
215     aTol=1.e-7;
216     aSC.Perform(aP3D, aTol);
217     aSt=aSC.State();
218     if (aSt==TopAbs_IN) {
219       myLSIN.Append(aV);
220     }
221     else if (aSt==TopAbs_OUT) {
222       myLSOUT.Append(aV);
223     }
224   }//for (i=iBeg; i<iEnd; ++i) {
225 }