Salome HOME
23239: [CEA 1739] Regression : crash trying to create mesh
[modules/smesh.git] / src / SMESH / SMESH_Mesh.hxx
index 2f77fa2e5050b9ea7e7755189fa6703bc95016b0..fc00be6a77c0a7070148737f9b5ec05d4bee565f 100644 (file)
-//  SMESH SMESH : implementaion of SMESH idl descriptions
+// Copyright (C) 2007-2015  CEA/DEN, EDF R&D, OPEN CASCADE
 //
-//  Copyright (C) 2003  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.opencascade.org/SALOME/ or email : webmaster.salome@opencascade.org 
+// 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, or (at your option) any later version.
 //
+// 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_Mesh.hxx
 //  Author : Paul RASCLE, EDF
 //  Module : SMESH
-//  $Header$
-
+//
 #ifndef _SMESH_MESH_HXX_
 #define _SMESH_MESH_HXX_
 
-#include "SMESHDS_Document.hxx"
-#include "SMESHDS_Mesh.hxx"
+#include "SMESH_SMESH.hxx"
+
+#include "SMDSAbs_ElementType.hxx"
 #include "SMESHDS_Command.hxx"
+#include "SMESHDS_Mesh.hxx"
+#include "SMESH_ComputeError.hxx"
+#include "SMESH_Controls.hxx"
 #include "SMESH_Hypothesis.hxx"
-#include "SMESH_subMesh.hxx"
-#include <SMDSAbs_ElementType.hxx>
+
 #include "Utils_SALOME_Exception.hxx"
 
-#include <TopExp.hxx>
-#include <TopExp_Explorer.hxx>
-#include <TopoDS.hxx>
-#include <TopoDS_Iterator.hxx>
-#include <TopoDS_Compound.hxx>
-#include <TopoDS_CompSolid.hxx>
-#include <TopoDS_Solid.hxx>
-#include <TopoDS_Shell.hxx>
-#include <TopoDS_Face.hxx>
-#include <TopoDS_Wire.hxx>
-#include <TopoDS_Edge.hxx>
-#include <TopoDS_Vertex.hxx>
 #include <TopoDS_Shape.hxx>
-#include <TopTools_IndexedMapOfShape.hxx>
 #include <TopTools_IndexedDataMapOfShapeListOfShape.hxx>
 
-#include <vector>
-#include <list>
 #include <map>
+#include <list>
 
-#if (__GNUC__>2)
-#include <string>
-#include <istream>
-#include <ostream>
-#else 
-#include <string>
-#include <istream.h>
-#include <ostream.h>
+
+#ifdef WIN32
+#pragma warning(disable:4251) // Warning DLL Interface ...
+#pragma warning(disable:4290) // Warning Exception ...
 #endif
 
 class SMESH_Gen;
+class SMESHDS_Document;
 class SMESH_Group;
 class TopTools_ListOfShape;
+class SMESH_subMesh;
+class SMESH_HypoFilter;
+class TopoDS_Solid;
 
-class SMESH_Mesh
-{
-  public:
-       SMESH_Mesh();
-       SMESH_Mesh(int localId, int studyId, SMESH_Gen * gen,
-               SMESHDS_Document * myDocument);
-
-       virtual ~SMESH_Mesh();
-
-        void ShapeToMesh(const TopoDS_Shape & aShape);
-        /*!
-         * consult DriverMED_R_SMESHDS_Mesh::ReadStatus for returned value
-         */
-        int MEDToMesh(const char* theFileName, const char* theMeshName);
-
-       SMESH_Hypothesis::Hypothesis_Status
-                AddHypothesis(const TopoDS_Shape & aSubShape, int anHypId)
-               throw(SALOME_Exception);
-
-       SMESH_Hypothesis::Hypothesis_Status
-                RemoveHypothesis(const TopoDS_Shape & aSubShape, int anHypId)
-               throw(SALOME_Exception);
+typedef std::list<int> TListOfInt;
+typedef std::list<TListOfInt> TListOfListOfInt;
 
-       const list <const SMESHDS_Hypothesis * >&
-               GetHypothesisList(const TopoDS_Shape & aSubShape) const
-               throw(SALOME_Exception);
-
-       const list<SMESHDS_Command*> & GetLog() throw(SALOME_Exception);
-
-//   const SMESHDS_ListOfAsciiString& GetLog()
-//     throw (SALOME_Exception);
-
-       void ClearLog() throw(SALOME_Exception);
+class SMESH_EXPORT SMESH_Mesh
+{
+ public:
+  SMESH_Mesh(int               theLocalId,
+             int               theStudyId,
+             SMESH_Gen*        theGen,
+             bool              theIsEmbeddedMode,
+             SMESHDS_Document* theDocument);
+
+  virtual ~SMESH_Mesh();
+
+  /*!
+   * \brief Set geometry to be meshed
+   */
+  void ShapeToMesh(const TopoDS_Shape & aShape);
+  /*!
+   * \brief Return geometry to be meshed. (It may be a PseudoShape()!)
+   */
+  TopoDS_Shape GetShapeToMesh() const;
+  /*!
+   * \brief Return true if there is a geometry to be meshed, not PseudoShape()
+   */
+  bool HasShapeToMesh() const { return _isShapeToMesh; }
+  /*!
+   * \brief Return diagonal size of bounding box of shape to mesh.
+   */
+  double GetShapeDiagonalSize() const;
+  /*!
+   * \brief Return diagonal size of bounding box of a shape.
+   */
+  static double GetShapeDiagonalSize(const TopoDS_Shape & aShape);
+  /*!
+   * \brief Return a solid which is returned by GetShapeToMesh() if
+   *        a real geometry to be meshed was not set
+   */
+  static const TopoDS_Solid& PseudoShape();
+
+  /*!
+   * \brief Load mesh from study file
+   */
+  void Load();
+  /*!
+   * \brief Remove all nodes and elements
+   */
+  void Clear();
+  /*!
+   * \brief Remove all nodes and elements of indicated shape
+   */
+  void ClearSubMesh(const int theShapeId);
+
+  /*!
+   * consult DriverMED_R_SMESHDS_Mesh::ReadStatus for returned value
+   */
+  int UNVToMesh(const char* theFileName);
+
+  int MEDToMesh(const char* theFileName, const char* theMeshName);
+  
+  int STLToMesh(const char* theFileName);
 
-       int GetId();
+  int CGNSToMesh(const char* theFileName, const int theMeshIndex, std::string& theMeshName);
+  
+  SMESH_ComputeErrorPtr GMFToMesh(const char* theFileName,
+                                  bool        theMakeRequiredGroups = true );
 
-       SMESHDS_Mesh * GetMeshDS();
+  SMESH_Hypothesis::Hypothesis_Status
+  AddHypothesis(const TopoDS_Shape & aSubShape, int anHypId, std::string* error=0)
+    throw(SALOME_Exception);
+  
+  SMESH_Hypothesis::Hypothesis_Status
+  RemoveHypothesis(const TopoDS_Shape & aSubShape, int anHypId)
+    throw(SALOME_Exception);
+  
+  const std::list <const SMESHDS_Hypothesis * >&
+  GetHypothesisList(const TopoDS_Shape & aSubShape) const
+    throw(SALOME_Exception);
+
+  const SMESH_Hypothesis * GetHypothesis(const TopoDS_Shape &    aSubShape,
+                                         const SMESH_HypoFilter& aFilter,
+                                         const bool              andAncestors,
+                                         TopoDS_Shape*           assignedTo=0) const;
+  
+  int GetHypotheses(const TopoDS_Shape &                     aSubShape,
+                    const SMESH_HypoFilter&                  aFilter,
+                    std::list< const SMESHDS_Hypothesis * >& aHypList,
+                    const bool                               andAncestors,
+                    std::list< TopoDS_Shape > *              assignedTo=0) const;
+
+  const SMESH_Hypothesis * GetHypothesis(const SMESH_subMesh *   aSubMesh,
+                                         const SMESH_HypoFilter& aFilter,
+                                         const bool              andAncestors,
+                                         TopoDS_Shape*           assignedTo=0) const;
+  
+  int GetHypotheses(const SMESH_subMesh *                    aSubMesh,
+                    const SMESH_HypoFilter&                  aFilter,
+                    std::list< const SMESHDS_Hypothesis * >& aHypList,
+                    const bool                               andAncestors,
+                    std::list< TopoDS_Shape > *              assignedTo=0) const;
 
-       SMESH_Gen *GetGen();
+  SMESH_Hypothesis * GetHypothesis(const int aHypID) const;
 
-       SMESH_subMesh *GetSubMesh(const TopoDS_Shape & aSubShape)
-               throw(SALOME_Exception);
+  const std::list<SMESHDS_Command*> & GetLog() throw(SALOME_Exception);
+  
+  void ClearLog() throw(SALOME_Exception);
+  
+  int GetId() const          { return _id; }
+  
+  bool MeshExists( int meshId ) const;
+  
+  SMESH_Mesh* FindMesh( int meshId ) const;
 
-       SMESH_subMesh *GetSubMeshContaining(const TopoDS_Shape & aSubShape)
-               throw(SALOME_Exception);
+  SMESHDS_Mesh * GetMeshDS() { return _myMeshDS; }
 
-       const list < SMESH_subMesh * >&
-               GetSubMeshUsingHypothesis(SMESHDS_Hypothesis * anHyp)
-               throw(SALOME_Exception);
+  const SMESHDS_Mesh * GetMeshDS() const { return _myMeshDS; }
+  
+  SMESH_Gen *GetGen()        { return _gen; }
 
-        bool IsUsedHypothesis(SMESHDS_Hypothesis * anHyp,
-                              const TopoDS_Shape & aSubShape);
-        // Return True if anHyp is used to mesh aSubShape
+  SMESH_subMesh *GetSubMesh(const TopoDS_Shape & aSubShape)
+    throw(SALOME_Exception);
+  
+  SMESH_subMesh *GetSubMeshContaining(const TopoDS_Shape & aSubShape) const
+    throw(SALOME_Exception);
+  
+  SMESH_subMesh *GetSubMeshContaining(const int aShapeID) const
+    throw(SALOME_Exception);
+  /*!
+   * \brief Return submeshes of groups containing the given subshape
+   */
+  std::list<SMESH_subMesh*> GetGroupSubMeshesContaining(const TopoDS_Shape & shape) const
+    throw(SALOME_Exception);
+  /*!
+   * \brief Say all submeshes that theChangedHyp has been modified
+   */
+  void NotifySubMeshesHypothesisModification(const SMESH_Hypothesis* theChangedHyp);
+
+  // const std::list < SMESH_subMesh * >&
+  // GetSubMeshUsingHypothesis(SMESHDS_Hypothesis * anHyp) throw(SALOME_Exception);
+  /*!
+   * \brief Return True if anHyp is used to mesh aSubShape
+   */
+  bool IsUsedHypothesis(SMESHDS_Hypothesis *  anHyp,
+                        const SMESH_subMesh * aSubMesh);
+  /*!
+   * \brief check if a hypothesis alowing notconform mesh is present
+   */
+  bool IsNotConformAllowed() const;
+  
+  bool IsMainShape(const TopoDS_Shape& theShape) const;
+  /*!
+   * \brief Return list of ancestors of theSubShape in the order
+   *        that lower dimention shapes come first
+   */
+  const TopTools_ListOfShape& GetAncestors(const TopoDS_Shape& theSubShape) const;
+
+  void SetAutoColor(bool theAutoColor) throw(SALOME_Exception);
+
+  bool GetAutoColor() throw(SALOME_Exception);
+
+  /*!
+   * \brief Set the flag meaning that the mesh has been edited "manually".
+   * It is to set to false after Clear() and to set to true by MeshEditor
+   */
+  void SetIsModified(bool isModified);
+
+  bool GetIsModified() const { return _isModified; }
+
+  /*!
+   * \brief Return true if the mesh has been edited since a total re-compute
+   *        and those modifications may prevent successful partial re-compute.
+   *        As a side effect reset _isModified flag if mesh is empty
+   */
+  bool HasModificationsToDiscard() const;
+
+  /*!
+   * \brief Return data map of descendant to ancestor shapes
+   */
+  typedef TopTools_IndexedDataMapOfShapeListOfShape TAncestorMap;
+  const TAncestorMap& GetAncestorMap() const { return _mapAncestors; }
+  /*!
+   * \brief Check group names for duplications.
+   *  Consider maximum group name length stored in MED file
+   */
+  bool HasDuplicatedGroupNamesMED();
+
+  void ExportMED(const char *        theFile, 
+                 const char*         theMeshName = NULL, 
+                 bool                theAutoGroups = true, 
+                 int                 theVersion = 0,
+                 const SMESHDS_Mesh* theMeshPart = 0,
+                 bool                theAutoDimension = false,
+                 bool                theAddODOnVertices = false)
+    throw(SALOME_Exception);
+
+  void ExportDAT(const char *        file,
+                 const SMESHDS_Mesh* meshPart = 0) throw(SALOME_Exception);
+  void ExportUNV(const char *        file,
+                 const SMESHDS_Mesh* meshPart = 0) throw(SALOME_Exception);
+  void ExportSTL(const char *        file,
+                 const bool          isascii,
+                 const SMESHDS_Mesh* meshPart = 0) throw(SALOME_Exception);
+  void ExportCGNS(const char *        file,
+                  const SMESHDS_Mesh* mesh,
+                  const char *        meshName = 0);
+  void ExportGMF(const char *        file,
+                 const SMESHDS_Mesh* mesh,
+                 bool                withRequiredGroups = true );
+  void ExportSAUV(const char *file, 
+                  const char* theMeshName = NULL, 
+                  bool theAutoGroups = true) throw(SALOME_Exception);
+
+  double GetComputeProgress() const;
+  
+  int NbNodes() const throw(SALOME_Exception);
+  int Nb0DElements() const throw(SALOME_Exception);
+  int NbBalls() const throw(SALOME_Exception);
+  
+  int NbEdges(SMDSAbs_ElementOrder order = ORDER_ANY) const throw(SALOME_Exception);
+  
+  int NbFaces(SMDSAbs_ElementOrder order = ORDER_ANY) const throw(SALOME_Exception);
+  int NbTriangles(SMDSAbs_ElementOrder order = ORDER_ANY) const throw(SALOME_Exception);
+  int NbQuadrangles(SMDSAbs_ElementOrder order = ORDER_ANY) const throw(SALOME_Exception);
+  int NbBiQuadQuadrangles() const throw(SALOME_Exception);
+  int NbBiQuadTriangles() const throw(SALOME_Exception);
+  int NbPolygons(SMDSAbs_ElementOrder order = ORDER_ANY) const throw(SALOME_Exception);
+  
+  int NbVolumes(SMDSAbs_ElementOrder order = ORDER_ANY) const throw(SALOME_Exception);
+  int NbTetras(SMDSAbs_ElementOrder order = ORDER_ANY) const throw(SALOME_Exception);
+  int NbHexas(SMDSAbs_ElementOrder order = ORDER_ANY) const throw(SALOME_Exception);
+  int NbTriQuadraticHexas() const throw(SALOME_Exception);
+  int NbPyramids(SMDSAbs_ElementOrder order = ORDER_ANY) const throw(SALOME_Exception);
+  int NbPrisms(SMDSAbs_ElementOrder order = ORDER_ANY) const throw(SALOME_Exception);
+  int NbHexagonalPrisms() const throw(SALOME_Exception);
+  int NbPolyhedrons() const throw(SALOME_Exception);
+  
+  int NbSubMesh() const throw(SALOME_Exception);
+  
+  int NbGroup() const { return _mapGroup.size(); }
 
-        bool IsNotConformAllowed() const;
-        // check if a hypothesis alowing notconform mesh is present
+  int NbMeshes() const; // nb meshes in the Study
 
-        bool IsMainShape(const TopoDS_Shape& theShape) const;
+  SMESH_Group* AddGroup (const SMDSAbs_ElementType theType,
+                         const char*               theName,
+                         int&                      theId,
+                         const TopoDS_Shape&       theShape=TopoDS_Shape(),
+                         const SMESH_PredicatePtr& thePredicate=SMESH_PredicatePtr());
 
-        const TopTools_ListOfShape& GetAncestors(const TopoDS_Shape& theSubShape);
-        // return list of ancestors of theSubShape in the order
-        // that lower dimention shapes come first.
-        
-       void ExportDAT(const char *file) throw(SALOME_Exception);
-       void ExportMED(const char *file, const char* theMeshName = NULL, bool theAutoGroups = true) throw(SALOME_Exception);
-       void ExportUNV(const char *file) throw(SALOME_Exception);
+  SMESH_Group* AddGroup (SMESHDS_GroupBase* groupDS) throw(SALOME_Exception);
 
-       int NbNodes() throw(SALOME_Exception);
+  typedef boost::shared_ptr< SMDS_Iterator<SMESH_Group*> > GroupIteratorPtr;
+  GroupIteratorPtr GetGroups() const;
+  
+  std::list<int> GetGroupIds() const;
+  
+  SMESH_Group* GetGroup (const int theGroupID);
 
-       int NbEdges() throw(SALOME_Exception);
+  bool RemoveGroup (const int theGroupID);
 
-       int NbFaces() throw(SALOME_Exception);
+  SMESH_Group* ConvertToStandalone ( int theGroupID );
 
-       int NbTriangles() throw(SALOME_Exception);
+  struct TCallUp // callback from SMESH to SMESH_I level
+  {
+    virtual void RemoveGroup (const int theGroupID)=0;
+    virtual void HypothesisModified ()=0;
+    virtual void Load ()=0;
+    virtual ~TCallUp() {}
+  };
+  void SetCallUp( TCallUp * upCaller );
 
-       int NbQuadrangles() throw(SALOME_Exception);
+  bool SynchronizeGroups();
 
-       int NbVolumes() throw(SALOME_Exception);
 
-       int NbTetras() throw(SALOME_Exception);
+  SMDSAbs_ElementType GetElementType( const int id, const bool iselem );
 
-       int NbHexas() throw(SALOME_Exception);
+  void ClearMeshOrder();
+  void SetMeshOrder(const TListOfListOfInt& theOrder );
+  const TListOfListOfInt& GetMeshOrder() const;
 
-       int NbPyramids() throw(SALOME_Exception);
+  // sort submeshes according to stored mesh order
+  bool SortByMeshOrder(std::vector<SMESH_subMesh*>& theListToSort) const;
 
-       int NbPrisms() throw(SALOME_Exception);
+  // return true if given order of sub-meshes is OK
+  bool IsOrderOK( const SMESH_subMesh* smBefore,
+                  const SMESH_subMesh* smAfter ) const;
 
-       int NbSubMesh() throw(SALOME_Exception);
+  std::ostream& Dump(ostream & save);
+  
+private:
 
-        int NbGroup() const { return _mapGroup.size(); }
+  void fillAncestorsMap(const TopoDS_Shape& theShape);
+  void getAncestorsSubMeshes(const TopoDS_Shape&            theSubShape,
+                             std::vector< SMESH_subMesh* >& theSubMeshes) const;
   
-        SMESH_Group* AddGroup (const SMDSAbs_ElementType theType,
-                               const char*               theName,
-                               int&                      theId);
+protected:
+  int                        _id;           // id given by creator (unique within the creator instance)
+  int                        _studyId;
+  int                        _groupId;      // id generator for group objects
+  int                        _nbSubShapes;  // initial nb of subshapes in the shape to mesh
+  bool                       _isShapeToMesh;// set to true when a shape is given (only once)
+  SMESHDS_Document *         _myDocument;
+  SMESHDS_Mesh *             _myMeshDS;
+  SMESH_Gen *                _gen;
+  std::map <int, SMESH_Group*> _mapGroup;
+
+  class SubMeshHolder;
+  SubMeshHolder*             _subMeshHolder;
   
-        SMESH_Group* GetGroup (const int theGroupID);
+  bool                       _isAutoColor;
+  bool                       _isModified; //!< modified since last total re-compute, issue 0020693
 
-        list<int> GetGroupIds();
+  double                     _shapeDiagonal; //!< diagonal size of bounding box of shape to mesh
   
-        void RemoveGroup (const int theGroupID);
+  TopTools_IndexedDataMapOfShapeListOfShape _mapAncestors;
 
-        ostream& Dump(ostream & save);
+  mutable std::vector<SMESH_subMesh*> _ancestorSubMeshes; // to speed up GetHypothes[ei]s()
 
-  private:
+  TListOfListOfInt           _mySubMeshOrder;
 
-       int _id;                                        // id given by creator (unique within the creator instance)
-       int _studyId;
-       int _idDoc;                                     // id given by SMESHDS_Document
-        int _groupId;                                   // id generator for group objects
-       bool _isShapeToMesh;                            // set to true when a shape is given (only once)
-       list<const SMESHDS_Hypothesis *> _subShapeHypothesisList;
-       list <SMESH_subMesh *> _subMeshesUsingHypothesisList;
-       SMESHDS_Document * _myDocument;
-       SMESHDS_Mesh * _myMeshDS;
-       TopTools_IndexedMapOfShape _subShapes;
-       map <int, SMESH_subMesh *> _mapSubMesh;
-       map <int, SMESH_Group *>   _mapGroup;
-       SMESH_Gen *_gen;
+  // Struct calling methods at CORBA API implementation level, used to
+  // 1) make an upper level (SMESH_I) be consistent with a lower one (SMESH)
+  // when group removal is invoked by hyp modification (issue 0020918)
+  // 2) to forget not loaded mesh data at hyp modification
+  TCallUp*                    _callUp;
 
-        TopTools_IndexedDataMapOfShapeListOfShape _mapAncestors;
+protected:
+  SMESH_Mesh();
+  SMESH_Mesh(const SMESH_Mesh&) {};
 };
 
 #endif