/src/contrib/geom-5.1.2.7/src/NMTDS/NMTDS_ShapesDataStructure.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:	NMTDS_ShapesDataStructure.cxx
 23// Created:	Mon Dec  1 10:21:04 2003
 24// Author:	Peter KURNEV
 25//		<pkv@irinox>
 26//
 27#include <NMTDS_ShapesDataStructure.ixx>
 28#include <TopoDS_Iterator.hxx>
 29#include <TopoDS_Shape.hxx>
 30#include <BooleanOperations_IndexedDataMapOfShapeAncestorsSuccessors.hxx>
 31#include <NMTDS_ListOfIndexedDataMapOfShapeAncestorsSuccessors.hxx>
 32#include <NMTDS_ListIteratorOfListOfIndexedDataMapOfShapeAncestorsSuccessors.hxx>
 33#include <BooleanOperations_ShapeAndInterferences.hxx>
 34#include <NMTDS_IndexRange.hxx>
 35//
 36#include <BooleanOperations_IndexedDataMapOfShapeAncestorsSuccessors.hxx>
 37#include <BooleanOperations_AncestorsSeqAndSuccessorsSeq.hxx>
 38#include <TColStd_MapOfInteger.hxx>
 39#include <NMTDS_Tools.hxx>
 40
 41
 42static
 43  void ComputeBoxExS(const Standard_Integer aIx,
 44		     const NMTDS_ShapesDataStructure* pDS,
 45		     Bnd_Box& aBoxEx);
 46static
 47  void GetAllSuccessorsS(const Standard_Integer nS,
 48			 const NMTDS_ShapesDataStructure* myDS,
 49			 TColStd_IndexedMapOfInteger& aMA);
 50
 51//===========================================================================
 52//function : NMTDS_ShapesDataStructure::NMTDS_ShapesDataStructure
 53//purpose  : 
 54//===========================================================================
 55  NMTDS_ShapesDataStructure::NMTDS_ShapesDataStructure()
 56:
 57  BooleanOperations_ShapesDataStructure()
 58{}
 59//===========================================================================
 60//function : SetCompositeShape
 61//purpose  : 
 62//===========================================================================
 63  void NMTDS_ShapesDataStructure::SetCompositeShape(const TopoDS_Shape& aS)
 64{
 65  myCompositeShape=aS;
 66}
 67//===========================================================================
 68//function : CompositeShape
 69//purpose  : 
 70//===========================================================================
 71  const TopoDS_Shape& NMTDS_ShapesDataStructure::CompositeShape()const
 72{
 73  return myCompositeShape;
 74}
 75//===========================================================================
 76//function : Ranges
 77//purpose  : 
 78//===========================================================================
 79  const NMTDS_CArray1OfIndexRange& NMTDS_ShapesDataStructure::Ranges()const
 80{
 81  return myRanges;
 82}
 83// Modified to Add new methods Thu Sep 14 14:35:18 2006 
 84// Contribution of Samtech www.samcef.com BEGIN
 85//===========================================================================
 86//function : FillMap
 87//purpose  : 
 88//===========================================================================
 89  void NMTDS_ShapesDataStructure::FillMap
 90  (const TopoDS_Shape& aS,
 91   BooleanOperations_IndexedDataMapOfShapeAncestorsSuccessors& aMSA,
 92   BooleanOperations_IndexedDataMapOfShapeAncestorsSuccessors& aMS) const
 93{
 94  Standard_Integer iX, i, j, aIndex, aNbSc, aNbS;
 95  BooleanOperations_AncestorsSeqAndSuccessorsSeq aAS;
 96  //
 97  aMSA.Add(aS, aAS);
 98  aMS.Add(aS, aAS);
 99  FillSubshapes(aS, aMSA, aMS);
100  //
101  aNbS=aMS.Extent();
102  for(i=1; i<=aNbS; ++i) {
103    TColStd_MapOfInteger aMFence;
104    //
105    const TopoDS_Shape& aSX=aMS.FindKey(i);
106    iX=aMSA.FindIndex(aSX);
107    const BooleanOperations_AncestorsSeqAndSuccessorsSeq& aAS1=aMSA(iX);
108    //
109    aNbSc=aAS1.NumberOfSuccessors();
110    for(j=1; j<=aNbSc; ++j) {
111      aIndex=aAS1.GetSuccessor(j);
112      if(aMFence.Add(aIndex)) {
113	BooleanOperations_AncestorsSeqAndSuccessorsSeq& aAS2=aMSA.ChangeFromIndex(aIndex);
114	aAS2.SetNewAncestor(iX);
115      }
116    }
117  }
118}
119//===========================================================================
120//function : FillSubshapes
121//purpose  : 
122//===========================================================================
123  void NMTDS_ShapesDataStructure::FillSubshapes
124  (const TopoDS_Shape& aS,
125   BooleanOperations_IndexedDataMapOfShapeAncestorsSuccessors& aMSA,
126   BooleanOperations_IndexedDataMapOfShapeAncestorsSuccessors& aMS) const
127{
128  Standard_Boolean bIsNewSubShape;
129  Standard_Integer aIndexSubShape, aIndex;
130  BooleanOperations_AncestorsSeqAndSuccessorsSeq aASx;
131  //
132  aIndex=aMSA.FindIndex(aS);
133  BooleanOperations_AncestorsSeqAndSuccessorsSeq& aAS=aMSA.ChangeFromIndex(aIndex);
134  //
135  TopoDS_Iterator anIt(aS, Standard_True);
136  for(; anIt.More(); anIt.Next()) {
137    const TopoDS_Shape& aSubShape = anIt.Value();
138    bIsNewSubShape = Standard_False;
139    if(!aMSA.Contains(aSubShape)) {
140      bIsNewSubShape=!bIsNewSubShape;
141      aIndexSubShape=aMSA.Add(aSubShape, aASx);
142      aMS.Add(aSubShape, aASx);
143    }
144    else {
145      aIndexSubShape=aMSA.FindIndex(aSubShape);
146    }
147    aAS.SetNewSuccessor(aIndexSubShape);
148    aAS.SetNewOrientation(aSubShape.Orientation());
149    //
150    if(bIsNewSubShape && (aSubShape.ShapeType() != TopAbs_VERTEX)) {
151      FillSubshapes(aSubShape, aMSA, aMS);
152    }
153  }
154}
155// Contribution of Samtech www.samcef.com END 
156//===========================================================================
157//function : Init
158//purpose  : 
159//===========================================================================
160  void NMTDS_ShapesDataStructure::Init()
161{
162  Standard_Integer i, j, aNbSx, aNbS, aShift, aNbRanges;
163  Standard_Integer iFirst, iLast;
164  NMTDS_ListOfIndexedDataMapOfShapeAncestorsSuccessors aLx;
165  NMTDS_ListIteratorOfListOfIndexedDataMapOfShapeAncestorsSuccessors aLit;
166  TopoDS_Iterator anIt;
167  BooleanOperations_IndexedDataMapOfShapeAncestorsSuccessors aMSA;
168  //
169  anIt.Initialize(myCompositeShape);
170  for (; anIt.More(); anIt.Next()) {
171    const TopoDS_Shape& aSx=anIt.Value(); 
172    BooleanOperations_IndexedDataMapOfShapeAncestorsSuccessors aMS;
173    //modified by NIZNHY-PKV Tue Feb 27 17:05:47 2007f
174    //FillMap(aSx, aMSA, aMS);
175    //aLx.Append(aMS);
176    //
177    if (!aMSA.Contains(aSx)) {
178      FillMap(aSx, aMSA, aMS);
179      aLx.Append(aMS);
180    }
181    //modified by NIZNHY-PKV Tue Feb 27 17:06:03 2007t
182  }
183  aNbS=aMSA.Extent(); 
184  //
185  // Fill myRanges
186  //modified by NIZNHY-PKV Tue Feb 27 17:10:07 2007f
187  i=aLx.Extent();
188  //modified by NIZNHY-PKV Tue Feb 27 17:10:10 2007t
189  myRanges.Resize(i);
190  aLit.Initialize(aLx);
191  for (i=1; aLit.More(); aLit.Next(), ++i) {
192    // Modified  Thu Sep 14 14:35:18 2006 
193    // Contribution of Samtech www.samcef.com BEGIN
194    const BooleanOperations_IndexedDataMapOfShapeAncestorsSuccessors& aMSx=aLit.Value();
195    aNbSx=aMSx.Extent();
196    // Contribution of Samtech www.samcef.com END
197    if (i==1) {
198      iFirst=1;
199      iLast=aNbSx;
200      myRanges(i).SetFirst(iFirst);
201      myRanges(i).SetLast(iLast);
202      continue;
203    }
204    iFirst=myRanges(i-1).Last()+1;
205    iLast=iFirst+aNbSx-1;
206    myRanges(i).SetFirst(iFirst);
207    myRanges(i).SetLast(iLast);
208  }
209  //
210  myNumberOfShapesOfTheObject=aNbS;
211  // Modified  Thu Sep 14 14:35:18 2006 
212  // Contribution of Samtech www.samcef.com BEGIN
213  //myNumberOfShapesOfTheTool=aNbS;
214  //myLength=3*aNbS;
215  myNumberOfShapesOfTheTool=0;
216  myLength=2*aNbS;
217  // Contribution of Samtech www.samcef.com END
218  //
219  // Allocate the whole Table
220  myListOfShapeAndInterferences = (BooleanOperations_PShapeAndInterferences)
221    Standard::Allocate(myLength*sizeof(BooleanOperations_ShapeAndInterferences));
222  // Modified  Thu Sep 14 14:35:18 2006 
223  // Contribution of Samtech www.samcef.com BEGIN
224  //
225  // Fill the table
226    
227  //aShift=0;
228  //for (i=0; i<2; ++i) {
229  //  if (i) {
230  //    aShift=aNbS;
231  //  }
232  //  for (j=1; j<=aNbS; ++j) {
233  //    const TopoDS_Shape& aSx=aMSA.FindKey(j);
234  //    const BooleanOperations_AncestorsSeqAndSuccessorsSeq& aASx=
235  //      aMSA.FindFromIndex(j);
236  //    InsertShapeAndAncestorsSuccessors(aSx, aASx, aShift);
237  //  }
238  //}
239  
240  aShift=0;
241  for (j=1; j<=aNbS; ++j) {
242    const TopoDS_Shape& aSx=aMSA.FindKey(j);
243    const BooleanOperations_AncestorsSeqAndSuccessorsSeq& aASx=aMSA.FindFromIndex(j);
244    InsertShapeAndAncestorsSuccessors(aSx, aASx, aShift);
245  }
246  // Contribution of Samtech www.samcef.com END
247  //
248  // myShapeIndexMap
249  myShapeIndexMap.Clear();
250  //
251  aNbRanges=myRanges.Extent();
252  for (i=1; i<=aNbRanges; ++i){
253    BooleanOperations_IndexedDataMapOfShapeInteger aSIM;
254    //
255    const NMTDS_IndexRange& aR=myRanges(i);
256    iFirst=aR.First();
257    iLast =aR.Last();
258    for (j=iFirst; j<=iLast; ++j) { 
259      const TopoDS_Shape& aS=Shape(j);
260      aSIM.Add(aS, j);
261    }
262    myShapeIndexMap.Add(i, aSIM);
263  }
264  //
265  // myRefEdges
266  iLast=myNumberOfShapesOfTheObject+myNumberOfShapesOfTheTool;
267  myRefEdges.Resize(iLast);
268
269  for (i=1; i<=iLast; ++i) {
270    const TopoDS_Shape& aS=Shape(i);
271    myRefEdges(i)=0;
272    if (aS.ShapeType()==TopAbs_EDGE) {
273      myNbEdges++;
274      myRefEdges(i)=myNbEdges;
275    }
276  }
277}
278//===========================================================================
279//function : ShapeRangeIndex
280//purpose  : 
281//===========================================================================
282  Standard_Integer NMTDS_ShapesDataStructure::ShapeRangeIndex(const Standard_Integer aId)const
283{
284  Standard_Boolean bFound;
285  Standard_Integer i, aNbR, aIdx, aNbS, aZero=0;
286  //
287  aNbS=myNumberOfShapesOfTheObject+myNumberOfShapesOfTheTool;
288  
289  aIdx=aId;
290  if (aIdx > aNbS || aIdx < 1){
291    return aZero;
292  }
293  //  
294  if (aIdx > myNumberOfShapesOfTheObject) {
295    aIdx-=myNumberOfShapesOfTheObject;
296  }
297  //
298  aNbR=myRanges.Extent();
299  for (i=1; i<=aNbR; ++i) {
300    const NMTDS_IndexRange& aRange=myRanges(i);
301    bFound=aRange.IsInRange(aIdx);
302    if (bFound) {
303     return i;
304    }
305  }
306  return aZero;
307}
308//===========================================================================
309//function : Rank
310//purpose  : 
311//===========================================================================
312  Standard_Integer NMTDS_ShapesDataStructure::Rank(const Standard_Integer aId)const
313{
314  Standard_Boolean bFound;
315  Standard_Integer i, aNbR, aNbS, aZero=0;
316  //
317  aNbS=myNumberOfShapesOfTheObject;
318  
319  if (aId > aNbS || aId < 1){
320    return aZero;
321  }
322  //  
323  aNbR=myRanges.Extent();
324  for (i=1; i<=aNbR; ++i) {
325    const NMTDS_IndexRange& aRange=myRanges(i);
326    bFound=aRange.IsInRange(aId);
327    if (bFound) {
328     return i;
329    }
330  }
331  return aZero;
332}
333//===========================================================================
334//function : ShapeIndex
335//purpose  : 
336//===========================================================================
337  Standard_Integer NMTDS_ShapesDataStructure::ShapeIndex(const TopoDS_Shape& aS,
338							 const Standard_Integer aRank)const
339{
340  Standard_Boolean bFound;
341  Standard_Integer aIndex=0;
342  //
343  bFound=myShapeIndexMap.Contains(aRank);
344  if (!bFound) {
345    return aIndex;
346  }
347  // Modified Thu Sep 14 14:35:18 2006 
348  // Contribution of Samtech www.samcef.com BEGIN
349  //
350  //const BooleanOperations_IndexedDataMapOfShapeInteger& aSIM=
351  //  myShapeIndexMap.FindFromKey(aRank);
352  //
353  //bFound=aSIM.Contains(aS);
354  //if (!bFound) {
355  //  return aIndex;
356  //}
357  //
358  //aIndex=aSIM.FindFromKey(aS);
359  //
360  Standard_Integer i, aNbRanks;
361  //
362  aNbRanks=myShapeIndexMap.Extent();
363  for (i=1; i<=aNbRanks; ++i){
364    const BooleanOperations_IndexedDataMapOfShapeInteger& aSIM=myShapeIndexMap.FindFromKey(i);
365    bFound=aSIM.Contains(aS);
366    if (bFound) {
367      aIndex=aSIM.FindFromKey(aS);
368      return aIndex;
369    }
370  }
371  // Contribution of Samtech www.samcef.com END
372  return aIndex;
373}
374
375//=======================================================================
376//function : ComputeBoxEx
377//purpose  : 
378//=======================================================================
379  void NMTDS_ShapesDataStructure::ComputeBoxEx (const Standard_Integer aIx,
380						Bnd_Box& aBoxEx)const
381{
382  ComputeBoxExS(aIx, this, aBoxEx);
383}
384//=======================================================================
385//function : GetAllSuccessors
386//purpose  : 
387//=======================================================================
388  void NMTDS_ShapesDataStructure::GetAllSuccessors(const Standard_Integer nS,
389						   TColStd_IndexedMapOfInteger& aMA)const
390{
391  GetAllSuccessorsS(nS, this, aMA);
392}
393//=======================================================================
394//function : GetAllSuccessorsS
395//purpose  : 
396//=======================================================================
397void GetAllSuccessorsS(const Standard_Integer nS,
398		       const NMTDS_ShapesDataStructure* myDS,
399		       TColStd_IndexedMapOfInteger& aMA)
400{
401  TopAbs_ShapeEnum aT;
402  Standard_Integer i, nSx, aNbSuccessors, *pSuccessors;
403  Standard_Address xSuccessors;
404  //
405  const TopoDS_Shape& aS=myDS->Shape(nS);
406  aT=aS.ShapeType();
407  if(NMTDS_Tools::HasBRep(aT)) {
408    aMA.Add(nS);
409    //
410    if (aT==TopAbs_VERTEX) {
411      return;
412    }
413  }
414  //
415  myDS->GetSuccessors(nS, xSuccessors, aNbSuccessors);
416  pSuccessors=(Standard_Integer*)xSuccessors;
417  for (i=0; i<aNbSuccessors; ++i) {
418    nSx=pSuccessors[i];
419    GetAllSuccessorsS(nSx, myDS, aMA);
420  }
421}
422//=======================================================================
423// function: ComputeBoxExS
424// purpose: 
425//=======================================================================
426  void ComputeBoxExS(const Standard_Integer aIx,
427		     const NMTDS_ShapesDataStructure* pDS,
428		     Bnd_Box& aBoxEx)
429{
430  Standard_Integer i, aNbS, iS;
431  //
432  const Bnd_Box& aBox=pDS->GetBoundingBox(aIx);
433  aBoxEx.Add(aBox);
434  //
435  aNbS=pDS->NumberOfSuccessors(aIx);
436  for (i=1; i<=aNbS; ++i) {
437    Bnd_Box aBoxS;
438    iS=pDS->GetSuccessor(aIx, i);
439    ComputeBoxExS(iS, pDS, aBoxS);
440    aBoxEx.Add(aBoxS);
441  }
442}