src/GEOM/GEOM_Gen_i.cc

   1 //  GEOM GEOM : implementaion of GEOM_Gen.idl and GEOM_Shape.idl
   2 //
   3 //  Copyright (C) 2003  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   4 //  CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   5 //
   6 //  This library is free software; you can redistribute it and/or
   7 //  modify it under the terms of the GNU Lesser General Public
   8 //  License as published by the Free Software Foundation; either
   9 //  version 2.1 of the License.
  10 //
  11 //  This library is distributed in the hope that it will be useful,
  12 //  but WITHOUT ANY WARRANTY; without even the implied warranty of
  13 //  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14 //  Lesser General Public License for more details.
  15 //
  16 //  You should have received a copy of the GNU Lesser General Public
  17 //  License along with this library; if not, write to the Free Software
  18 //  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  19 //
  20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  21 //
  22 //
  23 //
  24 //  File   : GEOM_GEN_i.cc file
  25 //  Author : Lucien PIGNOLONI
  26 //  Module : GEOM
  27 //  $Header$
  28
  29 using namespace std;
  30 #include "GEOM_Gen_i.hh"
  31
  32 #include "Partition_Spliter.hxx"
  33 #include "Archimede_VolumeSection.hxx"
  34 #include "Sketcher_Profile.hxx"
  35
  36 #include "Utils_CorbaException.hxx"
  37 #include "utilities.h"
  38
  39 #include <stdio.h>
  40
  41 // Cascade headers
  42
  43 #include <Standard_Failure.hxx>
  44
  45 #include <gp_Circ.hxx>
  46 #include <gp_Pln.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>
  60
  61 #include <BRepCheck_Analyzer.hxx>
  62 #if OCC_VERSION_MAJOR >= 5
  63 #include <BRepAlgo.hxx>
  64 #else
  65 #include <BRepAlgoAPI.hxx>
  66 #endif
  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>
  76
  77 #include <BRepBuilderAPI_MakeVertex.hxx>
  78 #include <BRepBuilderAPI_MakeEdge.hxx>
  79 #include <BRepBuilderAPI_MakeFace.hxx>
  80
  81 #include <BRepLib.hxx>
  82 #include <BRepBndLib.hxx>
  83 #include <Bnd_Box.hxx>
  84
  85 #include <BRepBuilderAPI_MakeShell.hxx>
  86 #include <BRepPrim_Builder.hxx>
  87 #include <BRepBuilderAPI_MakeSolid.hxx>
  88 #include <BRepClass3d_SolidClassifier.hxx>
  89
  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>
 104
 105 #include <GeomAPI_ProjectPointOnCurve.hxx>
 106
 107 #include <BRepGProp.hxx>
 108 #include <GProp_GProps.hxx>
 109 #include <Precision.hxx>
 110
 111 //VRV: OCC 4.0 migration
 112 #include <STEPControl_Reader.hxx>
 113 #include <IGESControl_Reader.hxx>
 114 //VRV: OCC 4.0 migration
 115
 116 #include <IFSelect_ReturnStatus.hxx>
 117 #include <TColStd_HSequenceOfTransient.hxx>
 118
 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
 126
 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>
 133
 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>
 143
 144 #include <IGESData_IGESEntity.hxx>
 145
 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>
 156
 157 #include "SALOMEDS_Tool.hxx"
 158 #include "GEOMDS_Commands.hxx"
 159 #include "GEOMDS_Explorer.hxx"
 160
 161 #include CORBA_SERVER_HEADER(SALOMEDS_Attributes)
 162
 163 #include "Utils_ExceptHandlers.hxx"
 164
 165 Standard_EXPORT static Standard_Boolean IsValid(const TopoDS_Shape& S) {
 166 #if OCC_VERSION_MAJOR >= 5
 167   return BRepAlgo::IsValid(S);
 168 #else
 169   return BRepAlgoAPI::IsValid(S);
 170 #endif
 171 }
 172
 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)
 183 {
 184   _thisObj = this ;
 185   _id = _poa->activate_object(_thisObj);
 186   // SCRUTE(this)
 187   name_service = new SALOME_NamingService(_orb);
 188   myOCAFApp    = new GEOMDS_Application();
 189   myStudyID    = -1;
 190   GetCurrentStudy(0);//for correct work of SuperVisor
 191 }
 192
 193
 194
 195 //============================================================================
 196 // function : ~GEOM_Gen_i()
 197 // purpose  : destructor
 198 //============================================================================
 199 GEOM_Gen_i::~GEOM_Gen_i() {
 200   delete name_service;
 201 }
 202
 203
 204 //============================================================================
 205 // function : IORToLocalPersistentID()
 206 // purpose  :
 207 //============================================================================
 208 char* GEOM_Gen_i::IORToLocalPersistentID(SALOMEDS::SObject_ptr theSObject,
 209                                          const char* IORString,
 210                                          CORBA::Boolean isMultiFile,
 211                                          CORBA::Boolean isASCII)
 212 {
 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();
 216   }
 217   return 0;
 218 }
 219
 220
 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)
 231 {
 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;
 236
 237   /* For a GEOM::GEOM_Shape the pers_ref is the Entry in the OCAF document */
 238   TCollection_ExtendedString MainIOR;
 239   TDF_Label Lab;
 240   TDF_Tool::Label(aDoc->GetData(), aPersRefString, Lab );
 241
 242   if (Lab.IsNull()) {
 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
 246   }
 247
 248   Handle(TNaming_NamedShape) NS;
 249   Lab.FindAttribute( TNaming_NamedShape::GetID(), NS );
 250   TopoDS_Shape S = TNaming_Tool::GetShape(NS);
 251
 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() ) ;
 255
 256   if ( GC.HasIOR(Lab) ) { /* shape already created/loaded */
 257     return 0 ;
 258   }
 259
 260   /******************* Dependent object (not a main shape) *********************/
 261   if( GC.IsDependentShape(Lab) ) {
 262
 263     TDF_Label mainLabel ;
 264     Standard_Boolean mainShapeOk = GC.GetMainShapeLabel(Lab, mainLabel) ;
 265
 266     /* Main object not yet loaded we load/create it */
 267     if( !GC.HasIOR(mainLabel) ) {
 268
 269       TCollection_AsciiString entry;
 270       TDF_Tool::Entry(mainLabel,entry);
 271       CORBA::String_var ent = CORBA::string_dup(entry.ToCString());
 272
 273       /* Create the main object recursively */
 274       MainIOR = LocalPersistentIDToIOR(theSObject, ent, isMultiFile, isASCII) ;
 275     } else {
 276       GC.ReturnNameIOR( mainLabel, MainIOR );
 277     }
 278
 279     result->MainName( TCollection_AsciiString(MainIOR).ToCString() ) ;
 280     result->IsMainShape(false) ;
 281     result->ShapeId(aPersRefString);
 282
 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);
 289
 290     GEOM::GEOM_Shape::ListOfSubShapeID_var ListOfID = new GEOM::GEOM_Shape::ListOfSubShapeID;
 291
 292     if(S.ShapeType() != TopAbs_COMPOUND) {
 293       /* to set the index of a unique sub shape (Explode All ONLY for the moment !) */
 294       ListOfID->length(1);
 295       int index = 1;
 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) ) {
 301             ListOfID[0] = index;
 302             break;
 303           }
 304           index++ ;
 305         }
 306       }
 307       result->Index(ListOfID) ;
 308       return result->Name();
 309     }
 310     else {
 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;
 315       TopoDS_Iterator it ;
 316       TopExp_Explorer exp ;
 317       TopAbs_ShapeEnum subType ;
 318
 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() ) ;
 323       }
 324
 325       ListOfID->length( SS.Length() ) ;
 326       int j, k ;  /* in TopTools_SequenceOfShape index start at 1 */
 327
 328       for( k=1; k<=SS.Length(); k++ ) {
 329         j = 1 ;
 330         for( exp.Init( mainTopo, subType ); exp.More(); exp.Next() ) {
 331           if( exp.Current().IsSame( SS.Value(k) ) ) {
 332             ListOfID[k-1] = j ;
 333           }
 334           j++ ;
 335         }
 336       }
 337       result->Index(ListOfID) ;
 338       return result->Name();
 339     }
 340
 341   }
 342   /******************* Independent object (not a sub shape) *********************/
 343   else {
 344     result->IsMainShape(true) ;
 345     if( !GC.AddIORNameAttribute(Lab, result->Name() ) )  {
 346       MESSAGE("in LocalPersistentIDToIOR, NAME/IOR attribute already exist." << endl ) ;
 347     }
 348     Handle(TNaming_NamedShape) NamedShape ;
 349     bool notTested = Lab.FindAttribute(TNaming_NamedShape::GetID(), NamedShape) ;
 350     result->ShapeId(aPersRefString);
 351     return result->Name();
 352   }
 353 }
 354
 355 //============================================================================
 356 // function : CanPublishInStudy
 357 // purpose  :
 358 //============================================================================
 359 bool GEOM_Gen_i::CanPublishInStudy(CORBA::Object_ptr theIOR)
 360 {
 361   GEOM::GEOM_Shape_var aShape = GEOM::GEOM_Shape::_narrow(theIOR);
 362   return !(aShape->_is_nil());
 363 }
 364
 365
 366 //============================================================================
 367 // function : PublishInStudy
 368 // purpose  :
 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)
 374 {
 375   Unexpect aCatch(SALOME_SalomeException);
 376   SALOMEDS::SObject_var aResultSO;
 377   if(CORBA::is_nil(theObject)) return aResultSO;
 378
 379   GEOM::GEOM_Shape_var aShape = GEOM::GEOM_Shape::_narrow(theObject);
 380   if(aShape->_is_nil()) return aResultSO;
 381
 382   if(theStudy->_is_nil()) return aResultSO;
 383
 384   SALOMEDS::GenericAttribute_var anAttr;
 385   SALOMEDS::StudyBuilder_var     aStudyBuilder = theStudy->NewBuilder();
 386
 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());
 396   }
 397   if (aFather->_is_nil()) return aResultSO;
 398
 399   if (CORBA::is_nil(theSObject)) {
 400     aResultSO = aStudyBuilder->NewObject(aFather);
 401   } else {
 402     if (!theSObject->ReferencedObject(aResultSO))
 403       THROW_SALOME_CORBA_EXCEPTION("Publish in study supervision graph error",SALOME::BAD_PARAM);
 404   }
 405   anAttr = aStudyBuilder->FindOrCreateAttribute(aResultSO, "AttributeIOR");
 406   SALOMEDS::AttributeIOR_var anIOR = SALOMEDS::AttributeIOR::_narrow(anAttr);
 407   anIOR->SetValue(aShape->Name());
 408
 409   anAttr = aStudyBuilder->FindOrCreateAttribute(aResultSO, "AttributePixMap");
 410   SALOMEDS::AttributePixMap_var aPixmap = SALOMEDS::AttributePixMap::_narrow(anAttr);
 411   TCollection_AsciiString aShapeName("Shape_");
 412
 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_";
 437   }
 438   if (strlen(theName) == 0) aShapeName += TCollection_AsciiString(aResultSO->Tag());
 439   else aShapeName = TCollection_AsciiString((char*)theName);
 440
 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());
 445
 446   //Add a reference to published object
 447 //    aStudyBuilder->Addreference(theObject, aResultSO);
 448   return aResultSO._retn();
 449 }
 450
 451
 452 //============================================================================
 453 // function : Save()
 454 // purpose  : save OCAF/Geom document
 455 //============================================================================
 456 SALOMEDS::TMPFile* GEOM_Gen_i::Save(SALOMEDS::SComponent_ptr theComponent,
 457                                     const char* theURL,
 458                                     bool isMultiFile) {
 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;
 464   aSeq->length(1);
 465   // Prepare a file name to open
 466   TCollection_AsciiString aNameWithExt("");
 467   if (isMultiFile)
 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();
 481 }
 482
 483 SALOMEDS::TMPFile* GEOM_Gen_i::SaveASCII(SALOMEDS::SComponent_ptr theComponent,
 484                                          const char* theURL,
 485                                          bool isMultiFile) {
 486   SALOMEDS::TMPFile_var aStreamFile = Save(theComponent, theURL, isMultiFile);
 487   return aStreamFile._retn();
 488 }
 489
 490
 491 CORBA::Boolean GEOM_Gen_i::Load(SALOMEDS::SComponent_ptr theComponent,
 492                                 const SALOMEDS::TMPFile& theStream,
 493                                 const char* theURL,
 494                                 bool isMultiFile) {
 495
 496   if (theStream.length() <= 9) {
 497     MESSAGE("The TMPFile is too short : " << theStream.length() << " bytes ");
 498     return false;
 499   }
 500
 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,
 505                                                                        aTmpDir.ToCString(),
 506                                                                        isMultiFile);
 507
 508   // Prepare a file name to open
 509   TCollection_AsciiString aNameWithExt("");
 510   if (isMultiFile)
 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;
 514
 515   // Open document
 516   if (myOCAFApp->Open(aFullName, myCurrentOCAFDoc) != CDF_RS_OK) return false;
 517
 518   // Remove the created file and tmp directory
 519   if (!isMultiFile) SALOMEDS_Tool::RemoveTemporaryFiles(aTmpDir.ToCString(), aSeq.in(), true);
 520
 521   SALOMEDS::Study_var Study = theComponent->GetStudy();
 522   TCollection_AsciiString name( Study->Name() );
 523
 524   int StudyID = Study->StudyId();
 525   myStudyIDToDoc.Bind( StudyID, myCurrentOCAFDoc );
 526   myStudyID = StudyID;
 527
 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());
 532
 533   return true;
 534 }
 535
 536 CORBA::Boolean GEOM_Gen_i::LoadASCII(SALOMEDS::SComponent_ptr theComponent,
 537                                      const SALOMEDS::TMPFile& theStream,
 538                                      const char* theURL,
 539                                      bool isMultiFile) {
 540   return Load(theComponent, theStream, theURL, isMultiFile);
 541 }
 542
 543 //  //============================================================================
 544 //  // function : Save()
 545 //  // purpose  : save OCAF/Geom document
 546 //  //============================================================================
 547 //  void GEOM_Gen_i::Save(const char *IORSComponent, const char *aUrlOfFile)
 548 //  {
 549
 550 //    TCollection_ExtendedString path((char*)aUrlOfFile);
 551 //    TCollection_ExtendedString pathWithExt = path + TCollection_ExtendedString(".sgd");
 552 //    myOCAFApp->SaveAs(myCurrentOCAFDoc,pathWithExt);
 553 //  }
 554
 555
 556 //  //============================================================================
 557 //  // function : Load()
 558 //  // purpose  : Load OCAF/Geom document
 559 //  //============================================================================
 560 //  void GEOM_Gen_i::Load(const char *IORSComponent, const char *aUrlOfFile)
 561 //  {
 562
 563 //    TCollection_ExtendedString path((char*)aUrlOfFile);
 564 //    TCollection_ExtendedString pathWithExt = path + TCollection_ExtendedString(".sgd");
 565
 566 //    myOCAFApp->Open(pathWithExt,myCurrentOCAFDoc);
 567
 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() );
 571
 572 //    int StudyID = Study->StudyId();
 573 //    myStudyIDToDoc.Bind( StudyID, myCurrentOCAFDoc );
 574 //    myStudyID = StudyID;
 575
 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());
 580
 581 //    return ;
 582 //  }
 583
 584
 585 //============================================================================
 586 // function : Close()
 587 // purpose  :
 588 //============================================================================
 589 void GEOM_Gen_i::Close(SALOMEDS::SComponent_ptr theComponent)
 590 {
 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");
 601 //      }
 602     myOCAFApp->Close(aDoc);
 603     myStudyIDToDoc.UnBind(anID); // remove document from GEOM documents data map
 604     }
 605 }
 606
 607 //============================================================================
 608 // function : CanCopy()
 609 // purpose  :
 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;
 618   return true;
 619 }
 620
 621 //============================================================================
 622 // function : CopyFrom()
 623 // purpose  :
 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;
 628
 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();
 635
 636   GetCurrentStudy(theObject->GetStudy()->StudyId());
 637
 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);
 647
 648   aStreamFile = new SALOMEDS::TMPFile(aSize, aSize, (CORBA::Octet*)aBuffer, 1);
 649
 650   // Assign an ID = 1 the the type GEOM_Shape
 651   theObjectID = 1;
 652
 653   // Return created TMPFile
 654   return aStreamFile._retn();
 655 }
 656
 657 //============================================================================
 658 // function : CanPaste()
 659 // purpose  :
 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;
 666   return true;
 667 }
 668
 669 //============================================================================
 670 // function : PasteInto()
 671 // purpose  :
 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();
 679
 680   // Retrieve a TopoDS_Shape from byte stream
 681   TopoDS_Shape aTopology;
 682   istrstream aStreamedBrep((char*) &theStream[0], theStream.length());
 683   BRep_Builder aBuilder;
 684   try {
 685     BRepTools::Read(aTopology, aStreamedBrep, aBuilder);
 686   } catch (Standard_Failure) {
 687     MESSAGE("GEOM_Gen_i::PasteInto exception");
 688     return false;
 689   }
 690
 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) ;
 696
 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());
 706
 707   // Return the created in the Study SObject
 708   return aNewSO._retn();
 709 }
 710
 711 //============================================================================
 712 // function : ComponentDataType()
 713 // purpose  :
 714 //============================================================================
 715 char* GEOM_Gen_i::ComponentDataType()
 716 {
 717   return CORBA::string_dup("GEOM");
 718 }
 719
 720 //============================================================================
 721 // function : register()
 722 // purpose  : register 'name' in 'name_service'
 723 //============================================================================
 724 void GEOM_Gen_i::register_name(char * name)
 725 {
 726   GEOM::GEOM_Gen_ptr g = GEOM::GEOM_Gen::_narrow(POA_GEOM::GEOM_Gen::_this());
 727   name_service->Register(g, name);
 728 }
 729
 730
 731
 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 )
 739 {
 740   int nbShapes = listShapes.length() ;
 741   if( nbShapes < 1)
 742     return 0 ;
 743
 744   for(int i=0; i<nbShapes; i++) {
 745     GEOM::GEOM_Shape_var aGeomShape = listShapes[i] ;
 746     TopoDS_Shape aShape = GetTopoShape(aGeomShape) ;
 747     SS.Append(aShape) ;
 748   }
 749   return nbShapes ;
 750 }
 751
 752
 753
 754
 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)
 761 {
 762   TopoDS_Shape tds ;
 763
 764   TDF_Label lab ;
 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) ;
 769
 770   if( !res ) {
 771     return tds ; /* a null shape is returned */
 772   }
 773   else {
 774     return TNaming_Tool::GetShape(NamedShape) ;
 775   }
 776 }
 777
 778
 779
 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) ;
 786   return ret ;
 787 }
 788
 789
 790
 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()) ;
 801   }
 802   return shapeIOR._retn() ;
 803 }
 804
 805
 806
 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)
 813 {
 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);
 817
 818   TCollection_AsciiString entry;
 819   TDF_Tool::Entry(Lab,entry);
 820   const char *ent = entry.ToCString() ;
 821   return ent ;
 822 }
 823
 824
 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,
 832                                                       const char *nameIor,
 833                                                       GEOM::GEOM_Shape_ptr mainshape_ptr,
 834                                                       Handle(TDocStd_Document) OCAFDoc )
 835 {
 836   GEOMDS_Commands GC(OCAFDoc->Main());
 837   /* add attributs S and nameIor in a new label */
 838
 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);
 842
 843   /* add attributs : S, nameIor and ref to main */
 844   TDF_Label Lab = GC.AddDependentShape(S, (char*)nameIor, mainRefLab);
 845
 846   TCollection_AsciiString entry;
 847   TDF_Tool::Entry(Lab, entry);
 848   const char *ent = entry.ToCString() ;
 849   return ent ;
 850 }
 851
 852
 853 //=================================================================================
 854 // function : InsertInLabelOneArgument()
 855 // purpose  :
 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)
 862 {
 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() ) ;
 870
 871   /* Create a new label */
 872   TDF_Label NewLab = Lab.NewChild();
 873   TCollection_ExtendedString Value("Arguments");
 874   TDataStd_Name::Set(NewLab,Value);
 875
 876   TDF_Label NewLab1 = NewLab.NewChild();
 877   TDF_Label RefLab;
 878   TDF_Tool::Label(OCAFDoc->GetData(), shape_ptr->ShapeId(), RefLab);
 879   TDF_Reference::Set(NewLab1, RefLab);
 880 }
 881
 882
 883 //=================================================================================
 884 // function : InsertInLabelMoreArguments()
 885 // purpose  :
 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)
 891 {
 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);
 898
 899   /* Create a new label */
 900   TDF_Label NewLab = Lab.NewChild();
 901   TCollection_ExtendedString Value("Arguments");
 902   TDataStd_Name::Set(NewLab, Value);
 903
 904   for (unsigned int ind = 0; ind < ListShapes.length(); ind++) {
 905
 906     TDF_Label NewLab1 = NewLab.NewChild();
 907     GEOM::GEOM_Shape_var aShape = GetIORFromString( ListShapes[ind] );
 908
 909     TDF_Label RefLab;
 910     TDF_Tool::Label(OCAFDoc->GetData(), aShape->ShapeId(), RefLab);
 911     TDF_Reference::Set(NewLab1, RefLab);
 912   }
 913   result->ShapeId(entry.ToCString());
 914 }
 915
 916
 917
 918 //=================================================================================
 919 // function: NbLabels()
 920 // purpose :
 921 //=================================================================================
 922 CORBA::Short GEOM_Gen_i::NbLabels()
 923 {
 924   TDF_ChildIterator ChildIterator(myCurrentOCAFDoc->Main());
 925   unsigned int i = 1;
 926
 927   while (ChildIterator.More()) {
 928     i++;
 929     ChildIterator.Next();
 930   }
 931   return i;
 932   //  return TDF_Tool::NbLabels( myCurrentOCAFDoc->Main() );
 933 }
 934
 935
 936
 937 //=================================================================================
 938 // function: GetCurrentStudy()
 939 // purpose : Finds or creates the geom/OCAF document corresponding to the index
 940 // 'StudyID'
 941 //=================================================================================
 942 void GEOM_Gen_i::GetCurrentStudy(CORBA::Long StudyID)
 943 {
 944   /* If StudyID is known we link myCurrentOCAFDoc to it */
 945   if (myStudyIDToDoc.IsBound(StudyID)) {
 946     myCurrentOCAFDoc =  Handle(TDocStd_Document)::DownCast(myStudyIDToDoc(StudyID));
 947   }
 948   /* Create a new OCAFDoc and link it to 'StudyID' argument */
 949   else {
 950     myOCAFApp->NewDocument("SALOME_GEOM",myCurrentOCAFDoc);
 951     myStudyIDToDoc.Bind(StudyID,myCurrentOCAFDoc);
 952   }
 953   myStudyID = StudyID;
 954 }
 955
 956
 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)
 962 {
 963 //   if ( tds.ShapeType() == TopAbs_COMPOUND ) {
 964 //     TopoDS_Iterator itr(tds);
 965 //     TopoDS_Shape res;
 966 //     int i = 0;
 967 //     while (itr.More()) {
 968 //       i++;
 969 //       res = itr.Value();
 970 //       itr.Next();
 971 //     }
 972
 973 //     if ( i == 1 )
 974 //       tds = res;
 975 //   }
 976
 977   GEOM::shape_type st = GEOM::shape_type(tds.ShapeType()) ; /* casting */
 978
 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;
 982   index->length(0);
 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());
 985
 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( "" );
 990   return shape;
 991 }
 992
 993 //=======================================================================
 994 //function : CreateSubObject
 995 //purpose  :
 996 //=======================================================================
 997
 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)
1001 {
1002   GEOM::shape_type st = GEOM::shape_type(SubShape.ShapeType()) ; /* casting */
1003
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());
1009
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  */
1019   const char *entry =
1020     InsertInLabelDependentShape(SubShape, shape->Name(), MainShape, myCurrentOCAFDoc) ;
1021   shape->ShapeId( entry ) ;
1022
1023   return shape;
1024 }
1025
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)
1035 {
1036   Unexpect aCatch(SALOME_SalomeException);
1037   BRepTools_Quilt Glue;
1038   aCompoundOfShells.Nullify() ;
1039
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 */
1044       Glue.Add(F);
1045     }
1046   }
1047
1048   /* Use specif method to calculate the compound of shells */
1049   aCompoundOfShells = Glue.Shells();
1050
1051   if( aCompoundOfShells.ShapeType() != TopAbs_COMPOUND ) {
1052     THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::SuppressFacesGlue() : not a GEOM::COMPOUND", SALOME::BAD_PARAM);
1053   }
1054
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() )
1058     numberOfShell++ ;
1059
1060   return numberOfShell ;
1061 }
1062
1063
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)
1070 {
1071   if( mainShape.IsNull() || subshape.IsNull() )
1072     return -1 ;
1073
1074   int index = 1;
1075   if (subshape.ShapeType() == TopAbs_COMPOUND)
1076     {
1077       TopoDS_Iterator it;
1078       TopTools_ListOfShape CL;
1079       CL.Append( mainShape );
1080       TopTools_ListIteratorOfListOfShape itC;
1081       for (itC.Initialize( CL ); itC.More(); itC.Next())
1082         {
1083           for (it.Initialize( itC.Value() );  it.More(); it.Next())
1084             {
1085               if ( it.Value().ShapeType() == TopAbs_COMPOUND)
1086                 {
1087                   if (it.Value().IsSame(subshape))
1088                     return index;
1089                   else
1090                     index++;
1091                   CL.Append( it.Value() );
1092                 }
1093             }
1094         }
1095     }
1096   else
1097     {
1098       TopExp_Explorer Exp ( mainShape,  subshape.ShapeType() );
1099       TopTools_MapOfShape M;
1100       while ( Exp.More() )
1101         {
1102           if ( M.Add(Exp.Current()) )
1103             {
1104               if ( Exp.Current().IsSame(subshape) )
1105                 return index;
1106               index++;
1107             }
1108           Exp.Next();
1109         }
1110     }
1111   return -1;
1112 }
1113
1114
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)
1123 {
1124
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);
1130   }
1131
1132   /* put faces of subShell in a Map of faces */
1133   int j = 0 ;
1134   TopTools_MapOfShape mapFaces ;
1135   for( TopExp_Explorer Exp1( subShell, TopAbs_FACE );  Exp1.More(); Exp1.Next() ) {
1136     mapFaces.Add(Exp1.Current() ) ;
1137     j++ ;
1138   }
1139
1140   if( j<1 )
1141     THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::IndexOfFacesOfSubShell() : no faces in sub shell", SALOME::BAD_PARAM);
1142
1143   /* Find index of each face of subshell in the main topology and put its index in ListOfID */
1144   int size = 0 ;
1145   for ( TopExp_Explorer Exp2(S, TopAbs_FACE); Exp2.More();  Exp2.Next() ) {
1146
1147     const TopoDS_Face& F = TopoDS::Face( Exp2.Current() ) ;
1148
1149     if( mapFaces.Contains(F) )  {
1150       int n = GetIndexTopology( F, S ) ;
1151       if( n<=0 ) {
1152         THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::IndexOfFacesOfSubShell() : no index found", SALOME::BAD_PARAM);
1153       }
1154       size++;
1155       ListOfID->length(size) ;
1156       ListOfID[size-1] = n ;
1157     }
1158   }
1159
1160   return ListOfID._retn() ;
1161 }
1162
1163
1164
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 )
1175 {
1176   if( L.length() < 1 || S.IsNull() ) {
1177     return false ;
1178   }
1179
1180   aMap.Clear() ;
1181   for( int k=0; k<L.length(); k++ ) {
1182     /* indices start at 1 in list L */
1183     int j = 1 ;
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()) )
1188         {
1189           if( L[k] == j ) {
1190             aMap.Add( exp.Current() ) ;
1191           }
1192           j++ ;
1193         }
1194     }
1195   }
1196   return true ;
1197 }
1198
1199
1200
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 )
1211 {
1212   if( L.length() < 1 || S.IsNull() ) {
1213     return false ;
1214   }
1215
1216   aSequenceOfShapes.Clear() ;
1217   for( int k=0; k<L.length(); k++ ) {
1218     /* indices start at 1 in list L */
1219     int j = 1 ;
1220     TopExp_Explorer exp ;
1221     for(  exp.Init( S, TopAbs_ShapeEnum(subShapeType) ); exp.More(); exp.Next() ) {
1222       if( L[k] == j ) {
1223         aSequenceOfShapes.Append( exp.Current() ) ;
1224       }
1225       j++ ;
1226     }
1227   }
1228   return true ;
1229 }
1230
1231
1232
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)
1241 {
1242   Unexpect aCatch(SALOME_SalomeException);
1243   GEOM::GEOM_Gen::ListOfGeomShapes_var listOfGeomShapes = new GEOM::GEOM_Gen::ListOfGeomShapes;
1244   listOfGeomShapes->length(0) ;
1245
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);
1249
1250   if( ListOfID.length() < 1 )
1251     THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::SuppressFaces() : empty ListOfID", SALOME::BAD_PARAM);
1252
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();
1257   }
1258
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);
1265   }
1266
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   */
1270   int i = 0 ;
1271   for( TopExp_Explorer exp(aCompoundOfShells, TopAbs_SHELL); exp.More(); exp.Next() ) {
1272
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);
1276     }
1277
1278     GEOM::GEOM_Shape::ListOfSubShapeID_var aList = new GEOM::GEOM_Shape::ListOfSubShapeID;
1279     aList = IndexOfFacesOfSubShell(mainShape, aShell) ;
1280
1281      if( aList->length() < 1 ) {
1282       THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::SuppressFaces() : aList is empty", SALOME::BAD_PARAM);
1283     }
1284
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 ) ;
1290       exp.More() ;
1291       aShellOrFace = exp.Current() ;
1292     }
1293     else {
1294       aShellOrFace = aShell ;
1295     }
1296
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);
1301     }
1302
1303     InsertInLabelOneArgument(mainShape, shape, aShellOrFace, result, myCurrentOCAFDoc) ;
1304     i++ ;
1305     listOfGeomShapes->length(i) ;
1306     listOfGeomShapes[i-1] = result ;
1307   }
1308
1309   return listOfGeomShapes._retn() ;
1310 }
1311
1312
1313
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)
1319 {
1320   if( aShape.IsNull() || SS.IsEmpty() )
1321     return false ;
1322   for( int i=1; i<=SS.Length(); i++) {
1323     if( SS.Value(i).IsSame(aShape) )
1324       return true ;
1325   }
1326   return false ;
1327 }
1328
1329
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 )
1337 {
1338   MS.Clear() ;
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() ) ;
1348       }
1349     }
1350   }
1351   return ;
1352 }
1353
1354
1355 //================================================================================
1356 // function : MapRemoveSequence()
1357 // purpose  : In term of shapes ST = MS - SSRemove
1358 //          :
1359 //================================================================================
1360 void GEOM_Gen_i::MapRemoveSequence( const TopTools_MapOfShape& MS,
1361                                     const TopTools_SequenceOfShape& SSRemove,
1362                                     TopTools_SequenceOfShape& ST )
1363 {
1364   ST.Clear() ;
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 ) ;
1370   }
1371   return ;
1372 }
1373
1374
1375
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)
1390 {
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    */
1396
1397   if( SU.IsEmpty() ) {
1398     return ;
1399   }
1400
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) ;
1407     TopoDS_Shape F ;
1408     TopTools_ListOfShape LF ;
1409     int ind = aMapEdgesFaces.FindIndex(E) ;
1410
1411     /* LF is the list of faces for an edge of SU : may be empty no matter */
1412     LF = aMapEdgesFaces.FindFromIndex(ind) ;
1413
1414     TopTools_ListIteratorOfListOfShape it ;
1415     for( it.Initialize(LF); it.More(); it.Next() ) {
1416       F = it.Value() ;
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' !
1423         }
1424       }
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);
1428
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) ;
1433         }
1434         else {
1435           if( !MSfacesSuppress.Contains(F) ) {
1436             MSfacesSuppress.Add(F) ;
1437             MSfacesTmp.Add(F) ;
1438           }
1439         }
1440       }
1441     }
1442   }
1443   /* Call recursively this routine */
1444   SuppressHoleSubRoutine( mainShape, aFace, SSedgesOfWire, aMapEdgesFaces, MSfacesTmp, MSfacesSuppress, withEndFace, endFace, MSwireEndEdges ) ;
1445 }
1446
1447
1448
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,
1458                                     const int index,
1459                                     TopoDS_Shape& tds )
1460
1461 {
1462   if (aShape.IsNull() || index < 1)
1463     return false ;
1464   /* Indices start at 1 */
1465   int j = 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 */
1471       if( index == j ) {
1472         tds =  exp.Current() ;
1473         return true ;
1474       }
1475       j++ ;
1476     }
1477   }
1478   return false ;
1479 }
1480
1481
1482
1483 //================================================================================
1484 // function : SuppressHolesInFaceOrShell()  Main method.
1485 // purpose  : Suppress holes identified by wires in a single face or shell
1486 //
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)
1491 {
1492   Unexpect aCatch(SALOME_SalomeException);
1493   GEOM::GEOM_Shape_var result;
1494
1495   if( ListIdWires.length() < 1 )
1496     THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHolesInFace : no holes selected", SALOME::BAD_PARAM);
1497
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);
1501
1502   /* Create a map of wires/holes to suppress */
1503   TopTools_MapOfShape MapHoles ;
1504   for ( int i = 0; i < ListIdWires.length(); i++ ) {
1505     TopoDS_Shape W ;
1506     if( !GetShapeFromIndex( tds, TopAbs_WIRE, ListIdWires[i], W ) )
1507       THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHolesInFace() : bad index ?", SALOME::BAD_PARAM);
1508     MapHoles.Add( W ) ;
1509   }
1510
1511   /* Test if argument is a face or shell */
1512   bool isFace ;
1513   if( tds.ShapeType() == TopAbs_FACE )
1514     isFace = true ;
1515   else if ( tds.ShapeType() == TopAbs_SHELL )
1516     isFace = false ;
1517   else
1518     THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHolesInFace() : not a face or a shell", SALOME::BAD_PARAM);
1519
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) ;
1532       }
1533     }
1534   }
1535
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) ;
1543   }
1544
1545   if( MapFacesToModify.IsEmpty() )
1546     THROW_SALOME_CORBA_EXCEPTION("Error : empty map of faces", SALOME::BAD_PARAM);
1547
1548   if( isFace && MapFacesToModify.Extent() != 1 )
1549     THROW_SALOME_CORBA_EXCEPTION("Incoherent", SALOME::BAD_PARAM);
1550
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) ;
1559     return result ;
1560   }
1561
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() ) ;
1566     TopoDS_Face resF ;
1567     if( !RebuildFaceRemovingHoles( FF, MapHoles, resF ) )
1568       THROW_SALOME_CORBA_EXCEPTION(" Problem shell : !RebuildFaceRemovingHoles()", SALOME::BAD_PARAM);
1569     MapFacesModified.Add(resF) ;
1570   }
1571
1572   /* Rebuild the shell with faces modified and non modified */
1573   TopoDS_Shell resultShell ;
1574   BRepPrim_Builder B;
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() ) ) ;
1581
1582   B.CompleteShell(resultShell) ;
1583
1584   if( resultShell.IsNull() )
1585     THROW_SALOME_CORBA_EXCEPTION("Null or not valid result Shell", SALOME::BAD_PARAM) ;
1586
1587   /* Creation of CORBA object : shell topology */
1588   result = CreateObject(resultShell);
1589   InsertInLabelOneArgument(tds, shapeFaceShell, resultShell, result, myCurrentOCAFDoc) ;
1590   return result ;
1591 }
1592
1593
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 )
1602 {
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);
1607
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 ) ;
1611
1612   if( F2.Orientation() != aFace.Orientation() )
1613     F2.Orientation( aFace.Orientation() ) ;
1614
1615   BRepBuilderAPI_MakeFace aBuilder( F2 ) ;
1616   bool foundAndKeepHoles = false ;
1617   TopExp_Explorer exp ;
1618
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);
1625
1626       resultFace = TopoDS::Face(aBuilder) ;
1627       foundAndKeepHoles = true ;
1628     }
1629   }
1630
1631   if( !foundAndKeepHoles )
1632     resultFace = F2 ;
1633   else
1634     resultFace = TopoDS::Face(aBuilder) ;
1635
1636   return true ;
1637 }
1638
1639
1640
1641
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)
1654 {
1655   Unexpect aCatch(SALOME_SalomeException);
1656   GEOM::GEOM_Shape_var result;
1657   TopoDS_Face aFace ;
1658   TopoDS_Wire aWire ;
1659   TopoDS_Face endFace ;
1660   bool withEndFace ;
1661   TopoDS_Shape tmp ;
1662
1663   /* Retrieve 'aShape' the initial main shape selection */
1664   const TopoDS_Shape aShape = GetTopoShape(shape);
1665
1666   if( !IsValid(aShape) )
1667     THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHole() : non valid main shape", SALOME::BAD_PARAM);
1668
1669   if( ListIdFace.length() != 1 || ListIdWire.length() != 1 )
1670     THROW_SALOME_CORBA_EXCEPTION("bad list", SALOME::BAD_PARAM);
1671
1672   /* Retrieve 'aFace' selection */
1673   if( !GetShapeFromIndex( aShape, TopAbs_FACE, ListIdFace[0], tmp ) ) {
1674     THROW_SALOME_CORBA_EXCEPTION("face not found", SALOME::BAD_PARAM);
1675   }
1676   else {
1677     aFace = TopoDS::Face(tmp) ;
1678   }
1679   if( !IsValid(aFace) )
1680     THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHole() : face shape not valid", SALOME::BAD_PARAM);
1681
1682    /* Retrieve 'aWire' selection : Warning : index of wire refers to the face ! */
1683   TopoDS_Shape aTmp ;
1684   if( !GetShapeFromIndex( aFace, TopAbs_WIRE, ListIdWire[0], aTmp ) ) {
1685     THROW_SALOME_CORBA_EXCEPTION("wire not found", SALOME::BAD_PARAM);
1686   }
1687   else {
1688     aWire = TopoDS::Wire(aTmp) ;
1689   }
1690   if( !IsValid(aWire) )
1691     THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHole() : bad wire" , SALOME::BAD_PARAM);
1692
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);
1697
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);
1701
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 ;
1706   }
1707   else {
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);
1711
1712     /* Test if 'endFace' as at least one hole */
1713     endFace = TopoDS::Face(aTemp) ;
1714
1715     TopExp_Explorer fExp ;
1716     int nbWires = 0 ;
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) )
1720         nbWires++ ;
1721     }
1722     if(nbWires > 1)
1723       withEndFace = true ; /* at least 2 wires : outer wire plus an hole or more */
1724     else
1725       THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SuppressHole() : end face selection ?", SALOME::BAD_PARAM);
1726   }
1727
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);
1735     }
1736     else {
1737       SSedgesOfWire.Append( wireExp.Current() ) ;
1738     }
1739   }
1740   if( SSedgesOfWire.Length() < 1 )
1741     THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SupressHole() : no edge(s) for aWire", SALOME::BAD_PARAM);
1742
1743   /* Retrieve face ancestors of all edges of 'aWire' but avoiding 'aFace' */
1744
1745   TopTools_IndexedDataMapOfShapeListOfShape aMapEdgesFaces;
1746   TopTools_MapIteratorOfMapOfShape anIt ;
1747   TopTools_MapOfShape MFSuppress ;
1748   TopTools_MapOfShape MFSuppressTmp ;
1749   bool wireOnFace = false ;
1750
1751   TopExp::MapShapesAndAncestors(aShape, TopAbs_EDGE, TopAbs_FACE, aMapEdgesFaces) ;
1752   for( int h=1; h<=SSedgesOfWire.Length(); h++ ) {
1753
1754     TopoDS_Shape anEdgeOfWire = SSedgesOfWire.Value(h) ;
1755     int ind = aMapEdgesFaces.FindIndex(anEdgeOfWire) ;
1756     if(ind < 1)
1757       THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SupressHole() : index of edge", SALOME::BAD_PARAM);
1758
1759     TopTools_ListOfShape LF;
1760     LF = aMapEdgesFaces.FindFromIndex(ind) ; /* Contains all faces ancestors of an edge of wire */
1761     if( LF.IsEmpty() )
1762       THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SupressHole() : no face for an edge", SALOME::BAD_PARAM);
1763
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) ;
1770         MFSuppress.Add(F) ;
1771       }
1772       else {
1773         wireOnFace = true ;
1774       }
1775     }
1776   }
1777
1778   if( !wireOnFace ) {
1779     THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SupressHole() : wire not on selected face", SALOME::BAD_PARAM);
1780   }
1781
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 ) ;
1785
1786   TopoDS_Wire endWire ;
1787   if( withEndFace ) {
1788
1789     if( MSwireEndEdges.Extent() < 1 )
1790       THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SupressHole() : MSwireEndEdges.Extent() < 1", SALOME::BAD_PARAM);
1791
1792     if( !FindCompareWireHoleOnFace( endFace, MSwireEndEdges, endWire ) )
1793       THROW_SALOME_CORBA_EXCEPTION("in GEOM_Gen_i::SupressHole() : no endWire found", SALOME::BAD_PARAM);
1794    }
1795
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);
1800
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);
1803
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 ;
1807
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() ) ;
1813     nbSolid++ ;
1814     if( nbSolid > 1 )
1815       THROW_SALOME_CORBA_EXCEPTION("Limitation : main shape contents more than one solid", SALOME::BAD_PARAM);
1816   }
1817
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() ) ;
1822     nbShell++ ;
1823     if( nbShell > 1 )
1824       THROW_SALOME_CORBA_EXCEPTION("Limitation : main shape contents more than one shell", SALOME::BAD_PARAM);
1825   }
1826
1827   /* No shells and no solids : can send a compound even for a single face, see GUI ! */
1828   if( nbShell == 0 ) {
1829     finalShape = resTds ;
1830   }
1831
1832   /* a shell */
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 ;
1838   }
1839
1840   /* a solid with a shell */
1841   if( nbSolid == 1 && nbShell == 1) {
1842     BRepBuilderAPI_MakeSolid B;
1843     B.Add(shellResult) ;
1844     if( !B.IsDone() )
1845       THROW_SALOME_CORBA_EXCEPTION("Error : !B.IsDone()", SALOME::BAD_PARAM);
1846     finalShape = B.Solid() ;
1847   }
1848
1849   result = CreateObject(finalShape);
1850   InsertInLabelOneArgument(aShape, shape, finalShape, result, myCurrentOCAFDoc) ;
1851   return result ;
1852 }
1853
1854
1855
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 )
1862 {
1863   resultShell.Nullify() ;
1864   BRepPrim_Builder B ;
1865   B.MakeShell(resultShell) ;
1866   TopExp_Explorer ex ;
1867   int i = 0 ;
1868   for( ex.Init( Comp, TopAbs_FACE); ex.More(); ex.Next() ) {
1869     TopoDS_Face F = TopoDS::Face( ex.Current() ) ;
1870     if( !IsValid(F) ) {
1871       return false ;
1872     }
1873     B.AddShellFace( resultShell, F ) ;
1874   }
1875   B.CompleteShell(resultShell) ;
1876   if( resultShell.IsNull() ) {
1877     return false ;
1878   }
1879   return true ;
1880 }
1881
1882
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 )
1893 {
1894   aFoundWire.Nullify() ;
1895
1896   if( F.IsNull() )
1897     return false ;
1898
1899   /* Get the outer wire of aFace */
1900   TopoDS_Wire outerW = BRepTools::OuterWire(F) ;
1901   if( outerW.IsNull() || !IsValid(outerW) ) {
1902     return false ;
1903   }
1904
1905   int nbEdges = MSwireEdges.Extent() ;
1906   if( nbEdges < 1 ) {
1907     return false ;
1908   }
1909
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() ) ;
1914     int i = 0 ;
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) ) {
1919           i++ ;
1920           if( i == nbEdges ) {
1921             aFoundWire = W ;
1922             return true ;
1923           }
1924         }
1925       }
1926     }
1927   }
1928   return false ;
1929 }
1930
1931
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)
1944 {
1945   Unexpect aCatch(SALOME_SalomeException);
1946   BRep_Builder B;
1947   TopExp_Explorer exp ;
1948   TopoDS_Face newFace ;
1949
1950   resultTds.Nullify() ;
1951   TopoDS_Compound Comp ;
1952   B.MakeCompound (Comp);
1953
1954   try  {
1955
1956     /* Explore all faces of 'aShape' to rebuild a compound */
1957     for ( exp.Init(aShape, TopAbs_FACE) ; exp.More(); exp.Next() ) {
1958
1959       TopoDS_Face F1 = TopoDS::Face( exp.Current() );
1960       /* Rebuild face(s) not suppressed */
1961       if( !MFSuppress.Contains(F1) ) {
1962
1963         if( F1.IsEqual( aFace ) ) {
1964           TopTools_MapOfShape MSwire ;
1965           MSwire.Add(aWire) ;
1966           if( !RebuildFaceRemovingHoles(aFace, MSwire, newFace) ) {
1967             return false ;
1968           }
1969           B.Add( Comp, newFace ) ;
1970         }
1971         else {
1972           /* For any other face not suppressed */
1973           B.Add( Comp, F1 ) ;
1974         }
1975       }
1976     }
1977   }
1978   catch(Standard_Failure) {
1979     THROW_SALOME_CORBA_EXCEPTION("in BuildShapeHoleNotTraversing() : Exception catched", SALOME::BAD_PARAM);
1980   }
1981   resultTds = Comp ;
1982   return true ;
1983 }
1984
1985
1986
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)
2003 {
2004   Unexpect aCatch(SALOME_SalomeException);
2005   BRep_Builder B;
2006   TopExp_Explorer exp ;
2007   TopoDS_Face newFace ;
2008
2009   resultTds.Nullify() ;
2010   TopoDS_Compound Comp ;
2011   B.MakeCompound (Comp);
2012
2013   /* Necessary to use general method */
2014   TopTools_MapOfShape MSwire1  ;
2015   MSwire1.Add(aWire) ;
2016   TopTools_MapOfShape MSwire2 ;
2017   MSwire2.Add(endWire) ;
2018
2019   try  {
2020
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() );
2024
2025       /* Rebuild face(s) not suppressed */
2026       if( !MFSuppress.Contains(F1) ) {
2027
2028         /* Rebuild 'aFace' */
2029         if( F1.IsEqual( aFace ) && !F1.IsEqual( endFace ) ) {
2030           if( !RebuildFaceRemovingHoles(aFace, MSwire1, newFace) ) {
2031             return false ;
2032           }
2033           B.Add(Comp, newFace) ;
2034         }
2035
2036         /* Rebuild 'endFace' */
2037         if( !F1.IsEqual( aFace ) && F1.IsEqual( endFace ) ) {
2038           if( !RebuildFaceRemovingHoles(endFace, MSwire2, newFace) ) {
2039             return false ;
2040           }
2041           B.Add(Comp, newFace) ;
2042         }
2043
2044         /* Hole in the same face : aFace = endFace */
2045         if( F1.IsEqual( aFace ) && F1.IsEqual( endFace ) ) {
2046           TopoDS_Face FF ;
2047           if( !RebuildFaceRemovingHoles(aFace, MSwire1, newFace) || !RebuildFaceRemovingHoles(newFace, MSwire2, FF) ) {
2048             return false ;
2049           }
2050           B.Add( Comp, FF ) ;
2051         }
2052
2053         /* For any other face not suppressed */
2054         if( !F1.IsEqual(aFace) && !F1.IsEqual( endFace ) ) {
2055           B.Add( Comp, F1 ) ;
2056         }
2057
2058       }
2059     }
2060   }
2061   catch(Standard_Failure) {
2062     THROW_SALOME_CORBA_EXCEPTION("in BuildShapeHoleTraversing() : Exception catched", SALOME::BAD_PARAM);
2063   }
2064   resultTds = Comp ;
2065   return true ;
2066 }
2067
2068 //=======================================================================
2069 //function : SortShapes
2070 //purpose  :
2071 //=======================================================================
2072
2073 static void SortShapes(TopTools_ListOfShape& SL)
2074 {
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;
2080
2081   // Computing of CentreOfMass
2082   Standard_Integer Index;
2083   GProp_GProps GPr;
2084   gp_Pnt GPoint;
2085   TopTools_ListIteratorOfListOfShape it(SL);
2086   for (Index=1;  it.More();  Index++)
2087   {
2088     TopoDS_Shape S = it.Value();
2089     SL.Remove( it ); // == it.Next()
2090     aShapes(Index) = S;
2091     OrderInd.SetValue (Index, Index);
2092     if (S.ShapeType() == TopAbs_VERTEX)
2093     {
2094       GPoint = BRep_Tool::Pnt( TopoDS::Vertex( S ));
2095       Length.SetValue( Index, (Standard_Real) S.Orientation());
2096     }
2097     else
2098     {
2099       BRepGProp::LinearProperties (S, GPr);
2100       GPoint = GPr.CentreOfMass();
2101       Length.SetValue( Index, GPr.Mass() );
2102     }
2103     MidXYZ.SetValue(Index,
2104                     GPoint.X()*999 + GPoint.Y()*99 + GPoint.Z()*0.9);
2105   }
2106   // Sorting
2107   Standard_Integer aTemp;
2108   Standard_Boolean exchange, Sort = Standard_True;
2109   while (Sort)
2110   {
2111     Sort = Standard_False;
2112     for (Index=1; Index < MaxShapes; Index++)
2113     {
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;
2119       else
2120         exchange = Standard_False;
2121       if (exchange)
2122       {
2123         aTemp = OrderInd(Index);
2124         OrderInd(Index) = OrderInd(Index+1);
2125         OrderInd(Index+1) = aTemp;
2126         Sort = Standard_True;
2127       }
2128     }
2129   }
2130   for (Index=1; Index <= MaxShapes; Index++)
2131     SL.Append( aShapes( OrderInd(Index) ));
2132 }
2133
2134 //================================================================================
2135 // function : SubShape()
2136 // purpose  : Method for GUI or TUI
2137 //================================================================================
2138
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)
2143 {
2144   Unexpect aCatch(SALOME_SalomeException);
2145   return SubShapesOne(shape, (TopAbs_ShapeEnum) ShapeType, ListOfID);
2146 }
2147
2148 //================================================================================
2149 // function : SubShapeSorted()
2150 // purpose  : Method for GUI or TUI
2151 //================================================================================
2152
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)
2157 {
2158   Unexpect aCatch(SALOME_SalomeException);
2159   return SubShapesOne(shape, (TopAbs_ShapeEnum) ShapeType, ListOfID, Standard_True);
2160 }
2161
2162 //================================================================================
2163 // function : SubShapesOne()
2164 // purpose  :
2165 //================================================================================
2166
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)
2172 {
2173   Unexpect aCatch(SALOME_SalomeException);
2174   GEOM::GEOM_Shape_var result;
2175   TopoDS_Shape mainShape;
2176   TopoDS_Shape mainTopo = GetTopoShape(shape);
2177
2178   bool main = false;
2179   while ( !main ) {
2180     if ( shape->IsMainShape() ) {
2181       mainShape = GetTopoShape(shape);
2182       main = true;
2183     } else
2184       shape = GetIORFromString( shape->MainName() );
2185   }
2186
2187   if(ListOfID.length() < 1) {
2188     THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::SubShape() : bad list of shapes",
2189                                  SALOME::BAD_PARAM);
2190   }
2191
2192   /* Create a sequence of all sub shapes */
2193   TopTools_ListOfShape listShape;
2194
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() );
2200   }
2201
2202   if (listShape.IsEmpty()) return shape;
2203
2204   if (Sort)
2205     SortShapes(listShape);
2206
2207   TopTools_SequenceOfShape SS;
2208   TopTools_ListIteratorOfListOfShape it (listShape);
2209   for (; it.More(); it.Next())
2210     SS.Append( it.Value() );
2211
2212   TopoDS_Shape SubShape;
2213   if (ListOfID.length() == 1)
2214     SubShape = SS.Value(ListOfID[0]);
2215   else
2216   {
2217     BRep_Builder B;
2218     TopoDS_Compound Comp;
2219     B.MakeCompound (Comp);
2220     unsigned int ind;
2221     for ( ind = 0; ind < ListOfID.length(); ind++ )
2222       B.Add( Comp, SS.Value(ListOfID[ind]) );
2223     SubShape = Comp;
2224   }
2225
2226   if ( !TNaming_Tool::HasLabel( myCurrentOCAFDoc->Main(), SubShape ) )
2227     result = CreateSubObject( SubShape, shape, ListOfID);
2228   else {
2229     int TransDef;
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);
2237       }
2238     }
2239   }
2240
2241   return result;
2242 }
2243
2244
2245 //================================================================================
2246 // function : SubShapeAll()
2247 // purpose  : Explode a shape in all sub shapes with a type (Method for TUI or GUI)
2248 //================================================================================
2249
2250 GEOM::GEOM_Gen::ListOfGeomShapes* GEOM_Gen_i::SubShapeAll(GEOM::GEOM_Shape_ptr shape,
2251                                                           CORBA::Short ShapeType)
2252      throw (SALOME::SALOME_Exception)
2253 {
2254   Unexpect aCatch(SALOME_SalomeException);
2255   return SubShapesAll(shape, (TopAbs_ShapeEnum) ShapeType);
2256 }
2257
2258 //================================================================================
2259 // function : SubShapeAllSorted()
2260 // purpose  : Explode a shape in all sub shapes with a type (Method for TUI or GUI)
2261 //================================================================================
2262
2263 GEOM::GEOM_Gen::ListOfGeomShapes* GEOM_Gen_i::SubShapeAllSorted(GEOM::GEOM_Shape_ptr shape,
2264                                                                 CORBA::Short ShapeType)
2265      throw (SALOME::SALOME_Exception)
2266 {
2267   Unexpect aCatch(SALOME_SalomeException);
2268   return SubShapesAll(shape, (TopAbs_ShapeEnum) ShapeType, Standard_True);
2269 }
2270
2271 //================================================================================
2272 // function : SubShapeAllSorted()
2273 // purpose  :
2274 //================================================================================
2275
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)
2280 {
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) ;
2285
2286   TopoDS_Shape mainTopo = GetTopoShape(shape);
2287   TopoDS_Shape mainShape;
2288   bool main = false;
2289   while ( !main ) {
2290     if ( shape->IsMainShape() ) {
2291       mainShape = GetTopoShape(shape);
2292       main = true;
2293     } else
2294       shape = GetIORFromString( shape->MainName() );
2295   }
2296
2297   if( mainTopo.IsNull() ) {
2298     THROW_SALOME_CORBA_EXCEPTION("In GEOM_Gen_i::SubShapeAll() : null main shape",
2299                                  SALOME::BAD_PARAM);
2300   }
2301
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) ;
2305
2306     // retrieve all subshapes
2307   TopTools_MapOfShape mapShape;
2308   TopTools_ListOfShape listShape;
2309
2310   if (mainTopo.ShapeType()==TopAbs_COMPOUND  &&  ShapeType==TopAbs_SHAPE)
2311   {
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() );
2316   }
2317   else
2318   {
2319     TopExp_Explorer exp ( mainTopo, ShapeType);
2320     for ( ; exp.More(); exp.Next() )
2321       if (mapShape.Add( exp.Current() ))
2322         listShape.Append( exp.Current() );
2323   }
2324
2325   if (Sort)
2326     SortShapes(listShape);
2327
2328   /* Create all sub shapes */
2329   int index;
2330   GEOM::GEOM_Shape_var result;
2331
2332   TopTools_ListIteratorOfListOfShape itSub (listShape);
2333   for (index = 1; itSub.More(); itSub.Next(), ++index)
2334   {
2335     const TopoDS_Shape& SubShape = itSub.Value();
2336     // check if SubShape is already in OCAFDS and ...
2337     if ( TNaming_Tool::HasLabel( myCurrentOCAFDoc->Main(), SubShape ) )
2338     {
2339       int TransDef;
2340       TDF_Label Lab = TNaming_Tool::Label( myCurrentOCAFDoc->Main(), SubShape, TransDef );
2341       Handle(TDataStd_Name) Att;
2342       if ( Lab.FindAttribute( TDataStd_Name::GetID(), Att ) )
2343       {
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 )
2348         {
2349           listOfGeomShapes->length(index) ;
2350           listOfGeomShapes[index-1] = result ;
2351           continue;
2352         }
2353       }
2354     }
2355
2356 //      if (Sort)
2357 //        ListOfID[0] = index;
2358 //      else
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 ;
2364   }
2365
2366   return listOfGeomShapes._retn() ;
2367 }
2368
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)
2377 {
2378   Unexpect aCatch(SALOME_SalomeException);
2379   GEOM::GEOM_Shape_var result;
2380   TopoDS_Shape shape ;
2381   TopoDS_Shape aShape1  ;
2382   TopoDS_Shape aShape2  ;
2383
2384   try {
2385     aShape1 = GetTopoShape(shape1) ;
2386     aShape2 = GetTopoShape(shape2) ;
2387
2388     if( aShape1.IsNull() || aShape2.IsNull() ) {
2389       THROW_SALOME_CORBA_EXCEPTION("Boolean aborted : argument shape is null", SALOME::BAD_PARAM);
2390     }
2391
2392     switch (operation)
2393       {
2394       case 1 :   /* Common */
2395         shape = BRepAlgoAPI_Common(aShape1, aShape2).Shape();
2396         break ;
2397       case 2 :   /* Cut */
2398         shape = BRepAlgoAPI_Cut(aShape1, aShape2).Shape();
2399         break ;
2400       case 3 :   /* Fuse */
2401         shape = BRepAlgoAPI_Fuse(aShape1, aShape2).Shape();
2402         break ;
2403       case 4 :   /* Section */
2404         shape = BRepAlgoAPI_Section(aShape1, aShape2).Shape();
2405         break ;
2406       default :
2407         MESSAGE("Boolean operation not known : " << operation ) ;
2408         return result ;
2409       }
2410   }
2411   catch(Standard_Failure) {
2412     THROW_SALOME_CORBA_EXCEPTION("Exception catched in boolean operation", SALOME::BAD_PARAM);
2413   }
2414
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);
2418   }
2419
2420   result = CreateObject(shape) ;
2421
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() ) ;
2429
2430   /* Create a new label */
2431   TDF_Label NewLab = Lab.NewChild();
2432   TCollection_ExtendedString Value("Arguments");
2433   TDataStd_Name::Set(NewLab,Value);
2434
2435   TDF_Label NewLab1 = NewLab.NewChild();
2436   TDF_Label RefLab;
2437   TDF_Tool::Label(myCurrentOCAFDoc->GetData(), shape1->ShapeId(), RefLab);
2438   TDF_Reference::Set(NewLab1, RefLab);
2439
2440   TDF_Label NewLab2 = NewLab.NewChild();
2441   TDF_Tool::Label(myCurrentOCAFDoc->GetData(), shape2->ShapeId(), RefLab);
2442   TDF_Reference::Set(NewLab2, RefLab);
2443   return result ;
2444 }
2445
2446
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)
2454 {
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);
2461   }
2462   TopoDS_Shape shape;
2463   try {
2464     shape  = BRepAlgoAPI_Fuse(aShape1, aShape2).Shape();
2465   }
2466   catch(Standard_Failure) {
2467     THROW_SALOME_CORBA_EXCEPTION("Exception catched in Fuse operation", SALOME::BAD_PARAM);
2468   }
2469
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);
2473   }
2474
2475   result = CreateObject(shape) ;
2476
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() ) ;
2484
2485   /* Create a new label */
2486   TDF_Label NewLab = Lab.NewChild();
2487   TCollection_ExtendedString Value("Arguments");
2488   TDataStd_Name::Set(NewLab,Value);
2489
2490   TDF_Label NewLab1 = NewLab.NewChild();
2491   TDF_Label RefLab;
2492   TDF_Tool::Label(myCurrentOCAFDoc->GetData(), shape1->ShapeId(), RefLab);
2493   TDF_Reference::Set(NewLab1, RefLab);
2494
2495   TDF_Label NewLab2 = NewLab.NewChild();
2496   TDF_Tool::Label(myCurrentOCAFDoc->GetData(), shape2->ShapeId(), RefLab);
2497   TDF_Reference::Set(NewLab2, RefLab);
2498   return result ;
2499 }
2500
2501
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,
2507                                       CORBA::Double y,
2508                                       CORBA::Double z,
2509                                       CORBA::Double vx,
2510                                       CORBA::Double vy,
2511                                       CORBA::Double vz)
2512 {
2513   GEOM::AxisStruct A ;
2514   A.x  = x ;   A.y  = y ;   A.z  = z ;
2515   A.vx = vx ;  A.vy = vy ;  A.vz = vz ;
2516   return A ;
2517 }
2518
2519
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,
2525                                         CORBA::Double y,
2526                                         CORBA::Double z)
2527 {
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" );
2532   return p ;
2533 }
2534
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)
2540 {
2541   GEOM::DirStruct d ;
2542   d.PS.x = p.x ;  d.PS.y = p.y ;  d.PS.z = p.z ;
2543   return d ;
2544 }
2545
2546 //=================================================================================
2547 // function : MakeBox()
2548 // purpose  : Create a box topology.
2549 //=================================================================================
2550 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeBox(CORBA::Double x1,
2551                                          CORBA::Double y1,
2552                                          CORBA::Double z1,
2553                                          CORBA::Double x2,
2554                                          CORBA::Double y2,
2555                                          CORBA::Double z2)
2556      throw (SALOME::SALOME_Exception)
2557 {
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 ;
2562   TopoDS_Shape tds ;
2563   try {
2564     tds = BRepPrimAPI_MakeBox(P1,P2).Shape();
2565   }
2566   catch(Standard_Failure) {
2567     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeBox", SALOME::BAD_PARAM);
2568   }
2569
2570   if (tds.IsNull()) {
2571     THROW_SALOME_CORBA_EXCEPTION("Make Box aborted : null shape", SALOME::BAD_PARAM);
2572   }
2573
2574   result = CreateObject(tds);
2575   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
2576   result->ShapeId(entry) ;
2577   return result;
2578 }
2579
2580
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)
2590 {
2591   Unexpect aCatch(SALOME_SalomeException);
2592   GEOM::GEOM_Shape_var result;
2593   TopoDS_Shape tds ;
2594   gp_Pnt p(pstruct.x, pstruct.y, pstruct.z) ;
2595   gp_Dir d(dstruct.PS.x, dstruct.PS.y, dstruct.PS.z) ;
2596   gp_Ax2 axis(p, d) ;
2597
2598   try {
2599     tds = BRepPrimAPI_MakeCylinder(axis, radius, height).Shape();
2600   }
2601   catch(Standard_Failure) {
2602     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeCylinder", SALOME::BAD_PARAM);
2603   }
2604
2605   if (tds.IsNull()) {
2606     THROW_SALOME_CORBA_EXCEPTION("Make Cylinder aborted", SALOME::BAD_PARAM);
2607   }
2608   result = CreateObject(tds);
2609   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
2610   result->ShapeId(entry);
2611   return result ;
2612 }
2613
2614 //================================================================================
2615 // function : MakeSphere()
2616 // purpose  : Make a sphere topology
2617 //================================================================================
2618 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeSphere(CORBA::Double x1,
2619                                       CORBA::Double y1,
2620                                       CORBA::Double z1,
2621                                       CORBA::Double radius)
2622   throw (SALOME::SALOME_Exception)
2623 {
2624   Unexpect aCatch(SALOME_SalomeException);
2625   GEOM::GEOM_Shape_var result ;
2626   TopoDS_Shape tds ;
2627   try {
2628     tds = BRepPrimAPI_MakeSphere(gp_Pnt(x1,y1,z1), radius).Shape();
2629   }
2630   catch(Standard_Failure) {
2631     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeSphere", SALOME::BAD_PARAM);
2632   }
2633
2634   if (tds.IsNull()) {
2635     THROW_SALOME_CORBA_EXCEPTION("Make Sphere aborted", SALOME::BAD_PARAM);
2636   }
2637   result = CreateObject(tds) ;
2638   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
2639   result->ShapeId(entry);
2640   return result;
2641 }
2642
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)
2652 {
2653   Unexpect aCatch(SALOME_SalomeException);
2654   GEOM::GEOM_Shape_var result;
2655   TopoDS_Shape tds ;
2656   gp_Pnt p(pstruct.x, pstruct.y, pstruct.z) ;
2657   gp_Dir d(dstruct.PS.x, dstruct.PS.y, dstruct.PS.z) ;
2658   gp_Ax2 axis(p, d) ;
2659
2660   try {
2661   tds = BRepPrimAPI_MakeTorus(axis, major_radius, minor_radius).Shape();
2662   }
2663   catch(Standard_Failure) {
2664     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeTorus", SALOME::BAD_PARAM);
2665   }
2666
2667   if (tds.IsNull()) {
2668     THROW_SALOME_CORBA_EXCEPTION("Make torus aborted", SALOME::BAD_PARAM);
2669   }
2670   result = CreateObject(tds);
2671   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
2672   result->ShapeId(entry);
2673   return result ;
2674 }
2675
2676
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)
2687 {
2688   Unexpect aCatch(SALOME_SalomeException);
2689   GEOM::GEOM_Shape_var result;
2690   TopoDS_Shape tds ;
2691   gp_Pnt p(pstruct.x, pstruct.y, pstruct.z) ;
2692   gp_Dir d(dstruct.PS.x, dstruct.PS.y, dstruct.PS.z) ;
2693   gp_Ax2 axis(p, d) ;
2694
2695   try {
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();
2699     }
2700     else {
2701       tds = BRepPrimAPI_MakeCone(axis, radius1, radius2, height).Shape();
2702     }
2703   }
2704   catch(Standard_Failure) {
2705     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeCone", SALOME::BAD_PARAM);
2706   }
2707
2708   if (tds.IsNull()) {
2709     THROW_SALOME_CORBA_EXCEPTION("Make Cone aborted", SALOME::BAD_PARAM);
2710   }
2711   result = CreateObject(tds);
2712   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
2713   result->ShapeId(entry);
2714   return result ;
2715 }
2716
2717
2718
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)
2726 {
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
2732   try {
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);    }
2736
2737     MESSAGE("ImportIGES : all Geometry Transfer" << endl ) ;
2738 #if OCC_VERSION_MAJOR >= 5
2739     aReader.ClearShapes();
2740     aReader.TransferRoots();
2741 #else
2742     aReader.Clear();
2743     aReader.TransferRoots(false);
2744 #endif
2745     MESSAGE("ImportIGES : count of shapes produced = " << aReader.NbShapes() << endl );
2746     TopoDS_Shape shape = aReader.OneShape();
2747
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);
2753       return result ;
2754     }
2755   }
2756   catch(Standard_Failure) {
2757     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ImportIGES()", SALOME::BAD_PARAM);
2758   }
2759
2760   THROW_SALOME_CORBA_EXCEPTION("Import IGES aborted : internal error", SALOME::BAD_PARAM);
2761 }
2762
2763
2764
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)
2772 {
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
2778
2779   TopoDS_Compound compound;
2780   BRep_Builder B;
2781   B.MakeCompound( compound );
2782
2783   try {
2784     IFSelect_ReturnStatus status = aReader.ReadFile((char*)filename);
2785
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);
2792
2793       for ( Standard_Integer n=1; n <= nbr; n++) {
2794
2795         Standard_Boolean ok = aReader.TransferRoot(n);
2796         /* Collecting resulting entities */
2797         Standard_Integer nbs = aReader.NbShapes();
2798         if (nbs == 0)
2799           THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ImportStep", SALOME::BAD_PARAM) ;
2800
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) ;
2805
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);
2811             return result ;
2812           }
2813           else {
2814             B.Add( compound, aShape ) ;
2815           }
2816         }
2817       }
2818
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);
2825      return result ;
2826     }
2827
2828   }
2829   catch(Standard_Failure) {
2830     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ImportStep", SALOME::BAD_PARAM);
2831   }
2832   return result ;
2833 }
2834
2835
2836 //==================================================================================
2837 // function : Partition()
2838 // purpose  :
2839 //==================================================================================
2840
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)
2848 {
2849   Unexpect aCatch(SALOME_SalomeException);
2850   GEOM::GEOM_Shape_var aResult;
2851   TopoDS_Shape tds ;
2852   //MESSAGE ("In Partition");
2853   try {
2854
2855     unsigned int ind, nbshapes = 0;
2856     nbshapes += ListShapes.length() + ListTools.length();
2857     nbshapes += ListKeepInside.length() + ListRemoveInside.length();
2858
2859     Partition_Spliter PS;
2860     TopTools_MapOfShape ShapesMap(nbshapes), ToolsMap(nbshapes);
2861
2862     // add object shapes that are in ListShapes;
2863     for ( ind = 0; ind < ListShapes.length(); ind++) {
2864
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);
2870       }
2871       if ( ShapesMap.Add( Shape ))
2872         PS.AddShape(Shape);
2873     }
2874
2875     // add tool shapes that are in ListTools and not in ListShapes;
2876     for (ind = 0; ind < ListTools.length(); ind++) {
2877
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);
2882       }
2883       if ( !ShapesMap.Contains( Shape ) && ToolsMap.Add( Shape ))
2884         PS.AddTool(Shape);
2885     }
2886
2887     // add shapes that are in ListKeepInside, as object shapes;
2888     for (ind = 0; ind < ListKeepInside.length(); ind++) {
2889
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);
2895       }
2896       if (!ToolsMap.Contains( Shape ) &&
2897           ShapesMap.Add( Shape ))
2898         PS.AddShape(Shape);
2899     }
2900
2901     // add shapes that are in ListRemoveInside, as object shapes;
2902     for (ind = 0; ind < ListRemoveInside.length(); ind++) {
2903
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);
2909       }
2910       if (!ToolsMap.Contains( Shape ) &&
2911           ShapesMap.Add( Shape ) )
2912         PS.AddShape(Shape);
2913     }
2914
2915     //MESSAGE ( "Partition::Compute() " );
2916     PS.Compute ((TopAbs_ShapeEnum) Limit);
2917     //MESSAGE ( "Partition::Compute() - END" );
2918
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 );
2924     }
2925
2926     // suppress result inside of shapes in RInsideMap
2927     for (ind = 0; ind < ListRemoveInside.length(); ind++) {
2928
2929       GEOM::GEOM_Shape_var aShape = GetIORFromString( ListRemoveInside[ind] );
2930       TopoDS_Shape Shape = GetTopoShape(aShape);
2931       PS.RemoveShapesInside( Shape );
2932     }
2933
2934     tds = PS.Shape();
2935
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);
2939     }
2940   }
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);
2944   }
2945
2946   aResult = CreateObject(tds) ;
2947
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() ) ;
2954
2955   // add arguments
2956
2957   /* Create a new label */
2958   TDF_Label NewLab = Lab.NewChild();
2959   TCollection_ExtendedString Value("Arguments");
2960   TDataStd_Name::Set(NewLab,Value);
2961
2962   // object shapes
2963   for (unsigned int ind = 0; ind < ListShapes.length(); ind++) {
2964     TDF_Label NewLab1 = NewLab.NewChild();
2965
2966     GEOM::GEOM_Shape_var aShape = GetIORFromString( ListShapes[ind] );
2967     Standard_CString anEntry = aShape->ShapeId();
2968     TDF_Label RefLab;
2969     TDF_Tool::Label(myCurrentOCAFDoc->GetData(), anEntry, RefLab);
2970     TDF_Reference::Set(NewLab1,RefLab);
2971   }
2972   // tool shapes
2973   for (unsigned int ind = 0; ind < ListTools.length(); ind++) {
2974     TDF_Label NewLab1 = NewLab.NewChild();
2975
2976     GEOM::GEOM_Shape_var aShape = GetIORFromString( ListTools[ind] );
2977     Standard_CString anEntry = aShape->ShapeId();
2978     TDF_Label RefLab;
2979     TDF_Tool::Label(myCurrentOCAFDoc->GetData(), anEntry, RefLab);
2980     TDF_Reference::Set(NewLab1,RefLab);
2981   }
2982   // limit shapes 1
2983   for (unsigned int ind = 0; ind < ListKeepInside.length(); ind++) {
2984     TDF_Label NewLab1 = NewLab.NewChild();
2985
2986     GEOM::GEOM_Shape_var aShape = GetIORFromString( ListKeepInside[ind] );
2987     Standard_CString anEntry = aShape->ShapeId();
2988     TDF_Label RefLab;
2989     TDF_Tool::Label(myCurrentOCAFDoc->GetData(), anEntry, RefLab);
2990     TDF_Reference::Set(NewLab1,RefLab);
2991   }
2992   // limit shapes 2
2993   for (unsigned int ind = 0; ind < ListRemoveInside.length(); ind++) {
2994     TDF_Label NewLab1 = NewLab.NewChild();
2995
2996     GEOM::GEOM_Shape_var aShape = GetIORFromString( ListRemoveInside[ind] );
2997     Standard_CString anEntry = aShape->ShapeId();
2998     TDF_Label RefLab;
2999     TDF_Tool::Label(myCurrentOCAFDoc->GetData(), anEntry, RefLab);
3000     TDF_Reference::Set(NewLab1,RefLab);
3001   }
3002
3003   return aResult;
3004 }
3005
3006
3007
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)
3019 {
3020   Unexpect aCatch(SALOME_SalomeException);
3021   GEOM::GEOM_Shape_var result ;
3022   TopoDS_Face tds ;
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);
3026   }
3027
3028   try {
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 ;
3035
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);
3041       }
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) ;
3046       i++ ;
3047     }
3048
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) ;
3054
3055     GeomFill_AppSurf App(mindeg, maxdeg, tol3d, tol2d, nbiter) ; /* user parameters */
3056     App.Perform(Line, Section) ;
3057
3058     if (!App.IsDone()) {
3059       THROW_SALOME_CORBA_EXCEPTION("Filling aborted : non valid shape result", SALOME::BAD_PARAM);
3060     }
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(),
3064                                                               App.SurfWeights(),
3065                                                               App.SurfUKnots(),
3066                                                               App.SurfVKnots(),
3067                                                               App.SurfUMults(),
3068                                                               App.SurfVMults(),
3069                                                               App.UDegree(),
3070                                                               App.VDegree());
3071
3072     if( GBS.IsNull() )  {
3073       THROW_SALOME_CORBA_EXCEPTION("Make Filling aborted", SALOME::BAD_PARAM);
3074     }
3075     tds  = BRepBuilderAPI_MakeFace(GBS) ;
3076   }
3077   catch(Standard_Failure) {
3078     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeFilling", SALOME::BAD_PARAM);
3079   }
3080
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);
3084   }
3085   else {
3086     result = CreateObject(tds) ;
3087     InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
3088   }
3089
3090   return result ;
3091 }
3092
3093
3094 //=================================================================================
3095 // function : MakeGlueFaces()
3096 // purpose  :
3097 //=================================================================================
3098
3099 TopoDS_Face GEOM_Gen_i::FindSameFace(const TopoDS_Shape& aShape,
3100                                      const TopoDS_Face& F,
3101                                      double tol3d)
3102 {
3103   TopoDS_Face aFace;
3104   bool isSame = false;
3105   for (TopExp_Explorer exf(aShape,TopAbs_FACE); exf.More(); exf.Next())
3106     {
3107       //MESSAGE("--- test a face");
3108       int nbFound = 0;
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())
3113         {
3114           const TopoDS_Vertex& V = TopoDS::Vertex(exp.Current());
3115           liste1.Append(V);
3116         }
3117       for (TopExp_Explorer exp(F,TopAbs_VERTEX); exp.More(); exp.Next())
3118         {
3119           const TopoDS_Vertex& V = TopoDS::Vertex(exp.Current());
3120           liste2.Append(V);
3121         }
3122       isSame = false;
3123       if (liste1.Extent() == liste2.Extent())
3124         {
3125           TopTools_ListIteratorOfListOfShape it1(liste1);
3126           isSame = true;
3127           for (; it1.More(); it1.Next())
3128             {
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())
3133                 {
3134                   gp_Pnt P2 = BRep_Tool::Pnt(TopoDS::Vertex(it2.Value()));
3135                   double d = P1.Distance(P2);
3136                   if (d < tol3d)
3137                     {
3138                       nbFound++;
3139                       //MESSAGE("    found Same Point : "<<nbFound<<" - "<<P1.X()<<" "<<P1.Y()<<" "<<P1.Z());
3140                       foundSamePoint = true;
3141                       break;
3142                     }
3143                 }
3144               isSame = isSame && foundSamePoint;
3145               if (! isSame) break; // a vertex does not correspond : not same face
3146             }
3147         }
3148       if (isSame)
3149         {
3150           //MESSAGE("    --- Found Same Face");
3151           break; // a face corresponding to F is found
3152         }
3153     }
3154   if (! isSame) aFace.Nullify(); // return null face
3155   return aFace;
3156 }
3157
3158 TopoDS_Edge GEOM_Gen_i::FindSameEdge(const TopoDS_Face& nf,
3159                                      TopoDS_Edge& Eold,
3160                                      double tol3d)
3161 {
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);
3167   TopoDS_Edge Enew;
3168   for (TopExp_Explorer ee(newFace,TopAbs_EDGE); ee.More(); ee.Next())
3169     {
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))
3180         {
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));
3184           break;
3185         }
3186       if ((dfl < tol3d) && (dlf <tol3d))
3187         {
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));
3191           break;
3192         }
3193     }
3194   return Enew;
3195 }
3196
3197 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeGlueFaces(GEOM::GEOM_Shape_ptr myShape,
3198                                                double tol3d)
3199   throw (SALOME::SALOME_Exception)
3200 {
3201
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)
3217
3218   Unexpect aCatch(SALOME_SalomeException);
3219   GEOM::GEOM_Shape_var result ;
3220   TopoDS_Shape tds ;
3221   TopoDS_Shape aShape = GetTopoShape(myShape) ;
3222   TopoDS_Compound C;
3223   BRep_Builder bu;
3224   bu.MakeCompound(C); // empty compound;
3225   TopTools_ListOfShape shellList;
3226   for (TopExp_Explorer exp(aShape,TopAbs_SHELL); exp.More(); exp.Next())
3227     {
3228       const TopoDS_Shell& S = TopoDS::Shell(exp.Current());
3229       shellList.Append(S);
3230     }
3231   TopTools_ListIteratorOfListOfShape its(shellList);
3232   if ( ! its.More())
3233     {
3234       THROW_SALOME_CORBA_EXCEPTION("glue aborted : no shell in shape", SALOME::BAD_PARAM);
3235     }
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)
3242     {
3243       //MESSAGE("more shells : "<< shellList.Extent());
3244       shellAdded = false;
3245       its.Initialize(shellList);
3246       for(; its.More(); its.Next())
3247         {
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())
3260             {
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())
3265                 {
3266                   //MESSAGE("--- face found");
3267                   newFaces.Append(newFace);
3268                   bu.Add(CFN, newFace); // common faces from new compound
3269                   oldFaces.Append(F);
3270                   for (TopExp_Explorer ee(F,TopAbs_EDGE);ee.More();ee.Next())
3271                     {
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);
3277                     }
3278                 }
3279               else
3280                 {
3281                   //MESSAGE("---");
3282                   addFaces.Append(F);
3283                   bu.Add(CFO, F); // not common faces from shell to add
3284                 }
3285             }
3286           if ( !newFaces.IsEmpty())
3287             {
3288               //MESSAGE("--- some faces found ---");
3289               shellAdded = true;
3290               BRepTools_Quilt glue;
3291               glue.Add(CFN);
3292               TopTools_ListIteratorOfListOfShape ito(oldEdges);
3293               TopTools_ListIteratorOfListOfShape itn(newEdges);
3294               for (; ito.More(); ito.Next())
3295                 {
3296                   //MESSAGE("--- bind");
3297                   glue.Bind(TopoDS::Edge(ito.Value()), TopoDS::Edge(itn.Value()));
3298                   itn.Next();
3299                 }
3300               glue.Add(CFO);
3301               TopoDS_Compound newc = TopoDS::Compound(glue.Shells());
3302               for (TopExp_Explorer exs(newc,TopAbs_SHELL); exs.More(); exs.Next())
3303                 {
3304                   TopoDS_Shell NS = TopoDS::Shell(exs.Current());
3305                   bu.Add(C, NS);
3306                 }
3307               shellList.Remove(its);
3308               //MESSAGE("--- remove shell from list");
3309               break;
3310             }
3311         }
3312     }
3313   //MESSAGE("---" << shellList.Extent() << " " << shellAdded);
3314
3315   TopExp_Explorer  exp(C,TopAbs_SHELL);
3316   Standard_Integer ish=0;
3317   TopoDS_Compound  Res;
3318   TopoDS_Solid     Sol;
3319   BRep_Builder     B;
3320   B.MakeCompound(Res);
3321   TopoDS_Shape theShape;
3322
3323   for (; exp.More(); exp.Next())
3324     {
3325       TopoDS_Shape Sh = exp.Current();
3326       B.MakeSolid(Sol);
3327       B.Add(Sol,Sh);
3328       BRepClass3d_SolidClassifier SC(Sol);
3329       SC.PerformInfinitePoint(1.E-6); // cf. BRepFill_Confusion() - BRepFill_Evolved.cxx
3330       if (SC.State() == TopAbs_IN)
3331         {
3332           B.MakeSolid(Sol);
3333           B.Add(Sol,Sh.Reversed());
3334         }
3335       B.Add(Res,Sol);
3336       ish++;
3337     }
3338   if (ish == 1) { theShape = Sol;}
3339   else          { theShape = Res;}
3340
3341   BRepLib::SameParameter(theShape, 1.E-5, Standard_True);
3342   tds = theShape;
3343   result = CreateObject(tds);
3344   InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
3345   //MESSAGE("---");
3346   return result;
3347 }
3348
3349 //=================================================================================
3350 // function : MakeSewing()
3351 // purpose  :
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)
3356 {
3357   Unexpect aCatch(SALOME_SalomeException);
3358   GEOM::GEOM_Shape_var result ;
3359   TopoDS_Shape tds ;
3360   BRepOffsetAPI_Sewing aMethod ;
3361
3362   try {
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);
3370       }
3371       aMethod.Add(Shape) ;
3372     }
3373
3374     aMethod.Perform() ;
3375     tds = aMethod.SewedShape() ;
3376     if( !IsValid(tds) ) {
3377       THROW_SALOME_CORBA_EXCEPTION("Make Sewing aborted : non valid shape", SALOME::BAD_PARAM);
3378     }
3379     if( tds.IsNull() ) {
3380       THROW_SALOME_CORBA_EXCEPTION("Make Sewing aborted : null shape", SALOME::BAD_PARAM);
3381     }
3382   }
3383   catch (Standard_Failure) {
3384     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeSewing", SALOME::BAD_PARAM);
3385   }
3386
3387   result = CreateObject(tds);
3388   InsertInLabelMoreArguments(tds, result, ListShapes, myCurrentOCAFDoc) ;
3389   return result;
3390 }
3391
3392 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeSewingShape( GEOM::GEOM_Shape_ptr aShape,
3393                                             CORBA::Double precision )
3394   throw (SALOME::SALOME_Exception)
3395 {
3396   Unexpect aCatch(SALOME_SalomeException);
3397   GEOM::GEOM_Shape_var result ;
3398   TopoDS_Shape tds, S ;
3399   BRepOffsetAPI_Sewing aMethod ;
3400
3401   try {
3402     S = GetTopoShape(aShape) ;
3403     if(S.IsNull() ) {
3404       THROW_SALOME_CORBA_EXCEPTION("In Sewing a Shape is null", SALOME::BAD_PARAM);
3405     }
3406
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());
3411       aMethod.Add(F) ;
3412     }
3413
3414     aMethod.Perform() ;
3415     tds = aMethod.SewedShape() ;
3416     if( !IsValid(tds) ) {
3417       THROW_SALOME_CORBA_EXCEPTION("Make Sewing aborted : non valid shape", SALOME::BAD_PARAM);
3418     }
3419   }
3420   catch (Standard_Failure) {
3421     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeSewing", SALOME::BAD_PARAM);
3422   }
3423
3424   result = CreateObject(tds);
3425   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3426   result->ShapeId(entry) ;
3427   return result;
3428 }
3429
3430 //==================================================================================
3431 // function : OrientationChange()
3432 // purpose  : Change the orientation of a new shape
3433 //          : TopAbs_FORWARD < -- > TopAbs_REVERSED
3434 //
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)
3439 {
3440   Unexpect aCatch(SALOME_SalomeException);
3441   GEOM::GEOM_Shape_var result ;
3442   BRep_Builder aBuilder;
3443
3444   TopoDS_Shape shape = GetTopoShape(aShape) ;
3445   if( shape.IsNull() ) {
3446     THROW_SALOME_CORBA_EXCEPTION("Shape is null", SALOME::BAD_PARAM);
3447   }
3448
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);
3454     }
3455
3456     if( tds.Orientation() == TopAbs_FORWARD)
3457       tds.Orientation(TopAbs_REVERSED) ;
3458     else
3459       tds.Orientation(TopAbs_FORWARD) ;
3460
3461     result = CreateObject(tds);
3462     InsertInLabelOneArgument(shape, aShape, tds, result, myCurrentOCAFDoc) ;
3463   }
3464   return result ;
3465 }
3466
3467
3468 //==================================================================================
3469 // function : GetReferencedObjects()
3470 // purpose  :
3471 //==================================================================================
3472 GEOM::GEOM_Gen::ListOfIOR* GEOM_Gen_i::GetReferencedObjects(GEOM::GEOM_Shape_ptr shape)
3473 {
3474   GEOM::GEOM_Gen::ListOfIOR_var aList = new GEOM::GEOM_Gen::ListOfIOR;
3475   aList->length(0);
3476
3477   if (shape->_is_nil()) return aList._retn();
3478
3479   Standard_CString entry = shape->ShapeId();
3480   TDF_Label Lab;
3481   TDF_Tool::Label(myCurrentOCAFDoc->GetData(), entry, Lab);
3482
3483   Handle(TDataStd_Name) Att;
3484   Lab.FindAttribute(TDataStd_Name::GetID(),Att);
3485
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")) ) {
3492
3493       TDF_ChildIterator ChildIterator1(L);
3494       unsigned int i = 0;
3495
3496       while (ChildIterator1.More()) {
3497         TDF_Label L = ChildIterator1.Value();
3498
3499         Handle(TDF_Reference) Ref;
3500         if (L.FindAttribute(TDF_Reference::GetID(),Ref)) {
3501           i++;
3502         }
3503         ChildIterator1.Next();
3504       }
3505       aList->length(i);
3506       i = 0;
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();
3513
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() );
3518           i++;
3519         }
3520
3521         ChildIterator2.Next();
3522       }
3523     }
3524   }
3525   return aList._retn();
3526 }
3527
3528 //==================================================================================
3529 // function : GetObjects()
3530 // purpose  :
3531 //==================================================================================
3532 GEOM::GEOM_Gen::ListOfIOR* GEOM_Gen_i::GetObjects(GEOM::GEOM_Shape_ptr shape)
3533 {
3534   GEOM::GEOM_Gen::ListOfIOR_var aList = new GEOM::GEOM_Gen::ListOfIOR;
3535   aList->length(0);
3536
3537   Standard_CString entry = shape->ShapeId();
3538   TDF_Label Lab;
3539   TDF_Tool::Label(myCurrentOCAFDoc->GetData(), entry, Lab);
3540
3541   Handle(TDataStd_Name) Att;
3542   Lab.FindAttribute(TDataStd_Name::GetID(),Att);
3543
3544   TDF_ChildIterator ChildIterator(Lab);
3545   unsigned int i = 0;
3546   while (ChildIterator.More()) {
3547     TDF_Label L = ChildIterator.Value();
3548     Handle(TDataStd_Name) Att;
3549     L.FindAttribute(TDataStd_Name::GetID(),Att);
3550
3551     if (!Att->Get().IsEqual(TCollection_ExtendedString("Arguments")) ) {
3552       i++;
3553     }
3554     ChildIterator.Next();
3555   }
3556
3557   aList->length(i);
3558   i = 0;
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);
3564
3565     if (!Att->Get().IsEqual(TCollection_ExtendedString("Arguments")) ) {
3566       TCollection_AsciiString nameIOR (Att->Get());
3567       aList[i] = CORBA::string_dup( nameIOR.ToCString() );
3568       i++;
3569     }
3570     ChildIterator1.Next();
3571   }
3572   return aList._retn();
3573 }
3574
3575
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)
3582 {
3583   Unexpect aCatch(SALOME_SalomeException);
3584   TopoDS_Shape tds ;
3585   GEOM::GEOM_Shape_var result ;
3586
3587   try {
3588     BRep_Builder aBuilder;
3589     char* aCopyfilename = strdup(filename);
3590     BRepTools::Read(tds, aCopyfilename, aBuilder) ;
3591     free(aCopyfilename);
3592     if (tds.IsNull()) {
3593       THROW_SALOME_CORBA_EXCEPTION("Import BRep aborted", SALOME::BAD_PARAM);
3594     }
3595   }
3596   catch(Standard_Failure) {
3597     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ImportBREP", SALOME::BAD_PARAM);
3598   }
3599
3600   result = CreateObject(tds) ;
3601   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3602   result->ShapeId(entry);
3603   return result;
3604 }
3605
3606
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)
3615 {
3616   Unexpect aCatch(SALOME_SalomeException);
3617   GEOM::GEOM_Shape_var result ;
3618   TopoDS_Shape tds ;
3619
3620   try {
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) ;
3626   }
3627   catch(Standard_Failure) {
3628     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakePlane", SALOME::BAD_PARAM);
3629   }
3630
3631   if (tds.IsNull()) {
3632     THROW_SALOME_CORBA_EXCEPTION("Make Plane aborted : null shape", SALOME::BAD_PARAM);
3633    }
3634
3635   result = CreateObject(tds) ;
3636   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3637   result->ShapeId(entry);
3638   return result ;
3639 }
3640
3641 //=================================================================================
3642 // function : MakeVertex()
3643 // purpose  : Create a Vertex topology.
3644 //=================================================================================
3645 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeVertex(CORBA::Double x,
3646                                       CORBA::Double y,
3647                                       CORBA::Double z)
3648   throw (SALOME::SALOME_Exception)
3649 {
3650   Unexpect aCatch(SALOME_SalomeException);
3651   GEOM::GEOM_Shape_var result ;
3652   gp_Pnt P(x,y,z);
3653   TopoDS_Shape tds = BRepBuilderAPI_MakeVertex(P).Shape();
3654   if (tds.IsNull()) {
3655     THROW_SALOME_CORBA_EXCEPTION("Make Vertex/Point aborted", SALOME::BAD_PARAM);
3656   }
3657   tds.Infinite(true);
3658   result = CreateObject(tds) ;
3659   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3660   result->ShapeId(entry);
3661   return result ;
3662 }
3663
3664
3665 //=================================================================================
3666 // function : MakeFace()
3667 // purpose  :
3668 //=================================================================================
3669 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeFace( GEOM::GEOM_Shape_ptr wire,
3670                                      CORBA::Boolean wantplanarface )
3671   throw (SALOME::SALOME_Exception)
3672 {
3673   Unexpect aCatch(SALOME_SalomeException);
3674   GEOM::GEOM_Shape_var result ;
3675   TopoDS_Shape aShape;
3676   TopoDS_Shape tds;
3677
3678   try {
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);
3682     }
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) ;
3688     }
3689     else {
3690       THROW_SALOME_CORBA_EXCEPTION("Null result in GEOM_Gen_i::MakeFace", SALOME::BAD_PARAM);
3691     }
3692   }
3693   catch (Standard_Failure) {
3694     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeFace", SALOME::BAD_PARAM);
3695   }
3696   return result ;
3697 }
3698
3699
3700 //=================================================================================
3701 // function : MakeFaces()
3702 // purpose  :
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)
3707 {
3708   GEOM::GEOM_Shape_var result;
3709
3710   try {
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);
3715     }
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);
3722       }
3723       else if(ListShapes.length() >= 2) {
3724         TopoDS_Compound C;
3725         BRep_Builder aBuilder;
3726         aBuilder.MakeCompound(C);
3727         BRepAlgo_FaceRestrictor FR;
3728
3729         TopAbs_Orientation OriF = FFace.Orientation();
3730         TopoDS_Shape aLocalS = FFace.Oriented(TopAbs_FORWARD);
3731         FR.Init(TopoDS::Face(aLocalS), Standard_False, Standard_True);
3732
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);
3738           }
3739           FR.Add(TopoDS::Wire(Shape));
3740         }
3741
3742         FR.Perform();
3743
3744         if(FR.IsDone()) {
3745           int k = 0;
3746           TopoDS_Shape aFace;
3747           for(; FR.More(); FR.Next()) {
3748             aFace = FR.Current().Oriented(OriF);
3749             aBuilder.Add(C, aFace);
3750             k++;
3751           }
3752           if(k == 1) {
3753             result = CreateObject(aFace);
3754             InsertInLabelMoreArguments(aFace, result, ListShapes, myCurrentOCAFDoc);
3755           }
3756           else {
3757             result = CreateObject(C);
3758             InsertInLabelMoreArguments(C, result, ListShapes, myCurrentOCAFDoc);
3759           }
3760         }
3761       }
3762     }
3763     else {
3764       THROW_SALOME_CORBA_EXCEPTION("Null result in GEOM_Gen_i::MakeFace", SALOME::BAD_PARAM);
3765     }
3766   }
3767   catch (Standard_Failure) {
3768     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeFace", SALOME::BAD_PARAM);
3769   }
3770   return result;
3771 }
3772
3773
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)
3780 {
3781   Unexpect aCatch(SALOME_SalomeException);
3782   GEOM::GEOM_Shape_var result ;
3783   BRepTools_Quilt Glue;
3784   TopoDS_Shape C;
3785
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);
3791     }
3792     Glue.Add(Shape) ;
3793   }
3794
3795   TopExp_Explorer exp(Glue.Shells(), TopAbs_SHELL);
3796   Standard_Integer ish = 0;
3797   for (; exp.More(); exp.Next()) {
3798     C = exp.Current();
3799     ish++;
3800   }
3801
3802   if (ish != 1)
3803     C = Glue.Shells();
3804
3805   if ( C.IsNull() ) {
3806     THROW_SALOME_CORBA_EXCEPTION("Null result : Shell operation aborted", SALOME::BAD_PARAM);
3807   }
3808   else {
3809     result = CreateObject(C) ;
3810     InsertInLabelMoreArguments(C, result, ListShapes, myCurrentOCAFDoc) ;
3811   }
3812   return result;
3813 }
3814
3815
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)
3822 {
3823   Unexpect aCatch(SALOME_SalomeException);
3824   GEOM::GEOM_Shape_var result ;
3825   Standard_Integer ish = 0;
3826   TopoDS_Compound  Res;
3827   TopoDS_Solid     Sol;
3828   BRep_Builder     B;
3829   TopoDS_Shape     Shape;
3830
3831   B.MakeCompound(Res);
3832
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) ;
3836     if( Sh.IsNull() ) {
3837        THROW_SALOME_CORBA_EXCEPTION("Solid aborted : null shape during operation", SALOME::BAD_PARAM);
3838     }
3839     B.MakeSolid(Sol);
3840     B.Add(Sol,Sh);
3841     BRepClass3d_SolidClassifier SC(Sol);
3842     SC.PerformInfinitePoint(Precision::Confusion());
3843     if (SC.State() == TopAbs_IN) {
3844       B.MakeSolid(Sol);
3845       B.Add(Sol,Sh.Reversed());
3846     }
3847     B.Add(Res,Sol);
3848     ish++;
3849   }
3850   if (ish == 1) { Shape = Sol;}
3851   else          { Shape = Res;}
3852
3853   if ( Shape.IsNull() ) {
3854     THROW_SALOME_CORBA_EXCEPTION("Null result : Solid operation aborted", SALOME::BAD_PARAM);
3855   }
3856   else {
3857     result = CreateObject(Shape) ;
3858     InsertInLabelMoreArguments(Shape, result, ListShapes, myCurrentOCAFDoc) ;
3859   }
3860   return result;
3861 }
3862
3863
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)
3871 {
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) ;
3876   TopoDS_Shape tds ;
3877
3878   try {
3879     tds = BRepBuilderAPI_MakeEdge(P1, P2).Shape();
3880   }
3881   catch(Standard_Failure) {
3882     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeLine", SALOME::BAD_PARAM);
3883   }
3884
3885   if ( tds.IsNull() ) {
3886     THROW_SALOME_CORBA_EXCEPTION("Make Line aborted : null shape", SALOME::BAD_PARAM);
3887   }
3888   else {
3889     tds.Infinite(true);
3890     result = CreateObject(tds) ;
3891     const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3892     result->ShapeId(entry);
3893   }
3894   return result ;
3895 }
3896
3897
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)
3905 {
3906   Unexpect aCatch(SALOME_SalomeException);
3907   GEOM::GEOM_Shape_var result  ;
3908   TopoDS_Shape tds ;
3909
3910   try {
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();
3914   }
3915   catch(Standard_Failure) {
3916     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeVector", SALOME::BAD_PARAM);
3917   }
3918
3919   if ( tds.IsNull() ) {
3920     THROW_SALOME_CORBA_EXCEPTION("Make Vector aborted : null shape", SALOME::BAD_PARAM);
3921   }
3922   else {
3923     result = CreateObject(tds) ;
3924     const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3925     result->ShapeId(entry);
3926   }
3927   return result ;
3928 }
3929
3930
3931 //================================================================================
3932 // function : MakeCircle()
3933 // purpose  :
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)
3939 {
3940   Unexpect aCatch(SALOME_SalomeException);
3941   GEOM::GEOM_Shape_var result;
3942   TopoDS_Shape tds ;
3943
3944   try {
3945     gp_Pnt p(pstruct.x, pstruct.y, pstruct.z) ;
3946     gp_Dir d(dstruct.PS.x, dstruct.PS.y, dstruct.PS.z) ;
3947     gp_Ax2 axis(p, d) ;
3948     gp_Circ circ( axis, radius);
3949     BRepBuilderAPI_MakeEdge MakeEdge( circ );
3950     tds = MakeEdge.Edge();
3951   }
3952   catch(Standard_Failure) {
3953     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeCircle", SALOME::BAD_PARAM);
3954   }
3955   if (tds.IsNull()) {
3956     THROW_SALOME_CORBA_EXCEPTION("Make Circle aborted", SALOME::BAD_PARAM);
3957   }
3958   result = CreateObject(tds);
3959   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3960   result->ShapeId(entry);
3961   return result ;
3962 }
3963
3964 //================================================================================
3965 // function : MakeEllipse()
3966 // purpose  :
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)
3973 {
3974   Unexpect aCatch(SALOME_SalomeException);
3975   GEOM::GEOM_Shape_var result;
3976   TopoDS_Shape tds ;
3977
3978   try {
3979     gp_Pnt p(pstruct.x, pstruct.y, pstruct.z) ;
3980     gp_Dir d(dstruct.PS.x, dstruct.PS.y, dstruct.PS.z) ;
3981
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();
3986
3987   }
3988   catch(Standard_Failure) {
3989     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeEllipse", SALOME::BAD_PARAM);
3990   }
3991   if (tds.IsNull()) {
3992     THROW_SALOME_CORBA_EXCEPTION("Make Ellipse aborted", SALOME::BAD_PARAM);
3993   }
3994   result = CreateObject(tds);
3995   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3996   result->ShapeId(entry);
3997   return result ;
3998 }
3999
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)
4008 {
4009   Unexpect aCatch(SALOME_SalomeException);
4010   GEOM::GEOM_Shape_var result;
4011   try {
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) ;
4015
4016     GC_MakeArcOfCircle arc( pI, pC, pE ) ;
4017     if( !arc.IsDone() ) {
4018       THROW_SALOME_CORBA_EXCEPTION("Arc not done", SALOME::BAD_PARAM);
4019     }
4020     BRepBuilderAPI_MakeEdge MakeEdge( arc );
4021     TopoDS_Shape tds = MakeEdge.Edge();
4022     if (tds.IsNull()) {
4023       THROW_SALOME_CORBA_EXCEPTION("Null result : arc not done", SALOME::BAD_PARAM);
4024     }
4025     else {
4026       result = CreateObject(tds);
4027       const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
4028       result->ShapeId(entry);
4029     }
4030   }
4031   catch(Standard_Failure) {
4032     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeArc", SALOME::BAD_PARAM);
4033   }
4034   return result ;
4035 }
4036
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)
4043 {
4044   GEOM::GEOM_Shape_var result ;
4045   TopoDS_Shape tds ;
4046   try {
4047     Sketcher_Profile aProfile (Cmd);
4048     if(aProfile.IsDone())
4049       tds = aProfile.GetShape();
4050   }
4051   catch(Standard_Failure) {
4052     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeSketcher", SALOME::BAD_PARAM);
4053   }
4054
4055   if (tds.IsNull()) {
4056     THROW_SALOME_CORBA_EXCEPTION("MakeSketcher aborted : null shape", SALOME::BAD_PARAM);
4057   }
4058   else {
4059     result = CreateObject(tds);
4060     const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
4061     result->ShapeId(entry) ;
4062   }
4063   return result;
4064 }
4065
4066
4067
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)
4074 {
4075   GEOM::GEOM_Shape_var result;
4076   TopoDS_Shape tds, Shape;
4077   TColgp_Array1OfPnt CurvePoints(1, ListShapes.length());
4078
4079   try {
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);
4085       }
4086       if(Shape.ShapeType() == TopAbs_VERTEX) {
4087         const gp_Pnt& P = BRep_Tool::Pnt(TopoDS::Vertex(Shape));
4088         CurvePoints.SetValue(i + 1, P);
4089       }
4090     }
4091     Handle(Geom_BezierCurve) GBC = new Geom_BezierCurve(CurvePoints);
4092     tds = BRepBuilderAPI_MakeEdge(GBC);
4093   }
4094   catch(Standard_Failure) {
4095     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeBezier", SALOME::BAD_PARAM);
4096   }
4097
4098   if( tds.IsNull() ) {
4099     THROW_SALOME_CORBA_EXCEPTION("Make Bezier operation aborted : null result", SALOME::BAD_PARAM);
4100   }
4101   else {
4102     result = CreateObject(tds);
4103     InsertInLabelMoreArguments(tds, result, ListShapes, myCurrentOCAFDoc);
4104   }
4105   return result;
4106 }
4107
4108
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)
4115 {
4116   GEOM::GEOM_Shape_var result;
4117   TopoDS_Shape tds, Shape;
4118   TColgp_Array1OfPnt CurvePoints(1, ListShapes.length());
4119
4120   try {
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);
4126       }
4127       if(Shape.ShapeType() == TopAbs_VERTEX) {
4128         const gp_Pnt& P = BRep_Tool::Pnt(TopoDS::Vertex(Shape));
4129         CurvePoints.SetValue(i + 1, P);
4130       }
4131     }
4132     GeomAPI_PointsToBSpline GBC(CurvePoints);
4133     tds = BRepBuilderAPI_MakeEdge(GBC);
4134   }
4135   catch(Standard_Failure) {
4136     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeBSpline", SALOME::BAD_PARAM);
4137   }
4138
4139   if( tds.IsNull() ) {
4140     THROW_SALOME_CORBA_EXCEPTION("Make BSpline operation aborted : null result", SALOME::BAD_PARAM);
4141   }
4142   else {
4143     result = CreateObject(tds) ;
4144     InsertInLabelMoreArguments(tds, result, ListShapes, myCurrentOCAFDoc) ;
4145   }
4146   return result;
4147 }
4148
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,
4154                                             CORBA::Double x,
4155                                             CORBA::Double y,
4156                                             CORBA::Double z)
4157   throw (SALOME::SALOME_Exception)
4158 {
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);
4164   }
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() ;
4170
4171   result = CreateObject(tds) ;
4172   if( CORBA::is_nil(result) ) {
4173     THROW_SALOME_CORBA_EXCEPTION("Translation aborted : null result", SALOME::BAD_PARAM);
4174   }
4175   InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4176   return result;
4177 }
4178
4179
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,
4186                                                    CORBA::Double step,
4187                                                    CORBA::Short nbtimes )
4188   throw (SALOME::SALOME_Exception)
4189 {
4190   Unexpect aCatch(SALOME_SalomeException);
4191   GEOM::GEOM_Shape_var result ;
4192   TopoDS_Shape tds ;
4193
4194   TopoDS_Shape aShape = GetTopoShape(myShape) ;
4195   if( aShape.IsNull() )
4196     THROW_SALOME_CORBA_EXCEPTION("MakeMultiTranslation1D aborted : null shape", SALOME::BAD_PARAM);
4197
4198   try {
4199     int i ;
4200     double DX, DY, DZ ;
4201     gp_Trsf theTransformation ;
4202     gp_Vec myVec ;
4203     gp_Vec Vec( dir.PS.x, dir.PS.y, dir.PS.z ) ;
4204     Vec.Normalize();
4205     TopoDS_Compound compound;
4206     BRep_Builder B;
4207     B.MakeCompound( compound );
4208
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() );
4217     }
4218     tds = compound ;
4219     result = CreateObject(tds) ;
4220   }
4221   catch (Standard_Failure) {
4222     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeMultiTranslation1D", SALOME::BAD_PARAM);
4223   }
4224
4225   if( CORBA::is_nil(result) ) {
4226     THROW_SALOME_CORBA_EXCEPTION("MakeMultiTranslation1D aborted : null result", SALOME::BAD_PARAM);
4227   }
4228   InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4229   return result;
4230 }
4231
4232
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)
4245 {
4246   Unexpect aCatch(SALOME_SalomeException);
4247   GEOM::GEOM_Shape_var result ;
4248   TopoDS_Shape tds ;
4249
4250   TopoDS_Shape aShape = GetTopoShape(myShape) ;
4251   if( aShape.IsNull() ) {
4252     THROW_SALOME_CORBA_EXCEPTION("MakeMultiTranslation2D aborted : null shape", SALOME::BAD_PARAM);
4253   }
4254
4255   try {
4256     int i, j ;
4257     double DX, DY, DZ ;
4258     gp_Trsf theTransformation ;
4259     gp_Vec myVec ;
4260     gp_Vec Vec1( dir1.PS.x, dir1.PS.y, dir1.PS.z ) ;
4261     Vec1.Normalize();
4262     gp_Vec Vec2( dir2.PS.x, dir2.PS.y, dir2.PS.z ) ;
4263     Vec2.Normalize();
4264     TopoDS_Compound compound;
4265     BRep_Builder B;
4266     B.MakeCompound( compound );
4267
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() );
4277       }
4278     }
4279     tds = compound ;
4280     result = CreateObject(tds) ;
4281   }
4282   catch(Standard_Failure) {
4283     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeMultiTranslation2D", SALOME::BAD_PARAM);
4284   }
4285
4286   if( CORBA::is_nil(result) ) {
4287     THROW_SALOME_CORBA_EXCEPTION("MakeMultiTranslation2D aborted : null result", SALOME::BAD_PARAM);
4288   }
4289   InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4290   return result;
4291 }
4292
4293
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)
4303 {
4304   Unexpect aCatch(SALOME_SalomeException);
4305   GEOM::GEOM_Shape_var result ;
4306   TopoDS_Shape tds ;
4307   TopoDS_Shape aShape = GetTopoShape(myShape) ;
4308   if( aShape.IsNull() ) {
4309     THROW_SALOME_CORBA_EXCEPTION("MakeMultiRotation1D aborted : null shape", SALOME::BAD_PARAM);
4310   }
4311
4312   try {
4313
4314     int i ;
4315     gp_Pnt P(loc.x, loc.y, loc.z) ;
4316     gp_Dir D(dir.PS.x, dir.PS.y, dir.PS.z) ;
4317     gp_Ax1 AX1(P, D) ;
4318
4319     double angle = 360.0/nbtimes ;
4320     gp_Trsf theTransformation ;
4321     TopoDS_Compound compound;
4322     BRep_Builder B;
4323     B.MakeCompound( compound );
4324
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() );
4329     }
4330     tds = compound ;
4331     result = CreateObject(tds) ;
4332   }
4333   catch(Standard_Failure) {
4334     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeMultiRotation1D", SALOME::BAD_PARAM);
4335   }
4336
4337   if( CORBA::is_nil(result) ) {
4338     THROW_SALOME_CORBA_EXCEPTION("MakeMultiRotation1D aborted : null result", SALOME::BAD_PARAM);
4339   }
4340   InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4341   return result;
4342 }
4343
4344
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,
4352                                                 CORBA::Double ang,
4353                                                 CORBA::Short nbtimes1,
4354                                                 CORBA::Double step,
4355                                                 CORBA::Short nbtimes2 )
4356   throw (SALOME::SALOME_Exception)
4357 {
4358   Unexpect aCatch(SALOME_SalomeException);
4359   GEOM::GEOM_Shape_var result ;
4360   TopoDS_Shape tds ;
4361   TopoDS_Shape aShape = GetTopoShape(myShape) ;
4362   if( aShape.IsNull() ) {
4363     THROW_SALOME_CORBA_EXCEPTION("MakeMultiRotation2D aborted : null shape", SALOME::BAD_PARAM);
4364   }
4365
4366   try {
4367
4368     int i, j ;
4369     double DX, DY, DZ ;
4370     gp_Pnt P(loc.x, loc.y, loc.z) ;
4371     gp_Dir D(dir.PS.x, dir.PS.y, dir.PS.z) ;
4372     gp_Ax1 AX1(P, D) ;
4373     gp_Trsf theTransformation1 ;
4374     gp_Trsf theTransformation2 ;
4375     gp_Pnt P1 ;
4376     GProp_GProps System ;
4377
4378     if ( aShape.ShapeType() == TopAbs_VERTEX) {
4379       P1 = BRep_Tool::Pnt(TopoDS::Vertex( aShape ));
4380     }
4381     else if ( aShape.ShapeType() == TopAbs_EDGE || aShape.ShapeType() == TopAbs_WIRE ) {
4382       BRepGProp::LinearProperties(aShape, System);
4383       P1 = System.CentreOfMass() ;
4384     }
4385     else if ( aShape.ShapeType() == TopAbs_FACE || aShape.ShapeType() == TopAbs_SHELL ) {
4386       BRepGProp::SurfaceProperties(aShape, System);
4387       P1 = System.CentreOfMass() ;
4388     }
4389     else {
4390       BRepGProp::VolumeProperties(aShape, System);
4391       P1 = System.CentreOfMass() ;
4392     }
4393
4394     Handle(Geom_Line) Line = new Geom_Line(AX1);
4395     gp_Pnt P2 = GeomAPI_ProjectPointOnCurve( P1, Line ) ;
4396
4397     if ( P1.IsEqual(P2, Precision::Confusion() ) )
4398       THROW_SALOME_CORBA_EXCEPTION("Points are confused", SALOME::BAD_PARAM);
4399
4400     gp_Vec Vec(P1.X()-P2.X(), P1.Y()-P2.Y(), P1.Z()-P2.Z()) ;
4401     Vec.Normalize();
4402
4403     gp_Vec myVec ;
4404     TopoDS_Compound compound;
4405     BRep_Builder B;
4406     B.MakeCompound( compound );
4407
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() );
4419       }
4420     }
4421     tds = compound ;
4422     result = CreateObject(tds) ;
4423   }
4424   catch(Standard_Failure) {
4425     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeMultiRotation2D", SALOME::BAD_PARAM);
4426   }
4427
4428   if( CORBA::is_nil(result) ) {
4429     THROW_SALOME_CORBA_EXCEPTION("MakeMultiRotation2D aborted : null result", SALOME::BAD_PARAM);
4430   }
4431   InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4432   return result;
4433 }
4434
4435
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)
4442 {
4443   Unexpect aCatch(SALOME_SalomeException);
4444   GEOM::GEOM_Shape_var result ;
4445   TopoDS_Shape tds ;
4446   TopoDS_Shape aShape = GetTopoShape(Shape) ;
4447   if(aShape.IsNull() ) {
4448     THROW_SALOME_CORBA_EXCEPTION("Copy aborted : null shape during operation", SALOME::BAD_PARAM);
4449   }
4450   BRepBuilderAPI_Copy Copy(aShape);
4451   if( Copy.IsDone() ) {
4452     tds = Copy.Shape();
4453     result = CreateObject(tds);
4454     InsertInLabelOneArgument(aShape, Shape, tds, result, myCurrentOCAFDoc) ;
4455   }
4456
4457   return result;
4458 }
4459
4460
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)
4468 {
4469   Unexpect aCatch(SALOME_SalomeException);
4470   GEOM::GEOM_Shape_var result ;
4471   TopoDS_Shape tds ;
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);
4476   }
4477
4478   try {
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) */
4484
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) ;
4490
4491     tds = myBRepTransformation.Shape() ;
4492     if(tds.IsNull() ) {
4493       THROW_SALOME_CORBA_EXCEPTION("Mirror aborted", SALOME::BAD_PARAM);
4494     }
4495   }
4496   catch(Standard_Failure) {
4497     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeMirrorByPlane", SALOME::BAD_PARAM);
4498   }
4499
4500   result = CreateObject(tds) ;
4501
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) ;
4508   return result ;
4509 }
4510
4511
4512
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)
4521 {
4522   Unexpect aCatch(SALOME_SalomeException);
4523   GEOM::GEOM_Shape_var result ;
4524   TopoDS_Shape tds ;
4525   TopoDS_Shape aShape = GetTopoShape(myShape) ;
4526   if( aShape.IsNull() ) {
4527     THROW_SALOME_CORBA_EXCEPTION("Rotation aborted : null shape during operation", SALOME::BAD_PARAM);
4528   }
4529
4530   try {
4531     gp_Pnt P(axis.x, axis.y, axis.z) ;
4532     gp_Dir D(axis.vx, axis.vy, axis.vz) ;
4533     gp_Ax1 AX(P, D) ;
4534
4535     gp_Trsf theTransformation ;
4536     theTransformation.SetRotation(AX, angle) ;
4537     BRepBuilderAPI_Transform myBRepTransformation(aShape, theTransformation, Standard_False) ;
4538     tds = myBRepTransformation.Shape() ;
4539   }
4540   catch(Standard_Failure) {
4541     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeRotation", SALOME::BAD_PARAM);
4542   }
4543
4544   if ( !tds.IsNull() ) {
4545     result = CreateObject(tds) ;
4546     InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4547   }
4548   return result ;
4549 }
4550
4551
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)
4560 {
4561   Unexpect aCatch(SALOME_SalomeException);
4562   GEOM::GEOM_Shape_var result ;
4563   TopoDS_Shape tds ;
4564   TopoDS_Shape aShape = GetTopoShape(myShape) ;
4565   if( aShape.IsNull() ) {
4566     THROW_SALOME_CORBA_EXCEPTION("Scale aborted : null shape during operation", SALOME::BAD_PARAM);
4567   }
4568
4569   try {
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() ;
4575   }
4576   catch(Standard_Failure) {
4577     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeScaleTransform", SALOME::BAD_PARAM);
4578   }
4579
4580   if ( !tds.IsNull() ) {
4581     result = CreateObject(tds) ;
4582     InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4583   }
4584   return result ;
4585 }
4586
4587
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)
4594 {
4595   Unexpect aCatch(SALOME_SalomeException);
4596   GEOM::GEOM_Shape_var result ;
4597   TopoDS_Compound C;
4598   BRep_Builder aBuilder;
4599   aBuilder.MakeCompound(C) ;
4600
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);
4606     }
4607     aBuilder.Add(C, Shape) ;
4608   }
4609
4610   if ( C.IsNull() ) {
4611     THROW_SALOME_CORBA_EXCEPTION("Null result : Compound operation aborted", SALOME::BAD_PARAM);
4612   }
4613   else {
4614     result = CreateObject(C) ;
4615     InsertInLabelMoreArguments(C, result, ListShapes, myCurrentOCAFDoc) ;
4616   }
4617   return result;
4618 }
4619
4620
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)
4628 {
4629   Unexpect aCatch(SALOME_SalomeException);
4630   GEOM::GEOM_Shape_var result  ;
4631   TopoDS_Shape tds ;
4632
4633   try {
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();
4637     if ( tds.IsNull() )
4638       THROW_SALOME_CORBA_EXCEPTION("MakeEdge aborted : null result", SALOME::BAD_PARAM);
4639   }
4640   catch (Standard_Failure) {
4641     THROW_SALOME_CORBA_EXCEPTION("Exception catched in MakeEdge", SALOME::BAD_PARAM);
4642   }
4643
4644   result = CreateObject(tds) ;
4645   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
4646   result->ShapeId(entry);
4647   return result ;
4648 }
4649
4650
4651
4652 //=================================================================================
4653 // function : MakeWire()
4654 // purpose  : Make a wire from a list containing one or more edges or wires that can
4655 // be connected
4656 //=================================================================================
4657 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeWire( const GEOM::GEOM_Gen::ListOfIOR& ListShapes )
4658   throw (SALOME::SALOME_Exception)
4659 {
4660   Unexpect aCatch(SALOME_SalomeException);
4661   GEOM::GEOM_Shape_var result ;
4662   BRepBuilderAPI_MakeWire MW ;
4663   TopoDS_Shape tds, Shape ;
4664
4665   try {
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);
4671       }
4672       if( Shape.ShapeType() == TopAbs_EDGE )
4673         MW.Add( TopoDS::Edge(Shape) ) ;
4674       if (Shape.ShapeType() == TopAbs_WIRE )
4675         MW.Add( TopoDS::Wire(Shape) ) ;
4676     }
4677     tds = MW  ;
4678
4679   }
4680   catch(Standard_Failure) {
4681     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeWire", SALOME::BAD_PARAM);
4682   }
4683
4684   if( tds.IsNull() ) {
4685     THROW_SALOME_CORBA_EXCEPTION("Make Wire operation aborted : null result", SALOME::BAD_PARAM);
4686   }
4687   else {
4688     result = CreateObject(tds) ;
4689     InsertInLabelMoreArguments(tds, result, ListShapes, myCurrentOCAFDoc) ;
4690   }
4691   return result;
4692 }
4693
4694
4695 //=================================================================================
4696 // function : MakeRevolution()
4697 // purpose  :
4698 //=================================================================================
4699 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeRevolution(GEOM::GEOM_Shape_ptr myShape,
4700                                           const GEOM::AxisStruct& axis,
4701                                           double angle)
4702   throw (SALOME::SALOME_Exception)
4703 {
4704   Unexpect aCatch(SALOME_SalomeException);
4705   GEOM::GEOM_Shape_var result ;
4706   TopoDS_Shape tds ;
4707   TopoDS_Shape aShape = GetTopoShape(myShape) ;
4708   if( aShape.IsNull() ) {
4709     THROW_SALOME_CORBA_EXCEPTION("Revolution aborted : null shape", SALOME::BAD_PARAM);
4710   }
4711   try {
4712     gp_Pnt P(axis.x, axis.y, axis.z) ;
4713     gp_Dir D(axis.vx, axis.vy, axis.vz);
4714     gp_Ax1 AX(P,D);
4715     tds = BRepPrimAPI_MakeRevol(aShape, AX, angle);
4716   }
4717   catch(Standard_Failure) {
4718     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeRevolution", SALOME::BAD_PARAM);
4719   }
4720
4721   if( tds.IsNull() ) {
4722     THROW_SALOME_CORBA_EXCEPTION("Revolution aborted", SALOME::BAD_PARAM);
4723   }
4724   result = CreateObject(tds) ;
4725   InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4726   return result ;
4727 }
4728
4729
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)
4737 {
4738   Unexpect aCatch(SALOME_SalomeException);
4739   GEOM::GEOM_Shape_var result ;
4740   TopoDS_Shape tds ;
4741   TopoDS_Wire aWire ;
4742   TopoDS_Shape pathTds = GetTopoShape(pathShape) ;
4743   TopoDS_Shape baseTds = GetTopoShape(baseShape) ;
4744
4745   if( baseTds.IsNull() || pathTds.IsNull() ) {
4746     THROW_SALOME_CORBA_EXCEPTION("MakePipe aborted : null shape argument", SALOME::BAD_PARAM);
4747   }
4748
4749   if( pathTds.ShapeType() == TopAbs_WIRE ) {
4750     aWire = TopoDS::Wire(pathTds) ;
4751   }
4752   else {
4753     if ( pathTds.ShapeType() == TopAbs_EDGE ) {
4754       TopoDS_Edge aEdge = TopoDS::Edge(pathTds) ;
4755       aWire = BRepBuilderAPI_MakeWire(aEdge);
4756     }
4757     else {
4758       THROW_SALOME_CORBA_EXCEPTION("MakePipe aborted : bad shape type", SALOME::BAD_PARAM);
4759     }
4760   }
4761
4762   try {
4763     tds = BRepOffsetAPI_MakePipe(aWire, baseTds) ;
4764   }
4765   catch(Standard_Failure) {
4766     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakePipe", SALOME::BAD_PARAM);
4767   }
4768
4769   if (  !IsValid(tds) ) {
4770     THROW_SALOME_CORBA_EXCEPTION("MakePipe aborted : non valid shape result", SALOME::BAD_PARAM);
4771   }
4772   else {
4773     result = CreateObject(tds) ;
4774
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) ;
4781   }
4782   return result ;
4783 }
4784
4785
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)
4794 {
4795   Unexpect aCatch(SALOME_SalomeException);
4796   GEOM::GEOM_Shape_var result ;
4797   TopoDS_Shape tds ;
4798   TopoDS_Shape aShape = GetTopoShape(myShape) ;
4799   if( aShape.IsNull() ) {
4800     THROW_SALOME_CORBA_EXCEPTION("Prism aborted : null shape operation", SALOME::BAD_PARAM);
4801   }
4802
4803   try {
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() ;
4806   }
4807   catch(Standard_Failure) {
4808     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakePipe", SALOME::BAD_PARAM);
4809   }
4810
4811   if ( tds.IsNull() ) {
4812     THROW_SALOME_CORBA_EXCEPTION("Prism aborted", SALOME::BAD_PARAM);
4813   }
4814   else {
4815     result = CreateObject(tds) ;
4816     InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4817   }
4818   return result ;
4819 }
4820
4821
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)
4828 {
4829   Unexpect aCatch(SALOME_SalomeException);
4830   GEOM::GEOM_Shape_var result ;
4831   TopoDS_Shape tds ;
4832   TopoDS_Shape shape = GetTopoShape(aShape) ;
4833   GProp_GProps System;
4834   gp_Pnt myCenterMass ;
4835
4836   if( shape.IsNull() ) {
4837     THROW_SALOME_CORBA_EXCEPTION("MakeCDG aborted : null shape argument", SALOME::BAD_PARAM);
4838   }
4839
4840   try {
4841     if ( shape.ShapeType() == TopAbs_VERTEX) {
4842       myCenterMass = BRep_Tool::Pnt(TopoDS::Vertex( shape ));
4843     }
4844     else if ( shape.ShapeType() == TopAbs_EDGE || shape.ShapeType() == TopAbs_WIRE ) {
4845       BRepGProp::LinearProperties(shape, System);
4846       myCenterMass = System.CentreOfMass() ;
4847     }
4848     else if ( shape.ShapeType() == TopAbs_FACE || shape.ShapeType() == TopAbs_SHELL ) {
4849       BRepGProp::SurfaceProperties(shape, System);
4850       myCenterMass = System.CentreOfMass() ;
4851     }
4852     else {
4853       BRepGProp::VolumeProperties(shape, System);
4854       myCenterMass = System.CentreOfMass() ;
4855     }
4856
4857     tds = BRepBuilderAPI_MakeVertex(myCenterMass).Shape() ;
4858   }
4859   catch(Standard_Failure) {
4860     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeCDG", SALOME::BAD_PARAM);
4861   }
4862
4863   if ( tds.IsNull() ) {
4864     THROW_SALOME_CORBA_EXCEPTION("Make CDG aborted : null shape result", SALOME::BAD_PARAM);
4865   }
4866   else {
4867     result = CreateObject(tds) ;
4868     InsertInLabelOneArgument(shape, aShape, tds, result, myCurrentOCAFDoc) ;
4869   }
4870   return result ;
4871 }
4872
4873
4874 //=================================================================================
4875 // function : Archimede()
4876 // purpose  :
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)
4883 {
4884   Unexpect aCatch(SALOME_SalomeException);
4885   GEOM::GEOM_Shape_var result;
4886
4887   double cste = -1;
4888   if (aWaterDensity != 0.)
4889     cste = aWeight/aWaterDensity;
4890   else
4891     THROW_SALOME_CORBA_EXCEPTION("Water density is null", SALOME::BAD_PARAM);
4892
4893   TopoDS_Shape shape = GetTopoShape(aShape) ;
4894   if( shape.IsNull() ) {
4895     THROW_SALOME_CORBA_EXCEPTION("Shape is null", SALOME::BAD_PARAM);
4896   }
4897
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);
4903
4904   gp_Dir Zdirection(0.0,0.0,1.0);
4905   VolumeSection VOL( shape, aMeshingDeflection);
4906   VOL.SetPlane(P);
4907   Handle (Geom_RectangularTrimmedSurface) SurfaceTrimmee;
4908
4909   if(Zdirection.IsEqual(direct,Precision::Angular()) == Standard_False) {
4910     VOL.MakeRotation(direct);
4911   }
4912
4913   VOL.CenterOfGravity();
4914   SurfaceTrimmee = VOL.TrimSurf();
4915   Standard_Real Cote = VOL.Archimede( cste, aMeshingDeflection );
4916
4917   if ( Cote == -1 ) {
4918     double Zmin,Zmax;
4919     VOL.getZ(Zmin,Zmax);
4920     double volume = VOL.CalculateVolume( Zmax ) * aWaterDensity;
4921
4922     char msg[100]="";
4923     sprintf(msg, "shape sinks to the bottom : Weigth max = %.1f", volume);
4924
4925     THROW_SALOME_CORBA_EXCEPTION(msg, SALOME::BAD_PARAM);
4926   }
4927
4928   SurfaceTrimmee=VOL.AjustePlan(SurfaceTrimmee,Cote,PosPlan);
4929   if(Zdirection.IsEqual(direct,Precision::Angular()) == Standard_False) {
4930     SurfaceTrimmee=VOL.InvMakeRotation(direct,SurfaceTrimmee);
4931   }
4932
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);
4936
4937   if (tirant.IsNull()) {
4938     THROW_SALOME_CORBA_EXCEPTION("Result is null", SALOME::BAD_PARAM);
4939   }
4940
4941   result = CreateObject(tirant);
4942   InsertInLabelOneArgument(shape, aShape, tirant, result, myCurrentOCAFDoc) ;
4943
4944   return result;
4945 }
4946
4947
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)
4957 {
4958   Unexpect aCatch(SALOME_SalomeException);
4959   GEOM::GEOM_Shape_var result;
4960   TopoDS_Shape tds ;
4961
4962   const TopoDS_Shape aShape = GetTopoShape(shape) ;
4963   if( aShape.IsNull() ) {
4964     THROW_SALOME_CORBA_EXCEPTION("Shape is null", SALOME::BAD_PARAM);
4965   }
4966
4967   BRepFilletAPI_MakeFillet fill(aShape);
4968
4969   try {
4970     /* case all */
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());
4975         fill.Add(E);
4976       }
4977       for (int i = 1;i<=fill.NbContours();i++) {
4978 #if OCC_VERSION_MAJOR >= 5
4979         fill.SetRadius(radius,i,i);
4980 #else
4981         fill.SetRadius(radius,i);
4982 #endif
4983       }
4984       tds = fill.Shape();
4985
4986     } else {
4987
4988       /* case selection */
4989       for ( unsigned int ind = 0; ind < ListOfID.length(); ind++ ) {
4990         TopoDS_Shape ss ;
4991         if( GetShapeFromIndex( aShape, (TopAbs_ShapeEnum)ShapeType, ListOfID[ind], ss ) ) {
4992           TopoDS_Edge E = TopoDS::Edge(ss) ;
4993           fill.Add( E );
4994         }
4995       }
4996       for (int i = 1;i<=fill.NbContours();i++) {
4997 #if OCC_VERSION_MAJOR >= 5
4998         fill.SetRadius(radius,i,i);
4999 #else
5000         fill.SetRadius(radius,i);
5001 #endif
5002       }
5003       tds = fill.Shape();
5004     }
5005   }
5006   catch(Standard_Failure) {
5007     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeFillet", SALOME::BAD_PARAM);
5008   }
5009
5010   if (tds.IsNull()) {
5011     THROW_SALOME_CORBA_EXCEPTION("Make Fillet aborted", SALOME::BAD_PARAM);
5012   }
5013   result = CreateObject(tds);
5014   InsertInLabelOneArgument(aShape, shape, tds, result, myCurrentOCAFDoc) ;
5015
5016   return result ;
5017 }
5018
5019
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,
5025                                               CORBA::Double d1,
5026                                               CORBA::Double d2,
5027                                               CORBA::Short ShapeType,
5028                                               const GEOM::GEOM_Shape::ListOfSubShapeID& ListOfID )
5029   throw (SALOME::SALOME_Exception)
5030 {
5031   Unexpect aCatch(SALOME_SalomeException);
5032   GEOM::GEOM_Shape_var result;
5033   TopoDS_Shape tds ;
5034
5035   const TopoDS_Shape aShape = GetTopoShape(shape) ;
5036   if( aShape.IsNull() ) {
5037     THROW_SALOME_CORBA_EXCEPTION("Shape is null", SALOME::BAD_PARAM);
5038   }
5039
5040   BRepFilletAPI_MakeChamfer MC(aShape);
5041
5042   try {
5043     /* case all */
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))
5051           MC.Add(d1,d2,E,F);
5052       }
5053       tds = MC.Shape();
5054
5055     } else {
5056
5057       /* case selection */
5058       for ( unsigned int ind = 0; ind < ListOfID.length(); ind++ ) {
5059         TopoDS_Shape ss ;
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))
5064             MC.Add(d1,d2,E,F);
5065         }
5066       }
5067       tds = MC.Shape();
5068     }
5069   }
5070   catch(Standard_Failure) {
5071     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeChamfer", SALOME::BAD_PARAM);
5072   }
5073
5074   if (tds.IsNull()) {
5075     THROW_SALOME_CORBA_EXCEPTION("Make Chamfer aborted", SALOME::BAD_PARAM);
5076   }
5077   result = CreateObject(tds);
5078   InsertInLabelOneArgument(aShape, shape, tds, result, myCurrentOCAFDoc) ;
5079
5080   return result ;
5081 }
5082
5083 //=================================================================================
5084 // function : CheckShape()
5085 // purpose  :
5086 //=================================================================================
5087 CORBA::Boolean GEOM_Gen_i::CheckShape(GEOM::GEOM_Shape_ptr shape)
5088   throw (SALOME::SALOME_Exception)
5089 {
5090   Unexpect aCatch(SALOME_SalomeException);
5091   TopoDS_Shape S = GetTopoShape(shape) ;
5092   if( S.IsNull() ) {
5093     THROW_SALOME_CORBA_EXCEPTION("Shape is null", SALOME::BAD_PARAM);
5094   }
5095
5096   BRepCheck_Analyzer ana(S,false);
5097   if (ana.IsValid())
5098     return 1;
5099
5100   return 0;
5101 }
5102
5103 //=================================================================================
5104 // function : MakePlacedBox()
5105 // purpose  :
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)
5110 {
5111   Unexpect aCatch(SALOME_SalomeException);
5112   GEOM::GEOM_Shape_var result ;
5113   TopoDS_Shape tds ;
5114
5115   CORBA::Double x2, y2, z2 ;
5116
5117   try {
5118     x2 = x1 + delta1 ;
5119     y2 = y1 + delta2 ;
5120     z2 = z1 + delta3 ;
5121
5122     gp_Pnt P1(x1,y1,z1);
5123     gp_Pnt P2(x2,y2,z2);
5124
5125     tds = BRepPrimAPI_MakeBox(P1,P2).Shape();
5126   }
5127   catch(Standard_Failure) {
5128     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeBox", SALOME::BAD_PARAM);
5129   }
5130
5131   if (tds.IsNull()) {
5132     THROW_SALOME_CORBA_EXCEPTION("Make Box aborted : null shape", SALOME::BAD_PARAM);
5133   }
5134
5135   result = CreateObject(tds);
5136   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
5137   result->ShapeId(entry) ;
5138
5139   return result;
5140 }
5141
5142 //=================================================================================
5143 // function : MakePanel()
5144 // purpose  :
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)
5150 {
5151   Unexpect aCatch(SALOME_SalomeException);
5152   GEOM::GEOM_Shape_var result ;
5153   TopoDS_Shape tds ;
5154   TopoDS_Shape aShape = GetTopoShape(shape) ;
5155   Bnd_Box B ;
5156   Standard_Real axmin,aymin,azmin,axmax,aymax,azmax ;
5157   GEOM::PointStruct pstruct1, pstruct2, pstruct3, pstruct4 ;
5158
5159   if(aShape.IsNull() ) {
5160     THROW_SALOME_CORBA_EXCEPTION("MakePanel aborted : null shape during operation", SALOME::BAD_PARAM);
5161   }
5162
5163   try {
5164     BRepBndLib::Add(aShape,B);
5165     B.Enlarge(10.);
5166     B.Get(axmin,aymin,azmin,axmax,aymax,azmax);
5167
5168     switch (directiontype)
5169       {
5170       case 1 :   /* X */
5171         pstruct1 = MakePointStruct( delta, aymin, azmin ) ;
5172         pstruct2 = MakePointStruct( delta, aymin, azmax ) ;
5173         pstruct3 = MakePointStruct( delta, aymax, azmax ) ;
5174         pstruct4 = MakePointStruct( delta, aymax, azmin ) ;
5175         break ;
5176       case 2 :   /* Y */
5177         pstruct1 = MakePointStruct( axmin, delta, azmin ) ;
5178         pstruct2 = MakePointStruct( axmin, delta, azmax ) ;
5179         pstruct3 = MakePointStruct( axmax, delta, azmax ) ;
5180         pstruct4 = MakePointStruct( axmax, delta, azmin ) ;
5181         break ;
5182       case 3 :   /* Z */
5183         pstruct1 = MakePointStruct( axmin, aymin, delta ) ;
5184         pstruct2 = MakePointStruct( axmin, aymax, delta ) ;
5185         pstruct3 = MakePointStruct( axmax, aymax, delta ) ;
5186         pstruct4 = MakePointStruct( axmax, aymin, delta ) ;
5187         break ;
5188       default :
5189         return result ;
5190       }
5191
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);
5196
5197     GEOM::GEOM_Gen::ListOfIOR_var aList = new GEOM::GEOM_Gen::ListOfIOR;
5198     aList->length(4);
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());
5203
5204     GEOM::GEOM_Shape_ptr aWire = MakeWire( aList );
5205     GEOM::GEOM_Shape_ptr aFace = MakeFace( aWire, true ) ;
5206     tds = GetTopoShape(aFace);
5207
5208   }
5209   catch(Standard_Failure) {
5210     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakePanel", SALOME::BAD_PARAM);
5211   }
5212
5213   if (tds.IsNull()) {
5214     THROW_SALOME_CORBA_EXCEPTION("Make PanelsPartition aborted : null shape", SALOME::BAD_PARAM);
5215   }
5216
5217   result = CreateObject(tds);
5218   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
5219   result->ShapeId(entry) ;
5220
5221   return result;
5222 }
5223
5224
5225 void GEOM_Gen_i::ExportIGES(const char* filename,GEOM::GEOM_Shape_ptr theShape)
5226   throw (SALOME::SALOME_Exception)
5227 {
5228   Unexpect aCatch(SALOME_SalomeException);
5229   if (theShape->_is_nil())
5230     {
5231       THROW_SALOME_CORBA_EXCEPTION("Export IGES aborted", SALOME::BAD_PARAM);
5232     }
5233   TopoDS_Shape tds = GetTopoShape(theShape);
5234   if (tds.IsNull())
5235     {
5236       THROW_SALOME_CORBA_EXCEPTION("Export IGES aborted", SALOME::BAD_PARAM);
5237     }
5238   try
5239     {
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
5245
5246       ICW.AddShape (tds);
5247       ICW.ComputeModel();
5248       char * aname = strdup(filename);
5249       Standard_Boolean result = ICW.Write( aname );
5250       free(aname);
5251     }
5252   catch(Standard_Failure)
5253     {
5254       THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ExportIGES", SALOME::BAD_PARAM);
5255     }
5256 }
5257
5258 void GEOM_Gen_i::ExportBREP(const char* filename,GEOM::GEOM_Shape_ptr theShape)
5259   throw (SALOME::SALOME_Exception)
5260 {
5261   Unexpect aCatch(SALOME_SalomeException);
5262   if (theShape->_is_nil())
5263     {
5264       THROW_SALOME_CORBA_EXCEPTION("Export BRep aborted", SALOME::BAD_PARAM);
5265     }
5266   TopoDS_Shape tds = GetTopoShape(theShape);
5267   if (tds.IsNull())
5268     {
5269       THROW_SALOME_CORBA_EXCEPTION("Export BRep aborted", SALOME::BAD_PARAM);
5270     }
5271   try
5272     {
5273       char * aname = strdup(filename);
5274       Standard_Boolean result = BRepTools::Write(tds,aname);
5275       free(aname);
5276     }
5277   catch(Standard_Failure)
5278     {
5279       THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ExportBREP", SALOME::BAD_PARAM);
5280     }
5281 }
5282
5283 void GEOM_Gen_i::ExportSTEP(const char* filename,GEOM::GEOM_Shape_ptr theShape)
5284   throw (SALOME::SALOME_Exception)
5285 {
5286   Unexpect aCatch(SALOME_SalomeException);
5287   if (theShape->_is_nil())
5288     {
5289       THROW_SALOME_CORBA_EXCEPTION("Export STEP aborted", SALOME::BAD_PARAM);
5290     }
5291   TopoDS_Shape tds = GetTopoShape(theShape);
5292   if (tds.IsNull())
5293     {
5294       THROW_SALOME_CORBA_EXCEPTION("Export STEP aborted", SALOME::BAD_PARAM);
5295     }
5296   try
5297     {
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 )
5304         {
5305           char * aname = strdup(filename);
5306           status = aWriter.Write( aname ) ;
5307           free(aname);
5308         }
5309     }
5310   catch(Standard_Failure)
5311     {
5312       THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ExportBREP", SALOME::BAD_PARAM);
5313     }
5314 }
5315
5316
5317 //=====================================================================================
5318 // EXPORTED METHODS
5319 //=====================================================================================
5320 extern "C"
5321 {
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)
5327   {
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() ;
5331   }
5332 }
5333