src/StdMeshers/StdMeshers_MaxElementVolume.cxx

   1 //  Copyright (C) 2007-2008  CEA/DEN, EDF R&D, OPEN CASCADE
   2 //
   3 //  Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   4 //  CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   5 //
   6 //  This library is free software; you can redistribute it and/or
   7 //  modify it under the terms of the GNU Lesser General Public
   8 //  License as published by the Free Software Foundation; either
   9 //  version 2.1 of the License.
  10 //
  11 //  This library is distributed in the hope that it will be useful,
  12 //  but WITHOUT ANY WARRANTY; without even the implied warranty of
  13 //  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14 //  Lesser General Public License for more details.
  15 //
  16 //  You should have received a copy of the GNU Lesser General Public
  17 //  License along with this library; if not, write to the Free Software
  18 //  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  19 //
  20 //  See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  21 //
  22 //  SMESH SMESH : implementaion of SMESH idl descriptions
  23 //  File   : StdMeshers_MaxElementVolume.cxx
  24 //           Moved here from SMESH_MaxElementVolume.cxx
  25 //  Author : Paul RASCLE, EDF
  26 //  Module : SMESH
  27 //  $Header$
  28 //
  29 #include "StdMeshers_MaxElementVolume.hxx"
  30
  31 #include "SMDS_MeshElement.hxx"
  32 #include "SMESHDS_SubMesh.hxx"
  33 #include "SMESH_ControlsDef.hxx"
  34 #include "SMESH_Mesh.hxx"
  35
  36 #include "utilities.h"
  37
  38 #include <TopExp.hxx>
  39 #include <TopExp_Explorer.hxx>
  40 #include <TopTools_IndexedMapOfShape.hxx>
  41
  42 using namespace std;
  43
  44 //=============================================================================
  45 /*!
  46  *
  47  */
  48 //=============================================================================
  49
  50 StdMeshers_MaxElementVolume::StdMeshers_MaxElementVolume(int hypId, int studyId, SMESH_Gen* gen)
  51   : SMESH_Hypothesis(hypId, studyId, gen)
  52 {
  53   _maxVolume = 1.;
  54   _name = "MaxElementVolume";
  55   _param_algo_dim = 3;
  56 }
  57
  58 //=============================================================================
  59 /*!
  60  *
  61  */
  62 //=============================================================================
  63
  64 StdMeshers_MaxElementVolume::~StdMeshers_MaxElementVolume()
  65 {
  66   MESSAGE("StdMeshers_MaxElementVolume::~StdMeshers_MaxElementVolume");
  67 }
  68
  69 //=============================================================================
  70 /*!
  71  *
  72  */
  73 //=============================================================================
  74
  75 void StdMeshers_MaxElementVolume::SetMaxVolume(double maxVolume)
  76   throw (SALOME_Exception)
  77 {
  78   double oldVolume = _maxVolume;
  79   if (maxVolume <= 0)
  80     throw SALOME_Exception(LOCALIZED("maxVolume must be positive"));
  81   _maxVolume = maxVolume;
  82   if (_maxVolume != oldVolume)
  83     NotifySubMeshesHypothesisModification();
  84 }
  85
  86 //=============================================================================
  87 /*!
  88  *
  89  */
  90 //=============================================================================
  91
  92 double StdMeshers_MaxElementVolume::GetMaxVolume() const
  93 {
  94   return _maxVolume;
  95 }
  96
  97 //=============================================================================
  98 /*!
  99  *
 100  */
 101 //=============================================================================
 102
 103 ostream & StdMeshers_MaxElementVolume::SaveTo(ostream & save)
 104 {
 105   save << this->_maxVolume;
 106   return save;
 107 }
 108
 109 //=============================================================================
 110 /*!
 111  *
 112  */
 113 //=============================================================================
 114
 115 istream & StdMeshers_MaxElementVolume::LoadFrom(istream & load)
 116 {
 117   bool isOK = true;
 118   double a;
 119   isOK = (load >> a);
 120   if (isOK)
 121     this->_maxVolume = a;
 122   else
 123     load.clear(ios::badbit | load.rdstate());
 124   return load;
 125 }
 126
 127 //=============================================================================
 128 /*!
 129  *
 130  */
 131 //=============================================================================
 132
 133 ostream & operator << (ostream & save, StdMeshers_MaxElementVolume & hyp)
 134 {
 135   return hyp.SaveTo( save );
 136 }
 137
 138 //=============================================================================
 139 /*!
 140  *
 141  */
 142 //=============================================================================
 143
 144 istream & operator >> (istream & load, StdMeshers_MaxElementVolume & hyp)
 145 {
 146   return hyp.LoadFrom( load );
 147 }
 148
 149
 150 //================================================================================
 151 /*!
 152  * \brief Initialize maximal area by the mesh built on the geometry
 153  * \param theMesh - the built mesh
 154  * \param theShape - the geometry of interest
 155  * \retval bool - true if parameter values have been successfully defined
 156  */
 157 //================================================================================
 158
 159 bool StdMeshers_MaxElementVolume::SetParametersByMesh(const SMESH_Mesh*   theMesh,
 160                                                       const TopoDS_Shape& theShape)
 161 {
 162   if ( !theMesh || theShape.IsNull() )
 163     return false;
 164
 165   _maxVolume = 0;
 166
 167   SMESH::Controls::Volume volumeControl;
 168
 169   TopTools_IndexedMapOfShape volMap;
 170   TopExp::MapShapes( theShape, TopAbs_SOLID, volMap );
 171   if ( volMap.IsEmpty() )
 172     TopExp::MapShapes( theShape, TopAbs_SHELL, volMap );
 173   if ( volMap.IsEmpty() )
 174     return false;
 175
 176   SMESHDS_Mesh* aMeshDS = const_cast< SMESH_Mesh* >( theMesh )->GetMeshDS();
 177
 178   for ( int iV = 1; iV <= volMap.Extent(); ++iV )
 179   {
 180     const TopoDS_Shape& S = volMap( iV );
 181     SMESHDS_SubMesh * subMesh = aMeshDS->MeshElements( S );
 182     if ( !subMesh && S.ShapeType() == TopAbs_SOLID ) {
 183       TopExp_Explorer shellExp( S, TopAbs_SHELL );
 184       if ( shellExp.More() )
 185         subMesh = aMeshDS->MeshElements( shellExp.Current() );
 186     }
 187     if ( !subMesh)
 188       return false;
 189     SMDS_ElemIteratorPtr vIt = subMesh->GetElements();
 190     while ( vIt->more() )
 191     {
 192       const SMDS_MeshElement* elem = vIt->next();
 193       if ( elem->GetType() == SMDSAbs_Volume ) {
 194         _maxVolume = max( _maxVolume, volumeControl.GetValue( elem->GetID() ));
 195       }
 196     }
 197   }
 198   return _maxVolume > 0;
 199 }