X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FSMDS%2FSMDS_MeshElement.cxx;h=48221c7c6861dda72d67072bf44434508399b404;hb=f513b3c72c0c9c1a0c68bae34813cb41ac6fdcd5;hp=8d4d9b75dc2e0bb913603c73ac5a5fefb72f9887;hpb=090aff07266d376ae028ae43434bdea7c0a0f9bb;p=modules%2Fsmesh.git diff --git a/src/SMDS/SMDS_MeshElement.cxx b/src/SMDS/SMDS_MeshElement.cxx index 8d4d9b75d..48221c7c6 100644 --- a/src/SMDS/SMDS_MeshElement.cxx +++ b/src/SMDS/SMDS_MeshElement.cxx @@ -1,194 +1,216 @@ -// SMESH SMDS : implementaion of Salome mesh data structure +// Copyright (C) 2007-2016 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, 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 // -// 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 +// SMESH SMDS : implementation of Salome mesh data structure +// #ifdef _MSC_VER #pragma warning(disable:4786) #endif #include "SMDS_MeshElement.hxx" -#include "SMDS_MeshNode.hxx" -#include "SMDS_MeshEdge.hxx" -#include "SMDS_MeshFace.hxx" -#include "SMDS_MeshVolume.hxx" + +#include "SMDS_Mesh.hxx" +#include "SMDS_ElementFactory.hxx" + #include "utilities.h" -using namespace std; +//================================================================================ +/*! + * \brief Constructor of a non-used element + */ +//================================================================================ -SMDS_MeshElement::SMDS_MeshElement(int ID):myID(ID) +SMDS_MeshElement::SMDS_MeshElement(): myHolder(0) { } -void SMDS_MeshElement::Print(ostream & OS) const +//================================================================================ +/*! + * \brief Check if a node is a medium node of a quadratic cell + */ +//================================================================================ + +bool SMDS_MeshElement::IsMediumNode(const SMDS_MeshNode* node) const { - OS << "dump of mesh element" << endl; + return !( GetNodeIndex( node ) < NbCornerNodes() ); } -ostream & operator <<(ostream & OS, const SMDS_MeshElement * ME) +//================================================================================ +/*! + * \brief Return true if index of node is valid (0 <= ind < NbNodes()) + * \param ind - node index + * \retval bool - index check result + */ +//================================================================================ + +bool SMDS_MeshElement::IsValidIndex(const int ind) const { - ME->Print(OS); - return OS; + return ( ind>-1 && indmore()) - { - it->next(); - nbnodes++; - } - return nbnodes; -} - -/////////////////////////////////////////////////////////////////////////////// -///Return the number of edges owned by or linked with the current element -/////////////////////////////////////////////////////////////////////////////// -int SMDS_MeshElement::NbEdges() const -{ - int nbedges=0; - SMDS_ElemIteratorPtr it=edgesIterator(); - while(it->more()) - { - it->next(); - nbedges++; - } - return nbedges; -} - -/////////////////////////////////////////////////////////////////////////////// -///Return the number of faces owned by or linked with the current element -/////////////////////////////////////////////////////////////////////////////// -int SMDS_MeshElement::NbFaces() const -{ - int nbfaces=0; - SMDS_ElemIteratorPtr it=facesIterator(); - while(it->more()) - { - it->next(); - nbfaces++; - } - return nbfaces; -} - -/////////////////////////////////////////////////////////////////////////////// -///Create an iterator which iterate on elements linked with the current element. -///@param type The of elements on which you want to iterate -///@return A smart pointer to iterator, you are not to take care of freeing memory -/////////////////////////////////////////////////////////////////////////////// -class SMDS_MeshElement_MyIterator:public SMDS_ElemIterator -{ - const SMDS_MeshElement * myElement; - bool myMore; - public: - SMDS_MeshElement_MyIterator(const SMDS_MeshElement * element): - myElement(element),myMore(true) {} - - bool more() - { - return myMore; - } - - const SMDS_MeshElement* next() - { - myMore=false; - return myElement; - } -}; -SMDS_ElemIteratorPtr SMDS_MeshElement:: - elementsIterator(SMDSAbs_ElementType type) const -{ - /** @todo Check that iterator in the child classes return elements - in the same order for each different implementation (i.e: SMDS_VolumeOfNodes - and SMDS_VolumeOfFaces */ - - if(type==GetType()) - return SMDS_ElemIteratorPtr(new SMDS_MeshElement_MyIterator(this)); - else - { - MESSAGE("Iterator not implemented"); - return SMDS_ElemIteratorPtr((SMDS_ElemIterator*)NULL); - } -} - -/////////////////////////////////////////////////////////////////////////////// -///Return the ID of the element -/////////////////////////////////////////////////////////////////////////////// +//================================================================================ +/*! + * \brief Return a valid corner node index, fixing the given one if necessary + * \param ind - node index + * \retval int - valid node index + */ +//================================================================================ + +int SMDS_MeshElement::WrappedIndex(const int ind) const +{ + if ( ind < 0 ) return NbCornerNodes() + ind % NbCornerNodes(); + if ( ind >= NbCornerNodes() ) return ind % NbCornerNodes(); + return ind; +} + +//================================================================================ +/*! + * \brief Check if a node belongs to the element + * \param node - the node to check + * \retval int - node index within the element, -1 if not found + */ +//================================================================================ + +int SMDS_MeshElement::GetNodeIndex( const SMDS_MeshNode* node ) const +{ + SMDS_ElemIteratorPtr nIt = nodesIterator(); + for ( int i = 0; nIt->more(); ++i ) + if ( nIt->next() == node ) + return i; + return -1; +} + +//================================================================================ +/*! + * \brief Return ID of an element + */ +//================================================================================ + int SMDS_MeshElement::GetID() const { - return myID; + return myHolder ? myHolder->GetID( this ) : -1; } -bool operator<(const SMDS_MeshElement& e1, const SMDS_MeshElement& e2) +//================================================================================ +/*! + * \brief Set ID of a shape this element was generated on + */ +//================================================================================ + +void SMDS_MeshElement::setShapeID( const int shapeID ) const { - if(e1.GetType()!=e2.GetType()) return false; - switch(e1.GetType()) - { - case SMDSAbs_Node: - return static_cast(e1) < - static_cast(e2); + const_cast( myHolder )->SetShapeID( this, shapeID ); +} + +//================================================================================ +/*! + * \brief Return ID of a shape this element was generated on + */ +//================================================================================ - case SMDSAbs_Edge: - return static_cast(e1) < - static_cast(e2); +int SMDS_MeshElement::GetShapeID() const +{ + return myHolder->GetShapeID( this ); +} - case SMDSAbs_Face: - return static_cast(e1) < - static_cast(e2); +//================================================================================ +/*! + * \brief Return VTK ID of this element + */ +//================================================================================ - case SMDSAbs_Volume: - return static_cast(e1) < - static_cast(e2); +int SMDS_MeshElement::GetVtkID() const +{ + return myHolder->GetVtkID( this ); +} - default : MESSAGE("Internal Error"); - } - return false; +//================================================================================ +/*! + * \brief Mark this element + */ +//================================================================================ + +void SMDS_MeshElement::setIsMarked( bool is ) const +{ + const_cast( myHolder )->SetIsMarked( this, is ); +} + +//================================================================================ +/*! + * \brief Check if this element is marked + */ +//================================================================================ + +bool SMDS_MeshElement::isMarked() const +{ + return myHolder->IsMarked( this ); +} + +//================================================================================ +/*! + * \brief Store VTK ID + */ +//================================================================================ + +void SMDS_MeshElement::setVtkID( const int vtkID ) +{ + myHolder->SetVTKID( this, vtkID ); +} + +//================================================================================ +/*! + * \brief Return the mesh this element belongs to + */ +//================================================================================ + +SMDS_Mesh* SMDS_MeshElement::GetMesh() const +{ + return const_cast( myHolder )->GetMesh(); +} + +//================================================================================ +/*! + * \brief Return a SMDS_UnstructuredGrid + */ +//================================================================================ + +SMDS_UnstructuredGrid* SMDS_MeshElement::getGrid() const +{ + return const_cast( myHolder )->GetMesh()->GetGrid(); +} + +//================================================================================ +/*! + * \brief Print self + */ +//================================================================================ + +void SMDS_MeshElement::Print(ostream & OS) const +{ + OS << "dump of mesh element" << endl; +} + +ostream & operator <<(ostream & OS, const SMDS_MeshElement * e) +{ + e->Print(OS); + return OS; }