smesh/src/SMESH_I/SMESH_Gen_i_1.cxx

// Copyright (C) 2007-2013  CEA/DEN, EDF R&D, OPEN CASCADE
//
// Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
// CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
//
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
// License as published by the Free Software Foundation; either
// version 2.1 of the License.
//
// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
//
// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
//
//  File      : SMESH_Gen_i_1.cxx
//  Created   : Thu Oct 21 17:24:06 2004
//  Author    : Edward AGAPOV (eap)
//  Module    : SMESH

#include "SMESH_Gen_i.hxx"

#include "SMESH_Algo_i.hxx"
#include "SMESH_Comment.hxx"
#include "SMESH_Group_i.hxx"
#include "SMESH_Hypothesis_i.hxx"
#include "SMESH_Mesh_i.hxx"
#include "SMESH_subMesh_i.hxx"

#include CORBA_CLIENT_HEADER(SALOME_ModuleCatalog)

#include <utilities.h>
#include <Utils_ExceptHandlers.hxx>
#include <SALOMEDS_wrap.hxx>
#include <SALOMEDS_Attributes_wrap.hxx>

#include <TCollection_AsciiString.hxx>
#include <TopoDS_Solid.hxx>

#ifdef _DEBUG_
static int MYDEBUG = 0;
//static int VARIABLE_DEBUG = 0;
#else
static int MYDEBUG = 0;
//static int VARIABLE_DEBUG = 0;
#endif

//=============================================================================
/*!
 *  Get...Tag [ static ]
 *
 *  Methods which determine SMESH data model structure
 */
//=============================================================================

long SMESH_Gen_i::GetHypothesisRootTag()
{
  return SMESH::Tag_HypothesisRoot;
}

long SMESH_Gen_i::GetAlgorithmsRootTag()
{
  return SMESH::Tag_AlgorithmsRoot;
}

long SMESH_Gen_i::GetRefOnShapeTag()
{
  return SMESH::Tag_RefOnShape;
}

long SMESH_Gen_i::GetRefOnAppliedHypothesisTag()
{
  return SMESH::Tag_RefOnAppliedHypothesis;
}

long SMESH_Gen_i::GetRefOnAppliedAlgorithmsTag()
{
  return SMESH::Tag_RefOnAppliedAlgorithms;
}

long SMESH_Gen_i::GetSubMeshOnVertexTag()
{
  return SMESH::Tag_SubMeshOnVertex;
}

long SMESH_Gen_i::GetSubMeshOnEdgeTag()
{
  return SMESH::Tag_SubMeshOnEdge;
}

long SMESH_Gen_i::GetSubMeshOnFaceTag()
{
  return SMESH::Tag_SubMeshOnFace;
}

long SMESH_Gen_i::GetSubMeshOnSolidTag()
{
  return SMESH::Tag_SubMeshOnSolid;
}

long SMESH_Gen_i::GetSubMeshOnCompoundTag()
{
  return SMESH::Tag_SubMeshOnCompound;
}

long SMESH_Gen_i::GetSubMeshOnWireTag()
{
  return SMESH::Tag_SubMeshOnWire;
}

long SMESH_Gen_i::GetSubMeshOnShellTag()
{
  return SMESH::Tag_SubMeshOnShell;
}

long SMESH_Gen_i::GetNodeGroupsTag()
{
  return SMESH::Tag_NodeGroups;
}

long SMESH_Gen_i::GetEdgeGroupsTag()
{
  return SMESH::Tag_EdgeGroups;
}

long SMESH_Gen_i::GetFaceGroupsTag()
{
  return SMESH::Tag_FaceGroups;
}

long SMESH_Gen_i::GetVolumeGroupsTag()
{
  return SMESH::Tag_VolumeGroups;
}

long SMESH_Gen_i::Get0DElementsGroupsTag()
{
  return SMESH::Tag_0DElementsGroups;
}

long SMESH_Gen_i::GetBallElementsGroupsTag()
{
  return SMESH::Tag_BallElementsGroups;
}

//=============================================================================
/*!
 *  SMESH_Gen_i::CanPublishInStudy
 *
 *  Returns true if object can be published in the study
 */
//=============================================================================

bool SMESH_Gen_i::CanPublishInStudy(CORBA::Object_ptr theIOR)
{
  if(MYDEBUG) MESSAGE("CanPublishInStudy - "<<!CORBA::is_nil(myCurrentStudy));
  if(CORBA::is_nil(myCurrentStudy))
    return false;
  
  SMESH::SMESH_Mesh_var aMesh       = SMESH::SMESH_Mesh::_narrow(theIOR);
  if( !aMesh->_is_nil() )
    return true;

  SMESH::SMESH_subMesh_var aSubMesh = SMESH::SMESH_subMesh::_narrow(theIOR);
  if( !aSubMesh->_is_nil() )
    return true;

  SMESH::SMESH_Hypothesis_var aHyp  = SMESH::SMESH_Hypothesis::_narrow(theIOR);
  if( !aHyp->_is_nil() )
    return true;

  SMESH::SMESH_GroupBase_var aGroup = SMESH::SMESH_GroupBase::_narrow(theIOR);
  if( !aGroup->_is_nil() )
    return true;

  if(MYDEBUG) MESSAGE("CanPublishInStudy--CANT");
  return false;
}

//=======================================================================
//function : ObjectToSObject
//purpose  : Put a result into a SALOMEDS::SObject_wrap or call UnRegister()!
//=======================================================================

SALOMEDS::SObject_ptr SMESH_Gen_i::ObjectToSObject(SALOMEDS::Study_ptr theStudy,
                                                   CORBA::Object_ptr   theObject)
{
  SALOMEDS::SObject_wrap aSO;
  if ( !CORBA::is_nil( theStudy ) && !CORBA::is_nil( theObject ))
  {
    CORBA::String_var objStr = SMESH_Gen_i::GetORB()->object_to_string( theObject );
    aSO = theStudy->FindObjectIOR( objStr.in() );
  }
  return aSO._retn();
}

//=======================================================================
//function : objectToServant
//purpose  : 
//=======================================================================

template<typename T> static inline T* objectToServant( CORBA::Object_ptr theIOR )
{
  return dynamic_cast<T*>( SMESH_Gen_i::GetServant( theIOR ).in() );
}

//=======================================================================
//function : ShapeToGeomObject
//purpose  : 
//=======================================================================

GEOM::GEOM_Object_ptr SMESH_Gen_i::ShapeToGeomObject (const TopoDS_Shape& theShape )
{
  GEOM::GEOM_Object_var aShapeObj;
  if ( !theShape.IsNull() ) {
    GEOM_Client* aClient = GetShapeReader();
    TCollection_AsciiString IOR;
    if ( aClient && aClient->Find( theShape, IOR ))
    {
      CORBA::Object_var obj = GetORB()->string_to_object( IOR.ToCString() );
      aShapeObj = GEOM::GEOM_Object::_narrow ( obj );
    }
  }
  return aShapeObj._retn();
}

//=======================================================================
//function : GeomObjectToShape
//purpose  : 
//=======================================================================

TopoDS_Shape SMESH_Gen_i::GeomObjectToShape(GEOM::GEOM_Object_ptr theGeomObject)
{
  TopoDS_Shape S;
  if ( !theGeomObject->_is_nil() ) {
    GEOM_Client* aClient = GetShapeReader();
    GEOM::GEOM_Gen_ptr aGeomEngine = GetGeomEngine();
    if ( aClient && !aGeomEngine->_is_nil () )
      S = aClient->GetShape( aGeomEngine, theGeomObject );
  }
  return S;
}

//=======================================================================
//function : publish
//purpose  : 
//=======================================================================

static SALOMEDS::SObject_ptr publish(SALOMEDS::Study_ptr   theStudy,
                                     CORBA::Object_ptr     theIOR,
                                     SALOMEDS::SObject_ptr theFatherObject,
                                     const int             theTag = 0,
                                     const char*           thePixMap = 0,
                                     const bool            theSelectable = true)
{
  SALOMEDS::SObject_wrap SO = SMESH_Gen_i::ObjectToSObject( theStudy, theIOR );
  SALOMEDS::StudyBuilder_var     aStudyBuilder = theStudy->NewBuilder();
  SALOMEDS::UseCaseBuilder_wrap useCaseBuilder = theStudy->GetUseCaseBuilder();
  SALOMEDS::SObject_wrap objAfter;
  if ( SO->_is_nil() )
  {
    if ( theTag == 0 ) {
      SO = aStudyBuilder->NewObject( theFatherObject );
    }
    else if ( !theFatherObject->FindSubObject( theTag, SO.inout() ))
    {
      SO = aStudyBuilder->NewObjectToTag( theFatherObject, theTag );

      // define the next tag after given one in the data tree to insert SObject
      std::string anEntry;
      int last2Pnt_pos = -1;
      int tagAfter = -1;
      CORBA::String_var entry;
      SALOMEDS::SObject_wrap curObj;
      SALOMEDS::UseCaseIterator_wrap anUseCaseIter = useCaseBuilder->GetUseCaseIterator(theFatherObject);
      for ( ; anUseCaseIter->More(); anUseCaseIter->Next() ) {
        curObj = anUseCaseIter->Value();
        entry = curObj->GetID();
        anEntry = entry.in();
        last2Pnt_pos = anEntry.rfind( ":" );
        tagAfter = atoi( anEntry.substr( last2Pnt_pos+1 ).c_str() );
        if ( tagAfter > theTag  ) {
          objAfter = curObj;
          break;
        }
      }
    }
  }

  SALOMEDS::GenericAttribute_wrap anAttr;
  if ( !CORBA::is_nil( theIOR )) {
    anAttr = aStudyBuilder->FindOrCreateAttribute( SO, "AttributeIOR" );
    CORBA::String_var objStr = SMESH_Gen_i::GetORB()->object_to_string( theIOR );
    SALOMEDS::AttributeIOR_wrap iorAttr = anAttr;
    iorAttr->SetValue( objStr.in() );
    // UnRegister() !!!
    SALOME::GenericObj_var genObj = SALOME::GenericObj::_narrow( theIOR );
    if ( !genObj->_is_nil() )
      genObj->UnRegister();
  }
  if ( thePixMap ) {
    anAttr  = aStudyBuilder->FindOrCreateAttribute( SO, "AttributePixMap" );
    SALOMEDS::AttributePixMap_wrap pm = anAttr;
    pm->SetPixMap( thePixMap );
  }
  if ( !theSelectable ) {
    anAttr = aStudyBuilder->FindOrCreateAttribute( SO, "AttributeSelectable" );
    SALOMEDS::AttributeSelectable_wrap selAttr = anAttr;
    selAttr->SetSelectable( false );
  }

  // add object to the use case tree
  // (to support tree representation customization and drag-n-drop)
  if ( !CORBA::is_nil( objAfter ) ) {
    useCaseBuilder->InsertBefore( SO, objAfter );    // insert at given tag
  } else if ( !useCaseBuilder->IsUseCaseNode( SO ) ) {
    useCaseBuilder->AppendTo( theFatherObject, SO ); // append to the end of list
  }

  return SO._retn();
}

//=======================================================================
//function : setName
//purpose  : 
//=======================================================================

void SMESH_Gen_i::SetName(SALOMEDS::SObject_ptr theSObject,
                          const char*           theName,
                          const char*           theDefaultName)
{
  if ( !theSObject->_is_nil() ) {
    SALOMEDS::Study_var               aStudy = theSObject->GetStudy();
    SALOMEDS::StudyBuilder_var aStudyBuilder = aStudy->NewBuilder();
    SALOMEDS::GenericAttribute_wrap   anAttr =
      aStudyBuilder->FindOrCreateAttribute( theSObject, "AttributeName" );
    SALOMEDS::AttributeName_wrap aNameAttr = anAttr;
    if ( theName && strlen( theName ) != 0 )
      aNameAttr->SetValue( theName );
    else {
      CORBA::String_var curName = aNameAttr->Value();
      if ( strlen( curName.in() ) == 0 ) {
        SMESH_Comment aName(theDefaultName);
        aNameAttr->SetValue( ( aName<< "_" << theSObject->Tag()).c_str() );
      }
    }
  }
}

//=======================================================================
//function : SetPixMap
//purpose  : 
//=======================================================================

void SMESH_Gen_i::SetPixMap(SALOMEDS::SObject_ptr theSObject,
                            const char*           thePixMap)
{
  if ( !theSObject->_is_nil() && thePixMap && strlen( thePixMap ))
  {
    SALOMEDS::Study_var               aStudy = theSObject->GetStudy();
    SALOMEDS::StudyBuilder_var aStudyBuilder = aStudy->NewBuilder();
    SALOMEDS::GenericAttribute_wrap anAttr =
      aStudyBuilder->FindOrCreateAttribute( theSObject, "AttributePixMap" );
    SALOMEDS::AttributePixMap_wrap aPMAttr = anAttr;
    aPMAttr->SetPixMap( thePixMap );
  }
}

//=======================================================================
//function : addReference
//purpose  : 
//=======================================================================

static void addReference (SALOMEDS::Study_ptr   theStudy,
                          SALOMEDS::SObject_ptr theSObject,
                          CORBA::Object_ptr     theToObject,
                          int                   theTag = 0)
{
  SALOMEDS::SObject_wrap aToObjSO = SMESH_Gen_i::ObjectToSObject( theStudy, theToObject );
  if ( !aToObjSO->_is_nil() && !theSObject->_is_nil() ) {
    SALOMEDS::StudyBuilder_var aStudyBuilder = theStudy->NewBuilder();
    SALOMEDS::SObject_wrap aReferenceSO;
    if ( !theTag ) {
      // check if the reference to theToObject already exists
      // and find a free label for the reference object
      bool isReferred = false;
      int tag = 1;
      SALOMEDS::ChildIterator_wrap anIter = theStudy->NewChildIterator( theSObject );
      for ( ; !isReferred && anIter->More(); anIter->Next(), ++tag ) {
        SALOMEDS::SObject_wrap curSO = anIter->Value();
        if ( curSO->ReferencedObject( aReferenceSO.inout() )) {
          CORBA::String_var refEntry = aReferenceSO->GetID();
          CORBA::String_var  toEntry = aToObjSO->GetID();
          if ( strcmp( refEntry, toEntry ) == 0 )
            isReferred = true;
        }
        else if ( !theTag ) {
          SALOMEDS::GenericAttribute_wrap anAttr;
          if ( !curSO->FindAttribute( anAttr.inout(), "AttributeIOR" ))
            theTag = tag;
        }
      }
      if ( isReferred )
        return;
      if ( !theTag )
        theTag = tag;
    }
    if ( !theSObject->FindSubObject( theTag, aReferenceSO.inout() ))
      aReferenceSO = aStudyBuilder->NewObjectToTag( theSObject, theTag );
    aStudyBuilder->Addreference( aReferenceSO, aToObjSO );
    // add reference to the use case tree
    // (to support tree representation customization and drag-n-drop)
    SALOMEDS::UseCaseBuilder_wrap useCaseBuilder = theStudy->GetUseCaseBuilder();
    useCaseBuilder->AppendTo( aReferenceSO->GetFather(), aReferenceSO );
  }
}

//=============================================================================
/*!
 *  SMESH_Gen_i::PublishInStudy
 *
 *  Publish object in the study
 */
//=============================================================================

SALOMEDS::SObject_ptr SMESH_Gen_i::PublishInStudy(SALOMEDS::Study_ptr   theStudy,
                                                  SALOMEDS::SObject_ptr theSObject,
                                                  CORBA::Object_ptr     theIOR,
                                                  const char*           theName)
     throw (SALOME::SALOME_Exception)
{
  Unexpect aCatch(SALOME_SalomeException);
  SALOMEDS::SObject_wrap aSO;
  if ( CORBA::is_nil( theStudy ) || CORBA::is_nil( theIOR ))
    return aSO._retn();
  if(MYDEBUG) MESSAGE("PublishInStudy");

  // Publishing a mesh
  SMESH::SMESH_Mesh_var aMesh = SMESH::SMESH_Mesh::_narrow( theIOR );
  if( !aMesh->_is_nil() )
    aSO = PublishMesh( theStudy, aMesh, theName );

  // Publishing a sub-mesh
  SMESH::SMESH_subMesh_var aSubMesh = SMESH::SMESH_subMesh::_narrow( theIOR );
  if( aSO->_is_nil() && !aSubMesh->_is_nil() ) {
    GEOM::GEOM_Object_var aShapeObject = aSubMesh->GetSubShape();
    aMesh = aSubMesh->GetFather();
    aSO = PublishSubMesh( theStudy, aMesh, aSubMesh, aShapeObject, theName );
  }

  // Publishing a hypothesis or algorithm
  SMESH::SMESH_Hypothesis_var aHyp = SMESH::SMESH_Hypothesis::_narrow( theIOR );
  if ( aSO->_is_nil() && !aHyp->_is_nil() )
    aSO = PublishHypothesis( theStudy, aHyp );

  // Publishing a group
  SMESH::SMESH_GroupBase_var aGroup = SMESH::SMESH_GroupBase::_narrow(theIOR);
  if ( aSO->_is_nil() && !aGroup->_is_nil() ) {
    GEOM::GEOM_Object_var aShapeObject;
    aMesh = aGroup->GetMesh();
    aSO = PublishGroup( theStudy, aMesh, aGroup, aShapeObject, theName );
  }
  if(MYDEBUG) MESSAGE("PublishInStudy_END");

  return aSO._retn();
}

//=======================================================================
//function : PublishComponent
//purpose  : 
//=======================================================================

SALOMEDS::SComponent_ptr SMESH_Gen_i::PublishComponent(SALOMEDS::Study_ptr theStudy)
{
  if ( CORBA::is_nil( theStudy ))
    return SALOMEDS::SComponent::_nil();
  if(MYDEBUG) MESSAGE("PublishComponent");

  SALOMEDS::StudyBuilder_var    aStudyBuilder  = theStudy->NewBuilder(); 
  SALOMEDS::UseCaseBuilder_wrap useCaseBuilder = theStudy->GetUseCaseBuilder();

  CORBA::String_var   compDataType = ComponentDataType();
  SALOMEDS::SComponent_wrap father = theStudy->FindComponent( compDataType.in() );
  if ( !CORBA::is_nil( father ) ) {
    // check that the component is added to the use case browser
    if ( !useCaseBuilder->IsUseCaseNode( father ) ) {
      useCaseBuilder->SetRootCurrent();
      useCaseBuilder->Append( father ); // component object is added as the top level item
    }
    return father._retn();
  }

  SALOME_ModuleCatalog::ModuleCatalog_var aCat =
    SALOME_ModuleCatalog::ModuleCatalog::_narrow( GetNS()->Resolve("/Kernel/ModulCatalog") );
  if ( CORBA::is_nil( aCat ) )
    return father._retn();

  SALOME_ModuleCatalog::Acomponent_var aComp = aCat->GetComponent( compDataType.in() );
  if ( CORBA::is_nil( aComp ) )
    return father._retn();

  SALOMEDS::GenericAttribute_wrap anAttr;
  SALOMEDS::AttributePixMap_wrap  aPixmap;

  father  = aStudyBuilder->NewComponent( compDataType.in() );
  aStudyBuilder->DefineComponentInstance( father, SMESH_Gen::_this() );
  anAttr  = aStudyBuilder->FindOrCreateAttribute( father, "AttributePixMap" );
  aPixmap = anAttr;
  aPixmap->SetPixMap( "ICON_OBJBROWSER_SMESH" );
  CORBA::String_var userName = aComp->componentusername();
  SetName( father, userName.in(), "MESH" );
  // add component to the use case tree
  // (to support tree representation customization and drag-n-drop)
  useCaseBuilder->SetRootCurrent();
  useCaseBuilder->Append( father ); // component object is added as the top level item
  if(MYDEBUG) MESSAGE("PublishComponent--END");

  return father._retn();
}

//=============================================================================
/*!
 *  findMaxChildTag [ static internal ]
 *
 *  Finds maximum child tag for the given object
 */
//=============================================================================

static long findMaxChildTag( SALOMEDS::SObject_ptr theSObject )
{
  long aTag = 0;
  if ( !theSObject->_is_nil() ) {
    SALOMEDS::Study_var aStudy = theSObject->GetStudy();
    if ( !aStudy->_is_nil() ) {
      SALOMEDS::ChildIterator_wrap anIter = aStudy->NewChildIterator( theSObject );
      for ( ; anIter->More(); anIter->Next() ) {
        SALOMEDS::SObject_wrap anSO = anIter->Value();
        long nTag = anSO->Tag();
        if ( nTag > aTag )
          aTag = nTag;
      }
    }
  }
  return aTag;
}

//=======================================================================
//function : PublishMesh
//purpose  : 
//=======================================================================

SALOMEDS::SObject_ptr SMESH_Gen_i::PublishMesh (SALOMEDS::Study_ptr   theStudy,
                                                SMESH::SMESH_Mesh_ptr theMesh,
                                                const char*           theName)
{
  if ( CORBA::is_nil( theStudy ) ||
       CORBA::is_nil( theMesh ))
    return SALOMEDS::SComponent::_nil();
  if(MYDEBUG) MESSAGE("PublishMesh--IN");

  // find or publish a mesh

  SALOMEDS::SObject_wrap aMeshSO = ObjectToSObject( theStudy, theMesh );
  if ( aMeshSO->_is_nil() )
  {
    SALOMEDS::SComponent_wrap father = PublishComponent( theStudy );
    if ( father->_is_nil() )
      return aMeshSO._retn();

    // Find correct free tag
    long aTag = findMaxChildTag( father.in() );
    if ( aTag <= GetAlgorithmsRootTag() )
      aTag = GetAlgorithmsRootTag() + 1;
    else
      aTag++;

    aMeshSO = publish (theStudy, theMesh, father, aTag, "ICON_SMESH_TREE_MESH_WARN" );
    if ( aMeshSO->_is_nil() )
      return aMeshSO._retn();
  }
  SetName( aMeshSO, theName, "Mesh" );

  // Add shape reference

  GEOM::GEOM_Object_var aShapeObject = theMesh->GetShapeToMesh();
  if ( !CORBA::is_nil( aShapeObject )) {
    addReference( theStudy, aMeshSO, aShapeObject, GetRefOnShapeTag() );

    // Publish global hypotheses

    SMESH::ListOfHypothesis_var hypList = theMesh->GetHypothesisList( aShapeObject );
    for ( int i = 0; i < hypList->length(); i++ )
    {
      SMESH::SMESH_Hypothesis_var aHyp = SMESH::SMESH_Hypothesis::_narrow( hypList[ i ]);
      SALOMEDS::SObject_wrap so = PublishHypothesis( theStudy, aHyp );
      AddHypothesisToShape( theStudy, theMesh, aShapeObject, aHyp );
    }
  }

  // Publish submeshes

  SMESH_Mesh_i* mesh_i = objectToServant<SMESH_Mesh_i>( theMesh );
  if ( !mesh_i )
    return aMeshSO._retn();
  map<int, SMESH_subMesh_i*>& subMap = mesh_i->_mapSubMesh_i;
  map<int, SMESH_subMesh_i*>::iterator subIt = subMap.begin();
  for ( ; subIt != subMap.end(); subIt++ ) {
    SMESH::SMESH_subMesh_ptr aSubMesh = (*subIt).second->_this();
    if ( !CORBA::is_nil( aSubMesh )) {
      aShapeObject = aSubMesh->GetSubShape();
      SALOMEDS::SObject_wrap( PublishSubMesh( theStudy, theMesh, aSubMesh, aShapeObject ));
    }
  }

  // Publish groups
  const map<int, SMESH::SMESH_GroupBase_ptr>& grMap = mesh_i->getGroups();
  map<int, SMESH::SMESH_GroupBase_ptr>::const_iterator it = grMap.begin();
  for ( ; it != grMap.end(); it++ )
  {
    SMESH::SMESH_GroupBase_ptr aGroup = (*it).second;
    if ( !aGroup->_is_nil() ) {
      GEOM::GEOM_Object_var aShapeObj;
      SMESH::SMESH_GroupOnGeom_var aGeomGroup = SMESH::SMESH_GroupOnGeom::_narrow( aGroup );
      if ( !aGeomGroup->_is_nil() )
        aShapeObj = aGeomGroup->GetShape();
      SALOMEDS::SObject_wrap( PublishGroup( theStudy, theMesh, aGroup, aShapeObj ));
    }
  }

  if(MYDEBUG) MESSAGE("PublishMesh_END");
  return aMeshSO._retn();
}

//=======================================================================
//function : PublishSubMesh
//purpose  : 
//=======================================================================

SALOMEDS::SObject_ptr SMESH_Gen_i::PublishSubMesh (SALOMEDS::Study_ptr      theStudy,
                                                   SMESH::SMESH_Mesh_ptr    theMesh,
                                                   SMESH::SMESH_subMesh_ptr theSubMesh,
                                                   GEOM::GEOM_Object_ptr    theShapeObject,
                                                   const char*              theName)
{
  if (theStudy->_is_nil() || theMesh->_is_nil() ||
      theSubMesh->_is_nil() || theShapeObject->_is_nil() )
    return SALOMEDS::SObject::_nil();

  SALOMEDS::SObject_wrap aSubMeshSO = ObjectToSObject( theStudy, theSubMesh );
  if ( aSubMeshSO->_is_nil() )
  {
    SALOMEDS::SObject_wrap aMeshSO = ObjectToSObject( theStudy, theMesh );
    if ( aMeshSO->_is_nil() ) {
      aMeshSO = PublishMesh( theStudy, theMesh );
      if ( aMeshSO->_is_nil())
        return SALOMEDS::SObject::_nil();
    }
    // Find submesh sub-tree tag
    long aRootTag;
    const char* aRootName = "";
    switch ( theShapeObject->GetShapeType() ) {
    case GEOM::VERTEX:
      aRootTag  = GetSubMeshOnVertexTag();
      aRootName = "SubMeshes on Vertex";
      break;
    case GEOM::EDGE:
      aRootTag  = GetSubMeshOnEdgeTag();
      aRootName = "SubMeshes on Edge";
      break;
    case GEOM::WIRE:
      aRootTag  = GetSubMeshOnWireTag();
      aRootName = "SubMeshes on Wire";
      break;
    case GEOM::FACE:
      aRootTag  = GetSubMeshOnFaceTag();
      aRootName = "SubMeshes on Face";    
      break;
    case GEOM::SHELL:
      aRootTag  = GetSubMeshOnShellTag();
      aRootName = "SubMeshes on Shell";   
      break;
    case GEOM::SOLID:
      aRootTag  = GetSubMeshOnSolidTag();
      aRootName = "SubMeshes on Solid";
      break;
    default:
      aRootTag  = GetSubMeshOnCompoundTag();
      aRootName = "SubMeshes on Compound";
      break;
    }

    // Find or create submesh root
    SALOMEDS::SObject_wrap aRootSO = publish (theStudy, CORBA::Object::_nil(),
                                             aMeshSO, aRootTag, 0, false );
    if ( aRootSO->_is_nil() )
      return aSubMeshSO._retn();

    SetName( aRootSO, aRootName );

    // Add new submesh to corresponding sub-tree
    SMESH::array_of_ElementType_var elemTypes = theSubMesh->GetTypes();
    const int isEmpty = ( elemTypes->length() == 0 );
    const char* pm[2] = { "ICON_SMESH_TREE_MESH", "ICON_SMESH_TREE_MESH_WARN" };
    aSubMeshSO = publish (theStudy, theSubMesh, aRootSO, 0, pm[isEmpty] );
    if ( aSubMeshSO->_is_nil() )
      return aSubMeshSO._retn();
  }
  SetName( aSubMeshSO, theName, "SubMesh" );

  // Add reference to theShapeObject

  addReference( theStudy, aSubMeshSO, theShapeObject, 1 );

  // Publish hypothesis

  SMESH::ListOfHypothesis_var hypList = theMesh->GetHypothesisList( theShapeObject );
  for ( int i = 0; i < hypList->length(); i++ ) {
    SMESH::SMESH_Hypothesis_var aHyp = SMESH::SMESH_Hypothesis::_narrow( hypList[ i ]);
    SALOMEDS::SObject_wrap so = PublishHypothesis( theStudy, aHyp );
    AddHypothesisToShape( theStudy, theMesh, theShapeObject, aHyp );
  }

  return aSubMeshSO._retn();
}

//=======================================================================
//function : PublishGroup
//purpose  : 
//=======================================================================

SALOMEDS::SObject_ptr SMESH_Gen_i::PublishGroup (SALOMEDS::Study_ptr    theStudy,
                                                 SMESH::SMESH_Mesh_ptr  theMesh,
                                                 SMESH::SMESH_GroupBase_ptr theGroup,
                                                 GEOM::GEOM_Object_ptr  theShapeObject,
                                                 const char*            theName)
{
  if (theStudy->_is_nil() || theMesh->_is_nil() || theGroup->_is_nil() )
    return SALOMEDS::SObject::_nil();

  SALOMEDS::SObject_wrap aGroupSO = ObjectToSObject( theStudy, theGroup );
  if ( aGroupSO->_is_nil() )
  {
    SALOMEDS::SObject_wrap aMeshSO = ObjectToSObject( theStudy, theMesh );
    if ( aMeshSO->_is_nil() ) {
      aMeshSO = PublishInStudy( theStudy, SALOMEDS::SObject::_nil(), theMesh, "");
      if ( aMeshSO->_is_nil())
        return SALOMEDS::SObject::_nil();
    }
    int aType = (int)theGroup->GetType();
    const char* aRootNames[] = {
      "Compound Groups", "Groups of Nodes", "Groups of Edges",
      "Groups of Faces", "Groups of Volumes", "Groups of 0D Elements",
      "Groups of Balls" };

    // Currently, groups with heterogenous content are not supported
    if ( aType != SMESH::ALL )
    {
      long aRootTag = GetNodeGroupsTag() + aType - 1;

      // Find or create groups root
      SALOMEDS::SObject_wrap aRootSO = publish (theStudy, CORBA::Object::_nil(),
                                               aMeshSO, aRootTag, 0, false );
      if ( aRootSO->_is_nil() ) return SALOMEDS::SObject::_nil();

      if ( aType < sizeof(aRootNames)/sizeof(char*) )
        SetName( aRootSO, aRootNames[aType] );

      // Add new group to corresponding sub-tree
      SMESH::array_of_ElementType_var elemTypes = theGroup->GetTypes();
      int isEmpty = ( elemTypes->length() == 0 );
      std::string pm[2] = { "ICON_SMESH_TREE_GROUP", "ICON_SMESH_TREE_MESH_WARN" };
      if ( SMESH::DownCast< SMESH_GroupOnFilter_i* > ( theGroup ))
      {
        pm[0] = "ICON_SMESH_TREE_GROUP_ON_FILTER";
      }
      else if ( SMESH::DownCast< SMESH_Group_i* > ( theGroup ))
      {
        SMESH::array_of_ElementType_var allElemTypes = theMesh->GetTypes();
        for ( size_t i =0; i < allElemTypes->length() && isEmpty; ++i )
          isEmpty = ( allElemTypes[i] != theGroup->GetType() );
      }
      aGroupSO = publish (theStudy, theGroup, aRootSO, 0, pm[isEmpty].c_str() );
    }
    if ( aGroupSO->_is_nil() )
      return aGroupSO._retn();
  }

  SetName( aGroupSO, theName, "Group" );

  //Add reference to geometry
  if ( !theShapeObject->_is_nil() )
    addReference( theStudy, aGroupSO, theShapeObject, 1 );

  return aGroupSO._retn();
}

//=======================================================================
//function : PublishHypothesis
//purpose  : 
//=======================================================================

SALOMEDS::SObject_ptr
  SMESH_Gen_i::PublishHypothesis (SALOMEDS::Study_ptr         theStudy,
                                  SMESH::SMESH_Hypothesis_ptr theHyp,
                                  const char*                 theName)
{
  if(MYDEBUG) MESSAGE("PublishHypothesis")
  if (theStudy->_is_nil() || theHyp->_is_nil())
    return SALOMEDS::SObject::_nil();

  CORBA::String_var hypType = theHyp->GetName();

  SALOMEDS::SObject_wrap aHypSO = ObjectToSObject( theStudy, theHyp );
  if ( aHypSO->_is_nil() )
  {
    SALOMEDS::SComponent_wrap father = PublishComponent( theStudy );
    if ( father->_is_nil() )
      return aHypSO._retn();

    //Find or Create Hypothesis root
    bool isAlgo = ( !SMESH::SMESH_Algo::_narrow( theHyp )->_is_nil() );
    int aRootTag = isAlgo ? GetAlgorithmsRootTag() : GetHypothesisRootTag();
    SALOMEDS::SObject_wrap aRootSO =
      publish (theStudy, CORBA::Object::_nil(),father, aRootTag,
               isAlgo ? "ICON_SMESH_TREE_ALGO" : "ICON_SMESH_TREE_HYPO", false);
    SetName( aRootSO, isAlgo ?  "Algorithms" : "Hypotheses" );

    // Add New Hypothesis
    string aPmName = isAlgo ? "ICON_SMESH_TREE_ALGO_" : "ICON_SMESH_TREE_HYPO_";
    aPmName += hypType.in();
    // prepend plugin name to pixmap name
    string pluginName = myHypCreatorMap[ hypType.in() ]->GetModuleName();
    if ( pluginName != "StdMeshers" )
      aPmName = pluginName + "::" + aPmName;
    aHypSO = publish( theStudy, theHyp, aRootSO, 0, aPmName.c_str() );
  }

  SetName( aHypSO, theName, hypType.in() );

  if(MYDEBUG) MESSAGE("PublishHypothesis--END")
  return aHypSO._retn();
}

//=======================================================================
//function : GetMeshOrSubmeshByShape
//purpose  : 
//=======================================================================

SALOMEDS::SObject_ptr
  SMESH_Gen_i::GetMeshOrSubmeshByShape (SALOMEDS::Study_ptr   theStudy,
                                        SMESH::SMESH_Mesh_ptr theMesh,
                                        GEOM::GEOM_Object_ptr theShape)
{
  if(MYDEBUG) MESSAGE("GetMeshOrSubmeshByShape")
  SALOMEDS::SObject_wrap aMeshOrSubMesh;
  if (theMesh->_is_nil() || ( theShape->_is_nil() && theMesh->HasShapeToMesh()))
    return aMeshOrSubMesh._retn();
  
  TopoDS_Shape aShape;
  if(theMesh->HasShapeToMesh())
    aShape = GeomObjectToShape( theShape );
  else
    aShape = SMESH_Mesh::PseudoShape();

  SMESH_Mesh_i* mesh_i = objectToServant<SMESH_Mesh_i>( theMesh );

  if ( !aShape.IsNull() && mesh_i && mesh_i->GetImpl().GetMeshDS() ) {
    SMESHDS_Mesh* meshDS = mesh_i->GetImpl().GetMeshDS();
    if ( aShape.IsSame( meshDS->ShapeToMesh() ))
      aMeshOrSubMesh = ObjectToSObject( theStudy, theMesh );
    else {
      int shapeID = meshDS->ShapeToIndex( aShape );
      SMESH::SMESH_subMesh_var aSubMesh = mesh_i->getSubMesh(shapeID);
      if ( !aSubMesh->_is_nil() )
        aMeshOrSubMesh = ObjectToSObject( theStudy, aSubMesh );
    }
  }
  if(MYDEBUG) MESSAGE("GetMeshOrSubmeshByShape--END")
  return aMeshOrSubMesh._retn();
}

//=======================================================================
//function : AddHypothesisToShape
//purpose  : 
//=======================================================================

bool SMESH_Gen_i::AddHypothesisToShape(SALOMEDS::Study_ptr         theStudy,
                                       SMESH::SMESH_Mesh_ptr       theMesh,
                                       GEOM::GEOM_Object_ptr       theShape,
                                       SMESH::SMESH_Hypothesis_ptr theHyp)
{
  if(MYDEBUG) MESSAGE("AddHypothesisToShape")
  if (theStudy->_is_nil() || theMesh->_is_nil() ||
      theHyp->_is_nil() || (theShape->_is_nil()
                            && theMesh->HasShapeToMesh()) )
    return false;

  SALOMEDS::SObject_wrap aMeshSO = ObjectToSObject( theStudy, theMesh );
  if ( aMeshSO->_is_nil() )
    aMeshSO = PublishMesh( theStudy, theMesh );
  SALOMEDS::SObject_wrap aHypSO = PublishHypothesis( theStudy, theHyp );
  if ( aMeshSO->_is_nil() || aHypSO->_is_nil())
    return false;

  // Find a mesh or submesh refering to theShape
  SALOMEDS::SObject_wrap aMeshOrSubMesh =
    GetMeshOrSubmeshByShape( theStudy, theMesh, theShape );
  if ( aMeshOrSubMesh->_is_nil() )
  {
    // publish submesh
    TopoDS_Shape aShape = GeomObjectToShape( theShape );
    SMESH_Mesh_i* mesh_i = objectToServant<SMESH_Mesh_i>( theMesh );
    if ( !aShape.IsNull() && mesh_i && mesh_i->GetImpl().GetMeshDS() ) {
      SMESHDS_Mesh* meshDS = mesh_i->GetImpl().GetMeshDS();
      int shapeID = meshDS->ShapeToIndex( aShape );
      SMESH::SMESH_subMesh_var aSubMesh = mesh_i->getSubMesh(shapeID);
      aMeshOrSubMesh = PublishSubMesh( theStudy, theMesh, aSubMesh, theShape );
    }
    if ( aMeshOrSubMesh->_is_nil() )
      return false;
  }

  //Find or Create Applied Hypothesis root
  bool aIsAlgo = !SMESH::SMESH_Algo::_narrow( theHyp )->_is_nil();
  SALOMEDS::SObject_wrap AHR =
    publish (theStudy, CORBA::Object::_nil(), aMeshOrSubMesh,
             aIsAlgo ? GetRefOnAppliedAlgorithmsTag() : GetRefOnAppliedHypothesisTag(),
             aIsAlgo ? "ICON_SMESH_TREE_ALGO" : "ICON_SMESH_TREE_HYPO", false);
  SetName( AHR, aIsAlgo ? "Applied algorithms" : "Applied hypotheses" );

  addReference( theStudy, AHR, theHyp );

  if(MYDEBUG) MESSAGE("AddHypothesisToShape--END")
  return true;
}

//=======================================================================
//function : RemoveHypothesisFromShape
//purpose  : 
//=======================================================================

bool SMESH_Gen_i::RemoveHypothesisFromShape(SALOMEDS::Study_ptr         theStudy,
                                            SMESH::SMESH_Mesh_ptr       theMesh,
                                            GEOM::GEOM_Object_ptr       theShape,
                                            SMESH::SMESH_Hypothesis_ptr theHyp)
{
  if (theStudy->_is_nil() || theMesh->_is_nil() ||
      theHyp->_is_nil() || (theShape->_is_nil()
                            && theMesh->HasShapeToMesh()))
    return false;

  SALOMEDS::SObject_wrap aHypSO = ObjectToSObject( theStudy, theHyp );
  if ( aHypSO->_is_nil() )
    return false;

  CORBA::String_var hypEntry = aHypSO->GetID();

  // Find a mesh or submesh refering to theShape
  SALOMEDS::SObject_wrap aMeshOrSubMesh =
    GetMeshOrSubmeshByShape( theStudy, theMesh, theShape );
  if ( aMeshOrSubMesh->_is_nil() )
    return false;

  // Find and remove a reference to aHypSO
  SALOMEDS::SObject_wrap aRef, anObj;
  SALOMEDS::ChildIterator_wrap it = theStudy->NewChildIterator( aMeshOrSubMesh );
  bool found = false;
  for ( it->InitEx( true ); ( it->More() && !found ); it->Next() ) {
    anObj = it->Value();
    if (anObj->ReferencedObject( aRef.inout() ))
    {
      CORBA::String_var refEntry = aRef->GetID();
      found = ( strcmp( refEntry, hypEntry ) == 0 );
    }
    if ( found )
    {
      SALOMEDS::StudyBuilder_var builder = theStudy->NewBuilder();
      builder->RemoveObject( anObj );
    }
  }

  return true;
}

//================================================================================
/*!
 * \brief Stores names of variables that WILL be passes as parameters when calling
 *        some method of a given object.
 *  \param [in] theObject - the object whose a method WILL be called with \a theParameters.
 *  \param [in] theParameters - a string contating parameters separated by ':'.
 */
//================================================================================

void SMESH_Gen_i::UpdateParameters(CORBA::Object_ptr theObject, const char* theParameters)
{
  SALOMEDS::Study_var aStudy = GetCurrentStudy();
  if ( aStudy->_is_nil() )
    return;

  // find variable names within theParameters

  myLastObj.clear();
  myLastParameters.clear();
  myLastParamIndex.clear(); /* vector holding indices of virables within the string
                               of all varibles used for theObject */ 
  int nbVars = 0;
  int pos = 0, prevPos = 0, len = strlen( theParameters );
  if ( len == 0 ) return;
  while ( pos <= len )
  {
    if ( pos == len || theParameters[pos] == ':' )
    {
      if ( prevPos < pos )
      {
        string val( theParameters + prevPos, theParameters + pos );
        if ( !aStudy->IsVariable( val.c_str() ))
          val.clear();
        myLastParameters.push_back( val );
        nbVars += (! myLastParameters.back().empty() );
      }
      else
      {
        myLastParameters.push_back("");
      }
      prevPos = pos+1;
    }
    ++pos;
  }

  if ( nbVars < 1 )
    return;

  // store
  // (1) variable names in the string of all varibles used for theObject and
  // (2) indices of found variables in myLastParamIndex.

  // remember theObject
  SALOMEDS::SObject_wrap aSObj =  ObjectToSObject(aStudy,theObject);
  if ( aSObj->_is_nil() )
    return;
  CORBA::String_var anObjEntry = aSObj->GetID();
  myLastObj = anObjEntry.in();

  // get a string of variable names
  SALOMEDS::StudyBuilder_var   aStudyBuilder = aStudy->NewBuilder();
  SALOMEDS::GenericAttribute_wrap     anAttr =
    aStudyBuilder->FindOrCreateAttribute( aSObj, "AttributeString" );
  SALOMEDS::AttributeString_wrap aStringAttr = anAttr;
  CORBA::String_var                  oldVars = aStringAttr->Value();
  std::string                         varStr = oldVars.in();

  // add new variables and find indices of variables
  for ( size_t i = 0; i < myLastParameters.size(); ++i )
  {
    int varIndex = -1;
    if ( !myLastParameters[i].empty() )
    {
      // find index of myLastParameters[i] in varStr
      int curIndex  = 0;
      bool varFound = false;
      size_t pos    = 0;
      // varStr can be "A|B::C;*=2|D"
      const std::string separators(":|;*=");
      while ( pos < varStr.size() && !varFound )
      {
        // skip separators
        while ( separators.find( varStr[ pos ]) != std::string::npos )
          if ( ++pos >= varStr.size() )
            break;
        // skip repetition number following '='
        if ( varStr[ pos-1 ] == '=' )
        {
          while ( '0' <= varStr[ pos ] && varStr[ pos ] <= '9' )
            ++pos;
          continue; // to skip next separator
        }
        // compare variable name
        if ( pos < varStr.size() )
        {
          varFound = ( varStr.compare( pos, myLastParameters[i].size(), myLastParameters[i] ) == 0 &&
                       // same string begining but is length same?
                       ( pos + myLastParameters[i].size() >= varStr.size() ||
                         separators.find( varStr[ pos+1 ]) != std::string::npos ));
          if ( varFound )
            varIndex = curIndex;
          else
            pos = varStr.find_first_of( separators, pos ); // goto the next separator
          ++curIndex;
        }
      }
      // add new variable
      if ( !varFound )
      {
        varStr += ":" + myLastParameters[i];
        varIndex = curIndex;
      }
    }
    myLastParamIndex.push_back( varIndex );
  }
  aStringAttr->SetValue( varStr.c_str() );
}

//================================================================================
/*!
 * \brief Return all variables used to create an object
 *  \param [in] theObjectEntry - an object entry in the current study
 *  \return std::vector< std::string > - all variable names (or values of removed variables)
 */
//================================================================================

std::vector< std::string > SMESH_Gen_i::GetAllParameters(const std::string& theObjectEntry) const
{
  std::vector< std::string > varNames;
  if ( myCurrentStudy->_is_nil() )
    return varNames;

  SALOMEDS::SObject_wrap aSObj = myCurrentStudy->FindObjectID( theObjectEntry.c_str() );
  if ( myCurrentStudy->_is_nil() )
    return varNames;

  // get a string of variable names
  SALOMEDS::StudyBuilder_var   aStudyBuilder = myCurrentStudy->NewBuilder();
  SALOMEDS::GenericAttribute_wrap     anAttr =
    aStudyBuilder->FindOrCreateAttribute( aSObj, "AttributeString" );
  SALOMEDS::AttributeString_wrap aStringAttr = anAttr;
  CORBA::String_var                  oldVars = aStringAttr->Value();
  std::string                         varStr = oldVars.in();

  // separate variables within varStr;
  // varStr can be "A|B::C;*=2|D"
  size_t pos = 0;
  const std::string separators(":|;*=");
  while ( pos < varStr.size() )
  {
    // skip separators
    pos = varStr.find_first_not_of( separators, pos );
    // while ( separators.find( varStr[ pos ]) != std::string::npos )
    //   if ( ++pos >= varStr.size() )
    //     break;
    // skip repetition number following '='
    if ( varStr[ pos-1 ] == '=' )
    {
      while ( '0' <= varStr[ pos ] && varStr[ pos ] <= '9' )
        ++pos;
      continue; // to skip next separator
    }
    // store variable name
    if ( pos < varStr.size() )
    {
      size_t pos2 = varStr.find_first_of( separators, pos );
      varNames.push_back( varStr.substr( pos, pos2 - pos));
      pos = pos2;
    }
  }
  return varNames;
}

//=======================================================================
//function : ParseParameters
//purpose  : Replace variables by their values
//=======================================================================
// char* SMESH_Gen_i::ParseParameters(const char* theParameters)
// {
//   //const char* aParameters = theParameters;
// //   const char* aParameters = CORBA::string_dup(theParameters);
//   TCollection_AsciiString anInputParams;
//   SALOMEDS::Study_var aStudy = GetCurrentStudy();
//   if( !aStudy->_is_nil() ) {
// //     SALOMEDS::ListOfListOfStrings_var aSections = aStudy->ParseVariables(theParameters);
// //     for(int j=0;j<aSections->length();j++) {
// //       SALOMEDS::ListOfStrings aVars= aSections[j];
// //       for(int i=0;i<aVars.length();i++ ) {
// //         anInputParams += aStudy->IsVariable(aVars[i].in()) ? 
// //           TCollection_AsciiString(aVars[i].in()) : TCollection_AsciiString("");
// //         if(i != aVars.length()-1)
// //           anInputParams+=":";
// //       }
// //       if(j!=aSections->length()-1)
// //         anInputParams+="|";
// //     }
//     TCollection_AsciiString paramStr( theParameters );
//     int beg = 0, end;
//     char sep, *pParams = (char*)paramStr.ToCString();
//     while ( beg < paramStr.Length() )
//     {
//       end = beg-1;
//       while ( ++end < paramStr.Length() )
//         if ( pParams[end] == ':' || pParams[end] == '|')
//           break;
//       if ( end < paramStr.Length())
//       {
//         sep = pParams[end];
//         pParams[end] = '\0';
//       }
//       if ( aStudy->IsVariable( pParams+beg ))
//         anInputParams += pParams+beg;
//       if ( end < paramStr.Length() )
//         anInputParams += sep;
//       else
//         break;
//       beg = end + 1;
//     }
//   }
//   return CORBA::string_dup(anInputParams.ToCString());
// }

//=======================================================================
//function : GetParameters
//purpose  : 
//=======================================================================

char* SMESH_Gen_i::GetParameters(CORBA::Object_ptr theObject)
{
  CORBA::String_var aResult("");

  SALOMEDS::SObject_wrap aSObj = ObjectToSObject( myCurrentStudy, theObject );
  if ( !aSObj->_is_nil() )
  {
    SALOMEDS::GenericAttribute_wrap attr;
    if ( aSObj->FindAttribute( attr.inout(), "AttributeString"))
    {
      SALOMEDS::AttributeString_wrap strAttr = attr;
      aResult = strAttr->Value();
    }
  }
  return aResult._retn();
}