Commit 4863aff3 authored by untereiner's avatar untereiner

some changes

parent 08553e53
This diff is collapsed.
......@@ -339,28 +339,28 @@ bool MeshTablesVolume<PFP>::importTet(const std::string& filename, std::vector<s
typename PFP::VEC3 B = position[verticesID[s1]];
typename PFP::VEC3 C = position[verticesID[s2]];
// 1 hexa ok avec cette partie
if (Geom::testOrientation3D<typename PFP::VEC3>(P,A,B,C) == Geom::UNDER)
{
unsigned int pt[8];
pt[0] = s0;
pt[1] = s1;
pt[2] = s3;
pt[3] = s2;
pt[4] = s4;
pt[5] = s5;
pt[6] = s7;
pt[7] = s6;
// // 1 hexa ok avec cette partie
// if (Geom::testOrientation3D<typename PFP::VEC3>(P,A,B,C) == Geom::UNDER)
// {
// unsigned int pt[8];
// pt[0] = s0;
// pt[1] = s1;
// pt[2] = s3;
// pt[3] = s2;
// pt[4] = s4;
// pt[5] = s5;
// pt[6] = s7;
// pt[7] = s6;
s0 = pt[0];
s1 = pt[1];
s2 = pt[2];
s3 = pt[3];
s4 = pt[4];
s5 = pt[5];
s6 = pt[6];
s7 = pt[7];
}
// s0 = pt[0];
// s1 = pt[1];
// s2 = pt[2];
// s3 = pt[3];
// s4 = pt[4];
// s5 = pt[5];
// s6 = pt[6];
// s7 = pt[7];
// }
m_emb.push_back(verticesID[s0]);
m_emb.push_back(verticesID[s1]);
......
......@@ -158,6 +158,20 @@ public:
return (*m_markVector)[d_index].testMark(m_mark) ;
}
template <unsigned int ORBIT>
inline bool isMarked(Cell<ORBIT> c) const
{
assert(m_map.template getMarkerSet<DART>(m_thread).testMark(m_mark));
m_map.foreach_dart_of_orbit(c, [&] (Dart d)
{
if(!(*m_markVector)[m_map.dartIndex(d)].testMark(m_mark))
return false;
}
, m_thread) ;
return true;
}
/**
* mark the darts of the given cell
*/
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment