/src/contrib/geom-5.1.2.7/src/GEOMAlgo/GEOMAlgo_GlueAnalyser.cpp

  1//  Copyright (C) 2007-2008  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.
 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// File:	GEOMAlgo_GlueDetector.cxx
 23// Created:	Wed Dec 15 11:08:09 2004
 24// Author:	Peter KURNEV
 25//		<pkv@irinox>
 26//
 27#include <GEOMAlgo_GlueAnalyser.ixx>
 28
 29#include <TopoDS.hxx>
 30#include <TopoDS_Shape.hxx>
 31#include <TopoDS_Face.hxx>
 32#include <TopoDS_Edge.hxx>
 33#include <TopoDS_Compound.hxx>
 34
 35#include <BRep_Builder.hxx>
 36
 37#include <TopExp.hxx>
 38
 39#include <TopTools_IndexedMapOfShape.hxx>
 40#include <TopTools_ListOfShape.hxx>
 41#include <TopTools_ListIteratorOfListOfShape.hxx>
 42#include <TopTools_IndexedDataMapOfShapeListOfShape.hxx>
 43#include <TopTools_DataMapIteratorOfDataMapOfShapeListOfShape.hxx>
 44
 45#include <GEOMAlgo_PassKeyShape.hxx>
 46#include <GEOMAlgo_IndexedDataMapOfPassKeyShapeListOfShape.hxx>
 47#include <GEOMAlgo_Tools.hxx>
 48#include <GEOMAlgo_CoupleOfShapes.hxx>
 49
 50#include <GEOMAlgo_Gluer.hxx>
 51#include <Bnd_HArray1OfBox.hxx>
 52#include <Bnd_BoundSortBox.hxx>
 53#include <TopoDS_Vertex.hxx>
 54#include <GEOMAlgo_IndexedDataMapOfIntegerShape.hxx>
 55#include <GEOMAlgo_IndexedDataMapOfShapeBox.hxx>
 56#include <Bnd_Box.hxx>
 57#include <TColStd_ListOfInteger.hxx>
 58#include <TopTools_MapOfShape.hxx>
 59#include <TColStd_ListIteratorOfListOfInteger.hxx>
 60#include <BRepBndLib.hxx>
 61
 62//=======================================================================
 63//function : 
 64//purpose  : 
 65//=======================================================================
 66  GEOMAlgo_GlueAnalyser::GEOMAlgo_GlueAnalyser()
 67:
 68  GEOMAlgo_Gluer()
 69{}
 70//=======================================================================
 71//function : ~
 72//purpose  : 
 73//=======================================================================
 74  GEOMAlgo_GlueAnalyser::~GEOMAlgo_GlueAnalyser()
 75{}
 76//=======================================================================
 77//function : HasSolidsToGlue
 78//purpose  : 
 79//=======================================================================
 80  Standard_Boolean GEOMAlgo_GlueAnalyser::HasSolidsToGlue()const
 81{
 82  return !mySolidsToGlue.IsEmpty();
 83}
 84//=======================================================================
 85//function : HasSolidsAlone
 86//purpose  : 
 87//=======================================================================
 88  Standard_Boolean GEOMAlgo_GlueAnalyser::HasSolidsAlone()const
 89{
 90  return !mySolidsAlone.IsEmpty();
 91}
 92//=======================================================================
 93//function : SolidsToGlue
 94//purpose  : 
 95//=======================================================================
 96  const GEOMAlgo_ListOfCoupleOfShapes& GEOMAlgo_GlueAnalyser::SolidsToGlue()const
 97{
 98  return mySolidsToGlue;
 99}
100//=======================================================================
101//function : SolidsAlone
102//purpose  : 
103//=======================================================================
104  const TopTools_ListOfShape& GEOMAlgo_GlueAnalyser::SolidsAlone()const
105{
106  return mySolidsAlone;
107}
108//=======================================================================
109//function : Perform
110//purpose  : 
111//=======================================================================
112  void GEOMAlgo_GlueAnalyser::Perform()
113{
114  myErrorStatus=0;
115  myWarningStatus=0;
116  //
117  mySolidsToGlue.Clear();
118  mySolidsAlone.Clear();
119  //
120  CheckData();
121  if (myErrorStatus) {
122    return;
123  }
124  //
125  InnerTolerance();
126  if (myErrorStatus) {
127    return;
128  }
129  //
130  DetectVertices();
131  if (myErrorStatus) {
132    return;
133  }
134  //
135  DetectEdges();
136  if (myErrorStatus) {
137    return;
138  }
139  //
140  DetectFaces();
141  if (myErrorStatus) {
142    return;
143  }
144  //
145  DetectSolids();
146  if (myErrorStatus) {
147    return;
148  }
149}
150//=======================================================================
151//function : DetectVertices
152//purpose  : 
153//=======================================================================
154  void GEOMAlgo_GlueAnalyser::DetectVertices()
155{
156  myErrorStatus=0;
157  //
158  Standard_Integer j, i, aNbV, aIndex, aNbVSD;
159  TColStd_ListIteratorOfListOfInteger aIt;
160  Handle(Bnd_HArray1OfBox) aHAB;
161  Bnd_BoundSortBox aBSB;
162  TopoDS_Shape aSTmp, aVF;
163  TopoDS_Vertex aVnew;
164  TopTools_IndexedMapOfShape aMV, aMVProcessed;
165  TopTools_ListIteratorOfListOfShape aItS;
166  TopTools_DataMapIteratorOfDataMapOfShapeListOfShape aItIm;
167  GEOMAlgo_IndexedDataMapOfIntegerShape aMIS;
168  GEOMAlgo_IndexedDataMapOfShapeBox aMSB;
169  //
170  TopExp::MapShapes(myShape, TopAbs_VERTEX, aMV);
171  aNbV=aMV.Extent();
172  if (!aNbV) {
173    myErrorStatus=2; // no vertices in source shape
174    return;
175  }
176  //
177  aHAB=new Bnd_HArray1OfBox(1, aNbV);
178  //
179  for (i=1; i<=aNbV; ++i) {
180    const TopoDS_Shape& aV=aMV(i);
181    Bnd_Box aBox;
182    //
183    aBox.SetGap(myTol); 
184    BRepBndLib::Add(aV, aBox);
185    aHAB->SetValue(i, aBox);
186    aMIS.Add(i, aV);
187    aMSB.Add(aV, aBox); 
188  }
189  //
190  aBSB.Initialize(aHAB);
191  //
192  for (i=1; i<=aNbV; ++i) {
193    const TopoDS_Shape& aV=aMV(i);
194    //
195    if (aMVProcessed.Contains(aV)) {
196      continue;
197    }
198    //
199    const Bnd_Box& aBoxV=aMSB.FindFromKey(aV);
200    const TColStd_ListOfInteger& aLI=aBSB.Compare(aBoxV);
201    aNbVSD=aLI.Extent();
202    if (!aNbVSD) {
203      myErrorStatus=3; // it must not be 
204      return;
205    }
206    //
207    // Images
208    TopTools_ListOfShape aLVSD;
209    //
210    aIt.Initialize(aLI);
211    for (j=0; aIt.More(); aIt.Next(), ++j) {
212      aIndex=aIt.Value();
213      const TopoDS_Shape& aVx=aMIS.FindFromKey(aIndex);
214      if(!j) {
215	aVF=aVx;
216      }
217      aLVSD.Append(aVx);
218      aMVProcessed.Add(aVx);
219    }
220    myImages.Bind(aVF, aLVSD);
221  }
222  // Origins
223  aItIm.Initialize(myImages);
224  for (; aItIm.More(); aItIm.Next()) {
225    const TopoDS_Shape& aV=aItIm.Key();
226    const TopTools_ListOfShape& aLVSD=aItIm.Value();
227    //
228    aItS.Initialize(aLVSD);
229    for (; aItS.More(); aItS.Next()) {
230      const TopoDS_Shape& aVSD=aItS.Value();
231      if (!myOrigins.IsBound(aVSD)) {
232	myOrigins.Bind(aVSD, aV);
233      }
234    }
235  }
236}
237//=======================================================================
238//function : DetectFaces
239//purpose  : 
240//=======================================================================
241  void GEOMAlgo_GlueAnalyser::DetectFaces()
242{
243  DetectShapes(TopAbs_FACE);
244}
245//=======================================================================
246//function : DetectEdges
247//purpose  : 
248//=======================================================================
249  void GEOMAlgo_GlueAnalyser::DetectEdges()
250{
251  DetectShapes(TopAbs_EDGE);
252}
253//=======================================================================
254//function : DetectShapes
255//purpose  : 
256//=======================================================================
257  void GEOMAlgo_GlueAnalyser::DetectShapes(const TopAbs_ShapeEnum aType)
258{
259  myErrorStatus=0;
260  //
261  Standard_Integer i, aNbF, aNbSDF, iErr;
262  TopoDS_Shape aNewShape;
263  TopTools_IndexedMapOfShape aMF;
264  TopTools_ListIteratorOfListOfShape aItS;
265  GEOMAlgo_PassKeyShape aPKF;
266  GEOMAlgo_IndexedDataMapOfPassKeyShapeListOfShape aMPKLF;
267  //
268  TopExp::MapShapes(myShape, aType, aMF);
269  //
270  aNbF=aMF.Extent();
271  for (i=1; i<=aNbF; ++i) {
272    const TopoDS_Shape& aS=aMF(i);
273    // 
274    //aPKF.Clear();//qft
275    if (aType==TopAbs_FACE) {
276      const TopoDS_Face& aF=TopoDS::Face(aS);
277      FacePassKey(aF, aPKF);
278    }
279    else if (aType==TopAbs_EDGE) {
280      const TopoDS_Edge& aE=TopoDS::Edge(aS);
281      EdgePassKey(aE, aPKF);
282    }
283    //
284    if (myErrorStatus) {
285      return;
286    }
287    //
288    if (aMPKLF.Contains(aPKF)) {
289      TopTools_ListOfShape& aLSDF=aMPKLF.ChangeFromKey(aPKF);
290      aLSDF.Append(aS);
291    }
292    else {
293      TopTools_ListOfShape aLSDF;
294      //
295      aLSDF.Append(aS);
296      aMPKLF.Add(aPKF, aLSDF);
297    }
298  }
299  // check geometric coincidence
300  if (myCheckGeometry) {
301    iErr=GEOMAlgo_Tools::RefineSDShapes(aMPKLF, myTol, myContext); //XX
302    if (iErr) {
303      myErrorStatus=200;
304      return;
305    }
306  }
307  //
308  // Images/Origins
309  aNbF=aMPKLF.Extent();
310  for (i=1; i<=aNbF; ++i) {
311    const TopTools_ListOfShape& aLSDF=aMPKLF(i);
312    aNbSDF=aLSDF.Extent();
313    if (!aNbSDF) {
314      myErrorStatus=4; // it must not be
315    }
316    //
317    const TopoDS_Shape& aS1=aLSDF.First();
318    aNewShape=aS1;
319    //
320    myImages.Bind(aNewShape, aLSDF);
321    // origins
322    aItS.Initialize(aLSDF);
323    for (; aItS.More(); aItS.Next()) {
324      const TopoDS_Shape& aFSD=aItS.Value();
325      if (!myOrigins.IsBound(aFSD)) {
326	myOrigins.Bind(aFSD, aNewShape);
327      }
328    }
329  }
330}
331//=======================================================================
332//function : DetectSolids
333//purpose  : 
334//=======================================================================
335  void GEOMAlgo_GlueAnalyser::DetectSolids()
336{
337  myErrorStatus=0;
338  //
339  Standard_Integer i, aNbF, aNbS, aNbC, aNbX;
340  TopoDS_Compound aCmp;
341  BRep_Builder aBB;
342  TopTools_IndexedDataMapOfShapeListOfShape aMFS;
343  TopTools_IndexedMapOfShape aMx, aMS;
344  TopTools_DataMapIteratorOfDataMapOfShapeListOfShape aItIm;
345  GEOMAlgo_CoupleOfShapes aCS;
346  //
347  GEOMAlgo_IndexedDataMapOfPassKeyShapeListOfShape aMPKLS;
348  GEOMAlgo_PassKeyShape aPKSx;
349  //
350  aBB.MakeCompound(aCmp);
351  //
352  TopExp::MapShapesAndAncestors(myShape, TopAbs_FACE, TopAbs_SOLID, aMFS);
353  //
354  aItIm.Initialize(myImages);
355  for (; aItIm.More(); aItIm.Next()) {
356    const TopoDS_Shape& aIm=aItIm.Key();
357    if (aIm.ShapeType()!=TopAbs_FACE) {
358      continue;
359    }
360    //
361    const TopTools_ListOfShape& aLF=aItIm.Value();
362    aNbF=aLF.Extent();
363    if (aNbF!=2) {
364      continue;
365    }
366    //
367    TopoDS_Shape aSx[2], aFx[2];
368    //
369    aFx[0]=aLF.First();
370    aFx[1]=aLF.Last();
371    for (i=0; i<2; ++i) {
372      if (!aMFS.Contains(aFx[i])) {
373	continue;// it must not be so
374      }
375      //
376      const TopTools_ListOfShape& aLS=aMFS.FindFromKey(aFx[i]);
377      aNbS=aLS.Extent();
378      if (aNbS!=1) {
379	continue;
380      }
381      aSx[i]=aLS.First();
382    }
383    //
384    if (aSx[0].IsNull() || aSx[1].IsNull()) {
385      continue;
386    }
387    //
388    //aPKSx.Clear();//qft
389    //qf
390    //aPKSx.SetIds(aSx[0], aSx[1]);
391    aPKSx.SetShapes(aSx[0], aSx[1]);
392    //qt
393    //
394    if (!aMPKLS.Contains(aPKSx)) {
395      TopTools_ListOfShape aLSx;
396      //
397      aLSx.Append(aSx[0]);
398      aLSx.Append(aSx[1]);
399      //
400      aMPKLS.Add(aPKSx, aLSx);
401    }
402  }
403  //
404  mySolidsToGlue.Clear();
405  mySolidsAlone.Clear();
406  
407  //
408  aNbC=aMPKLS.Extent();
409  if (!aNbC) {
410    return;
411  }
412  //
413  for (i=1; i<=aNbC; ++i) {
414    const TopTools_ListOfShape& aLSx=aMPKLS(i);
415    const TopoDS_Shape& aSx1=aLSx.First();
416    const TopoDS_Shape& aSx2=aLSx.Last();
417    aCS.SetShape1(aSx1);
418    aCS.SetShape2(aSx2);
419    mySolidsToGlue.Append(aCS);
420    //
421    if (!aMx.Contains(aSx1)) {
422      aBB.Add(aCmp, aSx1);
423      aMx.Add(aSx1);
424    }
425    if (!aMx.Contains(aSx2)) {
426      aBB.Add(aCmp, aSx2);
427      aMx.Add(aSx2);
428    }
429  }
430  myResult=aCmp;
431  //
432  // check alone solids
433  TopExp::MapShapes(myShape, TopAbs_SOLID, aMS);
434  //
435  aNbX=aMx.Extent();
436  for (i=1; i<=aNbX; ++i) {
437    const TopoDS_Shape& aSx=aMx(i);
438    if (!aMS.Contains(aSx)) {
439      mySolidsAlone.Append(aSx);
440    }
441  }
442}
443
444/*
445// A
446  // Make vertices
447  aMV.Clear();
448  aItIm.Initialize(myImages);
449  for (; aItIm.More(); aItIm.Next()) {
450    const TopoDS_Shape& aV=aItIm.Key();
451    aMV.Add(aV);
452    const TopTools_ListOfShape& aLVSD=aItIm.Value();
453    MakeVertex(aLVSD, aVnew);//ZZ
454    myImages.Bind(aVnew, aLVSD);
455  }
456  //
457  aNbV=aMV.Extent();
458  for (i=1; i<=aNbV; ++i) {
459    const TopoDS_Shape& aV=aMV(i);
460    myImages.UnBind(aV);
461  }
462  //
463  */