From 3ae392a03db81a84450302f83a7968f90d89394b Mon Sep 17 00:00:00 2001
From: Boris Sedji <sedji.boris@hotmail.com>
Date: Tue, 4 May 2010 12:11:25 +0000
Subject: [PATCH] *filter on level set physical tag
---
Solver/filters.h | 61 ++++++++++++++++++++++++++++++++++++++++++++++--
1 file changed, 59 insertions(+), 2 deletions(-)
diff --git a/Solver/filters.h b/Solver/filters.h
index 935e4fa662..8323bc06de 100644
--- a/Solver/filters.h
+++ b/Solver/filters.h
@@ -8,11 +8,13 @@
//
//
-#ifndef _XFEMFS_H_
-#define _XFEMFS_H_
+#ifndef _FILTERS_H_
+#define _FILTERS_H_
#include "simpleFunction.h"
#include "dofManager.h"
+#include "GModel.h"
+#include "groupOfElements.h"
class FilterNodeEnriched
{
@@ -64,5 +66,60 @@ class FilterNodeEnriched
// }
};
+class FilterElementsCutByLevelSet
+{
+
+ private :
+
+ std::set<int> _TagEnrichedVertex;
+ std::pair<int,int> _LevelSetEntity;
+ std::set<int> *_EnrichComp;
+
+ public :
+
+ FilterElementsCutByLevelSet(std::pair<int,int> LevelSetEntity , std::set<int> * EnrichComp)
+ {
+
+ _EnrichComp = EnrichComp;
+ _LevelSetEntity = LevelSetEntity;
+
+ // groupOfElements to get all the elements associate with the level set -- (work with *current GModel)
+ groupOfElements *LevelSetElements = new groupOfElements (_LevelSetEntity.first, _LevelSetEntity.second);
+
+ // tag enriched vertex determination
+ std::set<MElement*>::const_iterator it = LevelSetElements->begin();
+ for (; it != LevelSetElements->end(); it++)
+ {
+ MElement *e = *it;
+ if (e->getParent()) // if element got parents
+ {
+ for (int k = 0; k < e->getParent()->getNumVertices(); ++k)
+ { // for all vertices in the element parent
+ _TagEnrichedVertex.insert(e->getParent()->getVertex(k)->getNum());
+ }
+ }
+ }
+
+ }
+
+ virtual bool operator () (Dof & key) const
+ {
+ std::set<int>::iterator it1;
+ std::set<int>::iterator it2;
+ int i1,i2;
+ Dof::getTwoIntsFromType(key.getType(), i1,i2);
+ it2 = _EnrichComp->find(i1);
+ it1 = _TagEnrichedVertex.find(key.getEntity());
+ if (it1!=_TagEnrichedVertex.end() & it2 != _EnrichComp->end())
+ {
+ return true;
+ }
+ else return false;
+ }
+
+
+};
+
+
#endif
--
GitLab