Commit 789c24ee authored by Thomas's avatar Thomas
Browse files

Geometry : chgt de qq bornes

parent c021379e
......@@ -175,8 +175,8 @@ void ParticleCell2D<PFP>::faceState(const VEC3& current)
std::cout << "faceState" << d << std::endl;
#endif
// assert(std::isfinite(m_position[0]) && std::isfinite(m_position[1]) && std::isfinite(m_position[2]));
// assert(std::isfinite(current[0]) && std::isfinite(current[1]) && std::isfinite(current[2]));
assert(std::isfinite(m_position[0]) && std::isfinite(m_position[1]) && std::isfinite(m_position[2]));
assert(std::isfinite(current[0]) && std::isfinite(current[1]) && std::isfinite(current[2]));
// assert(Algo::Geometry::isPointInConvexFace2D<PFP>(m,d,m_positions,m_position,true));
Dart dd=d;
......
......@@ -163,7 +163,7 @@ bool arePointsEquals(const VEC3& point1, const VEC3& point2)
{
VEC3 v(point1 - point2);
#define PRECISION 1e-6
#define PRECISION 1e-3
return v.norm2() <= PRECISION ;
#undef PRECISION
}
......
......@@ -36,7 +36,7 @@ Orientation2D testOrientation2D(const VEC3& P, const VEC3& Pa, const VEC3& Pb)
typedef typename VEC3::DATA_TYPE T ;
// const T min = std::numeric_limits<T>::min()*T(100);
const T min = 0.00001;
const T min = 0.000001;
T wsof = (Pa[0]-P[0])*(P[1]-Pb[1])-(P[0]-Pb[0])*(Pa[1]-P[1]);
......
Supports Markdown
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