src/SMESH/SMESH_HypoFilter.cxx

   1 // Copyright (C) 2007-2013  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
  23 //  SMESH SMESH : implementaion of SMESH idl descriptions
  24 //  File   : SMESH_HypoFilter.cxx
  25 //  Module : SMESH
  26 //
  27 #include "SMESH_HypoFilter.hxx"
  28
  29 #include "SMESH_Gen.hxx"
  30 #include "SMESH_Hypothesis.hxx"
  31 #include "SMESH_MesherHelper.hxx"
  32 #include "SMESH_subMesh.hxx"
  33
  34 #include <TopExp_Explorer.hxx>
  35
  36 using namespace std;
  37
  38
  39 //=======================================================================
  40 //function : NamePredicate::Value
  41 //purpose  :
  42 //=======================================================================
  43
  44 bool SMESH_HypoFilter::NamePredicate::IsOk (const SMESH_Hypothesis* aHyp,
  45                                             const TopoDS_Shape&     /*aShape*/ ) const
  46 {
  47   return ( _name == aHyp->GetName() );
  48 }
  49
  50 //=======================================================================
  51 //function : TypePredicate::Value
  52 //purpose  :
  53 //=======================================================================
  54
  55 int SMESH_HypoFilter::TypePredicate::Value( const SMESH_Hypothesis* aHyp ) const
  56 {
  57   return aHyp->GetType();
  58 };
  59
  60 //=======================================================================
  61 //function : DimPredicate::Value
  62 //purpose  :
  63 //=======================================================================
  64
  65 int SMESH_HypoFilter::DimPredicate::Value( const SMESH_Hypothesis* aHyp ) const
  66 {
  67   return aHyp->GetDim();
  68 }
  69
  70 //=======================================================================
  71 //function : ApplicablePredicate::IsOk
  72 //purpose  :
  73 //=======================================================================
  74
  75 bool SMESH_HypoFilter::ApplicablePredicate::IsOk(const SMESH_Hypothesis* aHyp,
  76                                                  const TopoDS_Shape&     /*aShape*/) const
  77 {
  78   return SMESH_subMesh::IsApplicableHypotesis( aHyp, (TopAbs_ShapeEnum)_shapeType );
  79 };
  80
  81 //=======================================================================
  82 //function : IsAuxiliaryPredicate::IsOk
  83 //purpose  :
  84 //=======================================================================
  85
  86 bool SMESH_HypoFilter::IsAuxiliaryPredicate::IsOk(const SMESH_Hypothesis* aHyp,
  87                                                   const TopoDS_Shape&     /*aShape*/) const
  88 {
  89   return aHyp->IsAuxiliary();
  90 };
  91
  92 //=======================================================================
  93 //function : ApplicablePredicate::ApplicablePredicate
  94 //purpose  :
  95 //=======================================================================
  96
  97 SMESH_HypoFilter::ApplicablePredicate::ApplicablePredicate( const TopoDS_Shape& theShape )
  98 {
  99   _shapeType = ( theShape.IsNull() ? TopAbs_SHAPE : theShape.ShapeType());
 100 }
 101
 102 //=======================================================================
 103 //function : InstancePredicate::IsOk
 104 //purpose  :
 105 //=======================================================================
 106
 107 bool SMESH_HypoFilter::InstancePredicate::IsOk(const SMESH_Hypothesis* aHyp,
 108                                                const TopoDS_Shape&     /*aShape*/) const
 109 {
 110   return _hypo == aHyp;
 111 }
 112
 113 //=======================================================================
 114 //function : IsAssignedToPredicate::IsOk
 115 //purpose  :
 116 //=======================================================================
 117
 118 bool SMESH_HypoFilter::IsAssignedToPredicate::IsOk(const SMESH_Hypothesis* aHyp,
 119                                                    const TopoDS_Shape&     aShape) const
 120 {
 121   return ( !_mainShape.IsNull() && !aShape.IsNull() && _mainShape.IsSame( aShape ));
 122 }
 123
 124 //================================================================================
 125 /*!
 126  * \brief Finds shapes preferable over _shape due to sub-mesh order
 127  */
 128 //================================================================================
 129
 130 void SMESH_HypoFilter::IsMoreLocalThanPredicate::findPreferable()
 131 {
 132   const int shapeID = _mesh.GetMeshDS()->ShapeToIndex( _shape );
 133   const TListOfListOfInt& listOfShapeIDList = _mesh.GetMeshOrder();
 134   TListOfListOfInt::const_iterator listsIt = listOfShapeIDList.begin();
 135   for ( ; listsIt != listOfShapeIDList.end(); ++listsIt )
 136   {
 137     const TListOfInt& idList  = *listsIt;
 138     TListOfInt::const_iterator idIt =
 139       std::find( idList.begin(), idList.end(), shapeID );
 140     if ( idIt != idList.end() && *idIt != idList.front() )
 141     {
 142       for ( ; idIt != idList.end(); --idIt )
 143       {
 144         const TopoDS_Shape& shape = _mesh.GetMeshDS()->IndexToShape( *idIt );
 145         if ( !shape.IsNull())
 146           _preferableShapes.Add( shape );
 147       }
 148     }
 149   }
 150 }
 151
 152 //=======================================================================
 153 //function : IsMoreLocalThanPredicate::IsOk
 154 //purpose  :
 155 //=======================================================================
 156
 157 bool SMESH_HypoFilter::IsMoreLocalThanPredicate::IsOk(const SMESH_Hypothesis* aHyp,
 158                                                       const TopoDS_Shape&     aShape) const
 159 {
 160   if ( aShape.IsSame( _mesh.GetShapeToMesh() ))
 161     return false; // aHyp is global
 162
 163   if ( SMESH_MesherHelper::IsSubShape( aShape, /*mainShape=*/_shape ))
 164     return true;
 165
 166   if ( aShape.ShapeType() == TopAbs_COMPOUND &&
 167        !SMESH_MesherHelper::IsSubShape( _shape, /*mainShape=*/aShape)) // issue 0020963
 168   {
 169     for ( int type = TopAbs_SOLID; type < TopAbs_SHAPE; ++type )
 170       if ( aHyp->GetDim() == SMESH_Gen::GetShapeDim( TopAbs_ShapeEnum( type )))
 171         for ( TopExp_Explorer exp( aShape, TopAbs_ShapeEnum( type )); exp.More(); exp.Next())
 172           if ( SMESH_MesherHelper::IsSubShape( exp.Current(), /*mainShape=*/_shape ))
 173             return true;
 174   }
 175
 176   if ( _preferableShapes.Contains( aShape ))
 177     return true; // issue 21559, Mesh_6
 178
 179   return false;
 180 }
 181
 182 //=======================================================================
 183 //function : SMESH_HypoFilter
 184 //purpose  :
 185 //=======================================================================
 186
 187 SMESH_HypoFilter::SMESH_HypoFilter()
 188 {
 189 }
 190
 191 //=======================================================================
 192 //function : SMESH_HypoFilter
 193 //purpose  :
 194 //=======================================================================
 195
 196 SMESH_HypoFilter::SMESH_HypoFilter( SMESH_HypoPredicate* aPredicate, bool notNagate )
 197 {
 198   add( notNagate ? AND : AND_NOT, aPredicate );
 199 }
 200
 201 //=======================================================================
 202 //function : And
 203 //purpose  :
 204 //=======================================================================
 205
 206 SMESH_HypoFilter & SMESH_HypoFilter::And( SMESH_HypoPredicate* aPredicate )
 207 {
 208   add( AND, aPredicate );
 209   return *this;
 210 }
 211
 212 //=======================================================================
 213 //function : AndNot
 214 //purpose  :
 215 //=======================================================================
 216
 217 SMESH_HypoFilter & SMESH_HypoFilter::AndNot( SMESH_HypoPredicate* aPredicate )
 218 {
 219   add( AND_NOT, aPredicate );
 220   return *this;
 221 }
 222
 223 //=======================================================================
 224 //function : Or
 225 //purpose  :
 226 //=======================================================================
 227
 228 SMESH_HypoFilter & SMESH_HypoFilter::Or( SMESH_HypoPredicate* aPredicate )
 229 {
 230   add( OR, aPredicate );
 231   return *this;
 232 }
 233
 234 //=======================================================================
 235 //function : OrNot
 236 //purpose  : Return predicates
 237 //=======================================================================
 238
 239 SMESH_HypoFilter & SMESH_HypoFilter::OrNot( SMESH_HypoPredicate* aPredicate )
 240 {
 241   add( OR_NOT, aPredicate );
 242   return *this;
 243 }
 244
 245 //=======================================================================
 246 //function : Is
 247 //purpose  :
 248 //=======================================================================
 249
 250 SMESH_HypoPredicate* SMESH_HypoFilter::Is(const SMESH_Hypothesis* theHypo)
 251 {
 252   return new InstancePredicate( theHypo );
 253 }
 254
 255 //=======================================================================
 256 //function : IsAlgo
 257 //purpose  :
 258 //=======================================================================
 259
 260 SMESH_HypoPredicate* SMESH_HypoFilter::IsAlgo()
 261 {
 262   return new TypePredicate( MORE, SMESHDS_Hypothesis::PARAM_ALGO );
 263 }
 264
 265 //=======================================================================
 266 //function : IsAuxiliary
 267 //purpose  :
 268 //=======================================================================
 269
 270 SMESH_HypoPredicate* SMESH_HypoFilter::IsAuxiliary()
 271 {
 272   return new IsAuxiliaryPredicate();
 273 }
 274
 275
 276 //=======================================================================
 277 //function : IsGlobal
 278 //purpose  :
 279 //=======================================================================
 280
 281  SMESH_HypoPredicate* SMESH_HypoFilter::IsGlobal(const TopoDS_Shape& theMainShape)
 282 {
 283   return new IsAssignedToPredicate( theMainShape );
 284 }
 285
 286 //=======================================================================
 287 //function : IsAssignedTo
 288 //purpose  :
 289 //=======================================================================
 290
 291  SMESH_HypoPredicate* SMESH_HypoFilter::IsAssignedTo(const TopoDS_Shape& theShape)
 292 {
 293   return new IsAssignedToPredicate( theShape );
 294 }
 295
 296 //=======================================================================
 297 //function : HasName
 298 //purpose  :
 299 //=======================================================================
 300
 301 SMESH_HypoPredicate* SMESH_HypoFilter::HasName(const string & theName)
 302 {
 303   return new NamePredicate( theName );
 304 }
 305
 306 //=======================================================================
 307 //function : HasDim
 308 //purpose  :
 309 //=======================================================================
 310
 311 SMESH_HypoPredicate* SMESH_HypoFilter::HasDim(const int theDim)
 312 {
 313   return new DimPredicate( EQUAL, theDim );
 314 }
 315
 316 //=======================================================================
 317 //function : IsApplicableTo
 318 //purpose  :
 319 //=======================================================================
 320
 321 SMESH_HypoPredicate* SMESH_HypoFilter::IsApplicableTo(const TopoDS_Shape& theShape)
 322 {
 323   return new ApplicablePredicate( theShape );
 324 }
 325
 326 //=======================================================================
 327 //function : IsMoreLocalThan
 328 //purpose  :
 329 //=======================================================================
 330
 331 SMESH_HypoPredicate* SMESH_HypoFilter::IsMoreLocalThan(const TopoDS_Shape& theShape,
 332                                                        const SMESH_Mesh&   theMesh)
 333 {
 334   return new IsMoreLocalThanPredicate( theShape, theMesh );
 335 }
 336
 337 //=======================================================================
 338 //function : HasType
 339 //purpose  :
 340 //=======================================================================
 341
 342 SMESH_HypoPredicate* SMESH_HypoFilter::HasType(const int theHypType)
 343 {
 344   return new TypePredicate( EQUAL, theHypType );
 345 }
 346
 347 //=======================================================================
 348 //function : IsOk
 349 //purpose  :
 350 //=======================================================================
 351
 352 bool SMESH_HypoFilter::IsOk (const SMESH_Hypothesis* aHyp,
 353                              const TopoDS_Shape&     aShape) const
 354 {
 355   if ( myPredicates.empty() )
 356     return true;
 357
 358   bool ok = ( myPredicates.front()->_logical_op <= AND_NOT );
 359   list<SMESH_HypoPredicate*>::const_iterator pred = myPredicates.begin();
 360   for ( ; pred != myPredicates.end(); ++pred )
 361   {
 362     bool ok2 = (*pred)->IsOk( aHyp, aShape );
 363     switch ( (*pred)->_logical_op ) {
 364     case AND:     ok = ok && ok2; break;
 365     case AND_NOT: ok = ok && !ok2; break;
 366     case OR:      ok = ok || ok2; break;
 367     case OR_NOT:  ok = ok || !ok2; break;
 368     default:;
 369     }
 370   }
 371   return ok;
 372 }
 373
 374 //=======================================================================
 375 //function : Init
 376 //purpose  :
 377 //=======================================================================
 378
 379 SMESH_HypoFilter & SMESH_HypoFilter::Init  ( SMESH_HypoPredicate* aPredicate, bool notNagate )
 380 {
 381   list<SMESH_HypoPredicate*>::const_iterator pred = myPredicates.begin();
 382   for ( ; pred != myPredicates.end(); ++pred )
 383     delete *pred;
 384   myPredicates.clear();
 385
 386   add( notNagate ? AND : AND_NOT, aPredicate );
 387   return *this;
 388 }
 389
 390
 391 //=======================================================================
 392 //function : IsOk
 393 //purpose  :
 394 //=======================================================================
 395
 396 SMESH_HypoFilter::~SMESH_HypoFilter()
 397 {
 398   Init(0);
 399 }
 400
 401
 402