1 // Copyright (C) 2007-2011 CEA/DEN, EDF R&D, OPEN CASCADE
3 // Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
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.
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.
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
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
23 #include <Standard_Stream.hxx>
25 #include <GEOMImpl_GlueDriver.hxx>
26 #include <GEOMImpl_IGlue.hxx>
27 #include <GEOMImpl_Types.hxx>
29 #include <GEOM_Object.hxx>
30 #include <GEOM_Function.hxx>
32 //#include <GEOMAlgo_Gluer.hxx>
33 //#include "GEOMAlgo_Gluer1.hxx"
34 #include "GEOMAlgo_Gluer2.hxx"
35 #include "GEOMAlgo_ListIteratorOfListOfCoupleOfShapes.hxx"
36 #include "GEOMAlgo_CoupleOfShapes.hxx"
37 #include "GEOMAlgo_ListOfCoupleOfShapes.hxx"
39 #include "utilities.h"
41 #include <TDataStd_IntegerArray.hxx>
44 #include <TopoDS_Shape.hxx>
45 #include <TopTools_DataMapIteratorOfDataMapOfShapeListOfShape.hxx>
46 #include <TopTools_IndexedMapOfShape.hxx>
47 #include <TopTools_ListOfShape.hxx>
48 #include <TopTools_ListIteratorOfListOfShape.hxx>
50 #include <ShapeFix_Shape.hxx>
52 #include <Standard_NullObject.hxx>
53 #include <Standard_Failure.hxx>
55 #define MSG_BAD_TOLERANCE "Tolerance is too big"
56 #define MSG_BAD_ARG_SHAPE "Argument shape is not a compound of hexahedral solids"
58 //=======================================================================
59 //function : GEOMImpl_GlueDriver
61 //=======================================================================
62 GEOMImpl_GlueDriver::GEOMImpl_GlueDriver()
66 //=======================================================================
69 //=======================================================================
70 const Standard_GUID& GEOMImpl_GlueDriver::GetID()
72 static Standard_GUID aGlueDriver("FF1BBB63-5D14-4df2-980B-3A668264EA16");
76 //=======================================================================
77 //function : GlueFacesWithWarnings
79 //=======================================================================
81 TopoDS_Shape GEOMImpl_GlueDriver::GlueFacesWithWarnings (const TopoDS_Shape& theShape,
82 const Standard_Real theTolerance,
83 const Standard_Boolean doKeepNonSolids,
84 TCollection_AsciiString& theWarning) const
86 Standard_Integer iErr, iWrn;
88 GEOMAlgo_Gluer aGluer;
90 aGluer.SetShape(theShape);
91 aGluer.SetTolerance(theTolerance);
92 aGluer.SetCheckGeometry(Standard_True);
93 aGluer.SetKeepNonSolids(doKeepNonSolids);
97 iErr = aGluer.ErrorStatus();
101 Standard_Failure::Raise("No vertices found in source shape");
105 Standard_Failure::Raise(MSG_BAD_TOLERANCE " or " MSG_BAD_ARG_SHAPE);
108 Standard_Failure::Raise("Source shape is Null");
111 Standard_Failure::Raise("Result shape is Null");
114 Standard_Failure::Raise(MSG_BAD_TOLERANCE);
118 Standard_Failure::Raise(MSG_BAD_ARG_SHAPE);
121 Standard_Failure::Raise("Error occured during check of geometric coincidence");
125 // description of all errors see in GEOMAlgo_Gluer.cxx
126 TCollection_AsciiString aMsg ("Error in GEOMAlgo_Gluer with code ");
127 aMsg += TCollection_AsciiString(iErr);
128 Standard_Failure::Raise(aMsg.ToCString());
135 iWrn = aGluer.WarningStatus();
140 Standard_Integer nbAlone = aGluer.AloneShapes();
141 theWarning = TCollection_AsciiString(nbAlone);
142 theWarning += " solid(s) can not be glued by faces";
146 // description of all warnings see in GEOMAlgo_Gluer.cxx
147 theWarning = "Warning in GEOMAlgo_Gluer with code ";
148 theWarning += TCollection_AsciiString(iWrn);
153 aRes = aGluer.Result();
155 // SKL 18.01.2010 - patch for 20662
156 Handle(ShapeFix_Shape) aSfs = new ShapeFix_Shape(aRes);
157 aSfs->SetPrecision(Precision::Confusion());
159 aRes = aSfs->Shape();
161 // Fill history to be used by GetInPlace functionality
162 TopTools_IndexedMapOfShape aResIndices;
163 TopExp::MapShapes(aRes, aResIndices);
165 Handle(GEOM_Function) aFunction = GEOM_Function::GetFunction(Label());
167 // history for all argument shapes
168 TDF_LabelSequence aLabelSeq;
169 aFunction->GetDependency(aLabelSeq);
170 Standard_Integer nbArg = aLabelSeq.Length();
172 for (Standard_Integer iarg = 1; iarg <= nbArg; iarg++) {
174 TDF_Label anArgumentRefLabel = aLabelSeq.Value(iarg);
176 Handle(GEOM_Object) anArgumentObject = GEOM_Object::GetReferencedObject(anArgumentRefLabel);
177 TopoDS_Shape anArgumentShape = anArgumentObject->GetValue();
179 TopTools_IndexedMapOfShape anArgumentIndices;
180 TopExp::MapShapes(anArgumentShape, anArgumentIndices);
181 Standard_Integer nbArgumentEntities = anArgumentIndices.Extent();
183 // Find corresponding label in history
184 TDF_Label anArgumentHistoryLabel =
185 aFunction->GetArgumentHistoryEntry(anArgumentRefLabel, Standard_True);
187 for (Standard_Integer ie = 1; ie <= nbArgumentEntities; ie++) {
188 TopoDS_Shape anEntity = anArgumentIndices.FindKey(ie);
189 const TopTools_ListOfShape& aModified = aGluer.Modified(anEntity);
190 Standard_Integer nbModified = aModified.Extent();
192 if (nbModified > 0) {
193 TDF_Label aWhatHistoryLabel = anArgumentHistoryLabel.FindChild(ie, Standard_True);
194 Handle(TDataStd_IntegerArray) anAttr =
195 TDataStd_IntegerArray::Set(aWhatHistoryLabel, 1, nbModified);
197 TopTools_ListIteratorOfListOfShape itM (aModified);
198 for (int im = 1; itM.More(); itM.Next(), ++im) {
199 int id = aResIndices.FindIndex(itM.Value());
200 anAttr->SetValue(im, id);
209 //=======================================================================
210 //function : GlueFaces
212 //=======================================================================
213 TopoDS_Shape GEOMImpl_GlueDriver::GlueFaces (const TopoDS_Shape& theShape,
214 const Standard_Real theTolerance,
215 const Standard_Boolean doKeepNonSolids)
218 GEOMAlgo_Gluer aGluer;
220 aGluer.SetShape(theShape);
221 aGluer.SetTolerance(theTolerance);
222 aGluer.SetCheckGeometry(Standard_True);
223 aGluer.SetKeepNonSolids(doKeepNonSolids);
227 Standard_Integer iErr = aGluer.ErrorStatus();
231 Standard_Failure::Raise("No vertices found in source shape");
234 Standard_Failure::Raise("Source shape is Null");
237 Standard_Failure::Raise("Result shape is Null");
240 Standard_Failure::Raise("Error occured during check of geometric coincidence");
244 // description of all errors see in GEOMAlgo_Gluer.cxx
245 TCollection_AsciiString aMsg ("Error in GEOMAlgo_Gluer with code ");
246 aMsg += TCollection_AsciiString(iErr);
247 Standard_Failure::Raise(aMsg.ToCString());
254 Standard_Integer iWrn = aGluer.WarningStatus();
258 MESSAGE("Some shapes can not be glued by faces");
261 // description of all warnings see in GEOMAlgo_Gluer.cxx
262 MESSAGE("Warning in GEOMAlgo_Gluer with code " << iWrn);
267 aRes = aGluer.Result();
272 //=======================================================================
273 //function : GlueFacesByList
275 //=======================================================================
276 TopoDS_Shape GEOMImpl_GlueDriver::GlueFacesByList (const TopoDS_Shape& theShape,
277 const Standard_Real theTolerance,
278 const Standard_Boolean doKeepNonSolids,
279 const TopTools_MapOfShape& aFaces)
283 GEOMAlgo_Gluer1 aGluer;
284 GEOMAlgo_ListIteratorOfListOfCoupleOfShapes aItCS;
285 GEOMAlgo_CoupleOfShapes aCS;
286 GEOMAlgo_ListOfCoupleOfShapes aLCS;
288 aGluer.SetShape(theShape);
289 aGluer.SetTolerance(theTolerance);
290 aGluer.SetKeepNonSolids(doKeepNonSolids);
292 Standard_Integer iErr = aGluer.ErrorStatus();
293 if (iErr) return aRes;
295 TopTools_ListOfShape listShape;
296 const GEOMAlgo_ListOfCoupleOfShapes& aLCSG = aGluer.GluedFaces();
298 aItCS.Initialize(aLCSG);
299 for (; aItCS.More(); aItCS.Next()) {
300 const GEOMAlgo_CoupleOfShapes& aCSG = aItCS.Value();
301 const TopoDS_Shape& aF1 = aCSG.Shape1();
302 const TopoDS_Shape& aF2 = aCSG.Shape2();
303 if (aFaces.Contains(aF1) || aFaces.Contains(aF2))
305 aCS.SetShapes(aF1,aF2);
309 //cout<<"aLCS.Extent() = "<<aLCS.Extent()<<endl;
310 if (aLCS.Extent()>0) {
311 aGluer.SetFacesToUnglue(aLCS);
312 aGluer.UnglueFaces();
313 iErr = aGluer.ErrorStatus();
314 if (iErr) return aRes;
317 aRes = aGluer.Result();
323 //=======================================================================
324 //function : GlueFaces
326 //=======================================================================
327 TopoDS_Shape GEOMImpl_GlueDriver::GlueFaces (const TopoDS_Shape& theShape,
328 const Standard_Real theTolerance,
329 const Standard_Boolean doKeepNonSolids)
336 aGA.SetArgument(theShape);
337 aGA.SetTolerance(theTolerance);
338 aGA.SetKeepNonSolids(doKeepNonSolids);
340 // 2. Detect interferred shapes
342 Standard_Integer iErr = aGA.ErrorStatus();
346 Standard_Failure::Raise("GEOMAlgo_GlueDetector failed");
350 Standard_Failure::Raise("PerformImagesToWork failed");
354 // description of all errors see in GEOMAlgo_Gluer2.cxx
355 TCollection_AsciiString aMsg ("Error in GEOMAlgo_Gluer2 with code ");
356 aMsg += TCollection_AsciiString(iErr);
357 Standard_Failure::Raise(aMsg.ToCString());
364 // 3. Set shapes to glue. If the operator is absent, the whole gluing will be done
365 //aGA.SetShapesToGlue(aMSG);
369 iErr = aGA.ErrorStatus();
373 Standard_Failure::Raise("GEOMAlgo_GlueDetector failed");
377 Standard_Failure::Raise("PerformImagesToWork failed");
381 // description of all errors see in GEOMAlgo_Gluer2.cxx
382 TCollection_AsciiString aMsg ("Error in GEOMAlgo_Gluer2 with code ");
383 aMsg += TCollection_AsciiString(iErr);
384 Standard_Failure::Raise(aMsg.ToCString());
391 Standard_Integer iWrn = aGA.WarningStatus();
395 MESSAGE("No shapes to glue");
398 // description of all warnings see in GEOMAlgo_Gluer2.cxx
399 MESSAGE("Warning in GEOMAlgo_Gluer2 with code " << iWrn);
410 //=======================================================================
411 //function : GlueWithWarnings
413 //=======================================================================
414 TopoDS_Shape GEOMImpl_GlueDriver::GlueWithWarnings (const TopoDS_Shape& theShape,
415 const Standard_Real theTolerance,
416 const TopAbs_ShapeEnum theShapeType,
417 const Standard_Boolean doKeepNonSolids,
418 TCollection_AsciiString& theWarning) const
425 aGA.SetArgument(theShape);
426 aGA.SetTolerance(theTolerance);
427 aGA.SetKeepNonSolids(doKeepNonSolids);
429 // 2. Detect interferred shapes
431 Standard_Integer iErr = aGA.ErrorStatus();
435 Standard_Failure::Raise("GEOMAlgo_GlueDetector failed");
439 Standard_Failure::Raise("PerformImagesToWork failed");
443 // description of all errors see in GEOMAlgo_Gluer2.cxx
444 TCollection_AsciiString aMsg ("Error in GEOMAlgo_Gluer2 with code ");
445 aMsg += TCollection_AsciiString(iErr);
446 Standard_Failure::Raise(aMsg.ToCString());
453 if (theShapeType != TopAbs_FACE) {
454 // 3. Fill shapes to glue aMSG
455 TopTools_DataMapOfShapeListOfShape aMSG;
456 const TopTools_DataMapOfShapeListOfShape& aMSD = aGA.ShapesDetected();
457 TopTools_DataMapIteratorOfDataMapOfShapeListOfShape aItMSD;
458 aItMSD.Initialize(aMSD);
459 for (; aItMSD.More(); aItMSD.Next()) {
460 const TopoDS_Shape& aSx = aItMSD.Key();
461 const TopTools_ListOfShape& aLSD = aItMSD.Value();
462 if (aSx.ShapeType() == theShapeType) {
463 aMSG.Bind(aSx, aLSD);
467 // 4. Set shapes to glue. If the operator is absent, the whole gluing will be done
468 aGA.SetShapesToGlue(aMSG);
473 iErr = aGA.ErrorStatus();
477 Standard_Failure::Raise("GEOMAlgo_GlueDetector failed");
481 Standard_Failure::Raise("PerformImagesToWork failed");
485 // description of all errors see in GEOMAlgo_Gluer2.cxx
486 TCollection_AsciiString aMsg ("Error in GEOMAlgo_Gluer2 with code ");
487 aMsg += TCollection_AsciiString(iErr);
488 Standard_Failure::Raise(aMsg.ToCString());
495 Standard_Integer iWrn = aGA.WarningStatus();
499 theWarning = "No shapes to glue";
502 // description of all warnings see in GEOMAlgo_Gluer2.cxx
503 theWarning = "Warning in GEOMAlgo_Gluer2 with code ";
504 theWarning += TCollection_AsciiString(iWrn);
512 // 7. Fill history to be used by GetInPlace functionality
513 TopTools_IndexedMapOfShape aResIndices;
514 TopExp::MapShapes(aRes, aResIndices);
516 Handle(GEOM_Function) aFunction = GEOM_Function::GetFunction(Label());
518 // history for all argument shapes
519 TDF_LabelSequence aLabelSeq;
520 aFunction->GetDependency(aLabelSeq);
521 Standard_Integer nbArg = aLabelSeq.Length();
523 for (Standard_Integer iarg = 1; iarg <= nbArg; iarg++) {
525 TDF_Label anArgumentRefLabel = aLabelSeq.Value(iarg);
527 Handle(GEOM_Object) anArgumentObject = GEOM_Object::GetReferencedObject(anArgumentRefLabel);
528 TopoDS_Shape anArgumentShape = anArgumentObject->GetValue();
530 TopTools_IndexedMapOfShape anArgumentIndices;
531 TopExp::MapShapes(anArgumentShape, anArgumentIndices);
532 Standard_Integer nbArgumentEntities = anArgumentIndices.Extent();
534 // Find corresponding label in history
535 TDF_Label anArgumentHistoryLabel =
536 aFunction->GetArgumentHistoryEntry(anArgumentRefLabel, Standard_True);
538 for (Standard_Integer ie = 1; ie <= nbArgumentEntities; ie++) {
539 TopoDS_Shape anEntity = anArgumentIndices.FindKey(ie);
540 const TopTools_ListOfShape& aModified = aGA.Modified(anEntity);
541 Standard_Integer nbModified = aModified.Extent();
543 if (nbModified > 0) {
544 TDF_Label aWhatHistoryLabel = anArgumentHistoryLabel.FindChild(ie, Standard_True);
545 Handle(TDataStd_IntegerArray) anAttr =
546 TDataStd_IntegerArray::Set(aWhatHistoryLabel, 1, nbModified);
548 TopTools_ListIteratorOfListOfShape itM (aModified);
549 for (int im = 1; itM.More(); itM.Next(), ++im) {
550 int id = aResIndices.FindIndex(itM.Value());
551 anAttr->SetValue(im, id);
560 //=======================================================================
561 //function : GlueByList
563 //=======================================================================
564 TopoDS_Shape GEOMImpl_GlueDriver::GlueByList (const TopoDS_Shape& theShape,
565 const Standard_Real theTolerance,
566 const Standard_Boolean doKeepNonSolids,
567 const TopTools_MapOfShape& theShapesList,
568 const Standard_Boolean doGlueAllEdges)
575 aGA.SetArgument(theShape);
576 aGA.SetTolerance(theTolerance);
577 aGA.SetKeepNonSolids(doKeepNonSolids);
579 // 2. Detect interferred shapes
581 Standard_Integer iErr = aGA.ErrorStatus();
585 Standard_Failure::Raise("GEOMAlgo_GlueDetector failed");
589 Standard_Failure::Raise("PerformImagesToWork failed");
593 // description of all errors see in GEOMAlgo_Gluer2.cxx
594 TCollection_AsciiString aMsg ("Error in GEOMAlgo_Gluer2 with code ");
595 aMsg += TCollection_AsciiString(iErr);
596 Standard_Failure::Raise(aMsg.ToCString());
603 // 3. Fill shapes to glue aMSG
604 TopTools_DataMapOfShapeListOfShape aMSG;
605 const TopTools_DataMapOfShapeListOfShape& aMSD = aGA.ShapesDetected();
606 TopTools_DataMapIteratorOfDataMapOfShapeListOfShape aItMSD;
607 aItMSD.Initialize(aMSD);
608 for (; aItMSD.More(); aItMSD.Next()) {
609 const TopoDS_Shape& aSx = aItMSD.Key();
610 const TopTools_ListOfShape& aLSD = aItMSD.Value();
611 TopTools_ListIteratorOfListOfShape anItLSD (aLSD);
612 bool isToGlue = false;
613 if (doGlueAllEdges && aSx.ShapeType() == TopAbs_EDGE) {
617 for (; anItLSD.More() && !isToGlue; anItLSD.Next()) {
618 if (theShapesList.Contains(anItLSD.Value())) {
624 aMSG.Bind(aSx, aLSD);
628 // 4. Set shapes to glue. If the operator is absent, the whole gluing will be done
629 aGA.SetShapesToGlue(aMSG);
633 iErr = aGA.ErrorStatus();
637 Standard_Failure::Raise("GEOMAlgo_GlueDetector failed");
641 Standard_Failure::Raise("PerformImagesToWork failed");
645 // description of all errors see in GEOMAlgo_Gluer2.cxx
646 TCollection_AsciiString aMsg ("Error in GEOMAlgo_Gluer2 with code ");
647 aMsg += TCollection_AsciiString(iErr);
648 Standard_Failure::Raise(aMsg.ToCString());
655 Standard_Integer iWrn = aGA.WarningStatus();
659 MESSAGE("No shapes to glue");
662 // description of all warnings see in GEOMAlgo_Gluer2.cxx
663 MESSAGE("Warning in GEOMAlgo_Gluer2 with code " << iWrn);
674 //=======================================================================
677 //=======================================================================
678 Standard_Integer GEOMImpl_GlueDriver::Execute(TFunction_Logbook& log) const
680 if (Label().IsNull()) return 0;
681 Handle(GEOM_Function) aFunction = GEOM_Function::GetFunction(Label());
683 GEOMImpl_IGlue aCI (aFunction);
684 Standard_Integer aType = aFunction->GetType();
687 TCollection_AsciiString aWrn;
689 Handle(GEOM_Function) aRefBase = aCI.GetBase();
690 TopoDS_Shape aShapeBase = aRefBase->GetValue();
691 if (aShapeBase.IsNull()) {
692 Standard_NullObject::Raise("Shape for gluing is null");
695 Standard_Real tol3d = aCI.GetTolerance();
697 Standard_Boolean aKeepNonSolids = aCI.GetKeepNonSolids();
699 if (aType == GLUE_FACES) {
700 //aShape = GlueFacesWithWarnings(aShapeBase, tol3d, aKeepNonSolids, aWrn);
701 aShape = GlueWithWarnings(aShapeBase, tol3d, TopAbs_FACE, aKeepNonSolids, aWrn);
703 else if (aType == GLUE_EDGES) {
704 aShape = GlueWithWarnings(aShapeBase, tol3d, TopAbs_EDGE, aKeepNonSolids, aWrn);
706 else if (aType == GLUE_FACES_BY_LIST || aType == GLUE_EDGES_BY_LIST) {
707 Handle(TColStd_HSequenceOfTransient) SF = aCI.GetFaces();
708 TopTools_MapOfShape aFaces;
710 for (; i <= SF->Length(); i++) {
711 Handle(Standard_Transient) anItem = SF->Value(i);
714 Handle(GEOM_Function) aRefSh = Handle(GEOM_Function)::DownCast(anItem);
717 TopoDS_Shape aFace = aRefSh->GetValue();
720 if (!aFaces.Contains(aFace))
724 Standard_Boolean aGlueAllEdges = Standard_False;
725 if (aType == GLUE_FACES_BY_LIST)
726 aGlueAllEdges = aCI.GetGlueAllEdges();
728 //aShape = GlueFacesByList(aShapeBase, tol3d, aKeepNonSolids, aFaces);
729 aShape = GlueByList(aShapeBase, tol3d, aKeepNonSolids, aFaces, aGlueAllEdges);
732 if (aShape.IsNull()) return 0;
734 aFunction->SetValue(aShape);
736 log.SetTouched(Label());
738 if (!aWrn.IsEmpty()) {
739 Standard_Failure::Raise(aWrn.ToCString());
745 //=======================================================================
746 //function : GEOMImpl_GlueDriver_Type_
748 //=======================================================================
749 Standard_EXPORT Handle_Standard_Type& GEOMImpl_GlueDriver_Type_()
752 static Handle_Standard_Type aType1 = STANDARD_TYPE(TFunction_Driver);
753 if (aType1.IsNull()) aType1 = STANDARD_TYPE(TFunction_Driver);
754 static Handle_Standard_Type aType2 = STANDARD_TYPE(MMgt_TShared);
755 if (aType2.IsNull()) aType2 = STANDARD_TYPE(MMgt_TShared);
756 static Handle_Standard_Type aType3 = STANDARD_TYPE(Standard_Transient);
757 if (aType3.IsNull()) aType3 = STANDARD_TYPE(Standard_Transient);
760 static Handle_Standard_Transient _Ancestors[]= {aType1,aType2,aType3,NULL};
761 static Handle_Standard_Type _aType = new Standard_Type("GEOMImpl_GlueDriver",
762 sizeof(GEOMImpl_GlueDriver),
764 (Standard_Address)_Ancestors,
765 (Standard_Address)NULL);
770 //=======================================================================
771 //function : DownCast
773 //=======================================================================
774 const Handle(GEOMImpl_GlueDriver) Handle(GEOMImpl_GlueDriver)::DownCast(const Handle(Standard_Transient)& AnObject)
776 Handle(GEOMImpl_GlueDriver) _anOtherObject;
778 if (!AnObject.IsNull()) {
779 if (AnObject->IsKind(STANDARD_TYPE(GEOMImpl_GlueDriver))) {
780 _anOtherObject = Handle(GEOMImpl_GlueDriver)((Handle(GEOMImpl_GlueDriver)&)AnObject);
784 return _anOtherObject;