1 // GEOM GEOM : implementaion of GEOM_Gen.idl and GEOM_Shape.idl
3 // Copyright (C) 2003 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
24 // File : GEOM_GEN_i.cc file
25 // Author : Lucien PIGNOLONI
30 #include "GEOM_Gen_i.hh"
32 #include "Partition_Spliter.hxx"
33 #include "Archimede_VolumeSection.hxx"
34 #include "Sketcher_Profile.hxx"
36 #include "Utils_CorbaException.hxx"
37 #include "utilities.h"
43 #include <Standard_Failure.hxx>
45 #include <gp_Circ.hxx>
47 #include <gp_Elips.hxx>
48 #include <Geom_Plane.hxx>
49 #include <Geom_Line.hxx>
50 #include <Geom_BezierCurve.hxx>
51 #include <Geom_BSplineCurve.hxx>
52 #include <GeomFill_Line.hxx>
53 #include <GeomFill_AppSurf.hxx>
54 #include <GeomFill_SectionGenerator.hxx>
55 #include <Geom_BSplineSurface.hxx>
56 #include <Geom_TrimmedCurve.hxx>
57 #include <GC_MakeArcOfCircle.hxx>
58 #include <GeomAPI_PointsToBSpline.hxx>
59 #include <GC_Root.hxx>
61 #include <BRepCheck_Analyzer.hxx>
62 #if OCC_VERSION_MAJOR >= 5
63 #include <BRepAlgo.hxx>
65 #include <BRepAlgoAPI.hxx>
67 #include <BRepAlgo_FaceRestrictor.hxx>
68 #include <BRepAdaptor_Surface.hxx>
69 #include <BRepBuilderAPI_Copy.hxx>
70 #include <BRepAlgoAPI_Common.hxx>
71 #include <BRepAlgoAPI_Cut.hxx>
72 #include <BRepAlgoAPI_Fuse.hxx>
73 #include <BRepAlgoAPI_Section.hxx>
74 #include <BRepOffsetAPI_Sewing.hxx>
75 #include <BRepOffsetAPI_MakePipe.hxx>
77 #include <BRepBuilderAPI_MakeVertex.hxx>
78 #include <BRepBuilderAPI_MakeEdge.hxx>
79 #include <BRepBuilderAPI_MakeFace.hxx>
81 #include <BRepLib.hxx>
82 #include <BRepBndLib.hxx>
83 #include <Bnd_Box.hxx>
85 #include <BRepBuilderAPI_MakeShell.hxx>
86 #include <BRepPrim_Builder.hxx>
87 #include <BRepBuilderAPI_MakeSolid.hxx>
88 #include <BRepClass3d_SolidClassifier.hxx>
90 #include <BRepBuilderAPI_MakeWire.hxx>
91 #include <BRepBuilderAPI_Transform.hxx>
92 #include <BRepPrimAPI_MakeRevol.hxx>
93 #include <BRepPrimAPI_MakePrism.hxx>
94 #include <BRepPrimAPI_MakeTorus.hxx>
95 #include <BRepPrimAPI_MakeBox.hxx>
96 #include <BRepPrimAPI_MakeSphere.hxx>
97 #include <BRepPrimAPI_MakeCylinder.hxx>
98 #include <BRepPrimAPI_MakeCone.hxx>
99 #include <BRepFilletAPI_MakeFillet.hxx>
100 #include <BRepFilletAPI_MakeChamfer.hxx>
101 #include <BRepTools.hxx>
102 #include <BRepTools_Quilt.hxx>
103 #include <BRep_Tool.hxx>
105 #include <GeomAPI_ProjectPointOnCurve.hxx>
107 #include <BRepGProp.hxx>
108 #include <GProp_GProps.hxx>
109 #include <Precision.hxx>
111 //VRV: OCC 4.0 migration
112 #include <STEPControl_Reader.hxx>
113 #include <IGESControl_Reader.hxx>
114 //VRV: OCC 4.0 migration
116 #include <IFSelect_ReturnStatus.hxx>
117 #include <TColStd_HSequenceOfTransient.hxx>
119 //VRV: OCC 4.0 migration
120 #include <IGESControl_Writer.hxx>
121 #include <IGESControl_Controller.hxx>
122 #include <STEPControl_Writer.hxx>
123 #include <Interface_Static.hxx>
124 //#include <STEPControlStd_StepModelType.hxx>
125 //VRV: OCC 4.0 migration
127 #include <TopoDS_Shape.hxx>
128 #include <TopAbs.hxx>
129 #include <TopoDS_Wire.hxx>
130 #include <TopoDS_Edge.hxx>
131 #include <TopoDS_Compound.hxx>
132 #include <TopoDS_Solid.hxx>
134 #include <TopExp.hxx>
135 #include <TopExp_Explorer.hxx>
136 #include <TCollection_ExtendedString.hxx>
137 #include <TColgp_Array1OfPnt.hxx>
138 #include <TopoDS_Iterator.hxx>
139 #include <TopTools_MapOfShape.hxx>
140 #include <TopTools_MapIteratorOfMapOfShape.hxx>
141 #include <TopTools_ListIteratorOfListOfShape.hxx>
142 #include <TopTools_Array1OfShape.hxx>
144 #include <IGESData_IGESEntity.hxx>
146 #include <TDF_Tool.hxx>
147 #include <TDF_Label.hxx>
148 #include <TDataStd_Name.hxx>
149 #include <TDataStd_Comment.hxx>
150 #include <TDF_Reference.hxx>
151 #include <TDF_Data.hxx>
152 #include <TNaming_Builder.hxx>
153 #include <TNaming_NamedShape.hxx>
154 #include <TNaming_Tool.hxx>
155 // #include <TDocStd_Owner.hxx>
157 #include "SALOMEDS_Tool.hxx"
158 #include "GEOMDS_Commands.hxx"
159 #include "GEOMDS_Explorer.hxx"
161 #include CORBA_SERVER_HEADER(SALOMEDS_Attributes)
163 #include "Utils_ExceptHandlers.hxx"
165 Standard_EXPORT static Standard_Boolean IsValid(const TopoDS_Shape& S) {
166 #if OCC_VERSION_MAJOR >= 5
167 return BRepAlgo::IsValid(S);
169 return BRepAlgoAPI::IsValid(S);
173 //============================================================================
174 // function : GEOM_Gen_i()
175 // purpose : constructor to be called for servant creation.
176 //============================================================================
177 GEOM_Gen_i::GEOM_Gen_i(CORBA::ORB_ptr orb,
178 PortableServer::POA_ptr poa,
179 PortableServer::ObjectId * contId,
180 const char *instanceName,
181 const char *interfaceName) :
182 Engines_Component_i(orb, poa, contId, instanceName, interfaceName)
185 _id = _poa->activate_object(_thisObj);
187 name_service = new SALOME_NamingService(_orb);
188 myOCAFApp = new GEOMDS_Application();
190 GetCurrentStudy(0);//for correct work of SuperVisor
195 //============================================================================
196 // function : ~GEOM_Gen_i()
197 // purpose : destructor
198 //============================================================================
199 GEOM_Gen_i::~GEOM_Gen_i() {
204 //============================================================================
205 // function : IORToLocalPersistentID()
207 //============================================================================
208 char* GEOM_Gen_i::IORToLocalPersistentID(SALOMEDS::SObject_ptr theSObject,
209 const char* IORString,
210 CORBA::Boolean isMultiFile,
211 CORBA::Boolean isASCII)
213 GEOM::GEOM_Shape_var aShape = GEOM::GEOM_Shape::_narrow(_orb->string_to_object(IORString));
214 if (!CORBA::is_nil(aShape)) {
215 return aShape->ShapeId();
221 //============================================================================
222 // function : LocalPersistentIDToIOR()
223 // purpose : Create/Load CORBA object from a persistent ref (an entry)
224 // : Used when a study is loaded
225 // : The IOR (IORName) of object created is returned
226 //============================================================================
227 char* GEOM_Gen_i::LocalPersistentIDToIOR(SALOMEDS::SObject_ptr theSObject,
228 const char* aLocalPersistentID,
229 CORBA::Boolean isMultiFile,
230 CORBA::Boolean isASCII)
232 SALOMEDS::Study_var myStudy = theSObject->GetStudy();
233 GetCurrentStudy(myStudy->StudyId());
234 Handle(TDocStd_Document) aDoc = Handle(TDocStd_Document)::DownCast(myStudyIDToDoc(myStudy->StudyId()));
235 CORBA::String_var aPersRefString = aLocalPersistentID;
237 /* For a GEOM::GEOM_Shape the pers_ref is the Entry in the OCAF document */
238 TCollection_ExtendedString MainIOR;
240 TDF_Tool::Label(aDoc->GetData(), aPersRefString, Lab );
243 MESSAGE("Can not find label "<<aPersRefString<<" for study "<<myStudy->StudyId());
244 THROW_SALOME_CORBA_EXCEPTION("Incorrect GEOM data loaded",SALOME::BAD_PARAM);
245 //return ""; empty IORs raise new problem: "display" in popup crashes
248 Handle(TNaming_NamedShape) NS;
249 Lab.FindAttribute( TNaming_NamedShape::GetID(), NS );
250 TopoDS_Shape S = TNaming_Tool::GetShape(NS);
252 /* shapetype, index=0, topo, orb, shapetype, ismain=true and name are setted and modified later ? */
253 GEOM::GEOM_Shape_var result = CreateObject(S);
254 GEOMDS_Commands GC( aDoc->Main() ) ;
256 if ( GC.HasIOR(Lab) ) { /* shape already created/loaded */
260 /******************* Dependent object (not a main shape) *********************/
261 if( GC.IsDependentShape(Lab) ) {
263 TDF_Label mainLabel ;
264 Standard_Boolean mainShapeOk = GC.GetMainShapeLabel(Lab, mainLabel) ;
266 /* Main object not yet loaded we load/create it */
267 if( !GC.HasIOR(mainLabel) ) {
269 TCollection_AsciiString entry;
270 TDF_Tool::Entry(mainLabel,entry);
271 CORBA::String_var ent = CORBA::string_dup(entry.ToCString());
273 /* Create the main object recursively */
274 MainIOR = LocalPersistentIDToIOR(theSObject, ent, isMultiFile, isASCII) ;
276 GC.ReturnNameIOR( mainLabel, MainIOR );
279 result->MainName( TCollection_AsciiString(MainIOR).ToCString() ) ;
280 result->IsMainShape(false) ;
281 result->ShapeId(aPersRefString);
283 Handle(TDF_Reference) aRef;
284 Lab.FindAttribute( TDF_Reference::GetID(), aRef );
285 TDF_Label myL = aRef->Get() ;
286 Handle(TNaming_NamedShape) NN;
287 myL.FindAttribute( TNaming_NamedShape::GetID(), NN );
288 TopoDS_Shape mainTopo = TNaming_Tool::GetShape(NN);
290 GEOM::GEOM_Shape::ListOfSubShapeID_var ListOfID = new GEOM::GEOM_Shape::ListOfSubShapeID;
292 if(S.ShapeType() != TopAbs_COMPOUND) {
293 /* to set the index of a unique sub shape (Explode All ONLY for the moment !) */
296 TopTools_MapOfShape M;
297 TopExp_Explorer Exp ;
298 for( Exp.Init(mainTopo, TopAbs_ShapeEnum( result->ShapeType() )) ; Exp.More(); Exp.Next() ) {
299 if ( M.Add(Exp.Current()) ) {
300 if(Exp.Current().IsSame(S) ) {
307 result->Index(ListOfID) ;
308 return result->Name();
311 /* Here is a TopAbs_COMPOUND : we set the list/index for a compound : containing two or more sub shapes */
312 /* Warning : the Corba shape has a shapetype Compound : in GEOMDS_Client we have to retrieve the kind of */
313 /* subshapes contained in this compound ! */
314 TopTools_SequenceOfShape SS;
316 TopExp_Explorer exp ;
317 TopAbs_ShapeEnum subType ;
319 /* Set all sub shapes in a sequence of shapes */
320 for ( it.Initialize( S, true, true ) ; it.More(); it.Next() ) {
321 subType = it.Value().ShapeType() ;
322 SS.Append( it.Value() ) ;
325 ListOfID->length( SS.Length() ) ;
326 int j, k ; /* in TopTools_SequenceOfShape index start at 1 */
328 for( k=1; k<=SS.Length(); k++ ) {
330 for( exp.Init( mainTopo, subType ); exp.More(); exp.Next() ) {
331 if( exp.Current().IsSame( SS.Value(k) ) ) {
337 result->Index(ListOfID) ;
338 return result->Name();
342 /******************* Independent object (not a sub shape) *********************/
344 result->IsMainShape(true) ;
345 if( !GC.AddIORNameAttribute(Lab, result->Name() ) ) {
346 MESSAGE("in LocalPersistentIDToIOR, NAME/IOR attribute already exist." << endl ) ;
348 Handle(TNaming_NamedShape) NamedShape ;
349 bool notTested = Lab.FindAttribute(TNaming_NamedShape::GetID(), NamedShape) ;
350 result->ShapeId(aPersRefString);
351 return result->Name();
355 //============================================================================
356 // function : CanPublishInStudy
358 //============================================================================
359 bool GEOM_Gen_i::CanPublishInStudy(CORBA::Object_ptr theIOR)
361 GEOM::GEOM_Shape_var aShape = GEOM::GEOM_Shape::_narrow(theIOR);
362 return !(aShape->_is_nil());
366 //============================================================================
367 // function : PublishInStudy
369 //============================================================================
370 SALOMEDS::SObject_ptr GEOM_Gen_i::PublishInStudy(SALOMEDS::Study_ptr theStudy,
371 SALOMEDS::SObject_ptr theSObject,
372 CORBA::Object_ptr theObject,
373 const char* theName) throw (SALOME::SALOME_Exception)
375 Unexpect aCatch(SALOME_SalomeException);
376 SALOMEDS::SObject_var aResultSO;
377 if(CORBA::is_nil(theObject)) return aResultSO;
379 GEOM::GEOM_Shape_var aShape = GEOM::GEOM_Shape::_narrow(theObject);
380 if(aShape->_is_nil()) return aResultSO;
382 if(theStudy->_is_nil()) return aResultSO;
384 SALOMEDS::GenericAttribute_var anAttr;
385 SALOMEDS::StudyBuilder_var aStudyBuilder = theStudy->NewBuilder();
387 SALOMEDS::SComponent_var aFather = theStudy->FindComponent("GEOM");
388 if (aFather->_is_nil()) {
389 aFather = aStudyBuilder->NewComponent("GEOM");
390 anAttr = aStudyBuilder->FindOrCreateAttribute(aFather, "AttributeName");
391 SALOMEDS::AttributeName_var aName = SALOMEDS::AttributeName::_narrow(anAttr);
392 aName->SetValue("Geometry");
393 anAttr = aStudyBuilder->FindOrCreateAttribute(aFather, "AttributePixMap");
394 SALOMEDS::AttributePixMap::_narrow(anAttr)->SetPixMap("ICON_OBJBROWSER_Geometry");
395 aStudyBuilder->DefineComponentInstance(aFather, GEOM_Gen::_this());
397 if (aFather->_is_nil()) return aResultSO;
399 if (CORBA::is_nil(theSObject)) {
400 aResultSO = aStudyBuilder->NewObject(aFather);
402 if (!theSObject->ReferencedObject(aResultSO))
403 THROW_SALOME_CORBA_EXCEPTION("Publish in study supervision graph error",SALOME::BAD_PARAM);
405 anAttr = aStudyBuilder->FindOrCreateAttribute(aResultSO, "AttributeIOR");
406 SALOMEDS::AttributeIOR_var anIOR = SALOMEDS::AttributeIOR::_narrow(anAttr);
407 anIOR->SetValue(aShape->Name());
409 anAttr = aStudyBuilder->FindOrCreateAttribute(aResultSO, "AttributePixMap");
410 SALOMEDS::AttributePixMap_var aPixmap = SALOMEDS::AttributePixMap::_narrow(anAttr);
411 TCollection_AsciiString aShapeName("Shape_");
413 if ( aShape->ShapeType() == GEOM::COMPOUND ) {
414 aPixmap->SetPixMap( "ICON_OBJBROWSER_COMPOUND" );
415 aShapeName = "Compound_";
416 } else if ( aShape->ShapeType() == GEOM::COMPSOLID ) {
417 aPixmap->SetPixMap( "ICON_OBJBROWSER_COMPSOLID" );
418 aShapeName = "Compsolid_";
419 } else if ( aShape->ShapeType() == GEOM::SOLID ) {
420 aPixmap->SetPixMap( "ICON_OBJBROWSER_SOLID" );
421 aShapeName = "Solid_";
422 } else if ( aShape->ShapeType() == GEOM::SHELL ) {
423 aPixmap->SetPixMap( "ICON_OBJBROWSER_SHELL" );
424 aShapeName = "Shell_";
425 } else if ( aShape->ShapeType() == GEOM::FACE ) {
426 aPixmap->SetPixMap( "ICON_OBJBROWSER_FACE" );
427 aShapeName = "Face_";
428 } else if ( aShape->ShapeType() == GEOM::WIRE ) {
429 aPixmap->SetPixMap( "ICON_OBJBROWSER_WIRE" );
430 aShapeName = "Wire_";
431 } else if ( aShape->ShapeType() == GEOM::EDGE ) {
432 aPixmap->SetPixMap( "ICON_OBJBROWSER_EDGE" );
433 aShapeName = "Edge_";
434 } else if ( aShape->ShapeType() == GEOM::VERTEX ) {
435 aPixmap->SetPixMap( "ICON_OBJBROWSER_VERTEX" );
436 aShapeName = "Vertex_";
438 if (strlen(theName) == 0) aShapeName += TCollection_AsciiString(aResultSO->Tag());
439 else aShapeName = TCollection_AsciiString((char*)theName);
441 //Set a name of the added shape
442 anAttr = aStudyBuilder->FindOrCreateAttribute(aResultSO, "AttributeName");
443 SALOMEDS::AttributeName_var aNameAttrib = SALOMEDS::AttributeName::_narrow(anAttr);
444 aNameAttrib->SetValue(aShapeName.ToCString());
446 //Add a reference to published object
447 // aStudyBuilder->Addreference(theObject, aResultSO);
448 return aResultSO._retn();
452 //============================================================================
454 // purpose : save OCAF/Geom document
455 //============================================================================
456 SALOMEDS::TMPFile* GEOM_Gen_i::Save(SALOMEDS::SComponent_ptr theComponent,
459 SALOMEDS::TMPFile_var aStreamFile;
460 // Get a temporary directory to store a file
461 TCollection_AsciiString aTmpDir = (isMultiFile)?TCollection_AsciiString((char*)theURL):(char*)SALOMEDS_Tool::GetTmpDir().c_str();
462 // Create a list to store names of created files
463 SALOMEDS::ListOfFileNames_var aSeq = new SALOMEDS::ListOfFileNames;
465 // Prepare a file name to open
466 TCollection_AsciiString aNameWithExt("");
468 aNameWithExt = TCollection_AsciiString((char*)SALOMEDS_Tool::GetNameFromPath(theComponent->GetStudy()->URL()).c_str());
469 aNameWithExt += TCollection_AsciiString("_GEOM.sgd");
470 aSeq[0] = aNameWithExt.ToCString();
471 // Build a full file name of temporary file
472 TCollection_AsciiString aFullName = aTmpDir + aNameWithExt;
473 // Save GEOM component in this file
474 myOCAFApp->SaveAs(myCurrentOCAFDoc, aFullName);
475 // Conver a file to the byte stream
476 aStreamFile = SALOMEDS_Tool::PutFilesToStream(aTmpDir.ToCString(), aSeq.in(), isMultiFile);
477 // Remove the created file and tmp directory
478 if (!isMultiFile) SALOMEDS_Tool::RemoveTemporaryFiles(aTmpDir.ToCString(), aSeq.in(), true);
479 // Return the created byte stream
480 return aStreamFile._retn();
483 SALOMEDS::TMPFile* GEOM_Gen_i::SaveASCII(SALOMEDS::SComponent_ptr theComponent,
486 SALOMEDS::TMPFile_var aStreamFile = Save(theComponent, theURL, isMultiFile);
487 return aStreamFile._retn();
491 CORBA::Boolean GEOM_Gen_i::Load(SALOMEDS::SComponent_ptr theComponent,
492 const SALOMEDS::TMPFile& theStream,
496 if (theStream.length() <= 9) {
497 MESSAGE("The TMPFile is too short : " << theStream.length() << " bytes ");
501 // Get a temporary directory for a file
502 TCollection_AsciiString aTmpDir = isMultiFile?TCollection_AsciiString((char*)theURL):(char*)SALOMEDS_Tool::GetTmpDir().c_str();
503 // Conver the byte stream theStream to a file and place it in tmp directory
504 SALOMEDS::ListOfFileNames_var aSeq = SALOMEDS_Tool::PutStreamToFiles(theStream,
508 // Prepare a file name to open
509 TCollection_AsciiString aNameWithExt("");
511 aNameWithExt = TCollection_AsciiString((char*)SALOMEDS_Tool::GetNameFromPath(theComponent->GetStudy()->URL()).c_str());
512 aNameWithExt += TCollection_AsciiString("_GEOM.sgd");
513 TCollection_AsciiString aFullName = aTmpDir + aNameWithExt;
516 if (myOCAFApp->Open(aFullName, myCurrentOCAFDoc) != CDF_RS_OK) return false;
518 // Remove the created file and tmp directory
519 if (!isMultiFile) SALOMEDS_Tool::RemoveTemporaryFiles(aTmpDir.ToCString(), aSeq.in(), true);
521 SALOMEDS::Study_var Study = theComponent->GetStudy();
522 TCollection_AsciiString name( Study->Name() );
524 int StudyID = Study->StudyId();
525 myStudyIDToDoc.Bind( StudyID, myCurrentOCAFDoc );
528 /* We clear all IOR (nameIOR) attributes of all objects before reconstruction */
529 /* This information will be setted when each object is reconstructed */
530 GEOMDS_Commands GC( myCurrentOCAFDoc->Main() ) ;
531 GC.ClearAllIOR(myCurrentOCAFDoc->Main());
536 CORBA::Boolean GEOM_Gen_i::LoadASCII(SALOMEDS::SComponent_ptr theComponent,
537 const SALOMEDS::TMPFile& theStream,
540 return Load(theComponent, theStream, theURL, isMultiFile);
543 // //============================================================================
544 // // function : Save()
545 // // purpose : save OCAF/Geom document
546 // //============================================================================
547 // void GEOM_Gen_i::Save(const char *IORSComponent, const char *aUrlOfFile)
550 // TCollection_ExtendedString path((char*)aUrlOfFile);
551 // TCollection_ExtendedString pathWithExt = path + TCollection_ExtendedString(".sgd");
552 // myOCAFApp->SaveAs(myCurrentOCAFDoc,pathWithExt);
556 // //============================================================================
557 // // function : Load()
558 // // purpose : Load OCAF/Geom document
559 // //============================================================================
560 // void GEOM_Gen_i::Load(const char *IORSComponent, const char *aUrlOfFile)
563 // TCollection_ExtendedString path((char*)aUrlOfFile);
564 // TCollection_ExtendedString pathWithExt = path + TCollection_ExtendedString(".sgd");
566 // myOCAFApp->Open(pathWithExt,myCurrentOCAFDoc);
568 // SALOMEDS::SComponent_var SC = SALOMEDS::SComponent::_narrow(_orb->string_to_object(IORSComponent));
569 // SALOMEDS::Study_var Study = SC->GetStudy();
570 // TCollection_AsciiString name( Study->Name() );
572 // int StudyID = Study->StudyId();
573 // myStudyIDToDoc.Bind( StudyID, myCurrentOCAFDoc );
574 // myStudyID = StudyID;
576 // /* We clear all IOR (nameIOR) attributes of all objects before reconstruction */
577 // /* This information will be setted when each object is reconstructed */
578 // GEOMDS_Commands GC( myCurrentOCAFDoc->Main() ) ;
579 // GC.ClearAllIOR(myCurrentOCAFDoc->Main());
585 //============================================================================
586 // function : Close()
588 //============================================================================
589 void GEOM_Gen_i::Close(SALOMEDS::SComponent_ptr theComponent)
591 int anID = theComponent->GetStudy()->StudyId();
592 if (anID == myStudyID) GetCurrentStudy(0); // set default value of current study ID, if current is deleted
593 if (myStudyIDToDoc.IsBound(anID)) {
594 // close document in the application
595 // Handle(TDocStd_Owner) anOwner;
596 Handle(TDocStd_Document) aDoc = Handle(TDocStd_Document)::DownCast(myStudyIDToDoc.Find(anID));
597 // Handle(TDocStd_Document) anEmptyDoc;
598 // if (aDoc->Main().Root().FindAttribute(TDocStd_Owner::GetID(), anOwner)) {
599 // anOwner->SetDocument(anEmptyDoc);
600 // MESSAGE("********** Nullify owner of document");
602 myOCAFApp->Close(aDoc);
603 myStudyIDToDoc.UnBind(anID); // remove document from GEOM documents data map
607 //============================================================================
608 // function : CanCopy()
610 //============================================================================
611 CORBA::Boolean GEOM_Gen_i::CanCopy(SALOMEDS::SObject_ptr theObject) {
612 // Try to retrieve known by Geometry component GEOM_shape by given IOR
613 SALOMEDS::GenericAttribute_var anAttr;
614 if (!theObject->FindAttribute(anAttr, "AttributeIOR")) return false;
615 GEOM::GEOM_Shape_var aShape = GetIORFromString(SALOMEDS::AttributeIOR::_narrow(anAttr)->Value());
616 // If the object is null one it can't be copied: return false
617 if (aShape->_is_nil()) return false;
621 //============================================================================
622 // function : CopyFrom()
624 //============================================================================
625 SALOMEDS::TMPFile* GEOM_Gen_i::CopyFrom(SALOMEDS::SObject_ptr theObject, CORBA::Long& theObjectID) {
626 // Declare a sequence of the byte to store the copied object
627 SALOMEDS::TMPFile_var aStreamFile = new SALOMEDS::TMPFile;
629 // Try to get GEOM_Shape object by given SObject
630 SALOMEDS::GenericAttribute_var anAttr;
631 if (!theObject->FindAttribute(anAttr, "AttributeIOR")) return false;
632 GEOM::GEOM_Shape_var aShape = GetIORFromString(SALOMEDS::AttributeIOR::_narrow(anAttr)->Value());
633 // If the object is null one it can't be copied: return false
634 if (aShape->_is_nil()) return aStreamFile._retn();
636 GetCurrentStudy(theObject->GetStudy()->StudyId());
638 // Convert a TopoDS_Shape to a stream of bytes
639 TopoDS_Shape aTopology = GetTopoShape(aShape);
640 if (aTopology.IsNull()) return aStreamFile._retn();
641 ostrstream aStreamedShape;
642 BRepTools::Write(aTopology, aStreamedShape);
643 int aSize = aStreamedShape.pcount();
644 char* aBuffer = new char[aSize];
645 memcpy(aBuffer, aStreamedShape.str(), aSize);
646 aStreamedShape.rdbuf()->freeze(0);
648 aStreamFile = new SALOMEDS::TMPFile(aSize, aSize, (CORBA::Octet*)aBuffer, 1);
650 // Assign an ID = 1 the the type GEOM_Shape
653 // Return created TMPFile
654 return aStreamFile._retn();
657 //============================================================================
658 // function : CanPaste()
660 //============================================================================
661 CORBA::Boolean GEOM_Gen_i::CanPaste(const char* theComponentName, CORBA::Long theObjectID) {
662 // The Geometry component can paste only objects copied by Geometry component
663 // and with the object type = 1
664 // MESSAGE("********** GEOM_Gen_i::CanPaste ("<<theComponentName<<","<<theObjectID<<")");
665 if (strcmp(theComponentName, ComponentDataType()) != 0 || theObjectID != 1) return false;
669 //============================================================================
670 // function : PasteInto()
672 //============================================================================
673 SALOMEDS::SObject_ptr GEOM_Gen_i::PasteInto(const SALOMEDS::TMPFile& theStream,
674 CORBA::Long theObjectID,
675 SALOMEDS::SObject_ptr theObject) {
676 // Find the current Study and StudyBuilder
677 SALOMEDS::Study_var aStudy = theObject->GetStudy();
678 SALOMEDS::StudyBuilder_var aStudyBuilder = aStudy->NewBuilder();
680 // Retrieve a TopoDS_Shape from byte stream
681 TopoDS_Shape aTopology;
682 istrstream aStreamedBrep((char*) &theStream[0], theStream.length());
683 BRep_Builder aBuilder;
685 BRepTools::Read(aTopology, aStreamedBrep, aBuilder);
686 } catch (Standard_Failure) {
687 MESSAGE("GEOM_Gen_i::PasteInto exception");
691 // Create new object in Geometry component using retrieved topology
692 GEOM::GEOM_Shape_var aShape = CreateObject(aTopology);
693 GetCurrentStudy(aStudy->StudyId());
694 const char *anEntry = InsertInLabel(aTopology, aShape->Name(), myCurrentOCAFDoc) ;
695 aShape->ShapeId(anEntry) ;
697 // SObject of the created shape is theObject or new Child of Component if theObject == geom component
698 SALOMEDS::SObject_var aNewSO;
699 if (strcmp(theObject->GetFatherComponent()->GetID(),theObject->GetID()) == 0) {
700 aNewSO = aStudyBuilder->NewObject(theObject);
701 } else aNewSO = SALOMEDS::SObject::_duplicate(theObject);
702 // Add IORAttribute to the Study and set IOR of the created GEOM_Shape to it
703 SALOMEDS::GenericAttribute_var anAttr = aStudyBuilder->FindOrCreateAttribute(aNewSO, "AttributeIOR");
704 SALOMEDS::AttributeIOR_var anIOR = SALOMEDS::AttributeIOR::_narrow(anAttr);
705 anIOR->SetValue(aShape->Name());
707 // Return the created in the Study SObject
708 return aNewSO._retn();
711 //============================================================================
712 // function : ComponentDataType()
714 //============================================================================
715 char* GEOM_Gen_i::ComponentDataType()
717 return CORBA::string_dup("GEOM");
720 //============================================================================
721 // function : register()
722 // purpose : register 'name' in 'name_service'
723 //============================================================================
724 void GEOM_Gen_i::register_name(char * name)
726 GEOM::GEOM_Gen_ptr g = GEOM::GEOM_Gen::_narrow(POA_GEOM::GEOM_Gen::_this());
727 name_service->Register(g, name);
732 //================================================================================
733 // function : SequenceOfShapeFromListOfGeomShape()
734 // purpose : Define a sequence of shapes from 'listShapes' and return its length.
735 // : No control is made on shapes !
736 //================================================================================
737 int GEOM_Gen_i::SequenceOfShapeFromListOfGeomShape( const GEOM::GEOM_Gen::ListOfGeomShapes& listShapes,
738 TopTools_SequenceOfShape& SS )
740 int nbShapes = listShapes.length() ;
744 for(int i=0; i<nbShapes; i++) {
745 GEOM::GEOM_Shape_var aGeomShape = listShapes[i] ;
746 TopoDS_Shape aShape = GetTopoShape(aGeomShape) ;
755 //=================================================================================
756 // function : GetTopoShape()
757 // purpose : Returns a TopoDS_Shape from a GEOM::GEOM_Shape_ptr in 'myCurrentOCAFDoc'
758 // : A null shape is returned if not possible
759 //=================================================================================
760 TopoDS_Shape GEOM_Gen_i::GetTopoShape(GEOM::GEOM_Shape_ptr shape_ptr)
765 Handle(TDF_Data) D = myCurrentOCAFDoc->GetData() ;
766 TDF_Tool::Label( D, shape_ptr->ShapeId(), lab, true ) ;
767 Handle(TNaming_NamedShape) NamedShape ;
768 bool res = lab.FindAttribute(TNaming_NamedShape::GetID(), NamedShape) ;
771 return tds ; /* a null shape is returned */
774 return TNaming_Tool::GetShape(NamedShape) ;
780 //=================================================================================
781 // function : GetStringFromIOR()
782 // purpose : returns a string that represents a 'GEOM::GEOM_Shape_var'
783 //=================================================================================
784 const char* GEOM_Gen_i::GetStringFromIOR(GEOM::GEOM_Shape_var shapeIOR) {
785 const char * ret = _orb->object_to_string(shapeIOR) ;
791 //=================================================================================
792 // function : GetIORFromString()
793 // purpose : returns a 'GEOM::GEOM_Shape_var' from a string representing it
794 //=================================================================================
795 GEOM::GEOM_Shape_ptr GEOM_Gen_i::GetIORFromString(const char* stringIOR) {
796 GEOM::GEOM_Shape_var shapeIOR;
797 if(strcmp(stringIOR,"") != 0){
798 CORBA::Object_var anObject = _orb->string_to_object(stringIOR);
799 if(!CORBA::is_nil(anObject))
800 shapeIOR = GEOM::GEOM_Shape::_narrow(anObject.in()) ;
802 return shapeIOR._retn() ;
807 //==================================================================================
808 // function : InsertInLabel()
809 // purpose : Insert S = Shape and mystr = name in a new Label of Geom/OCAF document
810 // : and returns the corresponding OCAF entry
811 //==================================================================================
812 const char * GEOM_Gen_i::InsertInLabel(TopoDS_Shape S, const char *mystr, Handle(TDocStd_Document) OCAFDoc)
814 GEOMDS_Commands GC(OCAFDoc->Main());
815 /* add attributs S and mystr in a new label */
816 TDF_Label Lab = GC.AddShape (S, (char*)mystr);
818 TCollection_AsciiString entry;
819 TDF_Tool::Entry(Lab,entry);
820 const char *ent = entry.ToCString() ;
825 //==================================================================================
826 // function : InsertInLabelDependentShape()
827 // purpose : Insert S = Shape and its nameIor in a new Label of Geom/OCAF document
828 // : insert also a reference attribute (a label) to the main shape 'mainshap_ptr'.
829 // : and returns the corresponding OCAF entry of the new label.
830 //==================================================================================
831 const char * GEOM_Gen_i::InsertInLabelDependentShape( TopoDS_Shape S,
833 GEOM::GEOM_Shape_ptr mainshape_ptr,
834 Handle(TDocStd_Document) OCAFDoc )
836 GEOMDS_Commands GC(OCAFDoc->Main());
837 /* add attributs S and nameIor in a new label */
839 /* retrieve the label of the main shape in the document */
840 TDF_Label mainRefLab;
841 TDF_Tool::Label(OCAFDoc->GetData(), mainshape_ptr->ShapeId(), mainRefLab);
843 /* add attributs : S, nameIor and ref to main */
844 TDF_Label Lab = GC.AddDependentShape(S, (char*)nameIor, mainRefLab);
846 TCollection_AsciiString entry;
847 TDF_Tool::Entry(Lab, entry);
848 const char *ent = entry.ToCString() ;
853 //=================================================================================
854 // function : InsertInLabelOneArgument()
856 //=================================================================================
857 void GEOM_Gen_i::InsertInLabelOneArgument(TopoDS_Shape main_topo,
858 GEOM::GEOM_Shape_ptr shape_ptr,
859 TopoDS_Shape result_topo,
860 GEOM::GEOM_Shape_ptr result,
861 Handle(TDocStd_Document) OCAFDoc)
863 /* Put shape and name into geom/OCAF doc */
864 GEOMDS_Commands GC(OCAFDoc->Main());
865 /* Add attributs 'shape' and 'name_ior' in a new label */
866 TDF_Label Lab = GC.Generated( main_topo, result_topo, result->Name() );
867 TCollection_AsciiString entry;
868 TDF_Tool::Entry(Lab, entry);
869 result->ShapeId( entry.ToCString() ) ;
871 /* Create a new label */
872 TDF_Label NewLab = Lab.NewChild();
873 TCollection_ExtendedString Value("Arguments");
874 TDataStd_Name::Set(NewLab,Value);
876 TDF_Label NewLab1 = NewLab.NewChild();
878 TDF_Tool::Label(OCAFDoc->GetData(), shape_ptr->ShapeId(), RefLab);
879 TDF_Reference::Set(NewLab1, RefLab);
883 //=================================================================================
884 // function : InsertInLabelMoreArguments()
886 //=================================================================================
887 void GEOM_Gen_i::InsertInLabelMoreArguments(TopoDS_Shape main_topo,
888 GEOM::GEOM_Shape_ptr result,
889 const GEOM::GEOM_Gen::ListOfIOR& ListShapes,
890 Handle(TDocStd_Document) OCAFDoc)
892 /* Put shape and name into geom/OCAF doc */
893 GEOMDS_Commands GC(OCAFDoc->Main());
894 /* Add attributs TopoDS and name_ior in a new label */
895 TDF_Label Lab = GC.AddShape(main_topo, result->Name() );
896 TCollection_AsciiString entry;
897 TDF_Tool::Entry(Lab, entry);
899 /* Create a new label */
900 TDF_Label NewLab = Lab.NewChild();
901 TCollection_ExtendedString Value("Arguments");
902 TDataStd_Name::Set(NewLab, Value);
904 for (unsigned int ind = 0; ind < ListShapes.length(); ind++) {
906 TDF_Label NewLab1 = NewLab.NewChild();
907 GEOM::GEOM_Shape_var aShape = GetIORFromString( ListShapes[ind] );
910 TDF_Tool::Label(OCAFDoc->GetData(), aShape->ShapeId(), RefLab);
911 TDF_Reference::Set(NewLab1, RefLab);
913 result->ShapeId(entry.ToCString());
918 //=================================================================================
919 // function: NbLabels()
921 //=================================================================================
922 CORBA::Short GEOM_Gen_i::NbLabels()
924 TDF_ChildIterator ChildIterator(myCurrentOCAFDoc->Main());
927 while (ChildIterator.More()) {
929 ChildIterator.Next();
932 // return TDF_Tool::NbLabels( myCurrentOCAFDoc->Main() );
937 //=================================================================================
938 // function: GetCurrentStudy()
939 // purpose : Finds or creates the geom/OCAF document corresponding to the index
941 //=================================================================================
942 void GEOM_Gen_i::GetCurrentStudy(CORBA::Long StudyID)
944 /* If StudyID is known we link myCurrentOCAFDoc to it */
945 if (myStudyIDToDoc.IsBound(StudyID)) {
946 myCurrentOCAFDoc = Handle(TDocStd_Document)::DownCast(myStudyIDToDoc(StudyID));
948 /* Create a new OCAFDoc and link it to 'StudyID' argument */
950 myOCAFApp->NewDocument("SALOME_GEOM",myCurrentOCAFDoc);
951 myStudyIDToDoc.Bind(StudyID,myCurrentOCAFDoc);
957 //================================================================================
958 // function : CreateObject()
959 // purpose : private function to create a complete CORBA object and return it
960 //================================================================================
961 GEOM::GEOM_Shape_ptr GEOM_Gen_i::CreateObject(TopoDS_Shape& tds)
963 // if ( tds.ShapeType() == TopAbs_COMPOUND ) {
964 // TopoDS_Iterator itr(tds);
967 // while (itr.More()) {
969 // res = itr.Value();
977 GEOM::shape_type st = GEOM::shape_type(tds.ShapeType()) ; /* casting */
979 /* Create the CORBA servant holding the TopoDS_Shape */
980 GEOM::GEOM_Gen_ptr engine = POA_GEOM::GEOM_Gen::_this();
981 GEOM::GEOM_Shape::ListOfSubShapeID_var index = new GEOM::GEOM_Shape::ListOfSubShapeID;
983 GEOM_Shape_i * shape_servant = new GEOM_Shape_i(tds, _orb, engine, index, st, true);
984 GEOM::GEOM_Shape_var shape = GEOM::GEOM_Shape::_narrow(shape_servant->_this());
986 /* Create and set the name (IOR of shape converted into a string) */
987 string name_ior = _orb->object_to_string(shape) ;
988 shape->Name( name_ior.c_str() );
989 shape->NameType( "" );
993 //=======================================================================
994 //function : CreateSubObject
996 //=======================================================================
998 GEOM::GEOM_Shape_ptr GEOM_Gen_i::CreateSubObject(const TopoDS_Shape& SubShape,
999 const GEOM::GEOM_Shape_ptr MainShape,
1000 const GEOM::GEOM_Shape::ListOfSubShapeID& ListOfID)
1002 GEOM::shape_type st = GEOM::shape_type(SubShape.ShapeType()) ; /* casting */
1004 /* Create the CORBA servant holding the TopoDS_Shape */
1005 GEOM::GEOM_Gen_ptr engine = POA_GEOM::GEOM_Gen::_this();
1006 GEOM_Shape_i * shape_servant =
1007 new GEOM_Shape_i(SubShape, _orb, engine, ListOfID, st, false);
1008 GEOM::GEOM_Shape_var shape = GEOM::GEOM_Shape::_narrow(shape_servant->_this());
1010 /* Create and set the name (IOR of shape converted into a string) */
1011 string name_ior = _orb->object_to_string(shape) ;
1012 shape->Name( name_ior.c_str() );
1013 /* create and set the mainname (IOR of shape converted into a string) */
1014 const char *mainname_ior = _orb->object_to_string(MainShape) ;
1015 shape->MainName(mainname_ior);
1016 /* precaution : NameType will be set precisely in GUI */
1017 shape->NameType( "" );
1018 /* add 'SubShape' its 'nameIOR' and a reference to the main shape thanks to method below */
1020 InsertInLabelDependentShape(SubShape, shape->Name(), MainShape, myCurrentOCAFDoc) ;
1021 shape->ShapeId( entry ) ;
1026 //=======================================================================
1027 // function : SuppressFacesGlue()
1028 // purpose : Define a compound of shells after suppress of mapFaces in the
1029 // : shape S and return the number of shells of the compound.
1030 //=======================================================================
1031 int GEOM_Gen_i::SuppressFacesGlue( const TopoDS_Shape& S,
1032 const TopTools_MapOfShape& mapFaces,
1033 TopoDS_Shape& aCompoundOfShells )
1034 throw (SALOME::SALOME_Exception)
1036 Unexpect aCatch(SALOME_SalomeException);
1037 BRepTools_Quilt Glue;
1038 aCompoundOfShells.Nullify() ;
1040 for ( TopExp_Explorer exp( S, TopAbs_FACE); exp.More(); exp.Next() ) {
1041 const TopoDS_Face& F = TopoDS::Face(exp.Current());
1042 if ( !mapFaces.Contains(F) ) {
1043 /* this face must not to be suppressed */
1048 /* Use specif method to calculate the compound of shells */
1049 aCompoundOfShells = Glue.Shells();
1051 if( aCompoundOfShells.ShapeType() != TopAbs_COMPOUND ) {
1052 THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::SuppressFacesGlue() : not a GEOM::COMPOUND", SALOME::BAD_PARAM);
1055 /* explore Compound for verification and return the number of shells */
1056 int numberOfShell = 0 ;
1057 for ( TopExp_Explorer exp1( aCompoundOfShells, TopAbs_SHELL); exp1.More(); exp1.Next() )
1060 return numberOfShell ;
1064 //=====================================================================================
1065 // function : GetIndexTopology()
1066 // purpose : return the index of a sub shape in a shape (index starts at 1)
1067 // : Return -1 if not found
1068 //=====================================================================================
1069 int GEOM_Gen_i::GetIndexTopology(const TopoDS_Shape& subshape, const TopoDS_Shape& mainShape)
1071 if( mainShape.IsNull() || subshape.IsNull() )
1075 if (subshape.ShapeType() == TopAbs_COMPOUND)
1078 TopTools_ListOfShape CL;
1079 CL.Append( mainShape );
1080 TopTools_ListIteratorOfListOfShape itC;
1081 for (itC.Initialize( CL ); itC.More(); itC.Next())
1083 for (it.Initialize( itC.Value() ); it.More(); it.Next())
1085 if ( it.Value().ShapeType() == TopAbs_COMPOUND)
1087 if (it.Value().IsSame(subshape))
1091 CL.Append( it.Value() );
1098 TopExp_Explorer Exp ( mainShape, subshape.ShapeType() );
1099 TopTools_MapOfShape M;
1100 while ( Exp.More() )
1102 if ( M.Add(Exp.Current()) )
1104 if ( Exp.Current().IsSame(subshape) )
1115 //================================================================================
1116 // function : IndexOfFacesOfSubShell()
1117 // purpose : Return a list of indices corresponding to the faces of a 'subShell'
1118 // : in the main shape 'S'
1119 //================================================================================
1120 GEOM::GEOM_Shape::ListOfSubShapeID* GEOM_Gen_i::IndexOfFacesOfSubShell( const TopoDS_Shape& S,
1121 const TopoDS_Shape subShell )
1122 throw (SALOME::SALOME_Exception)
1125 Unexpect aCatch(SALOME_SalomeException);
1126 GEOM::GEOM_Shape::ListOfSubShapeID_var ListOfID = new GEOM::GEOM_Shape::ListOfSubShapeID;
1127 ListOfID->length(0) ;
1128 if( subShell.IsNull() || subShell.ShapeType() != TopAbs_SHELL ) {
1129 THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::IndexOfFacesOfSubShell() : null shape or not a GEOM::SHELL", SALOME::BAD_PARAM);
1132 /* put faces of subShell in a Map of faces */
1134 TopTools_MapOfShape mapFaces ;
1135 for( TopExp_Explorer Exp1( subShell, TopAbs_FACE ); Exp1.More(); Exp1.Next() ) {
1136 mapFaces.Add(Exp1.Current() ) ;
1141 THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::IndexOfFacesOfSubShell() : no faces in sub shell", SALOME::BAD_PARAM);
1143 /* Find index of each face of subshell in the main topology and put its index in ListOfID */
1145 for ( TopExp_Explorer Exp2(S, TopAbs_FACE); Exp2.More(); Exp2.Next() ) {
1147 const TopoDS_Face& F = TopoDS::Face( Exp2.Current() ) ;
1149 if( mapFaces.Contains(F) ) {
1150 int n = GetIndexTopology( F, S ) ;
1152 THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::IndexOfFacesOfSubShell() : no index found", SALOME::BAD_PARAM);
1155 ListOfID->length(size) ;
1156 ListOfID[size-1] = n ;
1160 return ListOfID._retn() ;
1165 //================================================================================
1166 // function : ListOfIDIntoMapOfShapes()
1167 // purpose : Define a MapOfShapes from a main topology 'S' a 'subShapeType'
1168 // : and a list of indices 'L'.
1169 // : Return true if 'aMap' is not empty
1170 //================================================================================
1171 bool GEOM_Gen_i::ListOfIDIntoMapOfShapes( const TopoDS_Shape& S,
1172 const GEOM::GEOM_Shape::ListOfSubShapeID& L,
1173 const int subShapeType,
1174 TopTools_MapOfShape& aMap )
1176 if( L.length() < 1 || S.IsNull() ) {
1181 for( int k=0; k<L.length(); k++ ) {
1182 /* indices start at 1 in list L */
1184 TopExp_Explorer exp ;
1185 TopTools_MapOfShape M;
1186 for( exp.Init( S, TopAbs_ShapeEnum(subShapeType) ); exp.More(); exp.Next() ) {
1187 if ( M.Add(exp.Current()) )
1190 aMap.Add( exp.Current() ) ;
1201 //================================================================================
1202 // function : ListOfIDIntoSequenceOfShapes()
1203 // purpose : Define 'aSequenceOfShapes' from a main topology 'S' a 'subShapeType'
1204 // : and a list of indices 'L'.
1205 // : Return true if 'aSequenceOfShapes' is not empty
1206 //================================================================================
1207 bool GEOM_Gen_i::ListOfIDIntoSequenceOfShapes( const TopoDS_Shape& S,
1208 const GEOM::GEOM_Shape::ListOfSubShapeID& L,
1209 const int subShapeType,
1210 TopTools_SequenceOfShape& aSequenceOfShapes )
1212 if( L.length() < 1 || S.IsNull() ) {
1216 aSequenceOfShapes.Clear() ;
1217 for( int k=0; k<L.length(); k++ ) {
1218 /* indices start at 1 in list L */
1220 TopExp_Explorer exp ;
1221 for( exp.Init( S, TopAbs_ShapeEnum(subShapeType) ); exp.More(); exp.Next() ) {
1223 aSequenceOfShapes.Append( exp.Current() ) ;
1233 //================================================================================
1234 // function : SuppressFaces()
1235 // purpose : Suppress faces contained in ListOfID from 'shape'.
1236 // : Return a list of Geom shapes each one is a main shape GEOM::FACE or GEOM::SHELL
1237 //================================================================================
1238 GEOM::GEOM_Gen::ListOfGeomShapes* GEOM_Gen_i::SuppressFaces( GEOM::GEOM_Shape_ptr shape,
1239 const GEOM::GEOM_Shape::ListOfSubShapeID& ListOfID )
1240 throw (SALOME::SALOME_Exception)
1242 Unexpect aCatch(SALOME_SalomeException);
1243 GEOM::GEOM_Gen::ListOfGeomShapes_var listOfGeomShapes = new GEOM::GEOM_Gen::ListOfGeomShapes;
1244 listOfGeomShapes->length(0) ;
1246 TopoDS_Shape mainShape = GetTopoShape(shape);
1247 if( mainShape.IsNull() )
1248 THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::SuppressFaces() : null argument shape", SALOME::BAD_PARAM);
1250 if( ListOfID.length() < 1 )
1251 THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::SuppressFaces() : empty ListOfID", SALOME::BAD_PARAM);
1253 /* Define 'mapFaces' a map of faces to be suppressed in mainShape */
1254 TopTools_MapOfShape mapFaces ;
1255 if( !ListOfIDIntoMapOfShapes(mainShape, ListOfID, TopAbs_FACE, mapFaces ) ) {
1256 return listOfGeomShapes._retn();
1259 /* Call algorithm to calculate a compound of shells resulting of face suppression */
1260 int numberOfShells = 0 ;
1261 TopoDS_Shape aCompoundOfShells ;
1262 numberOfShells = SuppressFacesGlue(mainShape, mapFaces, aCompoundOfShells) ;
1263 if(numberOfShells < 1) {
1264 THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::Suppressfaces() : no shells", SALOME::BAD_PARAM);
1267 /* Create a shell for each shell contained in 'aCompoundOfShells' and */
1268 /* put it in the list of GeomShapes to be returned. */
1269 /* But if the shell is composed of only a face we create a face and not a shell */
1271 for( TopExp_Explorer exp(aCompoundOfShells, TopAbs_SHELL); exp.More(); exp.Next() ) {
1273 const TopoDS_Shell& aShell = TopoDS::Shell( exp.Current() );
1274 if( aShell.IsNull() ) {
1275 THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::Suppressfaces() : null shell", SALOME::BAD_PARAM);
1278 GEOM::GEOM_Shape::ListOfSubShapeID_var aList = new GEOM::GEOM_Shape::ListOfSubShapeID;
1279 aList = IndexOfFacesOfSubShell(mainShape, aShell) ;
1281 if( aList->length() < 1 ) {
1282 THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::SuppressFaces() : aList is empty", SALOME::BAD_PARAM);
1285 TopoDS_Shape aShellOrFace ;
1286 /* Only a face into the shell : we create a single face instead of a shell : 'aList' is unchanged */
1287 if( aList->length() == 1 ) {
1288 TopExp_Explorer exp ;
1289 exp.Init( aShell, TopAbs_FACE ) ;
1291 aShellOrFace = exp.Current() ;
1294 aShellOrFace = aShell ;
1297 /* Create CORBA object */
1298 GEOM::GEOM_Shape_var result = CreateObject(aShellOrFace) ;
1299 if( CORBA::is_nil(result) ) {
1300 THROW_SALOME_CORBA_EXCEPTION("Suppress Faces aborted : null result", SALOME::BAD_PARAM);
1303 InsertInLabelOneArgument(mainShape, shape, aShellOrFace, result, myCurrentOCAFDoc) ;
1305 listOfGeomShapes->length(i) ;
1306 listOfGeomShapes[i-1] = result ;
1309 return listOfGeomShapes._retn() ;
1314 //================================================================================
1315 // function : IsShapeInSequence()
1316 // purpose : return true is aShape is in SS. The test method is 'IsSame()'
1317 //================================================================================
1318 bool GEOM_Gen_i::IsShapeInSequence(const TopTools_SequenceOfShape& SS, const TopoDS_Shape& aShape)
1320 if( aShape.IsNull() || SS.IsEmpty() )
1322 for( int i=1; i<=SS.Length(); i++) {
1323 if( SS.Value(i).IsSame(aShape) )
1330 //================================================================================
1331 // function : FreeEdgesFromMapOfFace()
1332 // purpose : Define MS a map of all edges of faces of 'MSfaces'
1333 // : All multiple edges are removed !
1334 //================================================================================
1335 void GEOM_Gen_i::FreeEdgesFromMapOfFace( const TopTools_MapOfShape& MSfaces,
1336 TopTools_MapOfShape& MS )
1339 TopTools_MapOfShape Multiple ;
1340 TopTools_MapIteratorOfMapOfShape it ;
1341 for( it.Initialize(MSfaces); it.More(); it.Next() ) {
1342 TopoDS_Shape aFace = it.Key() ;
1343 TopExp_Explorer exp ;
1344 for( exp.Init( aFace, TopAbs_EDGE); exp.More(); exp.Next() ) {
1345 if( !Multiple.Contains( exp.Current() ) && !MS.Add( exp.Current() ) ) {
1346 MS.Remove( exp.Current() ) ;
1347 Multiple.Add( exp.Current() ) ;
1355 //================================================================================
1356 // function : MapRemoveSequence()
1357 // purpose : In term of shapes ST = MS - SSRemove
1359 //================================================================================
1360 void GEOM_Gen_i::MapRemoveSequence( const TopTools_MapOfShape& MS,
1361 const TopTools_SequenceOfShape& SSRemove,
1362 TopTools_SequenceOfShape& ST )
1365 TopTools_MapIteratorOfMapOfShape it ;
1366 for( it.Initialize(MS); it.More(); it.Next() ) {
1367 TopoDS_Shape aShape = it.Key() ;
1368 if( !IsShapeInSequence( SSRemove, aShape ) )
1369 ST.Append( aShape ) ;
1376 //================================================================================
1377 // function : SuppressHoleSubRoutine()
1378 // purpose : Define recursively 'MSfacesSuppress' a list of faces to suppress in a hole
1379 //================================================================================
1380 void GEOM_Gen_i::SuppressHoleSubRoutine( const TopoDS_Shape& mainShape,
1381 const TopoDS_Face& aFace,
1382 const TopTools_SequenceOfShape& SSedgesOfWire,
1383 const TopTools_IndexedDataMapOfShapeListOfShape& aMapEdgesFaces,
1384 const TopTools_MapOfShape& MSfaces,
1385 TopTools_MapOfShape& MSfacesSuppress,
1386 const Standard_Boolean withEndFace,
1387 const TopoDS_Face& endFace,
1388 TopTools_MapOfShape& MSwireEndEdges )
1389 throw (SALOME::SALOME_Exception)
1391 Unexpect aCatch(SALOME_SalomeException);
1392 TopTools_MapOfShape MS ;
1393 TopTools_SequenceOfShape SU ;
1394 FreeEdgesFromMapOfFace(MSfaces, MS) ; /* MS = free edges of MSfaces */
1395 MapRemoveSequence(MS, SSedgesOfWire, SU) ; /* SU = MS - SSedgesOfWire */
1397 if( SU.IsEmpty() ) {
1401 /* Here SU contains new edges to find new faces to suppress */
1402 /* Define the list of faces of SU edges that aren't in faces of MSfacesSuppress in order to add into */
1403 /* For each edge we have a map of all its faces : it's in 'aMapEdgesFaces' */
1404 TopTools_MapOfShape MSfacesTmp ;
1405 for( int v=1; v<=SU.Length(); v++ ) {
1406 TopoDS_Shape E = SU.Value(v) ;
1408 TopTools_ListOfShape LF ;
1409 int ind = aMapEdgesFaces.FindIndex(E) ;
1411 /* LF is the list of faces for an edge of SU : may be empty no matter */
1412 LF = aMapEdgesFaces.FindFromIndex(ind) ;
1414 TopTools_ListIteratorOfListOfShape it ;
1415 for( it.Initialize(LF); it.More(); it.Next() ) {
1417 if( withEndFace == false ) {
1418 if( F.IsSame(aFace) )
1419 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHoleSubRoutine() : hole traversing or ?", SALOME::BAD_PARAM);
1420 if( !MSfacesSuppress.Contains(F) ) {
1421 MSfacesSuppress.Add(F) ;
1422 MSfacesTmp.Add(F) ; // Dont remove the 'if' !
1425 else { /* withEndFace == true */
1426 if( F.IsSame(aFace) && !F.IsSame(endFace) )
1427 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHoleSubRoutine() : hole traversing incoherent ?", SALOME::BAD_PARAM);
1429 if( F.IsSame(endFace) ) {
1430 /* We have reached endFace if selection was correct so we add */
1431 /* edge in a map to find later the corresponding endWire (hole */
1432 MSwireEndEdges.Add(E) ;
1435 if( !MSfacesSuppress.Contains(F) ) {
1436 MSfacesSuppress.Add(F) ;
1443 /* Call recursively this routine */
1444 SuppressHoleSubRoutine( mainShape, aFace, SSedgesOfWire, aMapEdgesFaces, MSfacesTmp, MSfacesSuppress, withEndFace, endFace, MSwireEndEdges ) ;
1449 //================================================================================
1450 // function : GetShapeFromIndex()
1451 // purpose : Find 'tds' a sub shape of 'aShape' according to 'aList' that contains
1452 // : a unique index !
1453 // : Warning : index must be setted with the same exploration logic !
1454 // : So 'index' is calculated with no shape doublons !
1455 //================================================================================
1456 bool GEOM_Gen_i::GetShapeFromIndex( const TopoDS_Shape& aShape,
1457 const TopAbs_ShapeEnum aType,
1462 if (aShape.IsNull() || index < 1)
1464 /* Indices start at 1 */
1466 bool found = false ;
1467 TopExp_Explorer exp ;
1468 TopTools_MapOfShape M;
1469 for( exp.Init( aShape, aType ); exp.More(); exp.Next() ) {
1470 if( M.Add(exp.Current()) ) { /* if not a doublon : we compare */
1472 tds = exp.Current() ;
1483 //================================================================================
1484 // function : SuppressHolesInFaceOrShell() Main method.
1485 // purpose : Suppress holes identified by wires in a single face or shell
1487 //================================================================================
1488 GEOM::GEOM_Shape_ptr GEOM_Gen_i::SuppressHolesInFaceOrShell( GEOM::GEOM_Shape_ptr shapeFaceShell,
1489 const GEOM::GEOM_Shape::ListOfSubShapeID& ListIdWires )
1490 throw (SALOME::SALOME_Exception)
1492 Unexpect aCatch(SALOME_SalomeException);
1493 GEOM::GEOM_Shape_var result;
1495 if( ListIdWires.length() < 1 )
1496 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHolesInFace : no holes selected", SALOME::BAD_PARAM);
1498 const TopoDS_Shape tds = GetTopoShape(shapeFaceShell) ;
1499 if( tds.IsNull() || !IsValid(tds) )
1500 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHolesInFace() : non valid main argument", SALOME::BAD_PARAM);
1502 /* Create a map of wires/holes to suppress */
1503 TopTools_MapOfShape MapHoles ;
1504 for ( int i = 0; i < ListIdWires.length(); i++ ) {
1506 if( !GetShapeFromIndex( tds, TopAbs_WIRE, ListIdWires[i], W ) )
1507 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHolesInFace() : bad index ?", SALOME::BAD_PARAM);
1511 /* Test if argument is a face or shell */
1513 if( tds.ShapeType() == TopAbs_FACE )
1515 else if ( tds.ShapeType() == TopAbs_SHELL )
1518 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHolesInFace() : not a face or a shell", SALOME::BAD_PARAM);
1520 /* Define two maps : all faces and faces to that will be modified */
1521 TopTools_MapOfShape MapFacesToModify ;
1522 TopTools_MapOfShape MapFacesAll ;
1523 TopExp_Explorer expF ;
1524 for( expF.Init( tds, TopAbs_FACE); expF.More(); expF.Next() ) {
1525 TopoDS_Face F = TopoDS::Face( expF.Current() ) ;
1526 MapFacesAll.Add(F) ;
1527 TopExp_Explorer expW ;
1528 for( expW.Init( F, TopAbs_WIRE); expW.More(); expW.Next() ) {
1529 TopoDS_Wire W = TopoDS::Wire( expW.Current() ) ;
1530 if( MapHoles.Contains(W) ) {
1531 MapFacesToModify.Add(F) ;
1536 /* Define faces not modified */
1537 TopTools_MapOfShape MapFacesNotModified ;
1538 TopTools_MapIteratorOfMapOfShape it ;
1539 for( it.Initialize(MapFacesAll); it.More(); it.Next() ) {
1540 TopoDS_Face FF = TopoDS::Face( it.Key() ) ;
1541 if( !MapFacesToModify.Contains(FF) )
1542 MapFacesNotModified.Add(FF) ;
1545 if( MapFacesToModify.IsEmpty() )
1546 THROW_SALOME_CORBA_EXCEPTION("Error : empty map of faces", SALOME::BAD_PARAM);
1548 if( isFace && MapFacesToModify.Extent() != 1 )
1549 THROW_SALOME_CORBA_EXCEPTION("Incoherent", SALOME::BAD_PARAM);
1551 /* Main argument is a face */
1552 if( isFace && MapFacesToModify.Extent() == 1 ) {
1553 TopoDS_Face resultFace ;
1554 if( !RebuildFaceRemovingHoles( TopoDS::Face(tds), MapHoles, resultFace ) )
1555 THROW_SALOME_CORBA_EXCEPTION(" Problem : !RebuildFaceRemovingHoles()", SALOME::BAD_PARAM);
1556 /* Creation of CORBA object : face topology */
1557 result = CreateObject(resultFace);
1558 InsertInLabelOneArgument(tds, shapeFaceShell, resultFace, result, myCurrentOCAFDoc) ;
1562 /* Main argument is a shell : rebuild faces modified */
1563 TopTools_MapOfShape MapFacesModified ;
1564 for( it.Initialize(MapFacesToModify); it.More(); it.Next() ) {
1565 TopoDS_Face FF = TopoDS::Face( it.Key() ) ;
1567 if( !RebuildFaceRemovingHoles( FF, MapHoles, resF ) )
1568 THROW_SALOME_CORBA_EXCEPTION(" Problem shell : !RebuildFaceRemovingHoles()", SALOME::BAD_PARAM);
1569 MapFacesModified.Add(resF) ;
1572 /* Rebuild the shell with faces modified and non modified */
1573 TopoDS_Shell resultShell ;
1575 B.MakeShell(resultShell) ;
1576 TopTools_MapIteratorOfMapOfShape it1 ;
1577 for( it1.Initialize(MapFacesModified); it1.More(); it1.Next() )
1578 B.AddShellFace( resultShell,TopoDS::Face( it1.Key() ) ) ;
1579 for( it1.Initialize(MapFacesNotModified); it1.More(); it1.Next() )
1580 B.AddShellFace( resultShell,TopoDS::Face( it1.Key() ) ) ;
1582 B.CompleteShell(resultShell) ;
1584 if( resultShell.IsNull() )
1585 THROW_SALOME_CORBA_EXCEPTION("Null or not valid result Shell", SALOME::BAD_PARAM) ;
1587 /* Creation of CORBA object : shell topology */
1588 result = CreateObject(resultShell);
1589 InsertInLabelOneArgument(tds, shapeFaceShell, resultShell, result, myCurrentOCAFDoc) ;
1594 //================================================================================
1595 // function : RebuildFaceRemovingHoles()
1596 // purpose : Rebuild a face removing holes that are in 'mapHoles'.
1597 // : NB : 'mapHoles' may content more holes than necessary
1598 //================================================================================
1599 bool GEOM_Gen_i::RebuildFaceRemovingHoles( const TopoDS_Face& aFace,
1600 const TopTools_MapOfShape& mapHoles,
1601 TopoDS_Shape& resultFace )
1603 /* Get the outer wire of the face 'aFace' */
1604 TopoDS_Wire outW = BRepTools::OuterWire( aFace ) ;
1605 if( outW.IsNull() || !IsValid(outW) )
1606 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHolesInFace : bad outer wire of 'aFace'", SALOME::BAD_PARAM);
1608 /* Rebuild a face avoiding holes in the map 'mapHoles' */
1609 Handle(Geom_Surface) Surface = BRep_Tool::Surface(aFace) ;
1610 TopoDS_Face F2 = BRepBuilderAPI_MakeFace( Surface, outW, true ) ;
1612 if( F2.Orientation() != aFace.Orientation() )
1613 F2.Orientation( aFace.Orientation() ) ;
1615 BRepBuilderAPI_MakeFace aBuilder( F2 ) ;
1616 bool foundAndKeepHoles = false ;
1617 TopExp_Explorer exp ;
1619 for( exp.Init( aFace, TopAbs_WIRE); exp.More(); exp.Next() ) {
1620 TopoDS_Wire hole = TopoDS::Wire( exp.Current() ) ;
1621 if( !mapHoles.Contains(hole) && !exp.Current().IsEqual(outW) ) {
1622 aBuilder.Add( hole) ;
1623 if( !aBuilder.IsDone() )
1624 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHolesInFace : builder problem !", SALOME::BAD_PARAM);
1626 resultFace = TopoDS::Face(aBuilder) ;
1627 foundAndKeepHoles = true ;
1631 if( !foundAndKeepHoles )
1634 resultFace = TopoDS::Face(aBuilder) ;
1642 //================================================================================
1643 // function : SuppressHole() Main method.
1644 // purpose : Suppress an hole identified by a wire in a face of shape
1645 // : ListIdFace contains a unique index of face in shape
1646 // : ListIdWire contains a unique index of wire in face !!!
1647 // : ListIdEndFace is used only when hole traverse.
1648 //================================================================================
1649 GEOM::GEOM_Shape_ptr GEOM_Gen_i::SuppressHole( GEOM::GEOM_Shape_ptr shape,
1650 const GEOM::GEOM_Shape::ListOfSubShapeID& ListIdFace,
1651 const GEOM::GEOM_Shape::ListOfSubShapeID& ListIdWire,
1652 const GEOM::GEOM_Shape::ListOfSubShapeID& ListIdEndFace )
1653 throw (SALOME::SALOME_Exception)
1655 Unexpect aCatch(SALOME_SalomeException);
1656 GEOM::GEOM_Shape_var result;
1659 TopoDS_Face endFace ;
1663 /* Retrieve 'aShape' the initial main shape selection */
1664 const TopoDS_Shape aShape = GetTopoShape(shape);
1666 if( !IsValid(aShape) )
1667 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHole() : non valid main shape", SALOME::BAD_PARAM);
1669 if( ListIdFace.length() != 1 || ListIdWire.length() != 1 )
1670 THROW_SALOME_CORBA_EXCEPTION("bad list", SALOME::BAD_PARAM);
1672 /* Retrieve 'aFace' selection */
1673 if( !GetShapeFromIndex( aShape, TopAbs_FACE, ListIdFace[0], tmp ) ) {
1674 THROW_SALOME_CORBA_EXCEPTION("face not found", SALOME::BAD_PARAM);
1677 aFace = TopoDS::Face(tmp) ;
1679 if( !IsValid(aFace) )
1680 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHole() : face shape not valid", SALOME::BAD_PARAM);
1682 /* Retrieve 'aWire' selection : Warning : index of wire refers to the face ! */
1684 if( !GetShapeFromIndex( aFace, TopAbs_WIRE, ListIdWire[0], aTmp ) ) {
1685 THROW_SALOME_CORBA_EXCEPTION("wire not found", SALOME::BAD_PARAM);
1688 aWire = TopoDS::Wire(aTmp) ;
1690 if( !IsValid(aWire) )
1691 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHole() : bad wire" , SALOME::BAD_PARAM);
1693 /* Get the outer wire of aFace */
1694 TopoDS_Wire outerW = BRepTools::OuterWire( aFace ) ;
1695 if( outerW.IsNull() || !IsValid(outerW) )
1696 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHole() : bad outer wire", SALOME::BAD_PARAM);
1698 /* Test bad user selection aWire */
1699 if( aWire.IsSame(outerW) )
1700 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHole() : outerW = aWire", SALOME::BAD_PARAM);
1702 /* Test if 'endFace' is used as argument and seems to be a valid one */
1703 /* NB : 'endFace' is optional and used when hole to suppress traverse 'aShape' */
1704 if( ListIdEndFace.length() == 0 ) {
1705 withEndFace = false ;
1708 TopoDS_Shape aTemp ;
1709 if( !GetShapeFromIndex( aShape, TopAbs_FACE, ListIdEndFace[0], aTemp ) || tmp.IsNull() || !IsValid(aTemp) )
1710 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHole() : non valid endFace", SALOME::BAD_PARAM);
1712 /* Test if 'endFace' as at least one hole */
1713 endFace = TopoDS::Face(aTemp) ;
1715 TopExp_Explorer fExp ;
1717 for( fExp.Init(endFace, TopAbs_WIRE); fExp.More(); fExp.Next() ) {
1718 TopoDS_Wire W = TopoDS::Wire( fExp.Current() ) ;
1719 if( !W.IsNull() && IsValid(W) )
1723 withEndFace = true ; /* at least 2 wires : outer wire plus an hole or more */
1725 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHole() : end face selection ?", SALOME::BAD_PARAM);
1728 /* Find edges of aWire and test if degenerated */
1729 TopTools_SequenceOfShape SSedgesOfWire ;
1730 TopExp_Explorer wireExp ;
1731 for( wireExp.Init(aWire, TopAbs_EDGE); wireExp.More(); wireExp.Next() ) {
1732 TopoDS_Edge E = TopoDS::Edge( wireExp.Current() ) ;
1733 if( E.IsNull() || BRep_Tool::Degenerated(E) ) {
1734 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SupressHole() : found bad edge", SALOME::BAD_PARAM);
1737 SSedgesOfWire.Append( wireExp.Current() ) ;
1740 if( SSedgesOfWire.Length() < 1 )
1741 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SupressHole() : no edge(s) for aWire", SALOME::BAD_PARAM);
1743 /* Retrieve face ancestors of all edges of 'aWire' but avoiding 'aFace' */
1745 TopTools_IndexedDataMapOfShapeListOfShape aMapEdgesFaces;
1746 TopTools_MapIteratorOfMapOfShape anIt ;
1747 TopTools_MapOfShape MFSuppress ;
1748 TopTools_MapOfShape MFSuppressTmp ;
1749 bool wireOnFace = false ;
1751 TopExp::MapShapesAndAncestors(aShape, TopAbs_EDGE, TopAbs_FACE, aMapEdgesFaces) ;
1752 for( int h=1; h<=SSedgesOfWire.Length(); h++ ) {
1754 TopoDS_Shape anEdgeOfWire = SSedgesOfWire.Value(h) ;
1755 int ind = aMapEdgesFaces.FindIndex(anEdgeOfWire) ;
1757 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SupressHole() : index of edge", SALOME::BAD_PARAM);
1759 TopTools_ListOfShape LF;
1760 LF = aMapEdgesFaces.FindFromIndex(ind) ; /* Contains all faces ancestors of an edge of wire */
1762 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SupressHole() : no face for an edge", SALOME::BAD_PARAM);
1764 /* Filter faces avoiding 'aFace' */
1765 TopTools_ListIteratorOfListOfShape it ;
1766 for( it.Initialize(LF); it.More(); it.Next() ) {
1767 TopoDS_Face F = TopoDS::Face( it.Value() ) ;
1768 if( !F.IsSame(aFace) ) {
1769 MFSuppressTmp.Add(F) ;
1779 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SupressHole() : wire not on selected face", SALOME::BAD_PARAM);
1782 /* Call routine to define faces to suppress and and optional endWire on endFace */
1783 TopTools_MapOfShape MSwireEndEdges ; /* will contain edges of final wire (hole) */
1784 SuppressHoleSubRoutine( aShape, aFace, SSedgesOfWire, aMapEdgesFaces, MFSuppressTmp, MFSuppress, withEndFace, endFace, MSwireEndEdges ) ;
1786 TopoDS_Wire endWire ;
1789 if( MSwireEndEdges.Extent() < 1 )
1790 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SupressHole() : MSwireEndEdges.Extent() < 1", SALOME::BAD_PARAM);
1792 if( !FindCompareWireHoleOnFace( endFace, MSwireEndEdges, endWire ) )
1793 THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SupressHole() : no endWire found", SALOME::BAD_PARAM);
1796 /* Build 'resTds' : a shape containing a compound of faces */
1797 TopoDS_Shape resTds;
1798 if( !withEndFace && !BuildShapeHoleNotTraversing( aShape, aFace, aWire, MFSuppress, resTds ) )
1799 THROW_SALOME_CORBA_EXCEPTION("Rebuild result shape has aborted", SALOME::BAD_PARAM);
1801 if( withEndFace && !BuildShapeHoleTraversing( aShape, aFace, aWire, MFSuppress, endFace, endWire, resTds ) )
1802 THROW_SALOME_CORBA_EXCEPTION("Rebuild result shape has aborted (end hole)", SALOME::BAD_PARAM);
1804 /* Reconstruction of final shape with 'resTds' : a compound of faces after hole suppressing */
1805 /* Actual limitation is : 'aShape' must not contain more than a solid or a shell ! */
1806 TopoDS_Shape finalShape ;
1808 TopExp_Explorer exp ;
1809 unsigned int nbSolid = 0 ;
1810 TopoDS_Solid aSolid ;
1811 for( exp.Init(aShape, TopAbs_SOLID); exp.More(); exp.Next() ) {
1812 aSolid = TopoDS::Solid( exp.Current() ) ;
1815 THROW_SALOME_CORBA_EXCEPTION("Limitation : main shape contents more than one solid", SALOME::BAD_PARAM);
1818 unsigned int nbShell = 0 ;
1819 TopoDS_Shell aShell ;
1820 for( exp.Init(aShape, TopAbs_SHELL); exp.More(); exp.Next() ) {
1821 aShell = TopoDS::Shell( exp.Current() ) ;
1824 THROW_SALOME_CORBA_EXCEPTION("Limitation : main shape contents more than one shell", SALOME::BAD_PARAM);
1827 /* No shells and no solids : can send a compound even for a single face, see GUI ! */
1828 if( nbShell == 0 ) {
1829 finalShape = resTds ;
1833 TopoDS_Shell shellResult ;
1834 if( nbShell == 1 ) {
1835 if ( !BuildShellWithFaceCompound( TopoDS::Compound(resTds), shellResult ) )
1836 THROW_SALOME_CORBA_EXCEPTION("Error after BuildShellWithFaceCompound()", SALOME::BAD_PARAM);
1837 finalShape = shellResult ;
1840 /* a solid with a shell */
1841 if( nbSolid == 1 && nbShell == 1) {
1842 BRepBuilderAPI_MakeSolid B;
1843 B.Add(shellResult) ;
1845 THROW_SALOME_CORBA_EXCEPTION("Error : !B.IsDone()", SALOME::BAD_PARAM);
1846 finalShape = B.Solid() ;
1849 result = CreateObject(finalShape);
1850 InsertInLabelOneArgument(aShape, shape, finalShape, result, myCurrentOCAFDoc) ;
1856 //================================================================================
1857 // function : BuildShellWithFaceCompound()
1858 // purpose : Build a shell with a compound of faces.
1859 //================================================================================
1860 bool GEOM_Gen_i::BuildShellWithFaceCompound( const TopoDS_Compound Comp,
1861 TopoDS_Shell& resultShell )
1863 resultShell.Nullify() ;
1864 BRepPrim_Builder B ;
1865 B.MakeShell(resultShell) ;
1866 TopExp_Explorer ex ;
1868 for( ex.Init( Comp, TopAbs_FACE); ex.More(); ex.Next() ) {
1869 TopoDS_Face F = TopoDS::Face( ex.Current() ) ;
1873 B.AddShellFace( resultShell, F ) ;
1875 B.CompleteShell(resultShell) ;
1876 if( resultShell.IsNull() ) {
1883 //================================================================================
1884 // function : FindCompareWireHoleOnFace()
1885 // purpose : Try to find a wire on 'aFace' which edges are same than those
1886 // : into 'MSwireEdges' map. 'aFoundWire' is defined and 'true' returned.
1887 // : 'MSwireEdges' represents generally an hole an 'aFace'.
1888 // : The outer wire of 'aFace' is avoided !
1889 //================================================================================
1890 bool GEOM_Gen_i::FindCompareWireHoleOnFace( const TopoDS_Face& F,
1891 const TopTools_MapOfShape& MSwireEdges,
1892 TopoDS_Wire& aFoundWire )
1894 aFoundWire.Nullify() ;
1899 /* Get the outer wire of aFace */
1900 TopoDS_Wire outerW = BRepTools::OuterWire(F) ;
1901 if( outerW.IsNull() || !IsValid(outerW) ) {
1905 int nbEdges = MSwireEdges.Extent() ;
1910 TopExp_Explorer exp1 ;
1911 TopExp_Explorer exp2 ;
1912 for ( exp1.Init(F, TopAbs_WIRE) ; exp1.More(); exp1.Next() ) {
1913 TopoDS_Wire W = TopoDS::Wire( exp1.Current() ) ;
1915 if( !W.IsSame(outerW) ) {
1916 for ( exp2.Init( W, TopAbs_EDGE) ; exp2.More(); exp2.Next() ) {
1917 TopoDS_Edge E = TopoDS::Edge( exp2.Current() ) ;
1918 if( MSwireEdges.Contains(E) ) {
1920 if( i == nbEdges ) {
1932 //================================================================================
1933 // function : BuildShapeHoleNotTraversing()
1934 // purpose : Define 'resultTds' a reconstruction of 'aShape' after modification
1935 // : on 'aFace' where 'aWire is removed' and suppression of faces 'MFSuppress'
1936 // : ( Used as a sub routine of SuppressHole() )
1937 //================================================================================
1938 bool GEOM_Gen_i::BuildShapeHoleNotTraversing( const TopoDS_Shape& aShape,
1939 const TopoDS_Face& aFace,
1940 const TopoDS_Wire& aWire,
1941 const TopTools_MapOfShape& MFSuppress,
1942 TopoDS_Shape& resultTds )
1943 throw (SALOME::SALOME_Exception)
1945 Unexpect aCatch(SALOME_SalomeException);
1947 TopExp_Explorer exp ;
1948 TopoDS_Face newFace ;
1950 resultTds.Nullify() ;
1951 TopoDS_Compound Comp ;
1952 B.MakeCompound (Comp);
1956 /* Explore all faces of 'aShape' to rebuild a compound */
1957 for ( exp.Init(aShape, TopAbs_FACE) ; exp.More(); exp.Next() ) {
1959 TopoDS_Face F1 = TopoDS::Face( exp.Current() );
1960 /* Rebuild face(s) not suppressed */
1961 if( !MFSuppress.Contains(F1) ) {
1963 if( F1.IsEqual( aFace ) ) {
1964 TopTools_MapOfShape MSwire ;
1966 if( !RebuildFaceRemovingHoles(aFace, MSwire, newFace) ) {
1969 B.Add( Comp, newFace ) ;
1972 /* For any other face not suppressed */
1978 catch(Standard_Failure) {
1979 THROW_SALOME_CORBA_EXCEPTION("in BuildShapeHoleNotTraversing() : Exception catched", SALOME::BAD_PARAM);
1987 //================================================================================
1988 // function : BuildShapeHoleTraversing()
1989 // purpose : Define 'resultTds' a reconstruction of 'aShape' after modifications.
1990 // : On 'aFace' 'aWire is removed'
1991 // : On 'endFace' 'endWire' is removed.
1992 // : Faces of 'MFSuppress' are removed.
1993 // : ( Used as a sub routine of SuppressHole() )
1994 //================================================================================
1995 bool GEOM_Gen_i::BuildShapeHoleTraversing( const TopoDS_Shape& aShape,
1996 const TopoDS_Face& aFace,
1997 const TopoDS_Wire& aWire,
1998 const TopTools_MapOfShape& MFSuppress,
1999 const TopoDS_Face& endFace,
2000 const TopoDS_Wire& endWire,
2001 TopoDS_Shape& resultTds )
2002 throw (SALOME::SALOME_Exception)
2004 Unexpect aCatch(SALOME_SalomeException);
2006 TopExp_Explorer exp ;
2007 TopoDS_Face newFace ;
2009 resultTds.Nullify() ;
2010 TopoDS_Compound Comp ;
2011 B.MakeCompound (Comp);
2013 /* Necessary to use general method */
2014 TopTools_MapOfShape MSwire1 ;
2015 MSwire1.Add(aWire) ;
2016 TopTools_MapOfShape MSwire2 ;
2017 MSwire2.Add(endWire) ;
2021 /* Explore all faces of 'aShape' to rebuild a compound */
2022 for ( exp.Init(aShape, TopAbs_FACE) ; exp.More(); exp.Next() ) {
2023 TopoDS_Face F1 = TopoDS::Face( exp.Current() );
2025 /* Rebuild face(s) not suppressed */
2026 if( !MFSuppress.Contains(F1) ) {
2028 /* Rebuild 'aFace' */
2029 if( F1.IsEqual( aFace ) && !F1.IsEqual( endFace ) ) {
2030 if( !RebuildFaceRemovingHoles(aFace, MSwire1, newFace) ) {
2033 B.Add(Comp, newFace) ;
2036 /* Rebuild 'endFace' */
2037 if( !F1.IsEqual( aFace ) && F1.IsEqual( endFace ) ) {
2038 if( !RebuildFaceRemovingHoles(endFace, MSwire2, newFace) ) {
2041 B.Add(Comp, newFace) ;
2044 /* Hole in the same face : aFace = endFace */
2045 if( F1.IsEqual( aFace ) && F1.IsEqual( endFace ) ) {
2047 if( !RebuildFaceRemovingHoles(aFace, MSwire1, newFace) || !RebuildFaceRemovingHoles(newFace, MSwire2, FF) ) {
2053 /* For any other face not suppressed */
2054 if( !F1.IsEqual(aFace) && !F1.IsEqual( endFace ) ) {
2061 catch(Standard_Failure) {
2062 THROW_SALOME_CORBA_EXCEPTION("in BuildShapeHoleTraversing() : Exception catched", SALOME::BAD_PARAM);
2068 //=======================================================================
2069 //function : SortShapes
2071 //=======================================================================
2073 static void SortShapes(TopTools_ListOfShape& SL)
2075 Standard_Integer MaxShapes = SL.Extent();
2076 TopTools_Array1OfShape aShapes (1,MaxShapes);
2077 TColStd_Array1OfInteger OrderInd(1,MaxShapes);
2078 TColStd_Array1OfReal MidXYZ (1,MaxShapes); //X,Y,Z;
2079 TColStd_Array1OfReal Length (1,MaxShapes); //X,Y,Z;
2081 // Computing of CentreOfMass
2082 Standard_Integer Index;
2085 TopTools_ListIteratorOfListOfShape it(SL);
2086 for (Index=1; it.More(); Index++)
2088 TopoDS_Shape S = it.Value();
2089 SL.Remove( it ); // == it.Next()
2091 OrderInd.SetValue (Index, Index);
2092 if (S.ShapeType() == TopAbs_VERTEX)
2094 GPoint = BRep_Tool::Pnt( TopoDS::Vertex( S ));
2095 Length.SetValue( Index, (Standard_Real) S.Orientation());
2099 BRepGProp::LinearProperties (S, GPr);
2100 GPoint = GPr.CentreOfMass();
2101 Length.SetValue( Index, GPr.Mass() );
2103 MidXYZ.SetValue(Index,
2104 GPoint.X()*999 + GPoint.Y()*99 + GPoint.Z()*0.9);
2107 Standard_Integer aTemp;
2108 Standard_Boolean exchange, Sort = Standard_True;
2111 Sort = Standard_False;
2112 for (Index=1; Index < MaxShapes; Index++)
2114 if (MidXYZ(OrderInd(Index)) > MidXYZ(OrderInd(Index+1)))
2115 exchange = Standard_True;
2116 else if (MidXYZ(OrderInd(Index)) == MidXYZ(OrderInd(Index+1)) &&
2117 Length(OrderInd(Index)) > Length(OrderInd(Index+1)) )
2118 exchange = Standard_True;
2120 exchange = Standard_False;
2123 aTemp = OrderInd(Index);
2124 OrderInd(Index) = OrderInd(Index+1);
2125 OrderInd(Index+1) = aTemp;
2126 Sort = Standard_True;
2130 for (Index=1; Index <= MaxShapes; Index++)
2131 SL.Append( aShapes( OrderInd(Index) ));
2134 //================================================================================
2135 // function : SubShape()
2136 // purpose : Method for GUI or TUI
2137 //================================================================================
2139 GEOM::GEOM_Shape_ptr GEOM_Gen_i::SubShape(GEOM::GEOM_Shape_ptr shape,
2140 CORBA::Short ShapeType,
2141 const GEOM::GEOM_Shape::ListOfSubShapeID& ListOfID)
2142 throw (SALOME::SALOME_Exception)
2144 Unexpect aCatch(SALOME_SalomeException);
2145 return SubShapesOne(shape, (TopAbs_ShapeEnum) ShapeType, ListOfID);
2148 //================================================================================
2149 // function : SubShapeSorted()
2150 // purpose : Method for GUI or TUI
2151 //================================================================================
2153 GEOM::GEOM_Shape_ptr GEOM_Gen_i::SubShapeSorted(GEOM::GEOM_Shape_ptr shape,
2154 CORBA::Short ShapeType,
2155 const GEOM::GEOM_Shape::ListOfSubShapeID& ListOfID)
2156 throw (SALOME::SALOME_Exception)
2158 Unexpect aCatch(SALOME_SalomeException);
2159 return SubShapesOne(shape, (TopAbs_ShapeEnum) ShapeType, ListOfID, Standard_True);
2162 //================================================================================
2163 // function : SubShapesOne()
2165 //================================================================================
2167 GEOM::GEOM_Shape_ptr GEOM_Gen_i::SubShapesOne( GEOM::GEOM_Shape_ptr shape,
2168 const TopAbs_ShapeEnum ShapeType,
2169 const GEOM::GEOM_Shape::ListOfSubShapeID& ListOfID,
2170 const Standard_Boolean Sort)
2171 throw (SALOME::SALOME_Exception)
2173 Unexpect aCatch(SALOME_SalomeException);
2174 GEOM::GEOM_Shape_var result;
2175 TopoDS_Shape mainShape;
2176 TopoDS_Shape mainTopo = GetTopoShape(shape);
2180 if ( shape->IsMainShape() ) {
2181 mainShape = GetTopoShape(shape);
2184 shape = GetIORFromString( shape->MainName() );
2187 if(ListOfID.length() < 1) {
2188 THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::SubShape() : bad list of shapes",
2192 /* Create a sequence of all sub shapes */
2193 TopTools_ListOfShape listShape;
2195 TopTools_MapOfShape mapShape;
2196 TopExp_Explorer exp ;
2197 for ( exp.Init( mainShape, TopAbs_ShapeEnum(ShapeType)) ; exp.More(); exp.Next() ) {
2198 if ( mapShape.Add( exp.Current() ) )
2199 listShape.Append ( exp.Current() );
2202 if (listShape.IsEmpty()) return shape;
2205 SortShapes(listShape);
2207 TopTools_SequenceOfShape SS;
2208 TopTools_ListIteratorOfListOfShape it (listShape);
2209 for (; it.More(); it.Next())
2210 SS.Append( it.Value() );
2212 TopoDS_Shape SubShape;
2213 if (ListOfID.length() == 1)
2214 SubShape = SS.Value(ListOfID[0]);
2218 TopoDS_Compound Comp;
2219 B.MakeCompound (Comp);
2221 for ( ind = 0; ind < ListOfID.length(); ind++ )
2222 B.Add( Comp, SS.Value(ListOfID[ind]) );
2226 if ( !TNaming_Tool::HasLabel( myCurrentOCAFDoc->Main(), SubShape ) )
2227 result = CreateSubObject( SubShape, shape, ListOfID);
2230 TDF_Label Lab = TNaming_Tool::Label( myCurrentOCAFDoc->Main(), SubShape, TransDef );
2231 Handle(TDataStd_Name) Att;
2232 if ( Lab.FindAttribute( TDataStd_Name::GetID(), Att ) ) {
2233 TCollection_AsciiString nameIOR( Att->Get() );
2234 result = GEOM::GEOM_Shape::_narrow(_orb->string_to_object( nameIOR.ToCString() ));
2235 if ( strcmp( result->MainName(), _orb->object_to_string(shape) ) != 0 ) {
2236 result = CreateSubObject( SubShape, shape, ListOfID);
2245 //================================================================================
2246 // function : SubShapeAll()
2247 // purpose : Explode a shape in all sub shapes with a type (Method for TUI or GUI)
2248 //================================================================================
2250 GEOM::GEOM_Gen::ListOfGeomShapes* GEOM_Gen_i::SubShapeAll(GEOM::GEOM_Shape_ptr shape,
2251 CORBA::Short ShapeType)
2252 throw (SALOME::SALOME_Exception)
2254 Unexpect aCatch(SALOME_SalomeException);
2255 return SubShapesAll(shape, (TopAbs_ShapeEnum) ShapeType);
2258 //================================================================================
2259 // function : SubShapeAllSorted()
2260 // purpose : Explode a shape in all sub shapes with a type (Method for TUI or GUI)
2261 //================================================================================
2263 GEOM::GEOM_Gen::ListOfGeomShapes* GEOM_Gen_i::SubShapeAllSorted(GEOM::GEOM_Shape_ptr shape,
2264 CORBA::Short ShapeType)
2265 throw (SALOME::SALOME_Exception)
2267 Unexpect aCatch(SALOME_SalomeException);
2268 return SubShapesAll(shape, (TopAbs_ShapeEnum) ShapeType, Standard_True);
2271 //================================================================================
2272 // function : SubShapeAllSorted()
2274 //================================================================================
2276 GEOM::GEOM_Gen::ListOfGeomShapes* GEOM_Gen_i::SubShapesAll(GEOM::GEOM_Shape_ptr shape,
2277 const TopAbs_ShapeEnum ShapeType,
2278 const Standard_Boolean Sort)
2279 throw (SALOME::SALOME_Exception)
2281 /* List of sub shapes returned */
2282 Unexpect aCatch(SALOME_SalomeException);
2283 GEOM::GEOM_Gen::ListOfGeomShapes_var listOfGeomShapes = new GEOM::GEOM_Gen::ListOfGeomShapes;
2284 listOfGeomShapes->length(0) ;
2286 TopoDS_Shape mainTopo = GetTopoShape(shape);
2287 TopoDS_Shape mainShape;
2290 if ( shape->IsMainShape() ) {
2291 mainShape = GetTopoShape(shape);
2294 shape = GetIORFromString( shape->MainName() );
2297 if( mainTopo.IsNull() ) {
2298 THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::SubShapeAll() : null main shape",
2302 /* List/index : field set for any sub shape */
2303 GEOM::GEOM_Shape::ListOfSubShapeID_var ListOfID = new GEOM::GEOM_Shape::ListOfSubShapeID;
2304 ListOfID->length(1) ;
2306 // retrieve all subshapes
2307 TopTools_MapOfShape mapShape;
2308 TopTools_ListOfShape listShape;
2310 if (mainTopo.ShapeType()==TopAbs_COMPOUND && ShapeType==TopAbs_SHAPE)
2312 TopoDS_Iterator It(mainTopo,Standard_True,Standard_True );
2313 for ( ; It.More(); It.Next() )
2314 if (mapShape.Add( It.Value() ))
2315 listShape.Append( It.Value() );
2319 TopExp_Explorer exp ( mainTopo, ShapeType);
2320 for ( ; exp.More(); exp.Next() )
2321 if (mapShape.Add( exp.Current() ))
2322 listShape.Append( exp.Current() );
2326 SortShapes(listShape);
2328 /* Create all sub shapes */
2330 GEOM::GEOM_Shape_var result;
2332 TopTools_ListIteratorOfListOfShape itSub (listShape);
2333 for (index = 1; itSub.More(); itSub.Next(), ++index)
2335 const TopoDS_Shape& SubShape = itSub.Value();
2336 // check if SubShape is already in OCAFDS and ...
2337 if ( TNaming_Tool::HasLabel( myCurrentOCAFDoc->Main(), SubShape ) )
2340 TDF_Label Lab = TNaming_Tool::Label( myCurrentOCAFDoc->Main(), SubShape, TransDef );
2341 Handle(TDataStd_Name) Att;
2342 if ( Lab.FindAttribute( TDataStd_Name::GetID(), Att ) )
2344 TCollection_AsciiString nameIOR( Att->Get() );
2345 result = GEOM::GEOM_Shape::_narrow(_orb->string_to_object( nameIOR.ToCString() ));
2346 // ... it is subshape of <shape>
2347 if ( strcmp( result->MainName(), _orb->object_to_string(shape) ) == 0 )
2349 listOfGeomShapes->length(index) ;
2350 listOfGeomShapes[index-1] = result ;
2357 // ListOfID[0] = index;
2359 ListOfID[0] = GetIndexTopology( SubShape, mainShape ) ;
2360 result = CreateSubObject( SubShape, shape, ListOfID);
2361 /* Add each sub shape in the list returned */
2362 listOfGeomShapes->length(index) ;
2363 listOfGeomShapes[index-1] = result ;
2366 return listOfGeomShapes._retn() ;
2369 //=================================================================================
2370 // function : MakeBoolean()
2371 // purpose : Boolean operation according to the type 'operation'
2372 //=================================================================================
2373 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeBoolean(GEOM::GEOM_Shape_ptr shape1,
2374 GEOM::GEOM_Shape_ptr shape2,
2375 CORBA::Long operation)
2376 throw (SALOME::SALOME_Exception)
2378 Unexpect aCatch(SALOME_SalomeException);
2379 GEOM::GEOM_Shape_var result;
2380 TopoDS_Shape shape ;
2381 TopoDS_Shape aShape1 ;
2382 TopoDS_Shape aShape2 ;
2385 aShape1 = GetTopoShape(shape1) ;
2386 aShape2 = GetTopoShape(shape2) ;
2388 if( aShape1.IsNull() || aShape2.IsNull() ) {
2389 THROW_SALOME_CORBA_EXCEPTION("Boolean aborted : argument shape is null", SALOME::BAD_PARAM);
2394 case 1 : /* Common */
2395 shape = BRepAlgoAPI_Common(aShape1, aShape2).Shape();
2398 shape = BRepAlgoAPI_Cut(aShape1, aShape2).Shape();
2401 shape = BRepAlgoAPI_Fuse(aShape1, aShape2).Shape();
2403 case 4 : /* Section */
2404 shape = BRepAlgoAPI_Section(aShape1, aShape2).Shape();
2407 MESSAGE("Boolean operation not known : " << operation ) ;
2411 catch(Standard_Failure) {
2412 THROW_SALOME_CORBA_EXCEPTION("Exception catched in boolean operation", SALOME::BAD_PARAM);
2415 /* We test the validity of resulting shape */
2416 if( !IsValid(shape) ) {
2417 THROW_SALOME_CORBA_EXCEPTION("Boolean aborted : non valid shape result", SALOME::BAD_PARAM);
2420 result = CreateObject(shape) ;
2422 /* put shape and name into geom/OCAF doc */
2423 GEOMDS_Commands GC(myCurrentOCAFDoc->Main());
2424 /* add attributs 'shape' and' name_io'r in a new label */
2425 TDF_Label Lab = GC.Generated(aShape1, shape, result->Name() );
2426 TCollection_AsciiString entry;
2427 TDF_Tool::Entry(Lab, entry);
2428 result->ShapeId( entry.ToCString() ) ;
2430 /* Create a new label */
2431 TDF_Label NewLab = Lab.NewChild();
2432 TCollection_ExtendedString Value("Arguments");
2433 TDataStd_Name::Set(NewLab,Value);
2435 TDF_Label NewLab1 = NewLab.NewChild();
2437 TDF_Tool::Label(myCurrentOCAFDoc->GetData(), shape1->ShapeId(), RefLab);
2438 TDF_Reference::Set(NewLab1, RefLab);
2440 TDF_Label NewLab2 = NewLab.NewChild();
2441 TDF_Tool::Label(myCurrentOCAFDoc->GetData(), shape2->ShapeId(), RefLab);
2442 TDF_Reference::Set(NewLab2, RefLab);
2447 //=================================================================================
2448 // function : MakeFuse()
2449 // purpose : Special purpose !
2450 //=================================================================================
2451 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeFuse(GEOM::GEOM_Shape_ptr shape1,
2452 GEOM::GEOM_Shape_ptr shape2)
2453 throw (SALOME::SALOME_Exception)
2455 Unexpect aCatch(SALOME_SalomeException);
2456 GEOM::GEOM_Shape_var result;
2457 TopoDS_Shape aShape1 = GetTopoShape(shape1) ;
2458 TopoDS_Shape aShape2 = GetTopoShape(shape2) ;
2459 if( aShape1.IsNull() || aShape2.IsNull() ) {
2460 THROW_SALOME_CORBA_EXCEPTION("Fuse aborted : shape in argument is null", SALOME::BAD_PARAM);
2464 shape = BRepAlgoAPI_Fuse(aShape1, aShape2).Shape();
2466 catch(Standard_Failure) {
2467 THROW_SALOME_CORBA_EXCEPTION("Exception catched in Fuse operation", SALOME::BAD_PARAM);
2470 /* We test the validity of resulting shape */
2471 if( !IsValid(shape) ) {
2472 THROW_SALOME_CORBA_EXCEPTION("Fuse aborted : non valid shape result", SALOME::BAD_PARAM);
2475 result = CreateObject(shape) ;
2477 /* put shape and name into geom/OCAF doc */
2478 GEOMDS_Commands GC(myCurrentOCAFDoc->Main());
2479 /* add attributs 'shape' and' name_io'r in a new label */
2480 TDF_Label Lab = GC.Generated(aShape1, shape, result->Name() );
2481 TCollection_AsciiString entry;
2482 TDF_Tool::Entry(Lab, entry);
2483 result->ShapeId( entry.ToCString() ) ;
2485 /* Create a new label */
2486 TDF_Label NewLab = Lab.NewChild();
2487 TCollection_ExtendedString Value("Arguments");
2488 TDataStd_Name::Set(NewLab,Value);
2490 TDF_Label NewLab1 = NewLab.NewChild();
2492 TDF_Tool::Label(myCurrentOCAFDoc->GetData(), shape1->ShapeId(), RefLab);
2493 TDF_Reference::Set(NewLab1, RefLab);
2495 TDF_Label NewLab2 = NewLab.NewChild();
2496 TDF_Tool::Label(myCurrentOCAFDoc->GetData(), shape2->ShapeId(), RefLab);
2497 TDF_Reference::Set(NewLab2, RefLab);
2502 //================================================================================
2503 // function : MakeAxisStruct()
2504 // purpose : Create a structure GEOM::AxisStruct (see IDL file)
2505 //================================================================================
2506 GEOM::AxisStruct GEOM_Gen_i::MakeAxisStruct(CORBA::Double x,
2513 GEOM::AxisStruct A ;
2514 A.x = x ; A.y = y ; A.z = z ;
2515 A.vx = vx ; A.vy = vy ; A.vz = vz ;
2520 //================================================================================
2521 // function : MakePointStruct()
2522 // purpose : Create a structure GEOM::PointStruct (see IDL file)
2523 //================================================================================
2524 GEOM::PointStruct GEOM_Gen_i::MakePointStruct(CORBA::Double x,
2528 beginService( "GEOM_Gen_i::MakePointStruct" );
2529 GEOM::PointStruct p ;
2530 p.x = x ; p.y = y ; p.z = z ;
2531 endService( "GEOM_Gen_i::MakePointStruct" );
2535 //================================================================================
2536 // function : MakeDirection()
2537 // purpose : Create a structure GEOM::DirStruct (see IDL file)
2538 //================================================================================
2539 GEOM::DirStruct GEOM_Gen_i::MakeDirection(const GEOM::PointStruct& p)
2542 d.PS.x = p.x ; d.PS.y = p.y ; d.PS.z = p.z ;
2546 //=================================================================================
2547 // function : MakeBox()
2548 // purpose : Create a box topology.
2549 //=================================================================================
2550 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeBox(CORBA::Double x1,
2556 throw (SALOME::SALOME_Exception)
2558 Unexpect aCatch(SALOME_SalomeException);
2559 gp_Pnt P1(x1,y1,z1);
2560 gp_Pnt P2(x2,y2,z2);
2561 GEOM::GEOM_Shape_var result ;
2564 tds = BRepPrimAPI_MakeBox(P1,P2).Shape();
2566 catch(Standard_Failure) {
2567 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeBox", SALOME::BAD_PARAM);
2571 THROW_SALOME_CORBA_EXCEPTION("Make Box aborted : null shape", SALOME::BAD_PARAM);
2574 result = CreateObject(tds);
2575 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
2576 result->ShapeId(entry) ;
2581 //================================================================================
2582 // function : MakeCylinder
2583 // purpose : Create a cylinder topology
2584 //================================================================================
2585 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeCylinder(const GEOM::PointStruct& pstruct,
2586 const GEOM::DirStruct& dstruct,
2587 CORBA::Double radius,
2588 CORBA::Double height)
2589 throw (SALOME::SALOME_Exception)
2591 Unexpect aCatch(SALOME_SalomeException);
2592 GEOM::GEOM_Shape_var result;
2594 gp_Pnt p(pstruct.x, pstruct.y, pstruct.z) ;
2595 gp_Dir d(dstruct.PS.x, dstruct.PS.y, dstruct.PS.z) ;
2599 tds = BRepPrimAPI_MakeCylinder(axis, radius, height).Shape();
2601 catch(Standard_Failure) {
2602 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeCylinder", SALOME::BAD_PARAM);
2606 THROW_SALOME_CORBA_EXCEPTION("Make Cylinder aborted", SALOME::BAD_PARAM);
2608 result = CreateObject(tds);
2609 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
2610 result->ShapeId(entry);
2614 //================================================================================
2615 // function : MakeSphere()
2616 // purpose : Make a sphere topology
2617 //================================================================================
2618 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeSphere(CORBA::Double x1,
2621 CORBA::Double radius)
2622 throw (SALOME::SALOME_Exception)
2624 Unexpect aCatch(SALOME_SalomeException);
2625 GEOM::GEOM_Shape_var result ;
2628 tds = BRepPrimAPI_MakeSphere(gp_Pnt(x1,y1,z1), radius).Shape();
2630 catch(Standard_Failure) {
2631 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeSphere", SALOME::BAD_PARAM);
2635 THROW_SALOME_CORBA_EXCEPTION("Make Sphere aborted", SALOME::BAD_PARAM);
2637 result = CreateObject(tds) ;
2638 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
2639 result->ShapeId(entry);
2643 //================================================================================
2644 // function : MakeTorus()
2645 // purpose : Create a torus topology
2646 //================================================================================
2647 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeTorus( const GEOM::PointStruct& pstruct,
2648 const GEOM::DirStruct& dstruct,
2649 CORBA::Double major_radius,
2650 CORBA::Double minor_radius )
2651 throw (SALOME::SALOME_Exception)
2653 Unexpect aCatch(SALOME_SalomeException);
2654 GEOM::GEOM_Shape_var result;
2656 gp_Pnt p(pstruct.x, pstruct.y, pstruct.z) ;
2657 gp_Dir d(dstruct.PS.x, dstruct.PS.y, dstruct.PS.z) ;
2661 tds = BRepPrimAPI_MakeTorus(axis, major_radius, minor_radius).Shape();
2663 catch(Standard_Failure) {
2664 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeTorus", SALOME::BAD_PARAM);
2668 THROW_SALOME_CORBA_EXCEPTION("Make torus aborted", SALOME::BAD_PARAM);
2670 result = CreateObject(tds);
2671 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
2672 result->ShapeId(entry);
2677 //================================================================================
2678 // function : MakeCone()
2679 // purpose : Create a cone topology
2680 //================================================================================
2681 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeCone(const GEOM::PointStruct& pstruct,
2682 const GEOM::DirStruct& dstruct,
2683 CORBA::Double radius1,
2684 CORBA::Double radius2,
2685 CORBA::Double height)
2686 throw (SALOME::SALOME_Exception)
2688 Unexpect aCatch(SALOME_SalomeException);
2689 GEOM::GEOM_Shape_var result;
2691 gp_Pnt p(pstruct.x, pstruct.y, pstruct.z) ;
2692 gp_Dir d(dstruct.PS.x, dstruct.PS.y, dstruct.PS.z) ;
2696 /* Cone doesn't work if same radius */
2697 if( fabs(radius1-radius2) <= Precision::Confusion() ) {
2698 tds = BRepPrimAPI_MakeCylinder(axis, (radius1+radius2)/2.0, height).Shape();
2701 tds = BRepPrimAPI_MakeCone(axis, radius1, radius2, height).Shape();
2704 catch(Standard_Failure) {
2705 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeCone", SALOME::BAD_PARAM);
2709 THROW_SALOME_CORBA_EXCEPTION("Make Cone aborted", SALOME::BAD_PARAM);
2711 result = CreateObject(tds);
2712 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
2713 result->ShapeId(entry);
2719 //==================================================================================
2720 // function : ImportIGES()
2721 // purpose : Import shape from an IGES (IGS) file
2722 // : LPN modified 7 mai 2002
2723 //==================================================================================
2724 GEOM::GEOM_Shape_ptr GEOM_Gen_i::ImportIGES(const char* filename)
2725 throw (SALOME::SALOME_Exception)
2727 Unexpect aCatch(SALOME_SalomeException);
2728 GEOM::GEOM_Shape_var result ;
2729 //VRV: OCC 4.0 migration
2730 IGESControl_Reader aReader;
2731 //VRV: OCC 4.0 migration
2733 IFSelect_ReturnStatus stat = aReader.ReadFile((char*)filename);
2734 if ( stat != IFSelect_RetDone ) {
2735 THROW_SALOME_CORBA_EXCEPTION("Error in reading import file", SALOME::BAD_PARAM); }
2737 MESSAGE("ImportIGES : all Geometry Transfer" << endl ) ;
2738 #if OCC_VERSION_MAJOR >= 5
2739 aReader.ClearShapes();
2740 aReader.TransferRoots();
2743 aReader.TransferRoots(false);
2745 MESSAGE("ImportIGES : count of shapes produced = " << aReader.NbShapes() << endl );
2746 TopoDS_Shape shape = aReader.OneShape();
2748 if ( !shape.IsNull() ) {
2749 /* Final CORBA object creation */
2750 result = CreateObject(shape) ;
2751 const char *entry = InsertInLabel( shape, result->Name(), myCurrentOCAFDoc ) ;
2752 result->ShapeId(entry);
2756 catch(Standard_Failure) {
2757 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ImportIGES()", SALOME::BAD_PARAM);
2760 THROW_SALOME_CORBA_EXCEPTION("Import IGES aborted : internal error", SALOME::BAD_PARAM);
2765 //==================================================================================
2766 // function : ImportSTEP()
2767 // purpose : Import shape from an STEP (stp) file
2768 // : 'result' is a compound of shapes if file contains more entities.
2769 //==================================================================================
2770 GEOM::GEOM_Shape_ptr GEOM_Gen_i::ImportSTEP(const char* filename)
2771 throw (SALOME::SALOME_Exception)
2773 Unexpect aCatch(SALOME_SalomeException);
2774 GEOM::GEOM_Shape_var result ;
2775 //VRV: OCC 4.0 migration
2776 STEPControl_Reader aReader;
2777 //VRV: OCC 4.0 migration
2779 TopoDS_Compound compound;
2781 B.MakeCompound( compound );
2784 IFSelect_ReturnStatus status = aReader.ReadFile((char*)filename);
2786 if (status == IFSelect_RetDone) {
2787 Standard_Boolean failsonly = Standard_False ;
2788 aReader.PrintCheckLoad (failsonly, IFSelect_ItemsByEntity);
2789 /* Root transfers */
2790 Standard_Integer nbr = aReader.NbRootsForTransfer();
2791 aReader.PrintCheckTransfer (failsonly, IFSelect_ItemsByEntity);
2793 for ( Standard_Integer n=1; n <= nbr; n++) {
2795 Standard_Boolean ok = aReader.TransferRoot(n);
2796 /* Collecting resulting entities */
2797 Standard_Integer nbs = aReader.NbShapes();
2799 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ImportStep", SALOME::BAD_PARAM) ;
2801 for ( Standard_Integer i=1; i<=nbs; i++ ) {
2802 TopoDS_Shape aShape = aReader.Shape(i);
2803 if ( aShape.IsNull() )
2804 THROW_SALOME_CORBA_EXCEPTION("Null shape in GEOM_Gen_i::ImportStep", SALOME::BAD_PARAM) ;
2806 /* For a single entity */
2807 if(nbr == 1 && nbs == 1) {
2808 result = CreateObject(aShape) ;
2809 const char *entry = InsertInLabel(aShape, result->Name(), myCurrentOCAFDoc) ;
2810 result->ShapeId(entry);
2814 B.Add( compound, aShape ) ;
2819 TopoDS_Shape tds = compound ;
2820 result = CreateObject(tds) ;
2821 if( CORBA::is_nil(result) )
2822 THROW_SALOME_CORBA_EXCEPTION("Translation aborted : null result", SALOME::BAD_PARAM);
2823 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
2824 result->ShapeId(entry);
2829 catch(Standard_Failure) {
2830 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ImportStep", SALOME::BAD_PARAM);
2836 //==================================================================================
2837 // function : Partition()
2839 //==================================================================================
2841 GEOM::GEOM_Shape_ptr
2842 GEOM_Gen_i::Partition(const GEOM::GEOM_Gen::ListOfIOR& ListShapes,
2843 const GEOM::GEOM_Gen::ListOfIOR& ListTools,
2844 const GEOM::GEOM_Gen::ListOfIOR& ListKeepInside,
2845 const GEOM::GEOM_Gen::ListOfIOR& ListRemoveInside,
2846 const CORBA::Short Limit)
2847 throw (SALOME::SALOME_Exception)
2849 Unexpect aCatch(SALOME_SalomeException);
2850 GEOM::GEOM_Shape_var aResult;
2852 //MESSAGE ("In Partition");
2855 unsigned int ind, nbshapes = 0;
2856 nbshapes += ListShapes.length() + ListTools.length();
2857 nbshapes += ListKeepInside.length() + ListRemoveInside.length();
2859 Partition_Spliter PS;
2860 TopTools_MapOfShape ShapesMap(nbshapes), ToolsMap(nbshapes);
2862 // add object shapes that are in ListShapes;
2863 for ( ind = 0; ind < ListShapes.length(); ind++) {
2865 GEOM::GEOM_Shape_var aShape = GetIORFromString( ListShapes[ind] );
2866 TopoDS_Shape Shape = GetTopoShape(aShape);
2867 if(Shape.IsNull() ) {
2868 //MESSAGE ( "In Partition a shape is null" );
2869 THROW_SALOME_CORBA_EXCEPTION("In Partition a shape is null", SALOME::BAD_PARAM);
2871 if ( ShapesMap.Add( Shape ))
2875 // add tool shapes that are in ListTools and not in ListShapes;
2876 for (ind = 0; ind < ListTools.length(); ind++) {
2878 GEOM::GEOM_Shape_var aShape = GetIORFromString( ListTools[ind] );
2879 TopoDS_Shape Shape = GetTopoShape(aShape);
2880 if(Shape.IsNull() ) {
2881 THROW_SALOME_CORBA_EXCEPTION("In Partition a shape is null", SALOME::BAD_PARAM);
2883 if ( !ShapesMap.Contains( Shape ) && ToolsMap.Add( Shape ))
2887 // add shapes that are in ListKeepInside, as object shapes;
2888 for (ind = 0; ind < ListKeepInside.length(); ind++) {
2890 GEOM::GEOM_Shape_var aShape = GetIORFromString( ListKeepInside[ind] );
2891 TopoDS_Shape Shape = GetTopoShape(aShape);
2892 if(Shape.IsNull() ) {
2893 //MESSAGE ( "In Partition a Keep Inside shape is null" );
2894 THROW_SALOME_CORBA_EXCEPTION("In Partition a shape is null", SALOME::BAD_PARAM);
2896 if (!ToolsMap.Contains( Shape ) &&
2897 ShapesMap.Add( Shape ))
2901 // add shapes that are in ListRemoveInside, as object shapes;
2902 for (ind = 0; ind < ListRemoveInside.length(); ind++) {
2904 GEOM::GEOM_Shape_var aShape = GetIORFromString( ListRemoveInside[ind] );
2905 TopoDS_Shape Shape = GetTopoShape(aShape);
2906 if(Shape.IsNull() ) {
2907 //MESSAGE ( "In Partition a Remove Inside shape is null" );
2908 THROW_SALOME_CORBA_EXCEPTION("In Partition a shape is null", SALOME::BAD_PARAM);
2910 if (!ToolsMap.Contains( Shape ) &&
2911 ShapesMap.Add( Shape ) )
2915 //MESSAGE ( "Partition::Compute() " );
2916 PS.Compute ((TopAbs_ShapeEnum) Limit);
2917 //MESSAGE ( "Partition::Compute() - END" );
2919 // suppress result outside of shapes in KInsideMap
2920 for (ind = 0; ind < ListKeepInside.length(); ind++) {
2921 GEOM::GEOM_Shape_var aShape = GetIORFromString( ListKeepInside[ind] );
2922 TopoDS_Shape Shape = GetTopoShape(aShape);
2923 PS.KeepShapesInside( Shape );
2926 // suppress result inside of shapes in RInsideMap
2927 for (ind = 0; ind < ListRemoveInside.length(); ind++) {
2929 GEOM::GEOM_Shape_var aShape = GetIORFromString( ListRemoveInside[ind] );
2930 TopoDS_Shape Shape = GetTopoShape(aShape);
2931 PS.RemoveShapesInside( Shape );
2936 if( !IsValid(tds) ) {
2937 //MESSAGE ( "In Partition: non valid shape result" );
2938 THROW_SALOME_CORBA_EXCEPTION("Partition aborted : non valid shape result", SALOME::BAD_PARAM);
2941 catch (Standard_Failure) {
2942 //MESSAGE ( "In Partition: Exception catched in GEOM_Gen_i::Partition()" );
2943 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::Partition", SALOME::BAD_PARAM);
2946 aResult = CreateObject(tds) ;
2948 /* add attributs S and mystr in a new label */
2949 GEOMDS_Commands GC(myCurrentOCAFDoc->Main());
2950 TDF_Label Lab = GC.AddShape(tds, aResult->Name() );
2951 TCollection_AsciiString entry;
2952 TDF_Tool::Entry(Lab,entry);
2953 aResult->ShapeId( entry.ToCString() ) ;
2957 /* Create a new label */
2958 TDF_Label NewLab = Lab.NewChild();
2959 TCollection_ExtendedString Value("Arguments");
2960 TDataStd_Name::Set(NewLab,Value);
2963 for (unsigned int ind = 0; ind < ListShapes.length(); ind++) {
2964 TDF_Label NewLab1 = NewLab.NewChild();
2966 GEOM::GEOM_Shape_var aShape = GetIORFromString( ListShapes[ind] );
2967 Standard_CString anEntry = aShape->ShapeId();
2969 TDF_Tool::Label(myCurrentOCAFDoc->GetData(), anEntry, RefLab);
2970 TDF_Reference::Set(NewLab1,RefLab);
2973 for (unsigned int ind = 0; ind < ListTools.length(); ind++) {
2974 TDF_Label NewLab1 = NewLab.NewChild();
2976 GEOM::GEOM_Shape_var aShape = GetIORFromString( ListTools[ind] );
2977 Standard_CString anEntry = aShape->ShapeId();
2979 TDF_Tool::Label(myCurrentOCAFDoc->GetData(), anEntry, RefLab);
2980 TDF_Reference::Set(NewLab1,RefLab);
2983 for (unsigned int ind = 0; ind < ListKeepInside.length(); ind++) {
2984 TDF_Label NewLab1 = NewLab.NewChild();
2986 GEOM::GEOM_Shape_var aShape = GetIORFromString( ListKeepInside[ind] );
2987 Standard_CString anEntry = aShape->ShapeId();
2989 TDF_Tool::Label(myCurrentOCAFDoc->GetData(), anEntry, RefLab);
2990 TDF_Reference::Set(NewLab1,RefLab);
2993 for (unsigned int ind = 0; ind < ListRemoveInside.length(); ind++) {
2994 TDF_Label NewLab1 = NewLab.NewChild();
2996 GEOM::GEOM_Shape_var aShape = GetIORFromString( ListRemoveInside[ind] );
2997 Standard_CString anEntry = aShape->ShapeId();
2999 TDF_Tool::Label(myCurrentOCAFDoc->GetData(), anEntry, RefLab);
3000 TDF_Reference::Set(NewLab1,RefLab);
3008 //==================================================================================
3009 // function : MakeFilling()
3010 // purpose : Create a surface from section curves filling
3011 //==================================================================================
3012 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeFilling(GEOM::GEOM_Shape_ptr myShape,
3013 CORBA::Short mindeg,
3014 CORBA::Short maxdeg,
3015 CORBA::Double tol3d,
3016 CORBA::Double tol2d,
3017 CORBA::Short nbiter)
3018 throw (SALOME::SALOME_Exception)
3020 Unexpect aCatch(SALOME_SalomeException);
3021 GEOM::GEOM_Shape_var result ;
3023 TopoDS_Shape aShape = GetTopoShape(myShape) ;
3024 if( aShape.IsNull() || aShape.ShapeType() != TopAbs_COMPOUND ) {
3025 THROW_SALOME_CORBA_EXCEPTION("MakeFilling aborted : null shape or not a compound", SALOME::BAD_PARAM);
3029 /* we verify the contents of the shape */
3030 TopExp_Explorer Ex ;
3031 TopoDS_Shape Scurrent ;
3032 Standard_Real First, Last ;
3033 Handle(Geom_Curve) C ;
3034 GeomFill_SectionGenerator Section ;
3036 Standard_Integer i = 0 ;
3037 for(Ex.Init(aShape, TopAbs_EDGE); Ex.More(); Ex.Next()) {
3038 Scurrent = Ex.Current() ;
3039 if( Scurrent.IsNull() || Scurrent.ShapeType() != TopAbs_EDGE) {
3040 THROW_SALOME_CORBA_EXCEPTION("Initial shape doesn't contain only edges !", SALOME::BAD_PARAM);
3042 C = BRep_Tool::Curve(TopoDS::Edge(Scurrent), First, Last);
3043 if (C.IsNull()) continue;
3044 C = new Geom_TrimmedCurve(C, First, Last);
3045 Section.AddCurve(C) ;
3049 /* a 'tolerance' is used to compare 2 knots : see GeomFill_Generator.cdl */
3050 /* We set 'tolerance' = tol3d */
3051 // Section.Perform( tol3d ) ; NRI */
3052 Section.Perform( Precision::Confusion() ) ;
3053 Handle(GeomFill_Line) Line = new GeomFill_Line(i) ;
3055 GeomFill_AppSurf App(mindeg, maxdeg, tol3d, tol2d, nbiter) ; /* user parameters */
3056 App.Perform(Line, Section) ;
3058 if (!App.IsDone()) {
3059 THROW_SALOME_CORBA_EXCEPTION("Filling aborted : non valid shape result", SALOME::BAD_PARAM);
3061 Standard_Integer UDegree, VDegree, NbUPoles, NbVPoles, NbUKnots, NbVKnots;
3062 App.SurfShape(UDegree, VDegree, NbUPoles, NbVPoles, NbUKnots, NbVKnots);
3063 Handle(Geom_BSplineSurface) GBS = new Geom_BSplineSurface(App.SurfPoles(),
3072 if( GBS.IsNull() ) {
3073 THROW_SALOME_CORBA_EXCEPTION("Make Filling aborted", SALOME::BAD_PARAM);
3075 tds = BRepBuilderAPI_MakeFace(GBS) ;
3077 catch(Standard_Failure) {
3078 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeFilling", SALOME::BAD_PARAM);
3081 /* We test the validity of resulting shape */
3082 if( !IsValid(tds) ) {
3083 THROW_SALOME_CORBA_EXCEPTION("Filling aborted : non valid shape result", SALOME::BAD_PARAM);
3086 result = CreateObject(tds) ;
3087 InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
3094 //=================================================================================
3095 // function : MakeGlueFaces()
3097 //=================================================================================
3099 TopoDS_Face GEOM_Gen_i::FindSameFace(const TopoDS_Shape& aShape,
3100 const TopoDS_Face& F,
3104 bool isSame = false;
3105 for (TopExp_Explorer exf(aShape,TopAbs_FACE); exf.More(); exf.Next())
3107 //MESSAGE("--- test a face");
3109 aFace = TopoDS::Face(exf.Current());
3110 TopTools_ListOfShape liste1;
3111 TopTools_ListOfShape liste2;
3112 for (TopExp_Explorer exp(aFace,TopAbs_VERTEX); exp.More(); exp.Next())
3114 const TopoDS_Vertex& V = TopoDS::Vertex(exp.Current());
3117 for (TopExp_Explorer exp(F,TopAbs_VERTEX); exp.More(); exp.Next())
3119 const TopoDS_Vertex& V = TopoDS::Vertex(exp.Current());
3123 if (liste1.Extent() == liste2.Extent())
3125 TopTools_ListIteratorOfListOfShape it1(liste1);
3127 for (; it1.More(); it1.Next())
3129 bool foundSamePoint = false;
3130 gp_Pnt P1 = BRep_Tool::Pnt(TopoDS::Vertex(it1.Value()));
3131 TopTools_ListIteratorOfListOfShape it2(liste2);
3132 for (it2; it2.More(); it2.Next())
3134 gp_Pnt P2 = BRep_Tool::Pnt(TopoDS::Vertex(it2.Value()));
3135 double d = P1.Distance(P2);
3139 //MESSAGE(" found Same Point : "<<nbFound<<" - "<<P1.X()<<" "<<P1.Y()<<" "<<P1.Z());
3140 foundSamePoint = true;
3144 isSame = isSame && foundSamePoint;
3145 if (! isSame) break; // a vertex does not correspond : not same face
3150 //MESSAGE(" --- Found Same Face");
3151 break; // a face corresponding to F is found
3154 if (! isSame) aFace.Nullify(); // return null face
3158 TopoDS_Edge GEOM_Gen_i::FindSameEdge(const TopoDS_Face& nf,
3162 TopoDS_Face newFace = TopoDS::Face(nf.Oriented(TopAbs_REVERSED));
3163 TopoDS_Vertex VFirst, VLast;
3164 TopExp::Vertices(Eold, VFirst, VLast);
3165 gp_Pnt Pf = BRep_Tool::Pnt(VFirst);
3166 gp_Pnt Pl = BRep_Tool::Pnt(VLast);
3168 for (TopExp_Explorer ee(newFace,TopAbs_EDGE); ee.More(); ee.Next())
3170 const TopoDS_Edge& E = TopoDS::Edge(ee.Current());
3171 TopoDS_Vertex VFn, VLn;
3172 TopExp::Vertices(E, VFn, VLn);
3173 gp_Pnt Pfn = BRep_Tool::Pnt(VFn);
3174 gp_Pnt Pln = BRep_Tool::Pnt(VLn);
3175 double dff = Pf.Distance(Pfn);
3176 double dfl = Pf.Distance(Pln);
3177 double dlf = Pl.Distance(Pfn);
3178 double dll = Pl.Distance(Pln);
3179 if ((dff < tol3d) && (dll <tol3d))
3181 //MESSAGE("--- edge forward " <<Pf.X()<<" "<<Pf.Y()<<" "<<Pf.Z()<<" "<<Pl.X()<<" "<<Pl.Y()<<" "<<Pl.Z());
3182 Enew = TopoDS::Edge(E.Oriented(TopAbs_FORWARD));
3183 Eold = TopoDS::Edge(Eold.Oriented(TopAbs_FORWARD));
3186 if ((dfl < tol3d) && (dlf <tol3d))
3188 //MESSAGE("--- edge reversed " <<Pf.X()<<" "<<Pf.Y()<<" "<<Pf.Z()<<" "<<Pl.X()<<" "<<Pl.Y()<<" "<<Pl.Z());
3189 Enew = TopoDS::Edge(E.Oriented(TopAbs_REVERSED));
3190 Eold = TopoDS::Edge(Eold.Oriented(TopAbs_FORWARD));
3197 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeGlueFaces(GEOM::GEOM_Shape_ptr myShape,
3199 throw (SALOME::SALOME_Exception)
3202 // prendre un premier shell dans la liste des shells
3203 // initialiser un compshell avec ce shell
3204 // tant qu'il reste des shells dans la liste
3205 // chercher un shell qui a des faces en commun avec le compshell
3206 // creer un BRepTools_Quilt
3207 // recenser les faces communes issues du compshell, les ajouter au quilt
3208 // recenser les faces restantes du shell a inclure, les ajouter au quilt
3209 // recenser les edges en double, a remplacer
3210 // pour chaque paire d'edge
3211 // tester l'orientation relative des aretes
3212 // bind dans le quilt de Eold.Forward et Enew.Forward (ou reverse)
3213 // recuperer le nouveau shell
3214 // l'incorporer dans le compshell
3215 // appliquer BRepTools_SameParameter au compshell
3216 // (rendre parametres 2D des edges identiques aux parametres 3D)
3218 Unexpect aCatch(SALOME_SalomeException);
3219 GEOM::GEOM_Shape_var result ;
3221 TopoDS_Shape aShape = GetTopoShape(myShape) ;
3224 bu.MakeCompound(C); // empty compound;
3225 TopTools_ListOfShape shellList;
3226 for (TopExp_Explorer exp(aShape,TopAbs_SHELL); exp.More(); exp.Next())
3228 const TopoDS_Shell& S = TopoDS::Shell(exp.Current());
3229 shellList.Append(S);
3231 TopTools_ListIteratorOfListOfShape its(shellList);
3234 THROW_SALOME_CORBA_EXCEPTION("glue aborted : no shell in shape", SALOME::BAD_PARAM);
3236 TopoDS_Shell S = TopoDS::Shell(its.Value());
3237 bu.Add(C, S); // add first shell to compound
3238 shellList.Remove(its);
3239 its.Initialize(shellList);
3240 bool shellAdded = true;
3241 while ((shellList.Extent() > 0) && shellAdded)
3243 //MESSAGE("more shells : "<< shellList.Extent());
3245 its.Initialize(shellList);
3246 for(; its.More(); its.Next())
3248 //MESSAGE("one more shell to try");
3249 TopTools_ListOfShape newFaces; // common faces from new compound
3250 TopTools_ListOfShape oldFaces; // common faces from shell to add
3251 TopTools_ListOfShape addFaces; // not common faces from shell to add
3252 TopTools_ListOfShape newEdges; // common edges from new compound
3253 TopTools_ListOfShape oldEdges; // common edges from face to add
3254 TopoDS_Compound CFN;
3255 TopoDS_Compound CFO;
3256 bu.MakeCompound(CFN); // empty compound for new faces
3257 bu.MakeCompound(CFO); // empty compound for old faces
3258 S = TopoDS::Shell(its.Value());
3259 for (TopExp_Explorer exp(S,TopAbs_FACE); exp.More(); exp.Next())
3261 //MESSAGE("--- try to find corresponding face in new compound");
3262 TopoDS_Face F = TopoDS::Face(exp.Current());
3263 TopoDS_Face newFace = FindSameFace(C,F,tol3d);
3264 if (! newFace.IsNull())
3266 //MESSAGE("--- face found");
3267 newFaces.Append(newFace);
3268 bu.Add(CFN, newFace); // common faces from new compound
3270 for (TopExp_Explorer ee(F,TopAbs_EDGE);ee.More();ee.Next())
3272 //MESSAGE("--- find edge pair");
3273 TopoDS_Edge Eold = TopoDS::Edge(ee.Current());
3274 const TopoDS_Edge& Enew = FindSameEdge(newFace, Eold, tol3d);
3275 oldEdges.Append(Eold);
3276 newEdges.Append(Enew);
3283 bu.Add(CFO, F); // not common faces from shell to add
3286 if ( !newFaces.IsEmpty())
3288 //MESSAGE("--- some faces found ---");
3290 BRepTools_Quilt glue;
3292 TopTools_ListIteratorOfListOfShape ito(oldEdges);
3293 TopTools_ListIteratorOfListOfShape itn(newEdges);
3294 for (; ito.More(); ito.Next())
3296 //MESSAGE("--- bind");
3297 glue.Bind(TopoDS::Edge(ito.Value()), TopoDS::Edge(itn.Value()));
3301 TopoDS_Compound newc = TopoDS::Compound(glue.Shells());
3302 for (TopExp_Explorer exs(newc,TopAbs_SHELL); exs.More(); exs.Next())
3304 TopoDS_Shell NS = TopoDS::Shell(exs.Current());
3307 shellList.Remove(its);
3308 //MESSAGE("--- remove shell from list");
3313 //MESSAGE("---" << shellList.Extent() << " " << shellAdded);
3315 TopExp_Explorer exp(C,TopAbs_SHELL);
3316 Standard_Integer ish=0;
3317 TopoDS_Compound Res;
3320 B.MakeCompound(Res);
3321 TopoDS_Shape theShape;
3323 for (; exp.More(); exp.Next())
3325 TopoDS_Shape Sh = exp.Current();
3328 BRepClass3d_SolidClassifier SC(Sol);
3329 SC.PerformInfinitePoint(1.E-6); // cf. BRepFill_Confusion() - BRepFill_Evolved.cxx
3330 if (SC.State() == TopAbs_IN)
3333 B.Add(Sol,Sh.Reversed());
3338 if (ish == 1) { theShape = Sol;}
3339 else { theShape = Res;}
3341 BRepLib::SameParameter(theShape, 1.E-5, Standard_True);
3343 result = CreateObject(tds);
3344 InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
3349 //=================================================================================
3350 // function : MakeSewing()
3352 //=================================================================================
3353 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeSewing( const GEOM::GEOM_Gen::ListOfIOR& ListShapes,
3354 CORBA::Double precision )
3355 throw (SALOME::SALOME_Exception)
3357 Unexpect aCatch(SALOME_SalomeException);
3358 GEOM::GEOM_Shape_var result ;
3360 BRepOffsetAPI_Sewing aMethod ;
3363 /* default OCC is 1.0e-06 */
3364 aMethod.Init(precision, Standard_False);
3365 for ( unsigned int i = 0; i < ListShapes.length(); i++) {
3366 GEOM::GEOM_Shape_var aShape = GetIORFromString( ListShapes[i] );
3367 TopoDS_Shape Shape = GetTopoShape(aShape) ;
3368 if( Shape.IsNull() ) {
3369 THROW_SALOME_CORBA_EXCEPTION("MakeSewing aborted : null shape during operation", SALOME::BAD_PARAM);
3371 aMethod.Add(Shape) ;
3375 tds = aMethod.SewedShape() ;
3376 if( !IsValid(tds) ) {
3377 THROW_SALOME_CORBA_EXCEPTION("Make Sewing aborted : non valid shape", SALOME::BAD_PARAM);
3379 if( tds.IsNull() ) {
3380 THROW_SALOME_CORBA_EXCEPTION("Make Sewing aborted : null shape", SALOME::BAD_PARAM);
3383 catch (Standard_Failure) {
3384 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeSewing", SALOME::BAD_PARAM);
3387 result = CreateObject(tds);
3388 InsertInLabelMoreArguments(tds, result, ListShapes, myCurrentOCAFDoc) ;
3392 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeSewingShape( GEOM::GEOM_Shape_ptr aShape,
3393 CORBA::Double precision )
3394 throw (SALOME::SALOME_Exception)
3396 Unexpect aCatch(SALOME_SalomeException);
3397 GEOM::GEOM_Shape_var result ;
3398 TopoDS_Shape tds, S ;
3399 BRepOffsetAPI_Sewing aMethod ;
3402 S = GetTopoShape(aShape) ;
3404 THROW_SALOME_CORBA_EXCEPTION("In Sewing a Shape is null", SALOME::BAD_PARAM);
3407 /* default OCC is 1.0e-06 */
3408 aMethod.Init(precision, Standard_False);
3409 for ( TopExp_Explorer exp( S, TopAbs_FACE); exp.More(); exp.Next() ) {
3410 const TopoDS_Face& F = TopoDS::Face(exp.Current());
3415 tds = aMethod.SewedShape() ;
3416 if( !IsValid(tds) ) {
3417 THROW_SALOME_CORBA_EXCEPTION("Make Sewing aborted : non valid shape", SALOME::BAD_PARAM);
3420 catch (Standard_Failure) {
3421 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeSewing", SALOME::BAD_PARAM);
3424 result = CreateObject(tds);
3425 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3426 result->ShapeId(entry) ;
3430 //==================================================================================
3431 // function : OrientationChange()
3432 // purpose : Change the orientation of a new shape
3433 // : TopAbs_FORWARD < -- > TopAbs_REVERSED
3435 // : WARNING : for the moment we make a new shape !
3436 //==================================================================================
3437 GEOM::GEOM_Shape_ptr GEOM_Gen_i::OrientationChange(GEOM::GEOM_Shape_ptr aShape)
3438 throw (SALOME::SALOME_Exception)
3440 Unexpect aCatch(SALOME_SalomeException);
3441 GEOM::GEOM_Shape_var result ;
3442 BRep_Builder aBuilder;
3444 TopoDS_Shape shape = GetTopoShape(aShape) ;
3445 if( shape.IsNull() ) {
3446 THROW_SALOME_CORBA_EXCEPTION("Shape is null", SALOME::BAD_PARAM);
3449 BRepBuilderAPI_Copy Copy(shape);
3450 if( Copy.IsDone() ) {
3451 TopoDS_Shape tds = Copy.Shape();
3452 if( tds.IsNull() ) {
3453 THROW_SALOME_CORBA_EXCEPTION("Orientation aborted : null shape", SALOME::BAD_PARAM);
3456 if( tds.Orientation() == TopAbs_FORWARD)
3457 tds.Orientation(TopAbs_REVERSED) ;
3459 tds.Orientation(TopAbs_FORWARD) ;
3461 result = CreateObject(tds);
3462 InsertInLabelOneArgument(shape, aShape, tds, result, myCurrentOCAFDoc) ;
3468 //==================================================================================
3469 // function : GetReferencedObjects()
3471 //==================================================================================
3472 GEOM::GEOM_Gen::ListOfIOR* GEOM_Gen_i::GetReferencedObjects(GEOM::GEOM_Shape_ptr shape)
3474 GEOM::GEOM_Gen::ListOfIOR_var aList = new GEOM::GEOM_Gen::ListOfIOR;
3477 if (shape->_is_nil()) return aList._retn();
3479 Standard_CString entry = shape->ShapeId();
3481 TDF_Tool::Label(myCurrentOCAFDoc->GetData(), entry, Lab);
3483 Handle(TDataStd_Name) Att;
3484 Lab.FindAttribute(TDataStd_Name::GetID(),Att);
3486 TDF_ChildIterator ChildIterator(Lab);
3487 if (ChildIterator.More()) {
3488 TDF_Label L = ChildIterator.Value();
3489 Handle(TDataStd_Name) Att;
3490 L.FindAttribute(TDataStd_Name::GetID(),Att);
3491 if (Att->Get().IsEqual(TCollection_ExtendedString("Arguments")) ) {
3493 TDF_ChildIterator ChildIterator1(L);
3496 while (ChildIterator1.More()) {
3497 TDF_Label L = ChildIterator1.Value();
3499 Handle(TDF_Reference) Ref;
3500 if (L.FindAttribute(TDF_Reference::GetID(),Ref)) {
3503 ChildIterator1.Next();
3507 TDF_ChildIterator ChildIterator2(L);
3508 while (ChildIterator2.More()) {
3509 TDF_Label L = ChildIterator2.Value();
3510 Handle(TDF_Reference) Ref;
3511 if (L.FindAttribute(TDF_Reference::GetID(),Ref)) {
3512 TDF_Label L = Ref->Get();
3514 Handle(TDataStd_Name) Att;
3515 L.FindAttribute(TDataStd_Name::GetID(),Att);
3516 TCollection_AsciiString nameIOR (Att->Get()) ;
3517 aList[i] = CORBA::string_dup( nameIOR.ToCString() );
3521 ChildIterator2.Next();
3525 return aList._retn();
3528 //==================================================================================
3529 // function : GetObjects()
3531 //==================================================================================
3532 GEOM::GEOM_Gen::ListOfIOR* GEOM_Gen_i::GetObjects(GEOM::GEOM_Shape_ptr shape)
3534 GEOM::GEOM_Gen::ListOfIOR_var aList = new GEOM::GEOM_Gen::ListOfIOR;
3537 Standard_CString entry = shape->ShapeId();
3539 TDF_Tool::Label(myCurrentOCAFDoc->GetData(), entry, Lab);
3541 Handle(TDataStd_Name) Att;
3542 Lab.FindAttribute(TDataStd_Name::GetID(),Att);
3544 TDF_ChildIterator ChildIterator(Lab);
3546 while (ChildIterator.More()) {
3547 TDF_Label L = ChildIterator.Value();
3548 Handle(TDataStd_Name) Att;
3549 L.FindAttribute(TDataStd_Name::GetID(),Att);
3551 if (!Att->Get().IsEqual(TCollection_ExtendedString("Arguments")) ) {
3554 ChildIterator.Next();
3559 TDF_ChildIterator ChildIterator1(Lab);
3560 while (ChildIterator1.More()) {
3561 TDF_Label L = ChildIterator1.Value();
3562 Handle(TDataStd_Name) Att;
3563 L.FindAttribute(TDataStd_Name::GetID(),Att);
3565 if (!Att->Get().IsEqual(TCollection_ExtendedString("Arguments")) ) {
3566 TCollection_AsciiString nameIOR (Att->Get());
3567 aList[i] = CORBA::string_dup( nameIOR.ToCString() );
3570 ChildIterator1.Next();
3572 return aList._retn();
3576 //==================================================================================
3577 // function : Import
3578 // purpose : Import shape from a BREP file
3579 //==================================================================================
3580 GEOM::GEOM_Shape_ptr GEOM_Gen_i::ImportBREP(const char* filename)
3581 throw (SALOME::SALOME_Exception)
3583 Unexpect aCatch(SALOME_SalomeException);
3585 GEOM::GEOM_Shape_var result ;
3588 BRep_Builder aBuilder;
3589 char* aCopyfilename = strdup(filename);
3590 BRepTools::Read(tds, aCopyfilename, aBuilder) ;
3591 free(aCopyfilename);
3593 THROW_SALOME_CORBA_EXCEPTION("Import BRep aborted", SALOME::BAD_PARAM);
3596 catch(Standard_Failure) {
3597 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ImportBREP", SALOME::BAD_PARAM);
3600 result = CreateObject(tds) ;
3601 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3602 result->ShapeId(entry);
3607 //================================================================================
3608 // function : MakePlane()
3609 // purpose : Make a plane topology (non infinite)
3610 //================================================================================
3611 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakePlane(const GEOM::PointStruct& pstruct,
3612 const GEOM::DirStruct& dstruct,
3613 CORBA::Double trimsize)
3614 throw (SALOME::SALOME_Exception)
3616 Unexpect aCatch(SALOME_SalomeException);
3617 GEOM::GEOM_Shape_var result ;
3621 gp_Pnt aPoint(pstruct.x, pstruct.y, pstruct.z) ;
3622 gp_Dir aDirection(dstruct.PS.x, dstruct.PS.y, dstruct.PS.z) ;
3623 /* we make a trimmed plane */
3624 gp_Pln gplane(aPoint, aDirection) ;
3625 tds = BRepBuilderAPI_MakeFace(gplane, -trimsize, +trimsize, -trimsize, +trimsize) ;
3627 catch(Standard_Failure) {
3628 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakePlane", SALOME::BAD_PARAM);
3632 THROW_SALOME_CORBA_EXCEPTION("Make Plane aborted : null shape", SALOME::BAD_PARAM);
3635 result = CreateObject(tds) ;
3636 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3637 result->ShapeId(entry);
3641 //=================================================================================
3642 // function : MakeVertex()
3643 // purpose : Create a Vertex topology.
3644 //=================================================================================
3645 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeVertex(CORBA::Double x,
3648 throw (SALOME::SALOME_Exception)
3650 Unexpect aCatch(SALOME_SalomeException);
3651 GEOM::GEOM_Shape_var result ;
3653 TopoDS_Shape tds = BRepBuilderAPI_MakeVertex(P).Shape();
3655 THROW_SALOME_CORBA_EXCEPTION("Make Vertex/Point aborted", SALOME::BAD_PARAM);
3658 result = CreateObject(tds) ;
3659 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3660 result->ShapeId(entry);
3665 //=================================================================================
3666 // function : MakeFace()
3668 //=================================================================================
3669 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeFace( GEOM::GEOM_Shape_ptr wire,
3670 CORBA::Boolean wantplanarface )
3671 throw (SALOME::SALOME_Exception)
3673 Unexpect aCatch(SALOME_SalomeException);
3674 GEOM::GEOM_Shape_var result ;
3675 TopoDS_Shape aShape;
3679 aShape = GetTopoShape(wire) ;
3680 if( aShape.IsNull() || aShape.ShapeType() != TopAbs_WIRE ) {
3681 THROW_SALOME_CORBA_EXCEPTION("MakeFace aborted : null or inappropriate shape", SALOME::BAD_PARAM);
3683 TopoDS_Wire W = TopoDS::Wire(aShape) ;
3684 tds = BRepBuilderAPI_MakeFace(W, wantplanarface).Shape() ;
3685 if( !tds.IsNull() ) {
3686 result = CreateObject(tds) ;
3687 InsertInLabelOneArgument(aShape, wire, tds, result, myCurrentOCAFDoc) ;
3690 THROW_SALOME_CORBA_EXCEPTION("Null result in GEOM_Gen_i::MakeFace", SALOME::BAD_PARAM);
3693 catch (Standard_Failure) {
3694 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeFace", SALOME::BAD_PARAM);
3700 //=================================================================================
3701 // function : MakeFaces()
3703 //=================================================================================
3704 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeFaces(const GEOM::GEOM_Gen::ListOfIOR& ListShapes,
3705 CORBA::Boolean wantplanarface)
3706 throw (SALOME::SALOME_Exception)
3708 GEOM::GEOM_Shape_var result;
3711 GEOM::GEOM_Shape_var aShape = GetIORFromString(ListShapes[0]);
3712 TopoDS_Shape Shape = GetTopoShape(aShape);
3713 if(Shape.IsNull() || Shape.ShapeType() != TopAbs_WIRE) {
3714 THROW_SALOME_CORBA_EXCEPTION("Shell aborted : null shape during operation", SALOME::BAD_PARAM);
3716 TopoDS_Wire W = TopoDS::Wire(Shape);
3717 TopoDS_Shape FFace = BRepBuilderAPI_MakeFace(W, wantplanarface).Shape();
3718 if(!FFace.IsNull()) {
3719 if(ListShapes.length() == 1) {
3720 result = CreateObject(FFace);
3721 InsertInLabelMoreArguments(FFace, result, ListShapes, myCurrentOCAFDoc);
3723 else if(ListShapes.length() >= 2) {
3725 BRep_Builder aBuilder;
3726 aBuilder.MakeCompound(C);
3727 BRepAlgo_FaceRestrictor FR;
3729 TopAbs_Orientation OriF = FFace.Orientation();
3730 TopoDS_Shape aLocalS = FFace.Oriented(TopAbs_FORWARD);
3731 FR.Init(TopoDS::Face(aLocalS), Standard_False, Standard_True);
3733 for(unsigned int i = 0; i < ListShapes.length(); i++) {
3734 GEOM::GEOM_Shape_var aShape = GetIORFromString(ListShapes[i]);
3735 TopoDS_Shape Shape = GetTopoShape(aShape);
3736 if(Shape.IsNull()) {
3737 THROW_SALOME_CORBA_EXCEPTION("Shell aborted : null shape during operation", SALOME::BAD_PARAM);
3739 FR.Add(TopoDS::Wire(Shape));
3747 for(; FR.More(); FR.Next()) {
3748 aFace = FR.Current().Oriented(OriF);
3749 aBuilder.Add(C, aFace);
3753 result = CreateObject(aFace);
3754 InsertInLabelMoreArguments(aFace, result, ListShapes, myCurrentOCAFDoc);
3757 result = CreateObject(C);
3758 InsertInLabelMoreArguments(C, result, ListShapes, myCurrentOCAFDoc);
3764 THROW_SALOME_CORBA_EXCEPTION("Null result in GEOM_Gen_i::MakeFace", SALOME::BAD_PARAM);
3767 catch (Standard_Failure) {
3768 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeFace", SALOME::BAD_PARAM);
3774 //=================================================================================
3775 // function : MakeShell()
3776 // purpose : Make a compound from a list containing one or more shapes
3777 //=================================================================================
3778 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeShell( const GEOM::GEOM_Gen::ListOfIOR& ListShapes )
3779 throw (SALOME::SALOME_Exception)
3781 Unexpect aCatch(SALOME_SalomeException);
3782 GEOM::GEOM_Shape_var result ;
3783 BRepTools_Quilt Glue;
3786 for ( unsigned int i = 0; i < ListShapes.length(); i++) {
3787 GEOM::GEOM_Shape_var aShape = GetIORFromString( ListShapes[i] );
3788 TopoDS_Shape Shape = GetTopoShape(aShape) ;
3789 if( Shape.IsNull() ) {
3790 THROW_SALOME_CORBA_EXCEPTION("Shell aborted : null shape during operation", SALOME::BAD_PARAM);
3795 TopExp_Explorer exp(Glue.Shells(), TopAbs_SHELL);
3796 Standard_Integer ish = 0;
3797 for (; exp.More(); exp.Next()) {
3806 THROW_SALOME_CORBA_EXCEPTION("Null result : Shell operation aborted", SALOME::BAD_PARAM);
3809 result = CreateObject(C) ;
3810 InsertInLabelMoreArguments(C, result, ListShapes, myCurrentOCAFDoc) ;
3816 //=================================================================================
3817 // function : MakeSolid()
3818 // purpose : Make a compound from a list containing one or more shapes
3819 //=================================================================================
3820 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeSolid( const GEOM::GEOM_Gen::ListOfIOR& ListShapes )
3821 throw (SALOME::SALOME_Exception)
3823 Unexpect aCatch(SALOME_SalomeException);
3824 GEOM::GEOM_Shape_var result ;
3825 Standard_Integer ish = 0;
3826 TopoDS_Compound Res;
3831 B.MakeCompound(Res);
3833 for ( unsigned int i = 0; i < ListShapes.length(); i++) {
3834 GEOM::GEOM_Shape_var aShape = GetIORFromString( ListShapes[i] );
3835 TopoDS_Shape Sh = GetTopoShape(aShape) ;
3837 THROW_SALOME_CORBA_EXCEPTION("Solid aborted : null shape during operation", SALOME::BAD_PARAM);
3841 BRepClass3d_SolidClassifier SC(Sol);
3842 SC.PerformInfinitePoint(Precision::Confusion());
3843 if (SC.State() == TopAbs_IN) {
3845 B.Add(Sol,Sh.Reversed());
3850 if (ish == 1) { Shape = Sol;}
3851 else { Shape = Res;}
3853 if ( Shape.IsNull() ) {
3854 THROW_SALOME_CORBA_EXCEPTION("Null result : Solid operation aborted", SALOME::BAD_PARAM);
3857 result = CreateObject(Shape) ;
3858 InsertInLabelMoreArguments(Shape, result, ListShapes, myCurrentOCAFDoc) ;
3864 //================================================================================
3865 // function : MakeLine
3866 // purpose : Make a Line topology
3867 //================================================================================
3868 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeLine(const GEOM::PointStruct& pstruct,
3869 const GEOM::DirStruct& dstruct)
3870 throw (SALOME::SALOME_Exception)
3872 Unexpect aCatch(SALOME_SalomeException);
3873 GEOM::GEOM_Shape_var result ;
3874 gp_Pnt P1(pstruct.x, pstruct.y, pstruct.z);
3875 gp_Pnt P2(dstruct.PS.x, dstruct.PS.y, dstruct.PS.z) ;
3879 tds = BRepBuilderAPI_MakeEdge(P1, P2).Shape();
3881 catch(Standard_Failure) {
3882 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeLine", SALOME::BAD_PARAM);
3885 if ( tds.IsNull() ) {
3886 THROW_SALOME_CORBA_EXCEPTION("Make Line aborted : null shape", SALOME::BAD_PARAM);
3890 result = CreateObject(tds) ;
3891 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3892 result->ShapeId(entry);
3898 //================================================================================
3899 // function : MakeVector()
3900 // purpose : Make a vector
3901 //================================================================================
3902 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeVector(const GEOM::PointStruct& pstruct1,
3903 const GEOM::PointStruct& pstruct2)
3904 throw (SALOME::SALOME_Exception)
3906 Unexpect aCatch(SALOME_SalomeException);
3907 GEOM::GEOM_Shape_var result ;
3911 gp_Pnt P1(pstruct1.x, pstruct1.y, pstruct1.z);
3912 gp_Pnt P2(pstruct2.x, pstruct2.y, pstruct2.z) ;
3913 tds = BRepBuilderAPI_MakeEdge(P1, P2).Shape();
3915 catch(Standard_Failure) {
3916 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeVector", SALOME::BAD_PARAM);
3919 if ( tds.IsNull() ) {
3920 THROW_SALOME_CORBA_EXCEPTION("Make Vector aborted : null shape", SALOME::BAD_PARAM);
3923 result = CreateObject(tds) ;
3924 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3925 result->ShapeId(entry);
3931 //================================================================================
3932 // function : MakeCircle()
3934 //================================================================================
3935 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeCircle(const GEOM::PointStruct& pstruct,
3936 const GEOM::DirStruct& dstruct,
3937 CORBA::Double radius)
3938 throw (SALOME::SALOME_Exception)
3940 Unexpect aCatch(SALOME_SalomeException);
3941 GEOM::GEOM_Shape_var result;
3945 gp_Pnt p(pstruct.x, pstruct.y, pstruct.z) ;
3946 gp_Dir d(dstruct.PS.x, dstruct.PS.y, dstruct.PS.z) ;
3948 gp_Circ circ( axis, radius);
3949 BRepBuilderAPI_MakeEdge MakeEdge( circ );
3950 tds = MakeEdge.Edge();
3952 catch(Standard_Failure) {
3953 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeCircle", SALOME::BAD_PARAM);
3956 THROW_SALOME_CORBA_EXCEPTION("Make Circle aborted", SALOME::BAD_PARAM);
3958 result = CreateObject(tds);
3959 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3960 result->ShapeId(entry);
3964 //================================================================================
3965 // function : MakeEllipse()
3967 //================================================================================
3968 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeEllipse( const GEOM::PointStruct& pstruct,
3969 const GEOM::DirStruct& dstruct,
3970 CORBA::Double radius_major,
3971 CORBA::Double radius_minor )
3972 throw (SALOME::SALOME_Exception)
3974 Unexpect aCatch(SALOME_SalomeException);
3975 GEOM::GEOM_Shape_var result;
3979 gp_Pnt p(pstruct.x, pstruct.y, pstruct.z) ;
3980 gp_Dir d(dstruct.PS.x, dstruct.PS.y, dstruct.PS.z) ;
3982 const gp_Ax2 axis(p, d) ;
3983 gp_Elips anEllipse( axis, radius_major, radius_minor ) ;
3984 BRepBuilderAPI_MakeEdge MakeEdge( anEllipse );
3985 tds = MakeEdge.Edge();
3988 catch(Standard_Failure) {
3989 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeEllipse", SALOME::BAD_PARAM);
3992 THROW_SALOME_CORBA_EXCEPTION("Make Ellipse aborted", SALOME::BAD_PARAM);
3994 result = CreateObject(tds);
3995 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3996 result->ShapeId(entry);
4000 //================================================================================
4001 // function : MakeArc()
4002 // purpose : make an arc of circle from pInit to pEnd and passing on pCircle
4003 //================================================================================
4004 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeArc(const GEOM::PointStruct& pInit,
4005 const GEOM::PointStruct& pCircle,
4006 const GEOM::PointStruct& pEnd)
4007 throw (SALOME::SALOME_Exception)
4009 Unexpect aCatch(SALOME_SalomeException);
4010 GEOM::GEOM_Shape_var result;
4012 gp_Pnt pI(pInit.x, pInit.y, pInit.z) ;
4013 gp_Pnt pC(pCircle.x, pCircle.y, pCircle.z) ;
4014 gp_Pnt pE(pEnd.x, pEnd.y, pEnd.z) ;
4016 GC_MakeArcOfCircle arc( pI, pC, pE ) ;
4017 if( !arc.IsDone() ) {
4018 THROW_SALOME_CORBA_EXCEPTION("Arc not done", SALOME::BAD_PARAM);
4020 BRepBuilderAPI_MakeEdge MakeEdge( arc );
4021 TopoDS_Shape tds = MakeEdge.Edge();
4023 THROW_SALOME_CORBA_EXCEPTION("Null result : arc not done", SALOME::BAD_PARAM);
4026 result = CreateObject(tds);
4027 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
4028 result->ShapeId(entry);
4031 catch(Standard_Failure) {
4032 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeArc", SALOME::BAD_PARAM);
4037 //=================================================================================
4038 // function : MakeSketcher()
4039 // purpose : Make a wire from a list containing many points
4040 //=================================================================================
4041 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeSketcher( const char *Cmd )
4042 throw (SALOME::SALOME_Exception)
4044 GEOM::GEOM_Shape_var result ;
4047 Sketcher_Profile aProfile (Cmd);
4048 if(aProfile.IsDone())
4049 tds = aProfile.GetShape();
4051 catch(Standard_Failure) {
4052 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeSketcher", SALOME::BAD_PARAM);
4056 THROW_SALOME_CORBA_EXCEPTION("MakeSketcher aborted : null shape", SALOME::BAD_PARAM);
4059 result = CreateObject(tds);
4060 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
4061 result->ShapeId(entry) ;
4068 //=================================================================================
4069 // function : MakeBezier()
4070 // purpose : Make a wire from a list containing many points
4071 //=================================================================================
4072 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeBezier( const GEOM::GEOM_Gen::ListOfIOR& ListShapes )
4073 throw (SALOME::SALOME_Exception)
4075 GEOM::GEOM_Shape_var result;
4076 TopoDS_Shape tds, Shape;
4077 TColgp_Array1OfPnt CurvePoints(1, ListShapes.length());
4080 for(unsigned int i = 0; i < ListShapes.length(); i++) {
4081 GEOM::GEOM_Shape_var aShape = GetIORFromString(ListShapes[i]);
4082 Shape = GetTopoShape(aShape);
4083 if(Shape.IsNull()) {
4084 THROW_SALOME_CORBA_EXCEPTION("MakeBezier aborted : null shape during operation", SALOME::BAD_PARAM);
4086 if(Shape.ShapeType() == TopAbs_VERTEX) {
4087 const gp_Pnt& P = BRep_Tool::Pnt(TopoDS::Vertex(Shape));
4088 CurvePoints.SetValue(i + 1, P);
4091 Handle(Geom_BezierCurve) GBC = new Geom_BezierCurve(CurvePoints);
4092 tds = BRepBuilderAPI_MakeEdge(GBC);
4094 catch(Standard_Failure) {
4095 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeBezier", SALOME::BAD_PARAM);
4098 if( tds.IsNull() ) {
4099 THROW_SALOME_CORBA_EXCEPTION("Make Bezier operation aborted : null result", SALOME::BAD_PARAM);
4102 result = CreateObject(tds);
4103 InsertInLabelMoreArguments(tds, result, ListShapes, myCurrentOCAFDoc);
4109 //=================================================================================
4110 // function : MakeInterpol()
4111 // purpose : Make a wire from a list containing many points
4112 //=================================================================================
4113 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeInterpol( const GEOM::GEOM_Gen::ListOfIOR& ListShapes )
4114 throw (SALOME::SALOME_Exception)
4116 GEOM::GEOM_Shape_var result;
4117 TopoDS_Shape tds, Shape;
4118 TColgp_Array1OfPnt CurvePoints(1, ListShapes.length());
4121 for(unsigned int i = 0; i < ListShapes.length(); i++) {
4122 GEOM::GEOM_Shape_var aShape = GetIORFromString(ListShapes[i]);
4123 Shape = GetTopoShape(aShape);
4124 if(Shape.IsNull()) {
4125 THROW_SALOME_CORBA_EXCEPTION("MakeBSpline aborted : null shape during operation", SALOME::BAD_PARAM);
4127 if(Shape.ShapeType() == TopAbs_VERTEX) {
4128 const gp_Pnt& P = BRep_Tool::Pnt(TopoDS::Vertex(Shape));
4129 CurvePoints.SetValue(i + 1, P);
4132 GeomAPI_PointsToBSpline GBC(CurvePoints);
4133 tds = BRepBuilderAPI_MakeEdge(GBC);
4135 catch(Standard_Failure) {
4136 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeBSpline", SALOME::BAD_PARAM);
4139 if( tds.IsNull() ) {
4140 THROW_SALOME_CORBA_EXCEPTION("Make BSpline operation aborted : null result", SALOME::BAD_PARAM);
4143 result = CreateObject(tds) ;
4144 InsertInLabelMoreArguments(tds, result, ListShapes, myCurrentOCAFDoc) ;
4149 //=================================================================================
4150 // function : MakeTranslation()
4151 // purpose : Translate a 3D shape
4152 //=================================================================================
4153 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeTranslation( GEOM::GEOM_Shape_ptr myShape,
4157 throw (SALOME::SALOME_Exception)
4159 Unexpect aCatch(SALOME_SalomeException);
4160 GEOM::GEOM_Shape_var result ;
4161 TopoDS_Shape aShape = GetTopoShape(myShape) ;
4162 if( aShape.IsNull() ) {
4163 THROW_SALOME_CORBA_EXCEPTION("Translation aborted : null shape", SALOME::BAD_PARAM);
4165 gp_Vec theVector(x,y,z) ;
4166 gp_Trsf theTransformation ;
4167 theTransformation.SetTranslation(theVector) ;
4168 BRepBuilderAPI_Transform myBRepTransformation(aShape, theTransformation, Standard_False) ;
4169 TopoDS_Shape tds = myBRepTransformation.Shape() ;
4171 result = CreateObject(tds) ;
4172 if( CORBA::is_nil(result) ) {
4173 THROW_SALOME_CORBA_EXCEPTION("Translation aborted : null result", SALOME::BAD_PARAM);
4175 InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4180 //=================================================================================
4181 // function : MakeMultiTranslation1D()
4182 // purpose : Multi-Translate a 3D shape
4183 //=================================================================================
4184 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeMultiTranslation1D( GEOM::GEOM_Shape_ptr myShape,
4185 const GEOM::DirStruct& dir,
4187 CORBA::Short nbtimes )
4188 throw (SALOME::SALOME_Exception)
4190 Unexpect aCatch(SALOME_SalomeException);
4191 GEOM::GEOM_Shape_var result ;
4194 TopoDS_Shape aShape = GetTopoShape(myShape) ;
4195 if( aShape.IsNull() )
4196 THROW_SALOME_CORBA_EXCEPTION("MakeMultiTranslation1D aborted : null shape", SALOME::BAD_PARAM);
4201 gp_Trsf theTransformation ;
4203 gp_Vec Vec( dir.PS.x, dir.PS.y, dir.PS.z ) ;
4205 TopoDS_Compound compound;
4207 B.MakeCompound( compound );
4209 for ( i = 0; i < nbtimes; i++ ) {
4210 DX = i * step * Vec.X() ;
4211 DY = i * step * Vec.Y() ;
4212 DZ = i * step * Vec.Z() ;
4213 myVec.SetCoord( DX, DY, DZ ) ;
4214 theTransformation.SetTranslation(myVec) ;
4215 BRepBuilderAPI_Transform myBRepTransformation(aShape, theTransformation, Standard_False) ;
4216 B.Add( compound, myBRepTransformation.Shape() );
4219 result = CreateObject(tds) ;
4221 catch (Standard_Failure) {
4222 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeMultiTranslation1D", SALOME::BAD_PARAM);
4225 if( CORBA::is_nil(result) ) {
4226 THROW_SALOME_CORBA_EXCEPTION("MakeMultiTranslation1D aborted : null result", SALOME::BAD_PARAM);
4228 InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4233 //=================================================================================
4234 // function : MakeMultiTranslation2D()
4235 // purpose : Multi-Translate a 3D shape
4236 //=================================================================================
4237 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeMultiTranslation2D( GEOM::GEOM_Shape_ptr myShape,
4238 const GEOM::DirStruct& dir1,
4239 CORBA::Double step1,
4240 CORBA::Short nbtimes1,
4241 const GEOM::DirStruct& dir2,
4242 CORBA::Double step2,
4243 CORBA::Short nbtimes2 )
4244 throw (SALOME::SALOME_Exception)
4246 Unexpect aCatch(SALOME_SalomeException);
4247 GEOM::GEOM_Shape_var result ;
4250 TopoDS_Shape aShape = GetTopoShape(myShape) ;
4251 if( aShape.IsNull() ) {
4252 THROW_SALOME_CORBA_EXCEPTION("MakeMultiTranslation2D aborted : null shape", SALOME::BAD_PARAM);
4258 gp_Trsf theTransformation ;
4260 gp_Vec Vec1( dir1.PS.x, dir1.PS.y, dir1.PS.z ) ;
4262 gp_Vec Vec2( dir2.PS.x, dir2.PS.y, dir2.PS.z ) ;
4264 TopoDS_Compound compound;
4266 B.MakeCompound( compound );
4268 for ( i = 0; i < nbtimes1; i++ ) {
4269 for ( j = 0; j < nbtimes2; j++ ) {
4270 DX = i * step1 * Vec1.X() + j * step2 * Vec2.X() ;
4271 DY = i * step1 * Vec1.Y() + j * step2 * Vec2.Y() ;
4272 DZ = i * step1 * Vec1.Z() + j * step2 * Vec2.Z() ;
4273 myVec.SetCoord( DX, DY, DZ ) ;
4274 theTransformation.SetTranslation(myVec) ;
4275 BRepBuilderAPI_Transform myBRepTransformation(aShape, theTransformation, Standard_False) ;
4276 B.Add( compound, myBRepTransformation.Shape() );
4280 result = CreateObject(tds) ;
4282 catch(Standard_Failure) {
4283 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeMultiTranslation2D", SALOME::BAD_PARAM);
4286 if( CORBA::is_nil(result) ) {
4287 THROW_SALOME_CORBA_EXCEPTION("MakeMultiTranslation2D aborted : null result", SALOME::BAD_PARAM);
4289 InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4294 //=================================================================================
4295 // function : MakeMultiRotation1D()
4296 // purpose : Multi-Rotate a 3D shape
4297 //=================================================================================
4298 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeMultiRotation1D( GEOM::GEOM_Shape_ptr myShape,
4299 const GEOM::DirStruct& dir,
4300 const GEOM::PointStruct& loc,
4301 CORBA::Short nbtimes)
4302 throw (SALOME::SALOME_Exception)
4304 Unexpect aCatch(SALOME_SalomeException);
4305 GEOM::GEOM_Shape_var result ;
4307 TopoDS_Shape aShape = GetTopoShape(myShape) ;
4308 if( aShape.IsNull() ) {
4309 THROW_SALOME_CORBA_EXCEPTION("MakeMultiRotation1D aborted : null shape", SALOME::BAD_PARAM);
4315 gp_Pnt P(loc.x, loc.y, loc.z) ;
4316 gp_Dir D(dir.PS.x, dir.PS.y, dir.PS.z) ;
4319 double angle = 360.0/nbtimes ;
4320 gp_Trsf theTransformation ;
4321 TopoDS_Compound compound;
4323 B.MakeCompound( compound );
4325 for ( i = 0; i < nbtimes; i++ ) {
4326 theTransformation.SetRotation(AX1, i*angle*PI180) ;
4327 BRepBuilderAPI_Transform myBRepTransformation(aShape, theTransformation, Standard_False) ;
4328 B.Add( compound, myBRepTransformation.Shape() );
4331 result = CreateObject(tds) ;
4333 catch(Standard_Failure) {
4334 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeMultiRotation1D", SALOME::BAD_PARAM);
4337 if( CORBA::is_nil(result) ) {
4338 THROW_SALOME_CORBA_EXCEPTION("MakeMultiRotation1D aborted : null result", SALOME::BAD_PARAM);
4340 InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4345 //=================================================================================
4346 // function : MakeMultiRotation2D()
4347 // purpose : Multi-Rotate a 3D shape
4348 //=================================================================================
4349 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeMultiRotation2D( GEOM::GEOM_Shape_ptr myShape,
4350 const GEOM::DirStruct& dir,
4351 const GEOM::PointStruct& loc,
4353 CORBA::Short nbtimes1,
4355 CORBA::Short nbtimes2 )
4356 throw (SALOME::SALOME_Exception)
4358 Unexpect aCatch(SALOME_SalomeException);
4359 GEOM::GEOM_Shape_var result ;
4361 TopoDS_Shape aShape = GetTopoShape(myShape) ;
4362 if( aShape.IsNull() ) {
4363 THROW_SALOME_CORBA_EXCEPTION("MakeMultiRotation2D aborted : null shape", SALOME::BAD_PARAM);
4370 gp_Pnt P(loc.x, loc.y, loc.z) ;
4371 gp_Dir D(dir.PS.x, dir.PS.y, dir.PS.z) ;
4373 gp_Trsf theTransformation1 ;
4374 gp_Trsf theTransformation2 ;
4376 GProp_GProps System ;
4378 if ( aShape.ShapeType() == TopAbs_VERTEX) {
4379 P1 = BRep_Tool::Pnt(TopoDS::Vertex( aShape ));
4381 else if ( aShape.ShapeType() == TopAbs_EDGE || aShape.ShapeType() == TopAbs_WIRE ) {
4382 BRepGProp::LinearProperties(aShape, System);
4383 P1 = System.CentreOfMass() ;
4385 else if ( aShape.ShapeType() == TopAbs_FACE || aShape.ShapeType() == TopAbs_SHELL ) {
4386 BRepGProp::SurfaceProperties(aShape, System);
4387 P1 = System.CentreOfMass() ;
4390 BRepGProp::VolumeProperties(aShape, System);
4391 P1 = System.CentreOfMass() ;
4394 Handle(Geom_Line) Line = new Geom_Line(AX1);
4395 gp_Pnt P2 = GeomAPI_ProjectPointOnCurve( P1, Line ) ;
4397 if ( P1.IsEqual(P2, Precision::Confusion() ) )
4398 THROW_SALOME_CORBA_EXCEPTION("Points are confused", SALOME::BAD_PARAM);
4400 gp_Vec Vec(P1.X()-P2.X(), P1.Y()-P2.Y(), P1.Z()-P2.Z()) ;
4404 TopoDS_Compound compound;
4406 B.MakeCompound( compound );
4408 for ( i = 0; i < nbtimes2; i++ ) {
4409 for ( j = 0; j < nbtimes1; j++ ) {
4410 DX = i * step * Vec.X() ;
4411 DY = i * step * Vec.Y() ;
4412 DZ = i * step * Vec.Z() ;
4413 myVec.SetCoord( DX, DY, DZ ) ;
4414 theTransformation1.SetTranslation(myVec) ;
4415 theTransformation2.SetRotation(AX1, j*ang*PI180) ;
4416 BRepBuilderAPI_Transform myBRepTransformation1(aShape, theTransformation1, Standard_False) ;
4417 BRepBuilderAPI_Transform myBRepTransformation2(myBRepTransformation1.Shape(), theTransformation2, Standard_False) ;
4418 B.Add( compound, myBRepTransformation2.Shape() );
4422 result = CreateObject(tds) ;
4424 catch(Standard_Failure) {
4425 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeMultiRotation2D", SALOME::BAD_PARAM);
4428 if( CORBA::is_nil(result) ) {
4429 THROW_SALOME_CORBA_EXCEPTION("MakeMultiRotation2D aborted : null result", SALOME::BAD_PARAM);
4431 InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4436 //=================================================================================
4437 // function : MakeCopy()
4438 // purpose : Copy a 3D shape
4439 //=================================================================================
4440 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeCopy( GEOM::GEOM_Shape_ptr Shape)
4441 throw (SALOME::SALOME_Exception)
4443 Unexpect aCatch(SALOME_SalomeException);
4444 GEOM::GEOM_Shape_var result ;
4446 TopoDS_Shape aShape = GetTopoShape(Shape) ;
4447 if(aShape.IsNull() ) {
4448 THROW_SALOME_CORBA_EXCEPTION("Copy aborted : null shape during operation", SALOME::BAD_PARAM);
4450 BRepBuilderAPI_Copy Copy(aShape);
4451 if( Copy.IsDone() ) {
4453 result = CreateObject(tds);
4454 InsertInLabelOneArgument(aShape, Shape, tds, result, myCurrentOCAFDoc) ;
4461 //=================================================================================
4462 // function : MakeMirrorByPlane()
4463 // purpose : build a shape by symmetry of 'myShape' with 'shapePlane' in argument
4464 //=================================================================================
4465 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeMirrorByPlane(GEOM::GEOM_Shape_ptr myShape,
4466 GEOM::GEOM_Shape_ptr shapePlane)
4467 throw (SALOME::SALOME_Exception)
4469 Unexpect aCatch(SALOME_SalomeException);
4470 GEOM::GEOM_Shape_var result ;
4472 TopoDS_Shape aShape = GetTopoShape(myShape) ;
4473 TopoDS_Shape aShapePlane = GetTopoShape(shapePlane) ;
4474 if( aShape.IsNull() || aShapePlane.IsNull() ) {
4475 THROW_SALOME_CORBA_EXCEPTION("Mirror aborted : null shape argument", SALOME::BAD_PARAM);
4479 Handle(Geom_Surface) surf = BRep_Tool::Surface(TopoDS::Face(aShapePlane)) ;
4480 Handle(Geom_Plane) myPlane = Handle(Geom_Plane)::DownCast(surf) ;
4481 const gp_Ax3 pos = myPlane->Position() ;
4482 const gp_Pnt loc = pos.Location() ; /* location of the plane */
4483 const gp_Dir dir = pos.Direction() ; /* Main direction of the plane (Z axis) */
4485 /* plane used for mirroring */
4486 gp_Ax2 pln(loc, dir) ;
4487 gp_Trsf theTransformation ;
4488 theTransformation.SetMirror(pln) ;
4489 BRepBuilderAPI_Transform myBRepTransformation(aShape, theTransformation, Standard_False) ;
4491 tds = myBRepTransformation.Shape() ;
4493 THROW_SALOME_CORBA_EXCEPTION("Mirror aborted", SALOME::BAD_PARAM);
4496 catch(Standard_Failure) {
4497 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeMirrorByPlane", SALOME::BAD_PARAM);
4500 result = CreateObject(tds) ;
4502 /* Insert arguments in ocaf */
4503 GEOM::GEOM_Gen::ListOfIOR_var ListShapes = new GEOM::GEOM_Gen::ListOfIOR;
4504 ListShapes->length(2);
4505 ListShapes[0] = GetStringFromIOR(GEOM::GEOM_Shape::_duplicate(myShape)) ;
4506 ListShapes[1] = GetStringFromIOR(GEOM::GEOM_Shape::_duplicate(shapePlane)) ;
4507 InsertInLabelMoreArguments(tds, result, ListShapes, myCurrentOCAFDoc) ;
4513 //=================================================================================
4514 // function : MakeRotation()
4515 // purpose : Rotation of a 3D shape around an axis
4516 //=================================================================================
4517 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeRotation( GEOM::GEOM_Shape_ptr myShape,
4518 const GEOM::AxisStruct& axis,
4519 CORBA::Double angle)
4520 throw (SALOME::SALOME_Exception)
4522 Unexpect aCatch(SALOME_SalomeException);
4523 GEOM::GEOM_Shape_var result ;
4525 TopoDS_Shape aShape = GetTopoShape(myShape) ;
4526 if( aShape.IsNull() ) {
4527 THROW_SALOME_CORBA_EXCEPTION("Rotation aborted : null shape during operation", SALOME::BAD_PARAM);
4531 gp_Pnt P(axis.x, axis.y, axis.z) ;
4532 gp_Dir D(axis.vx, axis.vy, axis.vz) ;
4535 gp_Trsf theTransformation ;
4536 theTransformation.SetRotation(AX, angle) ;
4537 BRepBuilderAPI_Transform myBRepTransformation(aShape, theTransformation, Standard_False) ;
4538 tds = myBRepTransformation.Shape() ;
4540 catch(Standard_Failure) {
4541 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeRotation", SALOME::BAD_PARAM);
4544 if ( !tds.IsNull() ) {
4545 result = CreateObject(tds) ;
4546 InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4552 //=================================================================================
4553 // function : MakeScaleTransform()
4554 // purpose : Make a shape multipling another by a scale factor
4555 //=================================================================================
4556 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeScaleTransform(GEOM::GEOM_Shape_ptr myShape,
4557 const GEOM::PointStruct& theCenterOfScale,
4558 CORBA::Double factor)
4559 throw (SALOME::SALOME_Exception)
4561 Unexpect aCatch(SALOME_SalomeException);
4562 GEOM::GEOM_Shape_var result ;
4564 TopoDS_Shape aShape = GetTopoShape(myShape) ;
4565 if( aShape.IsNull() ) {
4566 THROW_SALOME_CORBA_EXCEPTION("Scale aborted : null shape during operation", SALOME::BAD_PARAM);
4570 gp_Pnt Pcenter(theCenterOfScale.x, theCenterOfScale.y, theCenterOfScale.z) ;
4571 gp_Trsf theTransformation ;
4572 theTransformation.SetScale(Pcenter, factor) ;
4573 BRepBuilderAPI_Transform myBRepTransformation(aShape, theTransformation, Standard_False) ;
4574 tds = myBRepTransformation.Shape() ;
4576 catch(Standard_Failure) {
4577 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeScaleTransform", SALOME::BAD_PARAM);
4580 if ( !tds.IsNull() ) {
4581 result = CreateObject(tds) ;
4582 InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4588 //=================================================================================
4589 // function : MakeCompound()
4590 // purpose : Make a compound from a list containing one or more shapes
4591 //=================================================================================
4592 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeCompound( const GEOM::GEOM_Gen::ListOfIOR& ListShapes )
4593 throw (SALOME::SALOME_Exception)
4595 Unexpect aCatch(SALOME_SalomeException);
4596 GEOM::GEOM_Shape_var result ;
4598 BRep_Builder aBuilder;
4599 aBuilder.MakeCompound(C) ;
4601 for ( unsigned int i = 0; i < ListShapes.length(); i++) {
4602 GEOM::GEOM_Shape_var aShape = GetIORFromString( ListShapes[i] );
4603 TopoDS_Shape Shape = GetTopoShape(aShape) ;
4604 if( Shape.IsNull() ) {
4605 THROW_SALOME_CORBA_EXCEPTION("Compound aborted : null shape during operation", SALOME::BAD_PARAM);
4607 aBuilder.Add(C, Shape) ;
4611 THROW_SALOME_CORBA_EXCEPTION("Null result : Compound operation aborted", SALOME::BAD_PARAM);
4614 result = CreateObject(C) ;
4615 InsertInLabelMoreArguments(C, result, ListShapes, myCurrentOCAFDoc) ;
4621 //================================================================================
4622 // function : MakeEdge()
4623 // purpose : Make a linear edge with 2 points
4624 //================================================================================
4625 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeEdge(const GEOM::PointStruct& pstruct1,
4626 const GEOM::PointStruct& pstruct2)
4627 throw (SALOME::SALOME_Exception)
4629 Unexpect aCatch(SALOME_SalomeException);
4630 GEOM::GEOM_Shape_var result ;
4634 gp_Pnt P1(pstruct1.x, pstruct1.y, pstruct1.z);
4635 gp_Pnt P2(pstruct2.x, pstruct2.y, pstruct2.z) ;
4636 tds = BRepBuilderAPI_MakeEdge(P1, P2).Shape();
4638 THROW_SALOME_CORBA_EXCEPTION("MakeEdge aborted : null result", SALOME::BAD_PARAM);
4640 catch (Standard_Failure) {
4641 THROW_SALOME_CORBA_EXCEPTION("Exception catched in MakeEdge", SALOME::BAD_PARAM);
4644 result = CreateObject(tds) ;
4645 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
4646 result->ShapeId(entry);
4652 //=================================================================================
4653 // function : MakeWire()
4654 // purpose : Make a wire from a list containing one or more edges or wires that can
4656 //=================================================================================
4657 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeWire( const GEOM::GEOM_Gen::ListOfIOR& ListShapes )
4658 throw (SALOME::SALOME_Exception)
4660 Unexpect aCatch(SALOME_SalomeException);
4661 GEOM::GEOM_Shape_var result ;
4662 BRepBuilderAPI_MakeWire MW ;
4663 TopoDS_Shape tds, Shape ;
4666 for ( unsigned int i = 0; i < ListShapes.length(); i++) {
4667 GEOM::GEOM_Shape_var aShape = GetIORFromString( ListShapes[i] );
4668 Shape = GetTopoShape(aShape) ;
4669 if( Shape.IsNull() ) {
4670 THROW_SALOME_CORBA_EXCEPTION("MakeWire aborted : null shape during operation", SALOME::BAD_PARAM);
4672 if( Shape.ShapeType() == TopAbs_EDGE )
4673 MW.Add( TopoDS::Edge(Shape) ) ;
4674 if (Shape.ShapeType() == TopAbs_WIRE )
4675 MW.Add( TopoDS::Wire(Shape) ) ;
4680 catch(Standard_Failure) {
4681 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeWire", SALOME::BAD_PARAM);
4684 if( tds.IsNull() ) {
4685 THROW_SALOME_CORBA_EXCEPTION("Make Wire operation aborted : null result", SALOME::BAD_PARAM);
4688 result = CreateObject(tds) ;
4689 InsertInLabelMoreArguments(tds, result, ListShapes, myCurrentOCAFDoc) ;
4695 //=================================================================================
4696 // function : MakeRevolution()
4698 //=================================================================================
4699 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeRevolution(GEOM::GEOM_Shape_ptr myShape,
4700 const GEOM::AxisStruct& axis,
4702 throw (SALOME::SALOME_Exception)
4704 Unexpect aCatch(SALOME_SalomeException);
4705 GEOM::GEOM_Shape_var result ;
4707 TopoDS_Shape aShape = GetTopoShape(myShape) ;
4708 if( aShape.IsNull() ) {
4709 THROW_SALOME_CORBA_EXCEPTION("Revolution aborted : null shape", SALOME::BAD_PARAM);
4712 gp_Pnt P(axis.x, axis.y, axis.z) ;
4713 gp_Dir D(axis.vx, axis.vy, axis.vz);
4715 tds = BRepPrimAPI_MakeRevol(aShape, AX, angle);
4717 catch(Standard_Failure) {
4718 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeRevolution", SALOME::BAD_PARAM);
4721 if( tds.IsNull() ) {
4722 THROW_SALOME_CORBA_EXCEPTION("Revolution aborted", SALOME::BAD_PARAM);
4724 result = CreateObject(tds) ;
4725 InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4730 //=================================================================================
4731 // function : MakePipe()
4732 // purpose : Create a shape by sweeping a baseShape along a pathShape
4733 //=================================================================================
4734 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakePipe( GEOM::GEOM_Shape_ptr pathShape,
4735 GEOM::GEOM_Shape_ptr baseShape )
4736 throw (SALOME::SALOME_Exception)
4738 Unexpect aCatch(SALOME_SalomeException);
4739 GEOM::GEOM_Shape_var result ;
4742 TopoDS_Shape pathTds = GetTopoShape(pathShape) ;
4743 TopoDS_Shape baseTds = GetTopoShape(baseShape) ;
4745 if( baseTds.IsNull() || pathTds.IsNull() ) {
4746 THROW_SALOME_CORBA_EXCEPTION("MakePipe aborted : null shape argument", SALOME::BAD_PARAM);
4749 if( pathTds.ShapeType() == TopAbs_WIRE ) {
4750 aWire = TopoDS::Wire(pathTds) ;
4753 if ( pathTds.ShapeType() == TopAbs_EDGE ) {
4754 TopoDS_Edge aEdge = TopoDS::Edge(pathTds) ;
4755 aWire = BRepBuilderAPI_MakeWire(aEdge);
4758 THROW_SALOME_CORBA_EXCEPTION("MakePipe aborted : bad shape type", SALOME::BAD_PARAM);
4763 tds = BRepOffsetAPI_MakePipe(aWire, baseTds) ;
4765 catch(Standard_Failure) {
4766 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakePipe", SALOME::BAD_PARAM);
4769 if ( !IsValid(tds) ) {
4770 THROW_SALOME_CORBA_EXCEPTION("MakePipe aborted : non valid shape result", SALOME::BAD_PARAM);
4773 result = CreateObject(tds) ;
4775 /* Insert arguments in ocaf */
4776 GEOM::GEOM_Gen::ListOfIOR_var ListShapes = new GEOM::GEOM_Gen::ListOfIOR;
4777 ListShapes->length(2);
4778 ListShapes[0] = GetStringFromIOR(GEOM::GEOM_Shape::_duplicate(pathShape)) ;
4779 ListShapes[1] = GetStringFromIOR(GEOM::GEOM_Shape::_duplicate(baseShape)) ;
4780 InsertInLabelMoreArguments(tds, result, ListShapes, myCurrentOCAFDoc) ;
4786 //=================================================================================
4787 // function : MakePrism()
4788 // purpose : uses myShape as base and the vector P1 to P2
4789 //=================================================================================
4790 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakePrism( GEOM::GEOM_Shape_ptr myShape,
4791 const GEOM::PointStruct& P1,
4792 const GEOM::PointStruct& P2 )
4793 throw (SALOME::SALOME_Exception)
4795 Unexpect aCatch(SALOME_SalomeException);
4796 GEOM::GEOM_Shape_var result ;
4798 TopoDS_Shape aShape = GetTopoShape(myShape) ;
4799 if( aShape.IsNull() ) {
4800 THROW_SALOME_CORBA_EXCEPTION("Prism aborted : null shape operation", SALOME::BAD_PARAM);
4804 gp_Vec Vector (P2.x - P1.x, P2.y - P1.y, P2.z - P1.z) ;
4805 tds = BRepPrimAPI_MakePrism(aShape, Vector, Standard_False).Shape() ;
4807 catch(Standard_Failure) {
4808 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakePipe", SALOME::BAD_PARAM);
4811 if ( tds.IsNull() ) {
4812 THROW_SALOME_CORBA_EXCEPTION("Prism aborted", SALOME::BAD_PARAM);
4815 result = CreateObject(tds) ;
4816 InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4822 //=================================================================================
4823 // function : MakeCDG()
4824 // purpose : Create a CDG topology.
4825 //=================================================================================
4826 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeCDG(GEOM::GEOM_Shape_ptr aShape)
4827 throw (SALOME::SALOME_Exception)
4829 Unexpect aCatch(SALOME_SalomeException);
4830 GEOM::GEOM_Shape_var result ;
4832 TopoDS_Shape shape = GetTopoShape(aShape) ;
4833 GProp_GProps System;
4834 gp_Pnt myCenterMass ;
4836 if( shape.IsNull() ) {
4837 THROW_SALOME_CORBA_EXCEPTION("MakeCDG aborted : null shape argument", SALOME::BAD_PARAM);
4841 if ( shape.ShapeType() == TopAbs_VERTEX) {
4842 myCenterMass = BRep_Tool::Pnt(TopoDS::Vertex( shape ));
4844 else if ( shape.ShapeType() == TopAbs_EDGE || shape.ShapeType() == TopAbs_WIRE ) {
4845 BRepGProp::LinearProperties(shape, System);
4846 myCenterMass = System.CentreOfMass() ;
4848 else if ( shape.ShapeType() == TopAbs_FACE || shape.ShapeType() == TopAbs_SHELL ) {
4849 BRepGProp::SurfaceProperties(shape, System);
4850 myCenterMass = System.CentreOfMass() ;
4853 BRepGProp::VolumeProperties(shape, System);
4854 myCenterMass = System.CentreOfMass() ;
4857 tds = BRepBuilderAPI_MakeVertex(myCenterMass).Shape() ;
4859 catch(Standard_Failure) {
4860 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeCDG", SALOME::BAD_PARAM);
4863 if ( tds.IsNull() ) {
4864 THROW_SALOME_CORBA_EXCEPTION("Make CDG aborted : null shape result", SALOME::BAD_PARAM);
4867 result = CreateObject(tds) ;
4868 InsertInLabelOneArgument(shape, aShape, tds, result, myCurrentOCAFDoc) ;
4874 //=================================================================================
4875 // function : Archimede()
4877 //=================================================================================
4878 GEOM::GEOM_Shape_ptr GEOM_Gen_i::Archimede(GEOM::GEOM_Shape_ptr aShape,
4879 CORBA::Double aWeight,
4880 CORBA::Double aWaterDensity,
4881 CORBA::Double aMeshingDeflection)
4882 throw (SALOME::SALOME_Exception)
4884 Unexpect aCatch(SALOME_SalomeException);
4885 GEOM::GEOM_Shape_var result;
4888 if (aWaterDensity != 0.)
4889 cste = aWeight/aWaterDensity;
4891 THROW_SALOME_CORBA_EXCEPTION("Water density is null", SALOME::BAD_PARAM);
4893 TopoDS_Shape shape = GetTopoShape(aShape) ;
4894 if( shape.IsNull() ) {
4895 THROW_SALOME_CORBA_EXCEPTION("Shape is null", SALOME::BAD_PARAM);
4898 gp_Dir direct(0.0,0.0,1.0);
4899 gp_Pnt PosPlan(0.0,0.0,0.0);
4900 Geom_Plane PP (PosPlan,direct);
4901 Handle(Geom_Geometry) G = PP.Copy();
4902 Handle(Geom_Plane) P = Handle(Geom_Plane)::DownCast(G);
4904 gp_Dir Zdirection(0.0,0.0,1.0);
4905 VolumeSection VOL( shape, aMeshingDeflection);
4907 Handle (Geom_RectangularTrimmedSurface) SurfaceTrimmee;
4909 if(Zdirection.IsEqual(direct,Precision::Angular()) == Standard_False) {
4910 VOL.MakeRotation(direct);
4913 VOL.CenterOfGravity();
4914 SurfaceTrimmee = VOL.TrimSurf();
4915 Standard_Real Cote = VOL.Archimede( cste, aMeshingDeflection );
4919 VOL.getZ(Zmin,Zmax);
4920 double volume = VOL.CalculateVolume( Zmax ) * aWaterDensity;
4923 sprintf(msg, "shape sinks to the bottom : Weigth max = %.1f", volume);
4925 THROW_SALOME_CORBA_EXCEPTION(msg, SALOME::BAD_PARAM);
4928 SurfaceTrimmee=VOL.AjustePlan(SurfaceTrimmee,Cote,PosPlan);
4929 if(Zdirection.IsEqual(direct,Precision::Angular()) == Standard_False) {
4930 SurfaceTrimmee=VOL.InvMakeRotation(direct,SurfaceTrimmee);
4933 Standard_Real u1,u2,v1,v2;
4934 SurfaceTrimmee->Bounds(u1,u2,v1,v2);
4935 TopoDS_Face tirant = BRepBuilderAPI_MakeFace(SurfaceTrimmee, u1, u2, v1, v2);
4937 if (tirant.IsNull()) {
4938 THROW_SALOME_CORBA_EXCEPTION("Result is null", SALOME::BAD_PARAM);
4941 result = CreateObject(tirant);
4942 InsertInLabelOneArgument(shape, aShape, tirant, result, myCurrentOCAFDoc) ;
4948 //================================================================================
4949 // function : MakeFillet()
4950 // purpose : Create a cylinder topology
4951 //================================================================================
4952 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeFillet( GEOM::GEOM_Shape_ptr shape,
4953 CORBA::Double radius,
4954 CORBA::Short ShapeType,
4955 const GEOM::GEOM_Shape::ListOfSubShapeID& ListOfID )
4956 throw (SALOME::SALOME_Exception)
4958 Unexpect aCatch(SALOME_SalomeException);
4959 GEOM::GEOM_Shape_var result;
4962 const TopoDS_Shape aShape = GetTopoShape(shape) ;
4963 if( aShape.IsNull() ) {
4964 THROW_SALOME_CORBA_EXCEPTION("Shape is null", SALOME::BAD_PARAM);
4967 BRepFilletAPI_MakeFillet fill(aShape);
4971 if(ListOfID.length() == 0) {
4972 TopExp_Explorer Exp ( aShape, TopAbs_EDGE );
4973 for (Exp; Exp.More(); Exp.Next()) {
4974 TopoDS_Edge E =TopoDS::Edge(Exp.Current());
4977 for (int i = 1;i<=fill.NbContours();i++) {
4978 #if OCC_VERSION_MAJOR >= 5
4979 fill.SetRadius(radius,i,i);
4981 fill.SetRadius(radius,i);
4988 /* case selection */
4989 for ( unsigned int ind = 0; ind < ListOfID.length(); ind++ ) {
4991 if( GetShapeFromIndex( aShape, (TopAbs_ShapeEnum)ShapeType, ListOfID[ind], ss ) ) {
4992 TopoDS_Edge E = TopoDS::Edge(ss) ;
4996 for (int i = 1;i<=fill.NbContours();i++) {
4997 #if OCC_VERSION_MAJOR >= 5
4998 fill.SetRadius(radius,i,i);
5000 fill.SetRadius(radius,i);
5006 catch(Standard_Failure) {
5007 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeFillet", SALOME::BAD_PARAM);
5011 THROW_SALOME_CORBA_EXCEPTION("Make Fillet aborted", SALOME::BAD_PARAM);
5013 result = CreateObject(tds);
5014 InsertInLabelOneArgument(aShape, shape, tds, result, myCurrentOCAFDoc) ;
5020 //================================================================================
5021 // function : MakeChamfer
5022 // purpose : Create a Chamfer topology
5023 //================================================================================
5024 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeChamfer( GEOM::GEOM_Shape_ptr shape,
5027 CORBA::Short ShapeType,
5028 const GEOM::GEOM_Shape::ListOfSubShapeID& ListOfID )
5029 throw (SALOME::SALOME_Exception)
5031 Unexpect aCatch(SALOME_SalomeException);
5032 GEOM::GEOM_Shape_var result;
5035 const TopoDS_Shape aShape = GetTopoShape(shape) ;
5036 if( aShape.IsNull() ) {
5037 THROW_SALOME_CORBA_EXCEPTION("Shape is null", SALOME::BAD_PARAM);
5040 BRepFilletAPI_MakeChamfer MC(aShape);
5044 TopTools_IndexedDataMapOfShapeListOfShape M;
5045 TopExp::MapShapesAndAncestors(aShape,TopAbs_EDGE,TopAbs_FACE,M);
5046 if(ListOfID.length() == 0) {
5047 for (int i = 1;i<=M.Extent();i++) {
5048 TopoDS_Edge E = TopoDS::Edge(M.FindKey(i));
5049 TopoDS_Face F = TopoDS::Face(M.FindFromIndex(i).First());
5050 if (!BRepTools::IsReallyClosed(E, F) && !BRep_Tool::Degenerated(E))
5057 /* case selection */
5058 for ( unsigned int ind = 0; ind < ListOfID.length(); ind++ ) {
5060 if( GetShapeFromIndex( aShape, (TopAbs_ShapeEnum)ShapeType, ListOfID[ind], ss ) ) {
5061 TopoDS_Edge E = TopoDS::Edge( ss ) ;
5062 TopoDS_Face F = TopoDS::Face(M.FindFromKey(E).First());
5063 if (!BRepTools::IsReallyClosed(E, F) && !BRep_Tool::Degenerated(E))
5070 catch(Standard_Failure) {
5071 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeChamfer", SALOME::BAD_PARAM);
5075 THROW_SALOME_CORBA_EXCEPTION("Make Chamfer aborted", SALOME::BAD_PARAM);
5077 result = CreateObject(tds);
5078 InsertInLabelOneArgument(aShape, shape, tds, result, myCurrentOCAFDoc) ;
5083 //=================================================================================
5084 // function : CheckShape()
5086 //=================================================================================
5087 CORBA::Boolean GEOM_Gen_i::CheckShape(GEOM::GEOM_Shape_ptr shape)
5088 throw (SALOME::SALOME_Exception)
5090 Unexpect aCatch(SALOME_SalomeException);
5091 TopoDS_Shape S = GetTopoShape(shape) ;
5093 THROW_SALOME_CORBA_EXCEPTION("Shape is null", SALOME::BAD_PARAM);
5096 BRepCheck_Analyzer ana(S,false);
5103 //=================================================================================
5104 // function : MakePlacedBox()
5106 //=================================================================================
5107 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakePlacedBox(CORBA::Double x1, CORBA::Double y1, CORBA::Double z1,
5108 CORBA::Double delta1, CORBA::Double delta2, CORBA::Double delta3)
5109 throw (SALOME::SALOME_Exception)
5111 Unexpect aCatch(SALOME_SalomeException);
5112 GEOM::GEOM_Shape_var result ;
5115 CORBA::Double x2, y2, z2 ;
5122 gp_Pnt P1(x1,y1,z1);
5123 gp_Pnt P2(x2,y2,z2);
5125 tds = BRepPrimAPI_MakeBox(P1,P2).Shape();
5127 catch(Standard_Failure) {
5128 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeBox", SALOME::BAD_PARAM);
5132 THROW_SALOME_CORBA_EXCEPTION("Make Box aborted : null shape", SALOME::BAD_PARAM);
5135 result = CreateObject(tds);
5136 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
5137 result->ShapeId(entry) ;
5142 //=================================================================================
5143 // function : MakePanel()
5145 //=================================================================================
5146 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakePanel(GEOM::GEOM_Shape_ptr shape,
5147 CORBA::Short directiontype,
5148 CORBA::Double delta)
5149 throw (SALOME::SALOME_Exception)
5151 Unexpect aCatch(SALOME_SalomeException);
5152 GEOM::GEOM_Shape_var result ;
5154 TopoDS_Shape aShape = GetTopoShape(shape) ;
5156 Standard_Real axmin,aymin,azmin,axmax,aymax,azmax ;
5157 GEOM::PointStruct pstruct1, pstruct2, pstruct3, pstruct4 ;
5159 if(aShape.IsNull() ) {
5160 THROW_SALOME_CORBA_EXCEPTION("MakePanel aborted : null shape during operation", SALOME::BAD_PARAM);
5164 BRepBndLib::Add(aShape,B);
5166 B.Get(axmin,aymin,azmin,axmax,aymax,azmax);
5168 switch (directiontype)
5171 pstruct1 = MakePointStruct( delta, aymin, azmin ) ;
5172 pstruct2 = MakePointStruct( delta, aymin, azmax ) ;
5173 pstruct3 = MakePointStruct( delta, aymax, azmax ) ;
5174 pstruct4 = MakePointStruct( delta, aymax, azmin ) ;
5177 pstruct1 = MakePointStruct( axmin, delta, azmin ) ;
5178 pstruct2 = MakePointStruct( axmin, delta, azmax ) ;
5179 pstruct3 = MakePointStruct( axmax, delta, azmax ) ;
5180 pstruct4 = MakePointStruct( axmax, delta, azmin ) ;
5183 pstruct1 = MakePointStruct( axmin, aymin, delta ) ;
5184 pstruct2 = MakePointStruct( axmin, aymax, delta ) ;
5185 pstruct3 = MakePointStruct( axmax, aymax, delta ) ;
5186 pstruct4 = MakePointStruct( axmax, aymin, delta ) ;
5192 GEOM::GEOM_Shape_ptr Edge1 = MakeEdge(pstruct1, pstruct2);
5193 GEOM::GEOM_Shape_ptr Edge2 = MakeEdge(pstruct2, pstruct3);
5194 GEOM::GEOM_Shape_ptr Edge3 = MakeEdge(pstruct3, pstruct4);
5195 GEOM::GEOM_Shape_ptr Edge4 = MakeEdge(pstruct4, pstruct1);
5197 GEOM::GEOM_Gen::ListOfIOR_var aList = new GEOM::GEOM_Gen::ListOfIOR;
5199 aList[0]=CORBA::string_dup(Edge1->Name());
5200 aList[1]=CORBA::string_dup(Edge2->Name());
5201 aList[2]=CORBA::string_dup(Edge3->Name());
5202 aList[3]=CORBA::string_dup(Edge4->Name());
5204 GEOM::GEOM_Shape_ptr aWire = MakeWire( aList );
5205 GEOM::GEOM_Shape_ptr aFace = MakeFace( aWire, true ) ;
5206 tds = GetTopoShape(aFace);
5209 catch(Standard_Failure) {
5210 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakePanel", SALOME::BAD_PARAM);
5214 THROW_SALOME_CORBA_EXCEPTION("Make PanelsPartition aborted : null shape", SALOME::BAD_PARAM);
5217 result = CreateObject(tds);
5218 const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
5219 result->ShapeId(entry) ;
5225 void GEOM_Gen_i::ExportIGES(const char* filename,GEOM::GEOM_Shape_ptr theShape)
5226 throw (SALOME::SALOME_Exception)
5228 Unexpect aCatch(SALOME_SalomeException);
5229 if (theShape->_is_nil())
5231 THROW_SALOME_CORBA_EXCEPTION("Export IGES aborted", SALOME::BAD_PARAM);
5233 TopoDS_Shape tds = GetTopoShape(theShape);
5236 THROW_SALOME_CORBA_EXCEPTION("Export IGES aborted", SALOME::BAD_PARAM);
5240 //VRV: OCC 4.0 migration
5241 IGESControl_Controller::Init();
5242 IGESControl_Writer ICW (Interface_Static::CVal("XSTEP.iges.unit"),
5243 Interface_Static::IVal("XSTEP.iges.writebrep.mode"));
5244 //VRV: OCC 4.0 migration
5248 char * aname = strdup(filename);
5249 Standard_Boolean result = ICW.Write( aname );
5252 catch(Standard_Failure)
5254 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ExportIGES", SALOME::BAD_PARAM);
5258 void GEOM_Gen_i::ExportBREP(const char* filename,GEOM::GEOM_Shape_ptr theShape)
5259 throw (SALOME::SALOME_Exception)
5261 Unexpect aCatch(SALOME_SalomeException);
5262 if (theShape->_is_nil())
5264 THROW_SALOME_CORBA_EXCEPTION("Export BRep aborted", SALOME::BAD_PARAM);
5266 TopoDS_Shape tds = GetTopoShape(theShape);
5269 THROW_SALOME_CORBA_EXCEPTION("Export BRep aborted", SALOME::BAD_PARAM);
5273 char * aname = strdup(filename);
5274 Standard_Boolean result = BRepTools::Write(tds,aname);
5277 catch(Standard_Failure)
5279 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ExportBREP", SALOME::BAD_PARAM);
5283 void GEOM_Gen_i::ExportSTEP(const char* filename,GEOM::GEOM_Shape_ptr theShape)
5284 throw (SALOME::SALOME_Exception)
5286 Unexpect aCatch(SALOME_SalomeException);
5287 if (theShape->_is_nil())
5289 THROW_SALOME_CORBA_EXCEPTION("Export STEP aborted", SALOME::BAD_PARAM);
5291 TopoDS_Shape tds = GetTopoShape(theShape);
5294 THROW_SALOME_CORBA_EXCEPTION("Export STEP aborted", SALOME::BAD_PARAM);
5298 IFSelect_ReturnStatus status ;
5299 //VRV: OCC 4.0 migration
5300 STEPControl_Writer aWriter;
5301 status = aWriter.Transfer( tds, STEPControl_ManifoldSolidBrep ) ;
5302 //VRV: OCC 4.0 migration
5303 if ( status == IFSelect_RetDone )
5305 char * aname = strdup(filename);
5306 status = aWriter.Write( aname ) ;
5310 catch(Standard_Failure)
5312 THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ExportBREP", SALOME::BAD_PARAM);
5317 //=====================================================================================
5319 //=====================================================================================
5322 PortableServer::ObjectId * GEOMEngine_factory(CORBA::ORB_ptr orb,
5323 PortableServer::POA_ptr poa,
5324 PortableServer::ObjectId * contId,
5325 const char *instanceName,
5326 const char * interfaceName)
5328 GEOM_Gen_i * myGEOM_Gen_i = new GEOM_Gen_i(orb, poa, contId, instanceName, interfaceName);
5329 myGEOM_Gen_i->register_name("/myGEOM_Gen"); // NRI : 11/07/2002 : Add for Supervision example
5330 return myGEOM_Gen_i->getId() ;