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.opencascade.org/SALOME/ or email : webmaster.salome@opencascade.org
  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 CORBA::string_dup(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):SALOMEDS_Tool::GetTmpDir();
 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(SALOMEDS_Tool::GetNameFromPath(theComponent->GetStudy()->URL()));
 469   aNameWithExt += TCollection_AsciiString("_GEOM.sgd");
 470   aSeq[0] = CORBA::string_dup(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):SALOMEDS_Tool::GetTmpDir();
 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(SALOMEDS_Tool::GetNameFromPath(theComponent->GetStudy()->URL()));
 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 //        cout<<"********** Nullify owner of document"<<endl;
 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 //    cout<<"********** GEOM_Gen_i::CanPaste ("<<theComponentName<<","<<theObjectID<<")"<<endl;
 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 //      cout<<"GEOM_Gen_i::PasteInto exception"<<endl;
 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     //OCC 5.1.2 porting
2739 //     aReader.Clear();
2740 //     aReader.TransferRoots(false);
2741     aReader.ClearShapes();
2742     aReader.TransferRoots();
2743
2744     MESSAGE("ImportIGES : count of shapes produced = " << aReader.NbShapes() << endl );
2745     TopoDS_Shape shape = aReader.OneShape();
2746
2747     if ( !shape.IsNull() ) {
2748       /* Final CORBA object creation */
2749       result = CreateObject(shape) ;
2750       const char *entry = InsertInLabel( shape, result->Name(), myCurrentOCAFDoc ) ;
2751       result->ShapeId(entry);
2752       return result ;
2753     }
2754   }
2755   catch(Standard_Failure) {
2756     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ImportIGES()", SALOME::BAD_PARAM);
2757   }
2758
2759   THROW_SALOME_CORBA_EXCEPTION("Import IGES aborted : internal error", SALOME::BAD_PARAM);
2760 }
2761
2762
2763
2764 //==================================================================================
2765 // function : ImportSTEP()
2766 // purpose  : Import shape from an STEP (stp) file
2767 //          : 'result' is a compound of shapes if file contains more entities.
2768 //==================================================================================
2769 GEOM::GEOM_Shape_ptr GEOM_Gen_i::ImportSTEP(const char* filename)
2770   throw (SALOME::SALOME_Exception)
2771 {
2772   Unexpect aCatch(SALOME_SalomeException);
2773   GEOM::GEOM_Shape_var result ;
2774   //VRV: OCC 4.0 migration
2775   STEPControl_Reader aReader;
2776   //VRV: OCC 4.0 migration
2777
2778   TopoDS_Compound compound;
2779   BRep_Builder B;
2780   B.MakeCompound( compound );
2781
2782   try {
2783     IFSelect_ReturnStatus status = aReader.ReadFile((char*)filename);
2784
2785     if (status == IFSelect_RetDone) {
2786       Standard_Boolean failsonly = Standard_False ;
2787       aReader.PrintCheckLoad (failsonly, IFSelect_ItemsByEntity);
2788       /* Root transfers */
2789       Standard_Integer nbr = aReader.NbRootsForTransfer();
2790       aReader.PrintCheckTransfer (failsonly, IFSelect_ItemsByEntity);
2791
2792       for ( Standard_Integer n=1; n <= nbr; n++) {
2793
2794         Standard_Boolean ok = aReader.TransferRoot(n);
2795         /* Collecting resulting entities */
2796         Standard_Integer nbs = aReader.NbShapes();
2797         if (nbs == 0)
2798           THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ImportStep", SALOME::BAD_PARAM) ;
2799
2800         for ( Standard_Integer i=1; i<=nbs; i++ ) {
2801           TopoDS_Shape aShape = aReader.Shape(i);
2802           if ( aShape.IsNull() )
2803             THROW_SALOME_CORBA_EXCEPTION("Null shape in GEOM_Gen_i::ImportStep", SALOME::BAD_PARAM) ;
2804
2805           /* For a single entity */
2806           if(nbr == 1 && nbs == 1) {
2807             result = CreateObject(aShape) ;
2808             const char *entry = InsertInLabel(aShape, result->Name(), myCurrentOCAFDoc) ;
2809             result->ShapeId(entry);
2810             return result ;
2811           }
2812           else {
2813             B.Add( compound, aShape ) ;
2814           }
2815         }
2816       }
2817
2818      TopoDS_Shape tds = compound ;
2819      result = CreateObject(tds) ;
2820      if( CORBA::is_nil(result) )
2821        THROW_SALOME_CORBA_EXCEPTION("Translation aborted : null result", SALOME::BAD_PARAM);
2822      const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
2823      result->ShapeId(entry);
2824      return result ;
2825     }
2826
2827   }
2828   catch(Standard_Failure) {
2829     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ImportStep", SALOME::BAD_PARAM);
2830   }
2831   return result ;
2832 }
2833
2834
2835 //==================================================================================
2836 // function : Partition()
2837 // purpose  :
2838 //==================================================================================
2839
2840 GEOM::GEOM_Shape_ptr
2841   GEOM_Gen_i::Partition(const GEOM::GEOM_Gen::ListOfIOR& ListShapes,
2842                         const GEOM::GEOM_Gen::ListOfIOR& ListTools,
2843                         const GEOM::GEOM_Gen::ListOfIOR& ListKeepInside,
2844                         const GEOM::GEOM_Gen::ListOfIOR& ListRemoveInside,
2845                         const CORBA::Short               Limit)
2846 throw (SALOME::SALOME_Exception)
2847 {
2848   Unexpect aCatch(SALOME_SalomeException);
2849   GEOM::GEOM_Shape_var aResult;
2850   TopoDS_Shape tds ;
2851   //MESSAGE ("In Partition");
2852   try {
2853
2854     unsigned int ind, nbshapes = 0;
2855     nbshapes += ListShapes.length() + ListTools.length();
2856     nbshapes += ListKeepInside.length() + ListRemoveInside.length();
2857
2858     Partition_Spliter PS;
2859     TopTools_MapOfShape ShapesMap(nbshapes), ToolsMap(nbshapes);
2860
2861     // add object shapes that are in ListShapes;
2862     for ( ind = 0; ind < ListShapes.length(); ind++) {
2863
2864       GEOM::GEOM_Shape_var aShape = GetIORFromString( ListShapes[ind] );
2865       TopoDS_Shape Shape = GetTopoShape(aShape);
2866       if(Shape.IsNull() ) {
2867         //MESSAGE ( "In Partition a shape is null" );
2868         THROW_SALOME_CORBA_EXCEPTION("In Partition a shape is null", SALOME::BAD_PARAM);
2869       }
2870       if ( ShapesMap.Add( Shape ))
2871         PS.AddShape(Shape);
2872     }
2873
2874     // add tool shapes that are in ListTools and not in ListShapes;
2875     for (ind = 0; ind < ListTools.length(); ind++) {
2876
2877       GEOM::GEOM_Shape_var aShape = GetIORFromString( ListTools[ind] );
2878       TopoDS_Shape Shape = GetTopoShape(aShape);
2879       if(Shape.IsNull() ) {
2880         THROW_SALOME_CORBA_EXCEPTION("In Partition a shape is null", SALOME::BAD_PARAM);
2881       }
2882       if ( !ShapesMap.Contains( Shape ) && ToolsMap.Add( Shape ))
2883         PS.AddTool(Shape);
2884     }
2885
2886     // add shapes that are in ListKeepInside, as object shapes;
2887     for (ind = 0; ind < ListKeepInside.length(); ind++) {
2888
2889       GEOM::GEOM_Shape_var aShape = GetIORFromString( ListKeepInside[ind] );
2890       TopoDS_Shape Shape = GetTopoShape(aShape);
2891       if(Shape.IsNull() ) {
2892         //MESSAGE ( "In Partition a Keep Inside shape is null" );
2893         THROW_SALOME_CORBA_EXCEPTION("In Partition a shape is null", SALOME::BAD_PARAM);
2894       }
2895       if (!ToolsMap.Contains( Shape ) &&
2896           ShapesMap.Add( Shape ))
2897         PS.AddShape(Shape);
2898     }
2899
2900     // add shapes that are in ListRemoveInside, as object shapes;
2901     for (ind = 0; ind < ListRemoveInside.length(); ind++) {
2902
2903       GEOM::GEOM_Shape_var aShape = GetIORFromString( ListRemoveInside[ind] );
2904       TopoDS_Shape Shape = GetTopoShape(aShape);
2905       if(Shape.IsNull() ) {
2906         //MESSAGE ( "In Partition a Remove Inside shape is null" );
2907         THROW_SALOME_CORBA_EXCEPTION("In Partition a shape is null", SALOME::BAD_PARAM);
2908       }
2909       if (!ToolsMap.Contains( Shape ) &&
2910           ShapesMap.Add( Shape ) )
2911         PS.AddShape(Shape);
2912     }
2913
2914     //MESSAGE ( "Partition::Compute() " );
2915     PS.Compute ((TopAbs_ShapeEnum) Limit);
2916     //MESSAGE ( "Partition::Compute() - END" );
2917
2918     // suppress result outside of shapes in KInsideMap
2919     for (ind = 0; ind < ListKeepInside.length(); ind++) {
2920       GEOM::GEOM_Shape_var aShape = GetIORFromString( ListKeepInside[ind] );
2921       TopoDS_Shape Shape = GetTopoShape(aShape);
2922       PS.KeepShapesInside( Shape );
2923     }
2924
2925     // suppress result inside of shapes in RInsideMap
2926     for (ind = 0; ind < ListRemoveInside.length(); ind++) {
2927
2928       GEOM::GEOM_Shape_var aShape = GetIORFromString( ListRemoveInside[ind] );
2929       TopoDS_Shape Shape = GetTopoShape(aShape);
2930       PS.RemoveShapesInside( Shape );
2931     }
2932
2933     tds = PS.Shape();
2934
2935     if( !IsValid(tds) ) {
2936       //MESSAGE ( "In Partition: non valid shape result" );
2937       THROW_SALOME_CORBA_EXCEPTION("Partition aborted : non valid shape result", SALOME::BAD_PARAM);
2938     }
2939   }
2940   catch (Standard_Failure) {
2941     //MESSAGE ( "In Partition: Exception catched in GEOM_Gen_i::Partition()" );
2942     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::Partition", SALOME::BAD_PARAM);
2943   }
2944
2945   aResult = CreateObject(tds) ;
2946
2947   /* add attributs S and mystr in a new label */
2948   GEOMDS_Commands GC(myCurrentOCAFDoc->Main());
2949   TDF_Label Lab = GC.AddShape(tds, aResult->Name() );
2950   TCollection_AsciiString entry;
2951   TDF_Tool::Entry(Lab,entry);
2952   aResult->ShapeId( entry.ToCString() ) ;
2953
2954   // add arguments
2955
2956   /* Create a new label */
2957   TDF_Label NewLab = Lab.NewChild();
2958   TCollection_ExtendedString Value("Arguments");
2959   TDataStd_Name::Set(NewLab,Value);
2960
2961   // object shapes
2962   for (unsigned int ind = 0; ind < ListShapes.length(); ind++) {
2963     TDF_Label NewLab1 = NewLab.NewChild();
2964
2965     GEOM::GEOM_Shape_var aShape = GetIORFromString( ListShapes[ind] );
2966     Standard_CString anEntry = aShape->ShapeId();
2967     TDF_Label RefLab;
2968     TDF_Tool::Label(myCurrentOCAFDoc->GetData(), anEntry, RefLab);
2969     TDF_Reference::Set(NewLab1,RefLab);
2970   }
2971   // tool shapes
2972   for (unsigned int ind = 0; ind < ListTools.length(); ind++) {
2973     TDF_Label NewLab1 = NewLab.NewChild();
2974
2975     GEOM::GEOM_Shape_var aShape = GetIORFromString( ListTools[ind] );
2976     Standard_CString anEntry = aShape->ShapeId();
2977     TDF_Label RefLab;
2978     TDF_Tool::Label(myCurrentOCAFDoc->GetData(), anEntry, RefLab);
2979     TDF_Reference::Set(NewLab1,RefLab);
2980   }
2981   // limit shapes 1
2982   for (unsigned int ind = 0; ind < ListKeepInside.length(); ind++) {
2983     TDF_Label NewLab1 = NewLab.NewChild();
2984
2985     GEOM::GEOM_Shape_var aShape = GetIORFromString( ListKeepInside[ind] );
2986     Standard_CString anEntry = aShape->ShapeId();
2987     TDF_Label RefLab;
2988     TDF_Tool::Label(myCurrentOCAFDoc->GetData(), anEntry, RefLab);
2989     TDF_Reference::Set(NewLab1,RefLab);
2990   }
2991   // limit shapes 2
2992   for (unsigned int ind = 0; ind < ListRemoveInside.length(); ind++) {
2993     TDF_Label NewLab1 = NewLab.NewChild();
2994
2995     GEOM::GEOM_Shape_var aShape = GetIORFromString( ListRemoveInside[ind] );
2996     Standard_CString anEntry = aShape->ShapeId();
2997     TDF_Label RefLab;
2998     TDF_Tool::Label(myCurrentOCAFDoc->GetData(), anEntry, RefLab);
2999     TDF_Reference::Set(NewLab1,RefLab);
3000   }
3001
3002   return aResult;
3003 }
3004
3005
3006
3007 //==================================================================================
3008 // function : MakeFilling()
3009 // purpose  : Create a surface from section curves filling
3010 //==================================================================================
3011 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeFilling(GEOM::GEOM_Shape_ptr myShape,
3012                                        CORBA::Short mindeg,
3013                                        CORBA::Short maxdeg,
3014                                        CORBA::Double tol3d,
3015                                        CORBA::Double tol2d,
3016                                        CORBA::Short nbiter)
3017   throw (SALOME::SALOME_Exception)
3018 {
3019   Unexpect aCatch(SALOME_SalomeException);
3020   GEOM::GEOM_Shape_var result ;
3021   TopoDS_Face tds ;
3022   TopoDS_Shape aShape = GetTopoShape(myShape) ;
3023   if( aShape.IsNull() || aShape.ShapeType() != TopAbs_COMPOUND ) {
3024     THROW_SALOME_CORBA_EXCEPTION("MakeFilling aborted : null shape or not a compound", SALOME::BAD_PARAM);
3025   }
3026
3027   try {
3028     /* we verify the contents of the shape */
3029     TopExp_Explorer Ex ;
3030     TopoDS_Shape Scurrent ;
3031     Standard_Real First, Last ;
3032     Handle(Geom_Curve) C ;
3033     GeomFill_SectionGenerator Section ;
3034
3035     Standard_Integer i = 0 ;
3036     for(Ex.Init(aShape, TopAbs_EDGE); Ex.More(); Ex.Next()) {
3037       Scurrent = Ex.Current() ;
3038       if( Scurrent.IsNull() || Scurrent.ShapeType() != TopAbs_EDGE)  {
3039         THROW_SALOME_CORBA_EXCEPTION("Initial shape doesn't contain only edges !", SALOME::BAD_PARAM);
3040       }
3041       C = BRep_Tool::Curve(TopoDS::Edge(Scurrent), First, Last);
3042       if (C.IsNull()) continue;
3043       C = new Geom_TrimmedCurve(C, First, Last);
3044       Section.AddCurve(C) ;
3045       i++ ;
3046     }
3047
3048     /* a 'tolerance' is used to compare 2 knots : see GeomFill_Generator.cdl */
3049     /* We set 'tolerance' = tol3d                                            */
3050     // Section.Perform( tol3d ) ; NRI */
3051     Section.Perform( Precision::Confusion() ) ;
3052     Handle(GeomFill_Line) Line = new GeomFill_Line(i) ;
3053
3054     GeomFill_AppSurf App(mindeg, maxdeg, tol3d, tol2d, nbiter) ; /* user parameters */
3055     App.Perform(Line, Section) ;
3056
3057     if (!App.IsDone()) {
3058       THROW_SALOME_CORBA_EXCEPTION("Filling aborted : non valid shape result", SALOME::BAD_PARAM);
3059     }
3060     Standard_Integer UDegree, VDegree, NbUPoles, NbVPoles, NbUKnots, NbVKnots;
3061     App.SurfShape(UDegree, VDegree, NbUPoles, NbVPoles, NbUKnots, NbVKnots);
3062     Handle(Geom_BSplineSurface) GBS = new Geom_BSplineSurface(App.SurfPoles(),
3063                                                               App.SurfWeights(),
3064                                                               App.SurfUKnots(),
3065                                                               App.SurfVKnots(),
3066                                                               App.SurfUMults(),
3067                                                               App.SurfVMults(),
3068                                                               App.UDegree(),
3069                                                               App.VDegree());
3070
3071     if( GBS.IsNull() )  {
3072       THROW_SALOME_CORBA_EXCEPTION("Make Filling aborted", SALOME::BAD_PARAM);
3073     }
3074     tds  = BRepBuilderAPI_MakeFace(GBS) ;
3075   }
3076   catch(Standard_Failure) {
3077     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeFilling", SALOME::BAD_PARAM);
3078   }
3079
3080   /* We test the validity of resulting shape */
3081   if( !IsValid(tds) ) {
3082     THROW_SALOME_CORBA_EXCEPTION("Filling aborted : non valid shape result", SALOME::BAD_PARAM);
3083   }
3084   else {
3085     result = CreateObject(tds) ;
3086     InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
3087   }
3088
3089   return result ;
3090 }
3091
3092
3093 //=================================================================================
3094 // function : MakeGlueFaces()
3095 // purpose  :
3096 //=================================================================================
3097
3098 TopoDS_Face GEOM_Gen_i::FindSameFace(const TopoDS_Shape& aShape,
3099                                      const TopoDS_Face& F,
3100                                      double tol3d)
3101 {
3102   TopoDS_Face aFace;
3103   bool isSame = false;
3104   for (TopExp_Explorer exf(aShape,TopAbs_FACE); exf.More(); exf.Next())
3105     {
3106       //MESSAGE("--- test a face");
3107       int nbFound = 0;
3108       aFace = TopoDS::Face(exf.Current());
3109       TopTools_ListOfShape liste1;
3110       TopTools_ListOfShape liste2;
3111       for (TopExp_Explorer exp(aFace,TopAbs_VERTEX); exp.More(); exp.Next())
3112         {
3113           const TopoDS_Vertex& V = TopoDS::Vertex(exp.Current());
3114           liste1.Append(V);
3115         }
3116       for (TopExp_Explorer exp(F,TopAbs_VERTEX); exp.More(); exp.Next())
3117         {
3118           const TopoDS_Vertex& V = TopoDS::Vertex(exp.Current());
3119           liste2.Append(V);
3120         }
3121       isSame = false;
3122       if (liste1.Extent() == liste2.Extent())
3123         {
3124           TopTools_ListIteratorOfListOfShape it1(liste1);
3125           isSame = true;
3126           for (; it1.More(); it1.Next())
3127             {
3128               bool foundSamePoint = false;
3129               gp_Pnt P1 = BRep_Tool::Pnt(TopoDS::Vertex(it1.Value()));
3130               TopTools_ListIteratorOfListOfShape it2(liste2);
3131               for (it2; it2.More(); it2.Next())
3132                 {
3133                   gp_Pnt P2 = BRep_Tool::Pnt(TopoDS::Vertex(it2.Value()));
3134                   double d = P1.Distance(P2);
3135                   if (d < tol3d)
3136                     {
3137                       nbFound++;
3138                       //MESSAGE("    found Same Point : "<<nbFound<<" - "<<P1.X()<<" "<<P1.Y()<<" "<<P1.Z());
3139                       foundSamePoint = true;
3140                       break;
3141                     }
3142                 }
3143               isSame = isSame && foundSamePoint;
3144               if (! isSame) break; // a vertex does not correspond : not same face
3145             }
3146         }
3147       if (isSame)
3148         {
3149           //MESSAGE("    --- Found Same Face");
3150           break; // a face corresponding to F is found
3151         }
3152     }
3153   if (! isSame) aFace.Nullify(); // return null face
3154   return aFace;
3155 }
3156
3157 TopoDS_Edge GEOM_Gen_i::FindSameEdge(const TopoDS_Face& nf,
3158                                      TopoDS_Edge& Eold,
3159                                      double tol3d)
3160 {
3161   TopoDS_Face newFace = TopoDS::Face(nf.Oriented(TopAbs_REVERSED));
3162   TopoDS_Vertex VFirst, VLast;
3163   TopExp::Vertices(Eold, VFirst, VLast);
3164   gp_Pnt Pf = BRep_Tool::Pnt(VFirst);
3165   gp_Pnt Pl = BRep_Tool::Pnt(VLast);
3166   TopoDS_Edge Enew;
3167   for (TopExp_Explorer ee(newFace,TopAbs_EDGE); ee.More(); ee.Next())
3168     {
3169       const TopoDS_Edge& E = TopoDS::Edge(ee.Current());
3170       TopoDS_Vertex VFn, VLn;
3171       TopExp::Vertices(E, VFn, VLn);
3172       gp_Pnt Pfn = BRep_Tool::Pnt(VFn);
3173       gp_Pnt Pln = BRep_Tool::Pnt(VLn);
3174       double dff = Pf.Distance(Pfn);
3175       double dfl = Pf.Distance(Pln);
3176       double dlf = Pl.Distance(Pfn);
3177       double dll = Pl.Distance(Pln);
3178       if ((dff < tol3d) && (dll <tol3d))
3179         {
3180           //MESSAGE("--- edge forward " <<Pf.X()<<" "<<Pf.Y()<<" "<<Pf.Z()<<" "<<Pl.X()<<" "<<Pl.Y()<<" "<<Pl.Z());
3181           Enew = TopoDS::Edge(E.Oriented(TopAbs_FORWARD));
3182           Eold = TopoDS::Edge(Eold.Oriented(TopAbs_FORWARD));
3183           break;
3184         }
3185       if ((dfl < tol3d) && (dlf <tol3d))
3186         {
3187           //MESSAGE("--- edge reversed " <<Pf.X()<<" "<<Pf.Y()<<" "<<Pf.Z()<<" "<<Pl.X()<<" "<<Pl.Y()<<" "<<Pl.Z());
3188           Enew = TopoDS::Edge(E.Oriented(TopAbs_REVERSED));
3189           Eold = TopoDS::Edge(Eold.Oriented(TopAbs_FORWARD));
3190           break;
3191         }
3192     }
3193   return Enew;
3194 }
3195
3196 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeGlueFaces(GEOM::GEOM_Shape_ptr myShape,
3197                                                double tol3d)
3198   throw (SALOME::SALOME_Exception)
3199 {
3200
3201   // prendre un premier shell dans la liste des shells
3202   // initialiser un compshell avec ce shell
3203   // tant qu'il reste des shells dans la liste
3204   //   chercher un shell qui a des faces en  commun avec le compshell
3205   //   creer un BRepTools_Quilt
3206   //   recenser les faces communes issues du compshell, les ajouter au quilt
3207   //   recenser les faces restantes du shell a inclure, les ajouter au quilt
3208   //   recenser les edges en double, a remplacer
3209   //   pour chaque paire d'edge
3210   //     tester l'orientation relative des aretes
3211   //     bind dans le quilt de Eold.Forward et Enew.Forward (ou reverse)
3212   //   recuperer le nouveau shell
3213   // l'incorporer dans le compshell
3214   // appliquer BRepTools_SameParameter au compshell
3215   // (rendre parametres 2D des edges identiques aux parametres 3D)
3216
3217   Unexpect aCatch(SALOME_SalomeException);
3218   GEOM::GEOM_Shape_var result ;
3219   TopoDS_Shape tds ;
3220   TopoDS_Shape aShape = GetTopoShape(myShape) ;
3221   TopoDS_Compound C;
3222   BRep_Builder bu;
3223   bu.MakeCompound(C); // empty compound;
3224   TopTools_ListOfShape shellList;
3225   for (TopExp_Explorer exp(aShape,TopAbs_SHELL); exp.More(); exp.Next())
3226     {
3227       const TopoDS_Shell& S = TopoDS::Shell(exp.Current());
3228       shellList.Append(S);
3229     }
3230   TopTools_ListIteratorOfListOfShape its(shellList);
3231   if ( ! its.More())
3232     {
3233       THROW_SALOME_CORBA_EXCEPTION("glue aborted : no shell in shape", SALOME::BAD_PARAM);
3234     }
3235   TopoDS_Shell S = TopoDS::Shell(its.Value());
3236   bu.Add(C, S); // add first shell to compound
3237   shellList.Remove(its);
3238   its.Initialize(shellList);
3239   bool shellAdded = true;
3240   while ((shellList.Extent() > 0) && shellAdded)
3241     {
3242       //MESSAGE("more shells : "<< shellList.Extent());
3243       shellAdded = false;
3244       its.Initialize(shellList);
3245       for(; its.More(); its.Next())
3246         {
3247           //MESSAGE("one more shell to try");
3248           TopTools_ListOfShape newFaces; // common faces from new compound
3249           TopTools_ListOfShape oldFaces; // common faces from shell to add
3250           TopTools_ListOfShape addFaces; // not common faces from shell to add
3251           TopTools_ListOfShape newEdges; // common edges from new compound
3252           TopTools_ListOfShape oldEdges; // common edges from face to add
3253           TopoDS_Compound CFN;
3254           TopoDS_Compound CFO;
3255           bu.MakeCompound(CFN);       // empty compound for new faces
3256           bu.MakeCompound(CFO);       // empty compound for old faces
3257           S = TopoDS::Shell(its.Value());
3258           for (TopExp_Explorer exp(S,TopAbs_FACE); exp.More(); exp.Next())
3259             {
3260               //MESSAGE("--- try to find corresponding face in new compound");
3261               TopoDS_Face F = TopoDS::Face(exp.Current());
3262               TopoDS_Face newFace = FindSameFace(C,F,tol3d);
3263               if (! newFace.IsNull())
3264                 {
3265                   //MESSAGE("--- face found");
3266                   newFaces.Append(newFace);
3267                   bu.Add(CFN, newFace); // common faces from new compound
3268                   oldFaces.Append(F);
3269                   for (TopExp_Explorer ee(F,TopAbs_EDGE);ee.More();ee.Next())
3270                     {
3271                       //MESSAGE("--- find edge pair");
3272                       TopoDS_Edge Eold = TopoDS::Edge(ee.Current());
3273                       const TopoDS_Edge& Enew = FindSameEdge(newFace, Eold, tol3d);
3274                       oldEdges.Append(Eold);
3275                       newEdges.Append(Enew);
3276                     }
3277                 }
3278               else
3279                 {
3280                   //MESSAGE("---");
3281                   addFaces.Append(F);
3282                   bu.Add(CFO, F); // not common faces from shell to add
3283                 }
3284             }
3285           if ( !newFaces.IsEmpty())
3286             {
3287               //MESSAGE("--- some faces found ---");
3288               shellAdded = true;
3289               BRepTools_Quilt glue;
3290               glue.Add(CFN);
3291               TopTools_ListIteratorOfListOfShape ito(oldEdges);
3292               TopTools_ListIteratorOfListOfShape itn(newEdges);
3293               for (; ito.More(); ito.Next())
3294                 {
3295                   //MESSAGE("--- bind");
3296                   glue.Bind(TopoDS::Edge(ito.Value()), TopoDS::Edge(itn.Value()));
3297                   itn.Next();
3298                 }
3299               glue.Add(CFO);
3300               TopoDS_Compound newc = TopoDS::Compound(glue.Shells());
3301               for (TopExp_Explorer exs(newc,TopAbs_SHELL); exs.More(); exs.Next())
3302                 {
3303                   TopoDS_Shell NS = TopoDS::Shell(exs.Current());
3304                   bu.Add(C, NS);
3305                 }
3306               shellList.Remove(its);
3307               //MESSAGE("--- remove shell from list");
3308               break;
3309             }
3310         }
3311     }
3312   //MESSAGE("---" << shellList.Extent() << " " << shellAdded);
3313
3314   TopExp_Explorer  exp(C,TopAbs_SHELL);
3315   Standard_Integer ish=0;
3316   TopoDS_Compound  Res;
3317   TopoDS_Solid     Sol;
3318   BRep_Builder     B;
3319   B.MakeCompound(Res);
3320   TopoDS_Shape theShape;
3321
3322   for (; exp.More(); exp.Next())
3323     {
3324       TopoDS_Shape Sh = exp.Current();
3325       B.MakeSolid(Sol);
3326       B.Add(Sol,Sh);
3327       BRepClass3d_SolidClassifier SC(Sol);
3328       SC.PerformInfinitePoint(1.E-6); // cf. BRepFill_Confusion() - BRepFill_Evolved.cxx
3329       if (SC.State() == TopAbs_IN)
3330         {
3331           B.MakeSolid(Sol);
3332           B.Add(Sol,Sh.Reversed());
3333         }
3334       B.Add(Res,Sol);
3335       ish++;
3336     }
3337   if (ish == 1) { theShape = Sol;}
3338   else          { theShape = Res;}
3339
3340   BRepLib::SameParameter(theShape, 1.E-5, Standard_True);
3341   tds = theShape;
3342   result = CreateObject(tds);
3343   InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
3344   //MESSAGE("---");
3345   return result;
3346 }
3347
3348 //=================================================================================
3349 // function : MakeSewing()
3350 // purpose  :
3351 //=================================================================================
3352 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeSewing( const GEOM::GEOM_Gen::ListOfIOR& ListShapes,
3353                                        CORBA::Double precision )
3354   throw (SALOME::SALOME_Exception)
3355 {
3356   Unexpect aCatch(SALOME_SalomeException);
3357   GEOM::GEOM_Shape_var result ;
3358   TopoDS_Shape tds ;
3359   BRepOffsetAPI_Sewing aMethod ;
3360
3361   try {
3362     /* default OCC is 1.0e-06 */
3363     aMethod.Init(precision, Standard_False);
3364     for ( unsigned int i = 0; i < ListShapes.length(); i++) {
3365       GEOM::GEOM_Shape_var aShape = GetIORFromString( ListShapes[i] );
3366       TopoDS_Shape Shape = GetTopoShape(aShape) ;
3367       if( Shape.IsNull() ) {
3368         THROW_SALOME_CORBA_EXCEPTION("MakeSewing aborted : null shape during operation", SALOME::BAD_PARAM);
3369       }
3370       aMethod.Add(Shape) ;
3371     }
3372
3373     aMethod.Perform() ;
3374     tds = aMethod.SewedShape() ;
3375     if( !IsValid(tds) ) {
3376       THROW_SALOME_CORBA_EXCEPTION("Make Sewing aborted : non valid shape", SALOME::BAD_PARAM);
3377     }
3378     if( tds.IsNull() ) {
3379       THROW_SALOME_CORBA_EXCEPTION("Make Sewing aborted : null shape", SALOME::BAD_PARAM);
3380     }
3381   }
3382   catch (Standard_Failure) {
3383     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeSewing", SALOME::BAD_PARAM);
3384   }
3385
3386   result = CreateObject(tds);
3387   InsertInLabelMoreArguments(tds, result, ListShapes, myCurrentOCAFDoc) ;
3388   return result;
3389 }
3390
3391 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeSewingShape( GEOM::GEOM_Shape_ptr aShape,
3392                                             CORBA::Double precision )
3393   throw (SALOME::SALOME_Exception)
3394 {
3395   Unexpect aCatch(SALOME_SalomeException);
3396   GEOM::GEOM_Shape_var result ;
3397   TopoDS_Shape tds, S ;
3398   BRepOffsetAPI_Sewing aMethod ;
3399
3400   try {
3401     S = GetTopoShape(aShape) ;
3402     if(S.IsNull() ) {
3403       THROW_SALOME_CORBA_EXCEPTION("In Sewing a Shape is null", SALOME::BAD_PARAM);
3404     }
3405
3406     /* default OCC is 1.0e-06 */
3407     aMethod.Init(precision, Standard_False);
3408     for ( TopExp_Explorer exp( S, TopAbs_FACE); exp.More(); exp.Next() ) {
3409       const TopoDS_Face& F = TopoDS::Face(exp.Current());
3410       aMethod.Add(F) ;
3411     }
3412
3413     aMethod.Perform() ;
3414     tds = aMethod.SewedShape() ;
3415     if( !IsValid(tds) ) {
3416       THROW_SALOME_CORBA_EXCEPTION("Make Sewing aborted : non valid shape", SALOME::BAD_PARAM);
3417     }
3418   }
3419   catch (Standard_Failure) {
3420     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeSewing", SALOME::BAD_PARAM);
3421   }
3422
3423   result = CreateObject(tds);
3424   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3425   result->ShapeId(entry) ;
3426   return result;
3427 }
3428
3429 //==================================================================================
3430 // function : OrientationChange()
3431 // purpose  : Change the orientation of a new shape
3432 //          : TopAbs_FORWARD < -- > TopAbs_REVERSED
3433 //
3434 //          : WARNING : for the moment we make a new shape !
3435 //==================================================================================
3436 GEOM::GEOM_Shape_ptr GEOM_Gen_i::OrientationChange(GEOM::GEOM_Shape_ptr aShape)
3437   throw (SALOME::SALOME_Exception)
3438 {
3439   Unexpect aCatch(SALOME_SalomeException);
3440   GEOM::GEOM_Shape_var result ;
3441   BRep_Builder aBuilder;
3442
3443   TopoDS_Shape shape = GetTopoShape(aShape) ;
3444   if( shape.IsNull() ) {
3445     THROW_SALOME_CORBA_EXCEPTION("Shape is null", SALOME::BAD_PARAM);
3446   }
3447
3448   BRepBuilderAPI_Copy Copy(shape);
3449   if( Copy.IsDone() ) {
3450     TopoDS_Shape tds = Copy.Shape();
3451     if( tds.IsNull() ) {
3452       THROW_SALOME_CORBA_EXCEPTION("Orientation aborted : null shape", SALOME::BAD_PARAM);
3453     }
3454
3455     if( tds.Orientation() == TopAbs_FORWARD)
3456       tds.Orientation(TopAbs_REVERSED) ;
3457     else
3458       tds.Orientation(TopAbs_FORWARD) ;
3459
3460     result = CreateObject(tds);
3461     InsertInLabelOneArgument(shape, aShape, tds, result, myCurrentOCAFDoc) ;
3462   }
3463   return result ;
3464 }
3465
3466
3467 //==================================================================================
3468 // function : GetReferencedObjects()
3469 // purpose  :
3470 //==================================================================================
3471 GEOM::GEOM_Gen::ListOfIOR* GEOM_Gen_i::GetReferencedObjects(GEOM::GEOM_Shape_ptr shape)
3472 {
3473   GEOM::GEOM_Gen::ListOfIOR_var aList = new GEOM::GEOM_Gen::ListOfIOR;
3474   aList->length(0);
3475
3476   if (shape->_is_nil()) return aList._retn();
3477
3478   Standard_CString entry = shape->ShapeId();
3479   TDF_Label Lab;
3480   TDF_Tool::Label(myCurrentOCAFDoc->GetData(), entry, Lab);
3481
3482   Handle(TDataStd_Name) Att;
3483   Lab.FindAttribute(TDataStd_Name::GetID(),Att);
3484
3485   TDF_ChildIterator ChildIterator(Lab);
3486   if (ChildIterator.More()) {
3487     TDF_Label L = ChildIterator.Value();
3488     Handle(TDataStd_Name) Att;
3489     L.FindAttribute(TDataStd_Name::GetID(),Att);
3490     if (Att->Get().IsEqual(TCollection_ExtendedString("Arguments")) ) {
3491
3492       TDF_ChildIterator ChildIterator1(L);
3493       unsigned int i = 0;
3494
3495       while (ChildIterator1.More()) {
3496         TDF_Label L = ChildIterator1.Value();
3497
3498         Handle(TDF_Reference) Ref;
3499         if (L.FindAttribute(TDF_Reference::GetID(),Ref)) {
3500           i++;
3501         }
3502         ChildIterator1.Next();
3503       }
3504       aList->length(i);
3505       i = 0;
3506       TDF_ChildIterator ChildIterator2(L);
3507       while (ChildIterator2.More()) {
3508         TDF_Label L = ChildIterator2.Value();
3509         Handle(TDF_Reference) Ref;
3510         if (L.FindAttribute(TDF_Reference::GetID(),Ref)) {
3511           TDF_Label L = Ref->Get();
3512
3513           Handle(TDataStd_Name) Att;
3514           L.FindAttribute(TDataStd_Name::GetID(),Att);
3515           TCollection_AsciiString nameIOR (Att->Get()) ;
3516           aList[i] = CORBA::string_dup( nameIOR.ToCString() );
3517           i++;
3518         }
3519
3520         ChildIterator2.Next();
3521       }
3522     }
3523   }
3524   return aList._retn();
3525 }
3526
3527 //==================================================================================
3528 // function : GetObjects()
3529 // purpose  :
3530 //==================================================================================
3531 GEOM::GEOM_Gen::ListOfIOR* GEOM_Gen_i::GetObjects(GEOM::GEOM_Shape_ptr shape)
3532 {
3533   GEOM::GEOM_Gen::ListOfIOR_var aList = new GEOM::GEOM_Gen::ListOfIOR;
3534   aList->length(0);
3535
3536   Standard_CString entry = shape->ShapeId();
3537   TDF_Label Lab;
3538   TDF_Tool::Label(myCurrentOCAFDoc->GetData(), entry, Lab);
3539
3540   Handle(TDataStd_Name) Att;
3541   Lab.FindAttribute(TDataStd_Name::GetID(),Att);
3542
3543   TDF_ChildIterator ChildIterator(Lab);
3544   unsigned int i = 0;
3545   while (ChildIterator.More()) {
3546     TDF_Label L = ChildIterator.Value();
3547     Handle(TDataStd_Name) Att;
3548     L.FindAttribute(TDataStd_Name::GetID(),Att);
3549
3550     if (!Att->Get().IsEqual(TCollection_ExtendedString("Arguments")) ) {
3551       i++;
3552     }
3553     ChildIterator.Next();
3554   }
3555
3556   aList->length(i);
3557   i = 0;
3558   TDF_ChildIterator ChildIterator1(Lab);
3559   while (ChildIterator1.More()) {
3560     TDF_Label L = ChildIterator1.Value();
3561     Handle(TDataStd_Name) Att;
3562     L.FindAttribute(TDataStd_Name::GetID(),Att);
3563
3564     if (!Att->Get().IsEqual(TCollection_ExtendedString("Arguments")) ) {
3565       TCollection_AsciiString nameIOR (Att->Get());
3566       aList[i] = CORBA::string_dup( nameIOR.ToCString() );
3567       i++;
3568     }
3569     ChildIterator1.Next();
3570   }
3571   return aList._retn();
3572 }
3573
3574
3575 //==================================================================================
3576 // function : Import
3577 // purpose  : Import shape from a BREP file
3578 //==================================================================================
3579 GEOM::GEOM_Shape_ptr GEOM_Gen_i::ImportBREP(const char* filename)
3580   throw (SALOME::SALOME_Exception)
3581 {
3582   Unexpect aCatch(SALOME_SalomeException);
3583   TopoDS_Shape tds ;
3584   GEOM::GEOM_Shape_var result ;
3585
3586   try {
3587     BRep_Builder aBuilder;
3588     char* aCopyfilename = strdup(filename);
3589     BRepTools::Read(tds, aCopyfilename, aBuilder) ;
3590     free(aCopyfilename);
3591     if (tds.IsNull()) {
3592       THROW_SALOME_CORBA_EXCEPTION("Import BRep aborted", SALOME::BAD_PARAM);
3593     }
3594   }
3595   catch(Standard_Failure) {
3596     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ImportBREP", SALOME::BAD_PARAM);
3597   }
3598
3599   result = CreateObject(tds) ;
3600   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3601   result->ShapeId(entry);
3602   return result;
3603 }
3604
3605
3606 //================================================================================
3607 // function : MakePlane()
3608 // purpose  : Make a plane topology (non infinite)
3609 //================================================================================
3610 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakePlane(const GEOM::PointStruct& pstruct,
3611                                      const GEOM::DirStruct& dstruct,
3612                                      CORBA::Double trimsize)
3613   throw (SALOME::SALOME_Exception)
3614 {
3615   Unexpect aCatch(SALOME_SalomeException);
3616   GEOM::GEOM_Shape_var result ;
3617   TopoDS_Shape tds ;
3618
3619   try {
3620     gp_Pnt aPoint(pstruct.x, pstruct.y, pstruct.z) ;
3621     gp_Dir aDirection(dstruct.PS.x, dstruct.PS.y, dstruct.PS.z) ;
3622     /*  we make a trimmed plane */
3623     gp_Pln gplane(aPoint, aDirection) ;
3624     tds = BRepBuilderAPI_MakeFace(gplane, -trimsize, +trimsize, -trimsize, +trimsize) ;
3625   }
3626   catch(Standard_Failure) {
3627     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakePlane", SALOME::BAD_PARAM);
3628   }
3629
3630   if (tds.IsNull()) {
3631     THROW_SALOME_CORBA_EXCEPTION("Make Plane aborted : null shape", SALOME::BAD_PARAM);
3632    }
3633
3634   result = CreateObject(tds) ;
3635   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3636   result->ShapeId(entry);
3637   return result ;
3638 }
3639
3640 //=================================================================================
3641 // function : MakeVertex()
3642 // purpose  : Create a Vertex topology.
3643 //=================================================================================
3644 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeVertex(CORBA::Double x,
3645                                       CORBA::Double y,
3646                                       CORBA::Double z)
3647   throw (SALOME::SALOME_Exception)
3648 {
3649   Unexpect aCatch(SALOME_SalomeException);
3650   GEOM::GEOM_Shape_var result ;
3651   gp_Pnt P(x,y,z);
3652   TopoDS_Shape tds = BRepBuilderAPI_MakeVertex(P).Shape();
3653   if (tds.IsNull()) {
3654     THROW_SALOME_CORBA_EXCEPTION("Make Vertex/Point aborted", SALOME::BAD_PARAM);
3655   }
3656   tds.Infinite(true);
3657   result = CreateObject(tds) ;
3658   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3659   result->ShapeId(entry);
3660   return result ;
3661 }
3662
3663
3664 //=================================================================================
3665 // function : MakeFace()
3666 // purpose  :
3667 //=================================================================================
3668 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeFace( GEOM::GEOM_Shape_ptr wire,
3669                                      CORBA::Boolean wantplanarface )
3670   throw (SALOME::SALOME_Exception)
3671 {
3672   Unexpect aCatch(SALOME_SalomeException);
3673   GEOM::GEOM_Shape_var result ;
3674   TopoDS_Shape aShape;
3675   TopoDS_Shape tds;
3676
3677   try {
3678     aShape = GetTopoShape(wire) ;
3679     if( aShape.IsNull() || aShape.ShapeType() != TopAbs_WIRE ) {
3680       THROW_SALOME_CORBA_EXCEPTION("MakeFace aborted : null or inappropriate shape", SALOME::BAD_PARAM);
3681     }
3682     TopoDS_Wire W = TopoDS::Wire(aShape) ;
3683     tds = BRepBuilderAPI_MakeFace(W, wantplanarface).Shape() ;
3684     if( !tds.IsNull() ) {
3685       result = CreateObject(tds) ;
3686       InsertInLabelOneArgument(aShape, wire, tds, result, myCurrentOCAFDoc) ;
3687     }
3688     else {
3689       THROW_SALOME_CORBA_EXCEPTION("Null result in GEOM_Gen_i::MakeFace", SALOME::BAD_PARAM);
3690     }
3691   }
3692   catch (Standard_Failure) {
3693     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeFace", SALOME::BAD_PARAM);
3694   }
3695   return result ;
3696 }
3697
3698
3699 //=================================================================================
3700 // function : MakeFaces()
3701 // purpose  :
3702 //=================================================================================
3703 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeFaces(const GEOM::GEOM_Gen::ListOfIOR& ListShapes,
3704                                           CORBA::Boolean wantplanarface)
3705   throw (SALOME::SALOME_Exception)
3706 {
3707   GEOM::GEOM_Shape_var result;
3708
3709   try {
3710     GEOM::GEOM_Shape_var aShape = GetIORFromString(ListShapes[0]);
3711     TopoDS_Shape Shape = GetTopoShape(aShape);
3712     if(Shape.IsNull() || Shape.ShapeType() != TopAbs_WIRE) {
3713       THROW_SALOME_CORBA_EXCEPTION("Shell aborted : null shape during operation", SALOME::BAD_PARAM);
3714     }
3715     TopoDS_Wire W = TopoDS::Wire(Shape);
3716     TopoDS_Shape FFace = BRepBuilderAPI_MakeFace(W, wantplanarface).Shape();
3717     if(!FFace.IsNull()) {
3718       if(ListShapes.length() == 1) {
3719         result = CreateObject(FFace);
3720         InsertInLabelMoreArguments(FFace, result, ListShapes, myCurrentOCAFDoc);
3721       }
3722       else if(ListShapes.length() >= 2) {
3723         TopoDS_Compound C;
3724         BRep_Builder aBuilder;
3725         aBuilder.MakeCompound(C);
3726         BRepAlgo_FaceRestrictor FR;
3727
3728         TopAbs_Orientation OriF = FFace.Orientation();
3729         TopoDS_Shape aLocalS = FFace.Oriented(TopAbs_FORWARD);
3730         FR.Init(TopoDS::Face(aLocalS), Standard_False, Standard_True);
3731
3732         for(unsigned int i = 0; i < ListShapes.length(); i++) {
3733           GEOM::GEOM_Shape_var aShape = GetIORFromString(ListShapes[i]);
3734           TopoDS_Shape Shape = GetTopoShape(aShape);
3735           if(Shape.IsNull()) {
3736             THROW_SALOME_CORBA_EXCEPTION("Shell aborted : null shape during operation", SALOME::BAD_PARAM);
3737           }
3738           FR.Add(TopoDS::Wire(Shape));
3739         }
3740
3741         FR.Perform();
3742
3743         if(FR.IsDone()) {
3744           int k = 0;
3745           TopoDS_Shape aFace;
3746           for(; FR.More(); FR.Next()) {
3747             aFace = FR.Current().Oriented(OriF);
3748             aBuilder.Add(C, aFace);
3749             k++;
3750           }
3751           if(k == 1) {
3752             result = CreateObject(aFace);
3753             InsertInLabelMoreArguments(aFace, result, ListShapes, myCurrentOCAFDoc);
3754           }
3755           else {
3756             result = CreateObject(C);
3757             InsertInLabelMoreArguments(C, result, ListShapes, myCurrentOCAFDoc);
3758           }
3759         }
3760       }
3761     }
3762     else {
3763       THROW_SALOME_CORBA_EXCEPTION("Null result in GEOM_Gen_i::MakeFace", SALOME::BAD_PARAM);
3764     }
3765   }
3766   catch (Standard_Failure) {
3767     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeFace", SALOME::BAD_PARAM);
3768   }
3769   return result;
3770 }
3771
3772
3773 //=================================================================================
3774 // function : MakeShell()
3775 // purpose  : Make a compound from a list containing one or more shapes
3776 //=================================================================================
3777 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeShell( const GEOM::GEOM_Gen::ListOfIOR& ListShapes )
3778   throw (SALOME::SALOME_Exception)
3779 {
3780   Unexpect aCatch(SALOME_SalomeException);
3781   GEOM::GEOM_Shape_var result ;
3782   BRepTools_Quilt Glue;
3783   TopoDS_Shape C;
3784
3785   for ( unsigned int i = 0; i < ListShapes.length(); i++) {
3786     GEOM::GEOM_Shape_var aShape = GetIORFromString( ListShapes[i] );
3787     TopoDS_Shape Shape = GetTopoShape(aShape) ;
3788     if( Shape.IsNull() ) {
3789        THROW_SALOME_CORBA_EXCEPTION("Shell aborted : null shape during operation", SALOME::BAD_PARAM);
3790     }
3791     Glue.Add(Shape) ;
3792   }
3793
3794   TopExp_Explorer exp(Glue.Shells(), TopAbs_SHELL);
3795   Standard_Integer ish = 0;
3796   for (; exp.More(); exp.Next()) {
3797     C = exp.Current();
3798     ish++;
3799   }
3800
3801   if (ish != 1)
3802     C = Glue.Shells();
3803
3804   if ( C.IsNull() ) {
3805     THROW_SALOME_CORBA_EXCEPTION("Null result : Shell operation aborted", SALOME::BAD_PARAM);
3806   }
3807   else {
3808     result = CreateObject(C) ;
3809     InsertInLabelMoreArguments(C, result, ListShapes, myCurrentOCAFDoc) ;
3810   }
3811   return result;
3812 }
3813
3814
3815 //=================================================================================
3816 // function : MakeSolid()
3817 // purpose  : Make a compound from a list containing one or more shapes
3818 //=================================================================================
3819 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeSolid( const GEOM::GEOM_Gen::ListOfIOR& ListShapes )
3820   throw (SALOME::SALOME_Exception)
3821 {
3822   Unexpect aCatch(SALOME_SalomeException);
3823   GEOM::GEOM_Shape_var result ;
3824   Standard_Integer ish = 0;
3825   TopoDS_Compound  Res;
3826   TopoDS_Solid     Sol;
3827   BRep_Builder     B;
3828   TopoDS_Shape     Shape;
3829
3830   B.MakeCompound(Res);
3831
3832   for ( unsigned int i = 0; i < ListShapes.length(); i++) {
3833     GEOM::GEOM_Shape_var aShape = GetIORFromString( ListShapes[i] );
3834     TopoDS_Shape Sh = GetTopoShape(aShape) ;
3835     if( Sh.IsNull() ) {
3836        THROW_SALOME_CORBA_EXCEPTION("Solid aborted : null shape during operation", SALOME::BAD_PARAM);
3837     }
3838     B.MakeSolid(Sol);
3839     B.Add(Sol,Sh);
3840     BRepClass3d_SolidClassifier SC(Sol);
3841     SC.PerformInfinitePoint(Precision::Confusion());
3842     if (SC.State() == TopAbs_IN) {
3843       B.MakeSolid(Sol);
3844       B.Add(Sol,Sh.Reversed());
3845     }
3846     B.Add(Res,Sol);
3847     ish++;
3848   }
3849   if (ish == 1) { Shape = Sol;}
3850   else          { Shape = Res;}
3851
3852   if ( Shape.IsNull() ) {
3853     THROW_SALOME_CORBA_EXCEPTION("Null result : Solid operation aborted", SALOME::BAD_PARAM);
3854   }
3855   else {
3856     result = CreateObject(Shape) ;
3857     InsertInLabelMoreArguments(Shape, result, ListShapes, myCurrentOCAFDoc) ;
3858   }
3859   return result;
3860 }
3861
3862
3863 //================================================================================
3864 // function : MakeLine
3865 // purpose  : Make a Line topology
3866 //================================================================================
3867 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeLine(const GEOM::PointStruct& pstruct,
3868                                     const GEOM::DirStruct& dstruct)
3869   throw (SALOME::SALOME_Exception)
3870 {
3871   Unexpect aCatch(SALOME_SalomeException);
3872   GEOM::GEOM_Shape_var result  ;
3873   gp_Pnt P1(pstruct.x, pstruct.y, pstruct.z);
3874   gp_Pnt P2(dstruct.PS.x, dstruct.PS.y, dstruct.PS.z) ;
3875   TopoDS_Shape tds ;
3876
3877   try {
3878     tds = BRepBuilderAPI_MakeEdge(P1, P2).Shape();
3879   }
3880   catch(Standard_Failure) {
3881     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeLine", SALOME::BAD_PARAM);
3882   }
3883
3884   if ( tds.IsNull() ) {
3885     THROW_SALOME_CORBA_EXCEPTION("Make Line aborted : null shape", SALOME::BAD_PARAM);
3886   }
3887   else {
3888     tds.Infinite(true);
3889     result = CreateObject(tds) ;
3890     const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3891     result->ShapeId(entry);
3892   }
3893   return result ;
3894 }
3895
3896
3897 //================================================================================
3898 // function : MakeVector()
3899 // purpose  : Make a vector
3900 //================================================================================
3901 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeVector(const GEOM::PointStruct& pstruct1,
3902                                       const GEOM::PointStruct& pstruct2)
3903   throw (SALOME::SALOME_Exception)
3904 {
3905   Unexpect aCatch(SALOME_SalomeException);
3906   GEOM::GEOM_Shape_var result  ;
3907   TopoDS_Shape tds ;
3908
3909   try {
3910     gp_Pnt P1(pstruct1.x, pstruct1.y, pstruct1.z);
3911     gp_Pnt P2(pstruct2.x, pstruct2.y, pstruct2.z) ;
3912     tds = BRepBuilderAPI_MakeEdge(P1, P2).Shape();
3913   }
3914   catch(Standard_Failure) {
3915     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeVector", SALOME::BAD_PARAM);
3916   }
3917
3918   if ( tds.IsNull() ) {
3919     THROW_SALOME_CORBA_EXCEPTION("Make Vector aborted : null shape", SALOME::BAD_PARAM);
3920   }
3921   else {
3922     result = CreateObject(tds) ;
3923     const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3924     result->ShapeId(entry);
3925   }
3926   return result ;
3927 }
3928
3929
3930 //================================================================================
3931 // function : MakeCircle()
3932 // purpose  :
3933 //================================================================================
3934 GEOM::GEOM_Shape_ptr  GEOM_Gen_i::MakeCircle(const GEOM::PointStruct& pstruct,
3935                                        const GEOM::DirStruct& dstruct,
3936                                        CORBA::Double radius)
3937   throw (SALOME::SALOME_Exception)
3938 {
3939   Unexpect aCatch(SALOME_SalomeException);
3940   GEOM::GEOM_Shape_var result;
3941   TopoDS_Shape tds ;
3942
3943   try {
3944     gp_Pnt p(pstruct.x, pstruct.y, pstruct.z) ;
3945     gp_Dir d(dstruct.PS.x, dstruct.PS.y, dstruct.PS.z) ;
3946     gp_Ax2 axis(p, d) ;
3947     gp_Circ circ( axis, radius);
3948     BRepBuilderAPI_MakeEdge MakeEdge( circ );
3949     tds = MakeEdge.Edge();
3950   }
3951   catch(Standard_Failure) {
3952     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeCircle", SALOME::BAD_PARAM);
3953   }
3954   if (tds.IsNull()) {
3955     THROW_SALOME_CORBA_EXCEPTION("Make Circle aborted", SALOME::BAD_PARAM);
3956   }
3957   result = CreateObject(tds);
3958   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3959   result->ShapeId(entry);
3960   return result ;
3961 }
3962
3963 //================================================================================
3964 // function : MakeEllipse()
3965 // purpose  :
3966 //================================================================================
3967 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeEllipse( const GEOM::PointStruct& pstruct,
3968                                               const GEOM::DirStruct& dstruct,
3969                                               CORBA::Double radius_major,
3970                                               CORBA::Double radius_minor )
3971   throw (SALOME::SALOME_Exception)
3972 {
3973   Unexpect aCatch(SALOME_SalomeException);
3974   GEOM::GEOM_Shape_var result;
3975   TopoDS_Shape tds ;
3976
3977   try {
3978     gp_Pnt p(pstruct.x, pstruct.y, pstruct.z) ;
3979     gp_Dir d(dstruct.PS.x, dstruct.PS.y, dstruct.PS.z) ;
3980
3981     const gp_Ax2 axis(p, d) ;
3982     gp_Elips anEllipse( axis, radius_major, radius_minor ) ;
3983     BRepBuilderAPI_MakeEdge MakeEdge( anEllipse );
3984     tds = MakeEdge.Edge();
3985
3986   }
3987   catch(Standard_Failure) {
3988     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeEllipse", SALOME::BAD_PARAM);
3989   }
3990   if (tds.IsNull()) {
3991     THROW_SALOME_CORBA_EXCEPTION("Make Ellipse aborted", SALOME::BAD_PARAM);
3992   }
3993   result = CreateObject(tds);
3994   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
3995   result->ShapeId(entry);
3996   return result ;
3997 }
3998
3999 //================================================================================
4000 // function : MakeArc()
4001 // purpose  : make an arc of circle from pInit to pEnd and passing on pCircle
4002 //================================================================================
4003 GEOM::GEOM_Shape_ptr  GEOM_Gen_i::MakeArc(const GEOM::PointStruct& pInit,
4004                                     const GEOM::PointStruct& pCircle,
4005                                     const GEOM::PointStruct& pEnd)
4006   throw (SALOME::SALOME_Exception)
4007 {
4008   Unexpect aCatch(SALOME_SalomeException);
4009   GEOM::GEOM_Shape_var result;
4010   try {
4011     gp_Pnt pI(pInit.x, pInit.y, pInit.z) ;
4012     gp_Pnt pC(pCircle.x, pCircle.y, pCircle.z) ;
4013     gp_Pnt pE(pEnd.x, pEnd.y, pEnd.z) ;
4014
4015     GC_MakeArcOfCircle arc( pI, pC, pE ) ;
4016     if( !arc.IsDone() ) {
4017       THROW_SALOME_CORBA_EXCEPTION("Arc not done", SALOME::BAD_PARAM);
4018     }
4019     BRepBuilderAPI_MakeEdge MakeEdge( arc );
4020     TopoDS_Shape tds = MakeEdge.Edge();
4021     if (tds.IsNull()) {
4022       THROW_SALOME_CORBA_EXCEPTION("Null result : arc not done", SALOME::BAD_PARAM);
4023     }
4024     else {
4025       result = CreateObject(tds);
4026       const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
4027       result->ShapeId(entry);
4028     }
4029   }
4030   catch(Standard_Failure) {
4031     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeArc", SALOME::BAD_PARAM);
4032   }
4033   return result ;
4034 }
4035
4036 //=================================================================================
4037  // function : MakeSketcher()
4038  // purpose  : Make a wire from a list containing many points
4039  //=================================================================================
4040  GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeSketcher( const char *Cmd )
4041    throw (SALOME::SALOME_Exception)
4042  {
4043    GEOM::GEOM_Shape_var result ;
4044    TopoDS_Shape tds ;
4045    try {
4046      Sketcher_Profile aProfile (Cmd);
4047      if(aProfile.IsDone())
4048        tds = aProfile.GetShape();
4049    }
4050    catch(Standard_Failure) {
4051      THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeSketcher", SALOME::BAD_PARAM);
4052    }
4053
4054    if (tds.IsNull()) {
4055      THROW_SALOME_CORBA_EXCEPTION("MakeSketcher aborted : null shape", SALOME::BAD_PARAM);
4056    }
4057    else {
4058      result = CreateObject(tds);
4059      const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
4060      result->ShapeId(entry) ;
4061    }
4062    return result;
4063  }
4064
4065 //=================================================================================
4066 // function : MakeTranslation()
4067 // purpose  : Translate a 3D shape
4068 //=================================================================================
4069 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeTranslation( GEOM::GEOM_Shape_ptr myShape,
4070                                             CORBA::Double x,
4071                                             CORBA::Double y,
4072                                             CORBA::Double z)
4073   throw (SALOME::SALOME_Exception)
4074 {
4075   Unexpect aCatch(SALOME_SalomeException);
4076   GEOM::GEOM_Shape_var result ;
4077   TopoDS_Shape aShape = GetTopoShape(myShape) ;
4078   if( aShape.IsNull() ) {
4079     THROW_SALOME_CORBA_EXCEPTION("Translation aborted : null shape", SALOME::BAD_PARAM);
4080   }
4081   gp_Vec theVector(x,y,z) ;
4082   gp_Trsf theTransformation ;
4083   theTransformation.SetTranslation(theVector) ;
4084   BRepBuilderAPI_Transform myBRepTransformation(aShape, theTransformation, Standard_False) ;
4085   TopoDS_Shape tds = myBRepTransformation.Shape() ;
4086
4087   result = CreateObject(tds) ;
4088   if( CORBA::is_nil(result) ) {
4089     THROW_SALOME_CORBA_EXCEPTION("Translation aborted : null result", SALOME::BAD_PARAM);
4090   }
4091   InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4092   return result;
4093 }
4094
4095
4096 //=================================================================================
4097 // function : MakeMultiTranslation1D()
4098 // purpose  : Multi-Translate a 3D shape
4099 //=================================================================================
4100 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeMultiTranslation1D( GEOM::GEOM_Shape_ptr myShape,
4101                                                    const GEOM::DirStruct& dir,
4102                                                    CORBA::Double step,
4103                                                    CORBA::Short nbtimes )
4104   throw (SALOME::SALOME_Exception)
4105 {
4106   Unexpect aCatch(SALOME_SalomeException);
4107   GEOM::GEOM_Shape_var result ;
4108   TopoDS_Shape tds ;
4109
4110   TopoDS_Shape aShape = GetTopoShape(myShape) ;
4111   if( aShape.IsNull() )
4112     THROW_SALOME_CORBA_EXCEPTION("MakeMultiTranslation1D aborted : null shape", SALOME::BAD_PARAM);
4113
4114   try {
4115     int i ;
4116     double DX, DY, DZ ;
4117     gp_Trsf theTransformation ;
4118     gp_Vec myVec ;
4119     gp_Vec Vec( dir.PS.x, dir.PS.y, dir.PS.z ) ;
4120     Vec.Normalize();
4121     TopoDS_Compound compound;
4122     BRep_Builder B;
4123     B.MakeCompound( compound );
4124
4125     for ( i = 0; i < nbtimes; i++ ) {
4126       DX = i * step * Vec.X() ;
4127       DY = i * step * Vec.Y() ;
4128       DZ = i * step * Vec.Z() ;
4129       myVec.SetCoord( DX, DY, DZ ) ;
4130       theTransformation.SetTranslation(myVec) ;
4131       BRepBuilderAPI_Transform myBRepTransformation(aShape, theTransformation, Standard_False) ;
4132       B.Add( compound, myBRepTransformation.Shape() );
4133     }
4134     tds = compound ;
4135     result = CreateObject(tds) ;
4136   }
4137   catch (Standard_Failure) {
4138     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeMultiTranslation1D", SALOME::BAD_PARAM);
4139   }
4140
4141   if( CORBA::is_nil(result) ) {
4142     THROW_SALOME_CORBA_EXCEPTION("MakeMultiTranslation1D aborted : null result", SALOME::BAD_PARAM);
4143   }
4144   InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4145   return result;
4146 }
4147
4148
4149 //=================================================================================
4150 // function : MakeMultiTranslation2D()
4151 // purpose  : Multi-Translate a 3D shape
4152 //=================================================================================
4153 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeMultiTranslation2D( GEOM::GEOM_Shape_ptr myShape,
4154                                                    const GEOM::DirStruct& dir1,
4155                                                    CORBA::Double step1,
4156                                                    CORBA::Short nbtimes1,
4157                                                    const GEOM::DirStruct& dir2,
4158                                                    CORBA::Double step2,
4159                                                    CORBA::Short nbtimes2 )
4160   throw (SALOME::SALOME_Exception)
4161 {
4162   Unexpect aCatch(SALOME_SalomeException);
4163   GEOM::GEOM_Shape_var result ;
4164   TopoDS_Shape tds ;
4165
4166   TopoDS_Shape aShape = GetTopoShape(myShape) ;
4167   if( aShape.IsNull() ) {
4168     THROW_SALOME_CORBA_EXCEPTION("MakeMultiTranslation2D aborted : null shape", SALOME::BAD_PARAM);
4169   }
4170
4171   try {
4172     int i, j ;
4173     double DX, DY, DZ ;
4174     gp_Trsf theTransformation ;
4175     gp_Vec myVec ;
4176     gp_Vec Vec1( dir1.PS.x, dir1.PS.y, dir1.PS.z ) ;
4177     Vec1.Normalize();
4178     gp_Vec Vec2( dir2.PS.x, dir2.PS.y, dir2.PS.z ) ;
4179     Vec2.Normalize();
4180     TopoDS_Compound compound;
4181     BRep_Builder B;
4182     B.MakeCompound( compound );
4183
4184     for ( i = 0; i < nbtimes1; i++ ) {
4185       for ( j = 0; j < nbtimes2; j++ ) {
4186         DX = i * step1 * Vec1.X() + j * step2 * Vec2.X() ;
4187         DY = i * step1 * Vec1.Y() + j * step2 * Vec2.Y() ;
4188         DZ = i * step1 * Vec1.Z() + j * step2 * Vec2.Z() ;
4189         myVec.SetCoord( DX, DY, DZ ) ;
4190         theTransformation.SetTranslation(myVec) ;
4191         BRepBuilderAPI_Transform myBRepTransformation(aShape, theTransformation, Standard_False) ;
4192         B.Add( compound, myBRepTransformation.Shape() );
4193       }
4194     }
4195     tds = compound ;
4196     result = CreateObject(tds) ;
4197   }
4198   catch(Standard_Failure) {
4199     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeMultiTranslation2D", SALOME::BAD_PARAM);
4200   }
4201
4202   if( CORBA::is_nil(result) ) {
4203     THROW_SALOME_CORBA_EXCEPTION("MakeMultiTranslation2D aborted : null result", SALOME::BAD_PARAM);
4204   }
4205   InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4206   return result;
4207 }
4208
4209
4210 //=================================================================================
4211 // function : MakeMultiRotation1D()
4212 // purpose  : Multi-Rotate a 3D shape
4213 //=================================================================================
4214 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeMultiRotation1D( GEOM::GEOM_Shape_ptr myShape,
4215                                                 const GEOM::DirStruct& dir,
4216                                                 const GEOM::PointStruct& loc,
4217                                                 CORBA::Short nbtimes)
4218   throw (SALOME::SALOME_Exception)
4219 {
4220   Unexpect aCatch(SALOME_SalomeException);
4221   GEOM::GEOM_Shape_var result ;
4222   TopoDS_Shape tds ;
4223   TopoDS_Shape aShape = GetTopoShape(myShape) ;
4224   if( aShape.IsNull() ) {
4225     THROW_SALOME_CORBA_EXCEPTION("MakeMultiRotation1D aborted : null shape", SALOME::BAD_PARAM);
4226   }
4227
4228   try {
4229
4230     int i ;
4231     gp_Pnt P(loc.x, loc.y, loc.z) ;
4232     gp_Dir D(dir.PS.x, dir.PS.y, dir.PS.z) ;
4233     gp_Ax1 AX1(P, D) ;
4234
4235     double angle = 360.0/nbtimes ;
4236     gp_Trsf theTransformation ;
4237     TopoDS_Compound compound;
4238     BRep_Builder B;
4239     B.MakeCompound( compound );
4240
4241     for ( i = 0; i < nbtimes; i++ ) {
4242       theTransformation.SetRotation(AX1, i*angle*PI180) ;
4243       BRepBuilderAPI_Transform myBRepTransformation(aShape, theTransformation, Standard_False) ;
4244       B.Add( compound, myBRepTransformation.Shape() );
4245     }
4246     tds = compound ;
4247     result = CreateObject(tds) ;
4248   }
4249   catch(Standard_Failure) {
4250     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeMultiRotation1D", SALOME::BAD_PARAM);
4251   }
4252
4253   if( CORBA::is_nil(result) ) {
4254     THROW_SALOME_CORBA_EXCEPTION("MakeMultiRotation1D aborted : null result", SALOME::BAD_PARAM);
4255   }
4256   InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4257   return result;
4258 }
4259
4260
4261 //=================================================================================
4262 // function : MakeMultiRotation2D()
4263 // purpose  : Multi-Rotate a 3D shape
4264 //=================================================================================
4265 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeMultiRotation2D( GEOM::GEOM_Shape_ptr myShape,
4266                                                 const GEOM::DirStruct& dir,
4267                                                 const GEOM::PointStruct& loc,
4268                                                 CORBA::Double ang,
4269                                                 CORBA::Short nbtimes1,
4270                                                 CORBA::Double step,
4271                                                 CORBA::Short nbtimes2 )
4272   throw (SALOME::SALOME_Exception)
4273 {
4274   Unexpect aCatch(SALOME_SalomeException);
4275   GEOM::GEOM_Shape_var result ;
4276   TopoDS_Shape tds ;
4277   TopoDS_Shape aShape = GetTopoShape(myShape) ;
4278   if( aShape.IsNull() ) {
4279     THROW_SALOME_CORBA_EXCEPTION("MakeMultiRotation2D aborted : null shape", SALOME::BAD_PARAM);
4280   }
4281
4282   try {
4283
4284     int i, j ;
4285     double DX, DY, DZ ;
4286     gp_Pnt P(loc.x, loc.y, loc.z) ;
4287     gp_Dir D(dir.PS.x, dir.PS.y, dir.PS.z) ;
4288     gp_Ax1 AX1(P, D) ;
4289     gp_Trsf theTransformation1 ;
4290     gp_Trsf theTransformation2 ;
4291     gp_Pnt P1 ;
4292     GProp_GProps System ;
4293
4294     if ( aShape.ShapeType() == TopAbs_VERTEX) {
4295       P1 = BRep_Tool::Pnt(TopoDS::Vertex( aShape ));
4296     }
4297     else if ( aShape.ShapeType() == TopAbs_EDGE || aShape.ShapeType() == TopAbs_WIRE ) {
4298       BRepGProp::LinearProperties(aShape, System);
4299       P1 = System.CentreOfMass() ;
4300     }
4301     else if ( aShape.ShapeType() == TopAbs_FACE || aShape.ShapeType() == TopAbs_SHELL ) {
4302       BRepGProp::SurfaceProperties(aShape, System);
4303       P1 = System.CentreOfMass() ;
4304     }
4305     else {
4306       BRepGProp::VolumeProperties(aShape, System);
4307       P1 = System.CentreOfMass() ;
4308     }
4309
4310     Handle(Geom_Line) Line = new Geom_Line(AX1);
4311     gp_Pnt P2 = GeomAPI_ProjectPointOnCurve( P1, Line ) ;
4312
4313     if ( P1.IsEqual(P2, Precision::Confusion() ) )
4314       THROW_SALOME_CORBA_EXCEPTION("Points are confused", SALOME::BAD_PARAM);
4315
4316     gp_Vec Vec(P1.X()-P2.X(), P1.Y()-P2.Y(), P1.Z()-P2.Z()) ;
4317     Vec.Normalize();
4318
4319     gp_Vec myVec ;
4320     TopoDS_Compound compound;
4321     BRep_Builder B;
4322     B.MakeCompound( compound );
4323
4324     for ( i = 0; i < nbtimes2; i++ ) {
4325       for ( j = 0; j < nbtimes1; j++ ) {
4326         DX = i * step * Vec.X() ;
4327         DY = i * step * Vec.Y() ;
4328         DZ = i * step * Vec.Z() ;
4329         myVec.SetCoord( DX, DY, DZ ) ;
4330         theTransformation1.SetTranslation(myVec) ;
4331         theTransformation2.SetRotation(AX1, j*ang*PI180) ;
4332         BRepBuilderAPI_Transform myBRepTransformation1(aShape, theTransformation1, Standard_False) ;
4333         BRepBuilderAPI_Transform myBRepTransformation2(myBRepTransformation1.Shape(), theTransformation2, Standard_False) ;
4334         B.Add( compound, myBRepTransformation2.Shape() );
4335       }
4336     }
4337     tds = compound ;
4338     result = CreateObject(tds) ;
4339   }
4340   catch(Standard_Failure) {
4341     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeMultiRotation2D", SALOME::BAD_PARAM);
4342   }
4343
4344   if( CORBA::is_nil(result) ) {
4345     THROW_SALOME_CORBA_EXCEPTION("MakeMultiRotation2D aborted : null result", SALOME::BAD_PARAM);
4346   }
4347   InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4348   return result;
4349 }
4350
4351
4352 //=================================================================================
4353 // function : MakeCopy()
4354 // purpose  : Copy a 3D shape
4355 //=================================================================================
4356 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeCopy( GEOM::GEOM_Shape_ptr Shape)
4357   throw (SALOME::SALOME_Exception)
4358 {
4359   Unexpect aCatch(SALOME_SalomeException);
4360   GEOM::GEOM_Shape_var result ;
4361   TopoDS_Shape tds ;
4362   TopoDS_Shape aShape = GetTopoShape(Shape) ;
4363   if(aShape.IsNull() ) {
4364     THROW_SALOME_CORBA_EXCEPTION("Copy aborted : null shape during operation", SALOME::BAD_PARAM);
4365   }
4366   BRepBuilderAPI_Copy Copy(aShape);
4367   if( Copy.IsDone() ) {
4368     tds = Copy.Shape();
4369     result = CreateObject(tds);
4370     InsertInLabelOneArgument(aShape, Shape, tds, result, myCurrentOCAFDoc) ;
4371   }
4372
4373   return result;
4374 }
4375
4376
4377 //=================================================================================
4378 // function : MakeMirrorByPlane()
4379 // purpose  : build a shape by symmetry of 'myShape' with 'shapePlane' in argument
4380 //=================================================================================
4381 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeMirrorByPlane(GEOM::GEOM_Shape_ptr myShape,
4382                                              GEOM::GEOM_Shape_ptr shapePlane)
4383   throw (SALOME::SALOME_Exception)
4384 {
4385   Unexpect aCatch(SALOME_SalomeException);
4386   GEOM::GEOM_Shape_var result ;
4387   TopoDS_Shape tds ;
4388   TopoDS_Shape aShape      = GetTopoShape(myShape) ;
4389   TopoDS_Shape aShapePlane = GetTopoShape(shapePlane) ;
4390   if( aShape.IsNull()  || aShapePlane.IsNull() ) {
4391     THROW_SALOME_CORBA_EXCEPTION("Mirror aborted : null shape argument", SALOME::BAD_PARAM);
4392   }
4393
4394   try {
4395     Handle(Geom_Surface) surf = BRep_Tool::Surface(TopoDS::Face(aShapePlane)) ;
4396     Handle(Geom_Plane) myPlane = Handle(Geom_Plane)::DownCast(surf) ;
4397     const gp_Ax3 pos = myPlane->Position() ;
4398     const gp_Pnt loc = pos.Location() ;  /* location of the plane */
4399     const gp_Dir dir = pos.Direction() ; /* Main direction of the plane (Z axis) */
4400
4401     /* plane used for mirroring */
4402     gp_Ax2 pln(loc, dir) ;
4403     gp_Trsf theTransformation ;
4404     theTransformation.SetMirror(pln) ;
4405     BRepBuilderAPI_Transform myBRepTransformation(aShape, theTransformation, Standard_False) ;
4406
4407     tds = myBRepTransformation.Shape() ;
4408     if(tds.IsNull() ) {
4409       THROW_SALOME_CORBA_EXCEPTION("Mirror aborted", SALOME::BAD_PARAM);
4410     }
4411   }
4412   catch(Standard_Failure) {
4413     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeMirrorByPlane", SALOME::BAD_PARAM);
4414   }
4415
4416   result = CreateObject(tds) ;
4417
4418   /* Insert arguments in ocaf */
4419   GEOM::GEOM_Gen::ListOfIOR_var ListShapes = new GEOM::GEOM_Gen::ListOfIOR;
4420   ListShapes->length(2);
4421   ListShapes[0] = GetStringFromIOR(GEOM::GEOM_Shape::_duplicate(myShape)) ;
4422   ListShapes[1] = GetStringFromIOR(GEOM::GEOM_Shape::_duplicate(shapePlane)) ;
4423   InsertInLabelMoreArguments(tds, result, ListShapes, myCurrentOCAFDoc) ;
4424   return result ;
4425 }
4426
4427
4428
4429 //=================================================================================
4430 // function : MakeRotation()
4431 // purpose  : Rotation of a 3D shape around an axis
4432 //=================================================================================
4433 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeRotation( GEOM::GEOM_Shape_ptr myShape,
4434                                          const GEOM::AxisStruct& axis,
4435                                          CORBA::Double angle)
4436   throw (SALOME::SALOME_Exception)
4437 {
4438   Unexpect aCatch(SALOME_SalomeException);
4439   GEOM::GEOM_Shape_var result ;
4440   TopoDS_Shape tds ;
4441   TopoDS_Shape aShape = GetTopoShape(myShape) ;
4442   if( aShape.IsNull() ) {
4443     THROW_SALOME_CORBA_EXCEPTION("Rotation aborted : null shape during operation", SALOME::BAD_PARAM);
4444   }
4445
4446   try {
4447     gp_Pnt P(axis.x, axis.y, axis.z) ;
4448     gp_Dir D(axis.vx, axis.vy, axis.vz) ;
4449     gp_Ax1 AX(P, D) ;
4450
4451     gp_Trsf theTransformation ;
4452     theTransformation.SetRotation(AX, angle) ;
4453     BRepBuilderAPI_Transform myBRepTransformation(aShape, theTransformation, Standard_False) ;
4454     tds = myBRepTransformation.Shape() ;
4455   }
4456   catch(Standard_Failure) {
4457     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeRotation", SALOME::BAD_PARAM);
4458   }
4459
4460   if ( !tds.IsNull() ) {
4461     result = CreateObject(tds) ;
4462     InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4463   }
4464   return result ;
4465 }
4466
4467
4468 //=================================================================================
4469 // function : MakeScaleTransform()
4470 // purpose  : Make a shape multipling another by a scale factor
4471 //=================================================================================
4472 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeScaleTransform(GEOM::GEOM_Shape_ptr myShape,
4473                                               const GEOM::PointStruct& theCenterOfScale,
4474                                               CORBA::Double factor)
4475   throw (SALOME::SALOME_Exception)
4476 {
4477   Unexpect aCatch(SALOME_SalomeException);
4478   GEOM::GEOM_Shape_var result ;
4479   TopoDS_Shape tds ;
4480   TopoDS_Shape aShape = GetTopoShape(myShape) ;
4481   if( aShape.IsNull() ) {
4482     THROW_SALOME_CORBA_EXCEPTION("Scale aborted : null shape during operation", SALOME::BAD_PARAM);
4483   }
4484
4485   try {
4486     gp_Pnt Pcenter(theCenterOfScale.x, theCenterOfScale.y, theCenterOfScale.z) ;
4487     gp_Trsf theTransformation ;
4488     theTransformation.SetScale(Pcenter, factor) ;
4489     BRepBuilderAPI_Transform myBRepTransformation(aShape, theTransformation, Standard_False) ;
4490     tds = myBRepTransformation.Shape() ;
4491   }
4492   catch(Standard_Failure) {
4493     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeScaleTransform", SALOME::BAD_PARAM);
4494   }
4495
4496   if ( !tds.IsNull() ) {
4497     result = CreateObject(tds) ;
4498     InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4499   }
4500   return result ;
4501 }
4502
4503
4504 //=================================================================================
4505 // function : MakeCompound()
4506 // purpose  : Make a compound from a list containing one or more shapes
4507 //=================================================================================
4508 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeCompound( const GEOM::GEOM_Gen::ListOfIOR& ListShapes )
4509   throw (SALOME::SALOME_Exception)
4510 {
4511   Unexpect aCatch(SALOME_SalomeException);
4512   GEOM::GEOM_Shape_var result ;
4513   TopoDS_Compound C;
4514   BRep_Builder aBuilder;
4515   aBuilder.MakeCompound(C) ;
4516
4517   for ( unsigned int i = 0; i < ListShapes.length(); i++) {
4518     GEOM::GEOM_Shape_var aShape = GetIORFromString( ListShapes[i] );
4519     TopoDS_Shape Shape = GetTopoShape(aShape) ;
4520     if( Shape.IsNull() ) {
4521        THROW_SALOME_CORBA_EXCEPTION("Compound aborted : null shape during operation", SALOME::BAD_PARAM);
4522     }
4523     aBuilder.Add(C, Shape) ;
4524   }
4525
4526   if ( C.IsNull() ) {
4527     THROW_SALOME_CORBA_EXCEPTION("Null result : Compound operation aborted", SALOME::BAD_PARAM);
4528   }
4529   else {
4530     result = CreateObject(C) ;
4531     InsertInLabelMoreArguments(C, result, ListShapes, myCurrentOCAFDoc) ;
4532   }
4533   return result;
4534 }
4535
4536
4537 //================================================================================
4538 // function : MakeEdge()
4539 // purpose  : Make a linear edge with 2 points
4540 //================================================================================
4541 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeEdge(const GEOM::PointStruct& pstruct1,
4542                                     const GEOM::PointStruct& pstruct2)
4543   throw (SALOME::SALOME_Exception)
4544 {
4545   Unexpect aCatch(SALOME_SalomeException);
4546   GEOM::GEOM_Shape_var result  ;
4547   TopoDS_Shape tds ;
4548
4549   try {
4550     gp_Pnt P1(pstruct1.x, pstruct1.y, pstruct1.z);
4551     gp_Pnt P2(pstruct2.x, pstruct2.y, pstruct2.z) ;
4552     tds = BRepBuilderAPI_MakeEdge(P1, P2).Shape();
4553     if ( tds.IsNull() )
4554       THROW_SALOME_CORBA_EXCEPTION("MakeEdge aborted : null result", SALOME::BAD_PARAM);
4555   }
4556   catch (Standard_Failure) {
4557     THROW_SALOME_CORBA_EXCEPTION("Exception catched in MakeEdge", SALOME::BAD_PARAM);
4558   }
4559
4560   result = CreateObject(tds) ;
4561   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
4562   result->ShapeId(entry);
4563   return result ;
4564 }
4565
4566
4567
4568 //=================================================================================
4569 // function : MakeWire()
4570 // purpose  : Make a wire from a list containing one or more edges or wires that can
4571 // be connected
4572 //=================================================================================
4573 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeWire( const GEOM::GEOM_Gen::ListOfIOR& ListShapes )
4574   throw (SALOME::SALOME_Exception)
4575 {
4576   Unexpect aCatch(SALOME_SalomeException);
4577   GEOM::GEOM_Shape_var result ;
4578   BRepBuilderAPI_MakeWire MW ;
4579   TopoDS_Shape tds, Shape ;
4580
4581   try {
4582     for ( unsigned int i = 0; i < ListShapes.length(); i++) {
4583       GEOM::GEOM_Shape_var aShape = GetIORFromString( ListShapes[i] );
4584       Shape = GetTopoShape(aShape) ;
4585       if( Shape.IsNull() ) {
4586         THROW_SALOME_CORBA_EXCEPTION("MakeWire aborted : null shape during operation", SALOME::BAD_PARAM);
4587       }
4588       if( Shape.ShapeType() == TopAbs_EDGE )
4589         MW.Add( TopoDS::Edge(Shape) ) ;
4590       if (Shape.ShapeType() == TopAbs_WIRE )
4591         MW.Add( TopoDS::Wire(Shape) ) ;
4592     }
4593     tds = MW  ;
4594
4595   }
4596   catch(Standard_Failure) {
4597     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeWire", SALOME::BAD_PARAM);
4598   }
4599
4600   if( tds.IsNull() ) {
4601     THROW_SALOME_CORBA_EXCEPTION("Make Wire operation aborted : null result", SALOME::BAD_PARAM);
4602   }
4603   else {
4604     result = CreateObject(tds) ;
4605     InsertInLabelMoreArguments(tds, result, ListShapes, myCurrentOCAFDoc) ;
4606   }
4607   return result;
4608 }
4609
4610
4611 //=================================================================================
4612 // function : MakeRevolution()
4613 // purpose  :
4614 //=================================================================================
4615 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeRevolution(GEOM::GEOM_Shape_ptr myShape,
4616                                           const GEOM::AxisStruct& axis,
4617                                           double angle)
4618   throw (SALOME::SALOME_Exception)
4619 {
4620   Unexpect aCatch(SALOME_SalomeException);
4621   GEOM::GEOM_Shape_var result ;
4622   TopoDS_Shape tds ;
4623   TopoDS_Shape aShape = GetTopoShape(myShape) ;
4624   if( aShape.IsNull() ) {
4625     THROW_SALOME_CORBA_EXCEPTION("Revolution aborted : null shape", SALOME::BAD_PARAM);
4626   }
4627   try {
4628     gp_Pnt P(axis.x, axis.y, axis.z) ;
4629     gp_Dir D(axis.vx, axis.vy, axis.vz);
4630     gp_Ax1 AX(P,D);
4631     tds = BRepPrimAPI_MakeRevol(aShape, AX, angle);
4632   }
4633   catch(Standard_Failure) {
4634     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeRevolution", SALOME::BAD_PARAM);
4635   }
4636
4637   if( tds.IsNull() ) {
4638     THROW_SALOME_CORBA_EXCEPTION("Revolution aborted", SALOME::BAD_PARAM);
4639   }
4640   result = CreateObject(tds) ;
4641   InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4642   return result ;
4643 }
4644
4645
4646 //=================================================================================
4647 // function : MakePipe()
4648 // purpose  : Create a shape by sweeping a baseShape along a pathShape
4649 //=================================================================================
4650 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakePipe( GEOM::GEOM_Shape_ptr pathShape,
4651                                      GEOM::GEOM_Shape_ptr baseShape )
4652   throw (SALOME::SALOME_Exception)
4653 {
4654   Unexpect aCatch(SALOME_SalomeException);
4655   GEOM::GEOM_Shape_var result ;
4656   TopoDS_Shape tds ;
4657   TopoDS_Wire aWire ;
4658   TopoDS_Shape pathTds = GetTopoShape(pathShape) ;
4659   TopoDS_Shape baseTds = GetTopoShape(baseShape) ;
4660
4661   if( baseTds.IsNull() || pathTds.IsNull() ) {
4662     THROW_SALOME_CORBA_EXCEPTION("MakePipe aborted : null shape argument", SALOME::BAD_PARAM);
4663   }
4664
4665   if( pathTds.ShapeType() == TopAbs_WIRE ) {
4666     aWire = TopoDS::Wire(pathTds) ;
4667   }
4668   else {
4669     if ( pathTds.ShapeType() == TopAbs_EDGE ) {
4670       TopoDS_Edge aEdge = TopoDS::Edge(pathTds) ;
4671       aWire = BRepBuilderAPI_MakeWire(aEdge);
4672     }
4673     else {
4674       THROW_SALOME_CORBA_EXCEPTION("MakePipe aborted : bad shape type", SALOME::BAD_PARAM);
4675     }
4676   }
4677
4678   try {
4679     tds = BRepOffsetAPI_MakePipe(aWire, baseTds) ;
4680   }
4681   catch(Standard_Failure) {
4682     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakePipe", SALOME::BAD_PARAM);
4683   }
4684
4685   if (  !IsValid(tds) ) {
4686     THROW_SALOME_CORBA_EXCEPTION("MakePipe aborted : non valid shape result", SALOME::BAD_PARAM);
4687   }
4688   else {
4689     result = CreateObject(tds) ;
4690
4691     /* Insert arguments in ocaf */
4692     GEOM::GEOM_Gen::ListOfIOR_var ListShapes = new GEOM::GEOM_Gen::ListOfIOR;
4693     ListShapes->length(2);
4694     ListShapes[0] = GetStringFromIOR(GEOM::GEOM_Shape::_duplicate(pathShape)) ;
4695     ListShapes[1] = GetStringFromIOR(GEOM::GEOM_Shape::_duplicate(baseShape)) ;
4696     InsertInLabelMoreArguments(tds, result, ListShapes, myCurrentOCAFDoc) ;
4697   }
4698   return result ;
4699 }
4700
4701
4702 //=================================================================================
4703 // function : MakePrism()
4704 // purpose  : uses myShape as base and the vector P1 to P2
4705 //=================================================================================
4706 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakePrism( GEOM::GEOM_Shape_ptr myShape,
4707                                            const GEOM::PointStruct& P1,
4708                                            const GEOM::PointStruct& P2 )
4709   throw (SALOME::SALOME_Exception)
4710 {
4711   Unexpect aCatch(SALOME_SalomeException);
4712   GEOM::GEOM_Shape_var result ;
4713   TopoDS_Shape tds ;
4714   TopoDS_Shape aShape = GetTopoShape(myShape) ;
4715   if( aShape.IsNull() ) {
4716     THROW_SALOME_CORBA_EXCEPTION("Prism aborted : null shape operation", SALOME::BAD_PARAM);
4717   }
4718
4719   try {
4720     gp_Vec Vector (P2.x - P1.x, P2.y - P1.y, P2.z - P1.z) ;
4721     tds = BRepPrimAPI_MakePrism(aShape, Vector, Standard_False).Shape() ;
4722   }
4723   catch(Standard_Failure) {
4724     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakePipe", SALOME::BAD_PARAM);
4725   }
4726
4727   if ( tds.IsNull() ) {
4728     THROW_SALOME_CORBA_EXCEPTION("Prism aborted", SALOME::BAD_PARAM);
4729   }
4730   else {
4731     result = CreateObject(tds) ;
4732     InsertInLabelOneArgument(aShape, myShape, tds, result, myCurrentOCAFDoc) ;
4733   }
4734   return result ;
4735 }
4736
4737
4738 //=================================================================================
4739 // function : MakeCDG()
4740 // purpose  : Create a CDG topology.
4741 //=================================================================================
4742 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeCDG(GEOM::GEOM_Shape_ptr aShape)
4743   throw (SALOME::SALOME_Exception)
4744 {
4745   Unexpect aCatch(SALOME_SalomeException);
4746   GEOM::GEOM_Shape_var result ;
4747   TopoDS_Shape tds ;
4748   TopoDS_Shape shape = GetTopoShape(aShape) ;
4749   GProp_GProps System;
4750   gp_Pnt myCenterMass ;
4751
4752   if( shape.IsNull() ) {
4753     THROW_SALOME_CORBA_EXCEPTION("MakeCDG aborted : null shape argument", SALOME::BAD_PARAM);
4754   }
4755
4756   try {
4757     if ( shape.ShapeType() == TopAbs_VERTEX) {
4758       myCenterMass = BRep_Tool::Pnt(TopoDS::Vertex( shape ));
4759     }
4760     else if ( shape.ShapeType() == TopAbs_EDGE || shape.ShapeType() == TopAbs_WIRE ) {
4761       BRepGProp::LinearProperties(shape, System);
4762       myCenterMass = System.CentreOfMass() ;
4763     }
4764     else if ( shape.ShapeType() == TopAbs_FACE || shape.ShapeType() == TopAbs_SHELL ) {
4765       BRepGProp::SurfaceProperties(shape, System);
4766       myCenterMass = System.CentreOfMass() ;
4767     }
4768     else {
4769       BRepGProp::VolumeProperties(shape, System);
4770       myCenterMass = System.CentreOfMass() ;
4771     }
4772
4773     tds = BRepBuilderAPI_MakeVertex(myCenterMass).Shape() ;
4774   }
4775   catch(Standard_Failure) {
4776     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeCDG", SALOME::BAD_PARAM);
4777   }
4778
4779   if ( tds.IsNull() ) {
4780     THROW_SALOME_CORBA_EXCEPTION("Make CDG aborted : null shape result", SALOME::BAD_PARAM);
4781   }
4782   else {
4783     result = CreateObject(tds) ;
4784     InsertInLabelOneArgument(shape, aShape, tds, result, myCurrentOCAFDoc) ;
4785   }
4786   return result ;
4787 }
4788
4789
4790 //=================================================================================
4791 // function : Archimede()
4792 // purpose  :
4793 //=================================================================================
4794 GEOM::GEOM_Shape_ptr GEOM_Gen_i::Archimede(GEOM::GEOM_Shape_ptr aShape,
4795                                      CORBA::Double aWeight,
4796                                      CORBA::Double aWaterDensity,
4797                                      CORBA::Double aMeshingDeflection)
4798   throw (SALOME::SALOME_Exception)
4799 {
4800   Unexpect aCatch(SALOME_SalomeException);
4801   GEOM::GEOM_Shape_var result;
4802
4803   double cste = -1;
4804   if (aWaterDensity != 0.)
4805     cste = aWeight/aWaterDensity;
4806   else
4807     THROW_SALOME_CORBA_EXCEPTION("Water density is null", SALOME::BAD_PARAM);
4808
4809   TopoDS_Shape shape = GetTopoShape(aShape) ;
4810   if( shape.IsNull() ) {
4811     THROW_SALOME_CORBA_EXCEPTION("Shape is null", SALOME::BAD_PARAM);
4812   }
4813
4814   gp_Dir direct(0.0,0.0,1.0);
4815   gp_Pnt PosPlan(0.0,0.0,0.0);
4816   Geom_Plane PP (PosPlan,direct);
4817   Handle(Geom_Geometry) G = PP.Copy();
4818   Handle(Geom_Plane) P = Handle(Geom_Plane)::DownCast(G);
4819
4820   gp_Dir Zdirection(0.0,0.0,1.0);
4821   VolumeSection VOL( shape, aMeshingDeflection);
4822   VOL.SetPlane(P);
4823   Handle (Geom_RectangularTrimmedSurface) SurfaceTrimmee;
4824
4825   if(Zdirection.IsEqual(direct,Precision::Angular()) == Standard_False) {
4826     VOL.MakeRotation(direct);
4827   }
4828
4829   VOL.CenterOfGravity();
4830   SurfaceTrimmee = VOL.TrimSurf();
4831   Standard_Real Cote = VOL.Archimede( cste, aMeshingDeflection );
4832
4833   if ( Cote == -1 ) {
4834     double Zmin,Zmax;
4835     VOL.getZ(Zmin,Zmax);
4836     double volume = VOL.CalculateVolume( Zmax ) * aWaterDensity;
4837
4838     char msg[100]="";
4839     sprintf(msg, "shape sinks to the bottom : Weigth max = %.1f", volume);
4840
4841     THROW_SALOME_CORBA_EXCEPTION(msg, SALOME::BAD_PARAM);
4842   }
4843
4844   SurfaceTrimmee=VOL.AjustePlan(SurfaceTrimmee,Cote,PosPlan);
4845   if(Zdirection.IsEqual(direct,Precision::Angular()) == Standard_False) {
4846     SurfaceTrimmee=VOL.InvMakeRotation(direct,SurfaceTrimmee);
4847   }
4848
4849   Standard_Real u1,u2,v1,v2;
4850   SurfaceTrimmee->Bounds(u1,u2,v1,v2);
4851   TopoDS_Face tirant = BRepBuilderAPI_MakeFace(SurfaceTrimmee, u1, u2, v1, v2);
4852
4853   if (tirant.IsNull()) {
4854     THROW_SALOME_CORBA_EXCEPTION("Result is null", SALOME::BAD_PARAM);
4855   }
4856
4857   result = CreateObject(tirant);
4858   InsertInLabelOneArgument(shape, aShape, tirant, result, myCurrentOCAFDoc) ;
4859
4860   return result;
4861 }
4862
4863
4864 //================================================================================
4865 // function : MakeFillet()
4866 // purpose  : Create a cylinder topology
4867 //================================================================================
4868 GEOM::GEOM_Shape_ptr  GEOM_Gen_i::MakeFillet( GEOM::GEOM_Shape_ptr shape,
4869                                               CORBA::Double radius,
4870                                               CORBA::Short ShapeType,
4871                                               const GEOM::GEOM_Shape::ListOfSubShapeID& ListOfID )
4872   throw (SALOME::SALOME_Exception)
4873 {
4874   Unexpect aCatch(SALOME_SalomeException);
4875   GEOM::GEOM_Shape_var result;
4876   TopoDS_Shape tds ;
4877
4878   const TopoDS_Shape aShape = GetTopoShape(shape) ;
4879   if( aShape.IsNull() ) {
4880     THROW_SALOME_CORBA_EXCEPTION("Shape is null", SALOME::BAD_PARAM);
4881   }
4882
4883   BRepFilletAPI_MakeFillet fill(aShape);
4884
4885   try {
4886     /* case all */
4887     if(ListOfID.length() == 0) {
4888       TopExp_Explorer Exp ( aShape, TopAbs_EDGE );
4889       for (Exp; Exp.More(); Exp.Next()) {
4890         TopoDS_Edge E =TopoDS::Edge(Exp.Current());
4891         fill.Add(E);
4892       }
4893       for (int i = 1;i<=fill.NbContours();i++) {
4894         fill.SetRadius(radius,i,i);
4895       }
4896       tds = fill.Shape();
4897
4898     } else {
4899
4900       /* case selection */
4901       for ( unsigned int ind = 0; ind < ListOfID.length(); ind++ ) {
4902         TopoDS_Shape ss ;
4903         if( GetShapeFromIndex( aShape, (TopAbs_ShapeEnum)ShapeType, ListOfID[ind], ss ) ) {
4904           TopoDS_Edge E = TopoDS::Edge(ss) ;
4905           fill.Add( E );
4906         }
4907       }
4908       for (int i = 1;i<=fill.NbContours();i++) {
4909         fill.SetRadius(radius,i,i);
4910       }
4911       tds = fill.Shape();
4912     }
4913   }
4914   catch(Standard_Failure) {
4915     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeFillet", SALOME::BAD_PARAM);
4916   }
4917
4918   if (tds.IsNull()) {
4919     THROW_SALOME_CORBA_EXCEPTION("Make Fillet aborted", SALOME::BAD_PARAM);
4920   }
4921   result = CreateObject(tds);
4922   InsertInLabelOneArgument(aShape, shape, tds, result, myCurrentOCAFDoc) ;
4923
4924   return result ;
4925 }
4926
4927
4928 //================================================================================
4929 // function : MakeChamfer
4930 // purpose  : Create a Chamfer topology
4931 //================================================================================
4932 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakeChamfer( GEOM::GEOM_Shape_ptr shape,
4933                                               CORBA::Double d1,
4934                                               CORBA::Double d2,
4935                                               CORBA::Short ShapeType,
4936                                               const GEOM::GEOM_Shape::ListOfSubShapeID& ListOfID )
4937   throw (SALOME::SALOME_Exception)
4938 {
4939   Unexpect aCatch(SALOME_SalomeException);
4940   GEOM::GEOM_Shape_var result;
4941   TopoDS_Shape tds ;
4942
4943   const TopoDS_Shape aShape = GetTopoShape(shape) ;
4944   if( aShape.IsNull() ) {
4945     THROW_SALOME_CORBA_EXCEPTION("Shape is null", SALOME::BAD_PARAM);
4946   }
4947
4948   BRepFilletAPI_MakeChamfer MC(aShape);
4949
4950   try {
4951     /* case all */
4952     TopTools_IndexedDataMapOfShapeListOfShape M;
4953     TopExp::MapShapesAndAncestors(aShape,TopAbs_EDGE,TopAbs_FACE,M);
4954     if(ListOfID.length() == 0) {
4955       for (int i = 1;i<=M.Extent();i++) {
4956         TopoDS_Edge E = TopoDS::Edge(M.FindKey(i));
4957         TopoDS_Face F = TopoDS::Face(M.FindFromIndex(i).First());
4958         if (!BRepTools::IsReallyClosed(E, F) && !BRep_Tool::Degenerated(E))
4959           MC.Add(d1,d2,E,F);
4960       }
4961       tds = MC.Shape();
4962
4963     } else {
4964
4965       /* case selection */
4966       for ( unsigned int ind = 0; ind < ListOfID.length(); ind++ ) {
4967         TopoDS_Shape ss ;
4968         if( GetShapeFromIndex( aShape, (TopAbs_ShapeEnum)ShapeType, ListOfID[ind], ss ) ) {
4969           TopoDS_Edge E = TopoDS::Edge( ss ) ;
4970           TopoDS_Face F = TopoDS::Face(M.FindFromKey(E).First());
4971           if (!BRepTools::IsReallyClosed(E, F) && !BRep_Tool::Degenerated(E))
4972             MC.Add(d1,d2,E,F);
4973         }
4974       }
4975       tds = MC.Shape();
4976     }
4977   }
4978   catch(Standard_Failure) {
4979     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeChamfer", SALOME::BAD_PARAM);
4980   }
4981
4982   if (tds.IsNull()) {
4983     THROW_SALOME_CORBA_EXCEPTION("Make Chamfer aborted", SALOME::BAD_PARAM);
4984   }
4985   result = CreateObject(tds);
4986   InsertInLabelOneArgument(aShape, shape, tds, result, myCurrentOCAFDoc) ;
4987
4988   return result ;
4989 }
4990
4991 //=================================================================================
4992 // function : CheckShape()
4993 // purpose  :
4994 //=================================================================================
4995 CORBA::Boolean GEOM_Gen_i::CheckShape(GEOM::GEOM_Shape_ptr shape)
4996   throw (SALOME::SALOME_Exception)
4997 {
4998   Unexpect aCatch(SALOME_SalomeException);
4999   TopoDS_Shape S = GetTopoShape(shape) ;
5000   if( S.IsNull() ) {
5001     THROW_SALOME_CORBA_EXCEPTION("Shape is null", SALOME::BAD_PARAM);
5002   }
5003
5004   BRepCheck_Analyzer ana(S,false);
5005   if (ana.IsValid())
5006     return 1;
5007
5008   return 0;
5009 }
5010
5011 //=================================================================================
5012 // function : MakePlacedBox()
5013 // purpose  :
5014 //=================================================================================
5015 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakePlacedBox(CORBA::Double x1,  CORBA::Double y1,  CORBA::Double z1,
5016                                          CORBA::Double delta1, CORBA::Double delta2, CORBA::Double delta3)
5017   throw (SALOME::SALOME_Exception)
5018 {
5019   Unexpect aCatch(SALOME_SalomeException);
5020   GEOM::GEOM_Shape_var result ;
5021   TopoDS_Shape tds ;
5022
5023   CORBA::Double x2, y2, z2 ;
5024
5025   try {
5026     x2 = x1 + delta1 ;
5027     y2 = y1 + delta2 ;
5028     z2 = z1 + delta3 ;
5029
5030     gp_Pnt P1(x1,y1,z1);
5031     gp_Pnt P2(x2,y2,z2);
5032
5033     tds = BRepPrimAPI_MakeBox(P1,P2).Shape();
5034   }
5035   catch(Standard_Failure) {
5036     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakeBox", SALOME::BAD_PARAM);
5037   }
5038
5039   if (tds.IsNull()) {
5040     THROW_SALOME_CORBA_EXCEPTION("Make Box aborted : null shape", SALOME::BAD_PARAM);
5041   }
5042
5043   result = CreateObject(tds);
5044   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
5045   result->ShapeId(entry) ;
5046
5047   return result;
5048 }
5049
5050 //=================================================================================
5051 // function : MakePanel()
5052 // purpose  :
5053 //=================================================================================
5054 GEOM::GEOM_Shape_ptr GEOM_Gen_i::MakePanel(GEOM::GEOM_Shape_ptr shape,
5055                                      CORBA::Short directiontype,
5056                                      CORBA::Double delta)
5057   throw (SALOME::SALOME_Exception)
5058 {
5059   Unexpect aCatch(SALOME_SalomeException);
5060   GEOM::GEOM_Shape_var result ;
5061   TopoDS_Shape tds ;
5062   TopoDS_Shape aShape = GetTopoShape(shape) ;
5063   Bnd_Box B ;
5064   Standard_Real axmin,aymin,azmin,axmax,aymax,azmax ;
5065   GEOM::PointStruct pstruct1, pstruct2, pstruct3, pstruct4 ;
5066
5067   if(aShape.IsNull() ) {
5068     THROW_SALOME_CORBA_EXCEPTION("MakePanel aborted : null shape during operation", SALOME::BAD_PARAM);
5069   }
5070
5071   try {
5072     BRepBndLib::Add(aShape,B);
5073     B.Enlarge(10.);
5074     B.Get(axmin,aymin,azmin,axmax,aymax,azmax);
5075
5076     switch (directiontype)
5077       {
5078       case 1 :   /* X */
5079         pstruct1 = MakePointStruct( delta, aymin, azmin ) ;
5080         pstruct2 = MakePointStruct( delta, aymin, azmax ) ;
5081         pstruct3 = MakePointStruct( delta, aymax, azmax ) ;
5082         pstruct4 = MakePointStruct( delta, aymax, azmin ) ;
5083         break ;
5084       case 2 :   /* Y */
5085         pstruct1 = MakePointStruct( axmin, delta, azmin ) ;
5086         pstruct2 = MakePointStruct( axmin, delta, azmax ) ;
5087         pstruct3 = MakePointStruct( axmax, delta, azmax ) ;
5088         pstruct4 = MakePointStruct( axmax, delta, azmin ) ;
5089         break ;
5090       case 3 :   /* Z */
5091         pstruct1 = MakePointStruct( axmin, aymin, delta ) ;
5092         pstruct2 = MakePointStruct( axmin, aymax, delta ) ;
5093         pstruct3 = MakePointStruct( axmax, aymax, delta ) ;
5094         pstruct4 = MakePointStruct( axmax, aymin, delta ) ;
5095         break ;
5096       default :
5097         return result ;
5098       }
5099
5100     GEOM::GEOM_Shape_ptr Edge1 = MakeEdge(pstruct1, pstruct2);
5101     GEOM::GEOM_Shape_ptr Edge2 = MakeEdge(pstruct2, pstruct3);
5102     GEOM::GEOM_Shape_ptr Edge3 = MakeEdge(pstruct3, pstruct4);
5103     GEOM::GEOM_Shape_ptr Edge4 = MakeEdge(pstruct4, pstruct1);
5104
5105     GEOM::GEOM_Gen::ListOfIOR_var aList = new GEOM::GEOM_Gen::ListOfIOR;
5106     aList->length(4);
5107     aList[0]=CORBA::string_dup(Edge1->Name());
5108     aList[1]=CORBA::string_dup(Edge2->Name());
5109     aList[2]=CORBA::string_dup(Edge3->Name());
5110     aList[3]=CORBA::string_dup(Edge4->Name());
5111
5112     GEOM::GEOM_Shape_ptr aWire = MakeWire( aList );
5113     GEOM::GEOM_Shape_ptr aFace = MakeFace( aWire, true ) ;
5114     tds = GetTopoShape(aFace);
5115
5116   }
5117   catch(Standard_Failure) {
5118     THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::MakePanel", SALOME::BAD_PARAM);
5119   }
5120
5121   if (tds.IsNull()) {
5122     THROW_SALOME_CORBA_EXCEPTION("Make PanelsPartition aborted : null shape", SALOME::BAD_PARAM);
5123   }
5124
5125   result = CreateObject(tds);
5126   const char *entry = InsertInLabel(tds, result->Name(), myCurrentOCAFDoc) ;
5127   result->ShapeId(entry) ;
5128
5129   return result;
5130 }
5131
5132
5133 void GEOM_Gen_i::ExportIGES(const char* filename,GEOM::GEOM_Shape_ptr theShape)
5134   throw (SALOME::SALOME_Exception)
5135 {
5136   Unexpect aCatch(SALOME_SalomeException);
5137   if (theShape->_is_nil())
5138     {
5139       THROW_SALOME_CORBA_EXCEPTION("Export IGES aborted", SALOME::BAD_PARAM);
5140     }
5141   TopoDS_Shape tds = GetTopoShape(theShape);
5142   if (tds.IsNull())
5143     {
5144       THROW_SALOME_CORBA_EXCEPTION("Export IGES aborted", SALOME::BAD_PARAM);
5145     }
5146   try
5147     {
5148       //VRV: OCC 4.0 migration
5149       IGESControl_Controller::Init();
5150       IGESControl_Writer ICW (Interface_Static::CVal("XSTEP.iges.unit"),
5151                                    Interface_Static::IVal("XSTEP.iges.writebrep.mode"));
5152       //VRV: OCC 4.0 migration
5153
5154       ICW.AddShape (tds);
5155       ICW.ComputeModel();
5156       char * aname = strdup(filename);
5157       Standard_Boolean result = ICW.Write( aname );
5158       free(aname);
5159     }
5160   catch(Standard_Failure)
5161     {
5162       THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ExportIGES", SALOME::BAD_PARAM);
5163     }
5164 }
5165
5166 void GEOM_Gen_i::ExportBREP(const char* filename,GEOM::GEOM_Shape_ptr theShape)
5167   throw (SALOME::SALOME_Exception)
5168 {
5169   Unexpect aCatch(SALOME_SalomeException);
5170   if (theShape->_is_nil())
5171     {
5172       THROW_SALOME_CORBA_EXCEPTION("Export BRep aborted", SALOME::BAD_PARAM);
5173     }
5174   TopoDS_Shape tds = GetTopoShape(theShape);
5175   if (tds.IsNull())
5176     {
5177       THROW_SALOME_CORBA_EXCEPTION("Export BRep aborted", SALOME::BAD_PARAM);
5178     }
5179   try
5180     {
5181       char * aname = strdup(filename);
5182       Standard_Boolean result = BRepTools::Write(tds,aname);
5183       free(aname);
5184     }
5185   catch(Standard_Failure)
5186     {
5187       THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ExportBREP", SALOME::BAD_PARAM);
5188     }
5189 }
5190
5191 void GEOM_Gen_i::ExportSTEP(const char* filename,GEOM::GEOM_Shape_ptr theShape)
5192   throw (SALOME::SALOME_Exception)
5193 {
5194   Unexpect aCatch(SALOME_SalomeException);
5195   if (theShape->_is_nil())
5196     {
5197       THROW_SALOME_CORBA_EXCEPTION("Export STEP aborted", SALOME::BAD_PARAM);
5198     }
5199   TopoDS_Shape tds = GetTopoShape(theShape);
5200   if (tds.IsNull())
5201     {
5202       THROW_SALOME_CORBA_EXCEPTION("Export STEP aborted", SALOME::BAD_PARAM);
5203     }
5204   try
5205     {
5206       IFSelect_ReturnStatus status ;
5207       //VRV: OCC 4.0 migration
5208       STEPControl_Writer aWriter;
5209       status = aWriter.Transfer( tds, STEPControl_ManifoldSolidBrep ) ;
5210       //VRV: OCC 4.0 migration
5211       if ( status == IFSelect_RetDone )
5212         {
5213           char * aname = strdup(filename);
5214           status = aWriter.Write( aname ) ;
5215           free(aname);
5216         }
5217     }
5218   catch(Standard_Failure)
5219     {
5220       THROW_SALOME_CORBA_EXCEPTION("Exception catched in GEOM_Gen_i::ExportBREP", SALOME::BAD_PARAM);
5221     }
5222 }
5223
5224
5225 //=====================================================================================
5226 // EXPORTED METHODS
5227 //=====================================================================================
5228 extern "C"
5229 {
5230   PortableServer::ObjectId * GEOMEngine_factory(CORBA::ORB_ptr orb,
5231                                                 PortableServer::POA_ptr poa,
5232                                                 PortableServer::ObjectId * contId,
5233                                                 const char *instanceName,
5234                                                 const char * interfaceName)
5235   {
5236    GEOM_Gen_i * myGEOM_Gen_i = new GEOM_Gen_i(orb, poa, contId, instanceName, interfaceName);
5237    myGEOM_Gen_i->register_name("/myGEOM_Gen"); // NRI : 11/07/2002 : Add for Supervision example
5238    return myGEOM_Gen_i->getId() ;
5239   }
5240 }
5241