From 55c49fd8d4662157a3e64ddbaf64a623d3dd577b Mon Sep 17 00:00:00 2001
From: Karsten Rink <karsten.rink@ufz.de>
Date: Wed, 24 Sep 2014 14:59:46 +0200
Subject: [PATCH] added optional eps for 3d calculation
---
GeoLib/AnalyticalGeometry.cpp | 16 ++++++++++------
GeoLib/AnalyticalGeometry.h | 2 +-
MeshLib/Elements/Element.h | 11 ++++++++---
MeshLib/Elements/Face.h | 8 --------
MeshLib/Elements/TemplateHex-impl.h | 14 +++++++-------
MeshLib/Elements/TemplateHex.h | 9 +++++++--
MeshLib/Elements/TemplateLine-impl.h | 4 ++--
MeshLib/Elements/TemplateLine.h | 9 +++++++--
MeshLib/Elements/TemplatePrism-impl.h | 8 ++++----
MeshLib/Elements/TemplatePrism.h | 9 +++++++--
MeshLib/Elements/TemplatePyramid-impl.h | 6 +++---
MeshLib/Elements/TemplatePyramid.h | 9 +++++++--
MeshLib/Elements/TemplateQuad-impl.h | 6 ++++--
MeshLib/Elements/TemplateQuad.h | 11 ++++-------
MeshLib/Elements/TemplateTet-impl.h | 4 ++--
MeshLib/Elements/TemplateTet.h | 9 +++++++--
MeshLib/Elements/TemplateTri-impl.h | 2 +-
MeshLib/Elements/TemplateTri.h | 11 ++++-------
.../Mesh2MeshPropertyInterpolation.cpp | 2 +-
19 files changed, 86 insertions(+), 64 deletions(-)
diff --git a/GeoLib/AnalyticalGeometry.cpp b/GeoLib/AnalyticalGeometry.cpp
index 00ae380c172..e4d8aae3bc8 100644
--- a/GeoLib/AnalyticalGeometry.cpp
+++ b/GeoLib/AnalyticalGeometry.cpp
@@ -243,24 +243,28 @@ bool isPointInTriangle(GeoLib::Point const& q,
}
bool isPointInTetrahedron(GeoLib::Point const& p, GeoLib::Point const& a, GeoLib::Point const& b,
- GeoLib::Point const& c, GeoLib::Point const& d)
+ GeoLib::Point const& c, GeoLib::Point const& d, double eps)
{
double const d0 (orient3d(d,a,b,c));
// if tetrahedron is not coplanar
- if (d0 != 0)
+ if (std::abs(d0) > std::numeric_limits<double>::epsilon())
{
bool const d0_sign (d0>0);
// if p is on the same side of bcd as a
- if (d0_sign != (orient3d(d, p, b, c)>=0))
+ double const d1 (orient3d(d, p, b, c));
+ if (!(d0_sign == (d1>=0) || std::abs(d1) < eps))
return false;
// if p is on the same side of acd as b
- if (d0_sign != (orient3d(d, a, p, c)>=0))
+ double const d2 (orient3d(d, a, p, c));
+ if (!(d0_sign == (d2>=0) || std::abs(d2) < eps))
return false;
// if p is on the same side of abd as c
- if (d0_sign != (orient3d(d, a, b, p)>=0))
+ double const d3 (orient3d(d, a, b, p));
+ if (!(d0_sign == (d3>=0) || std::abs(d3) < eps))
return false;
// if p is on the same side of abc as d
- if (d0_sign != (orient3d(p, a, b, c)>=0))
+ double const d4 (orient3d(p, a, b, c));
+ if (!(d0_sign == (d4>=0) || std::abs(d4) < eps))
return false;
return true;
}
diff --git a/GeoLib/AnalyticalGeometry.h b/GeoLib/AnalyticalGeometry.h
index 3c11356b5e9..a8a8cb204ef 100644
--- a/GeoLib/AnalyticalGeometry.h
+++ b/GeoLib/AnalyticalGeometry.h
@@ -139,7 +139,7 @@ bool isPointInTriangle(GeoLib::Point const& p,
* @return true if the test point p is not located outside of abcd (i.e. inside or on a plane/edge).
*/
bool isPointInTetrahedron(GeoLib::Point const& p, GeoLib::Point const& a, GeoLib::Point const& b,
- GeoLib::Point const& c, GeoLib::Point const& d);
+ GeoLib::Point const& c, GeoLib::Point const& d, double eps = std::numeric_limits<double>::epsilon());
/**
* test for intersections of the line segments of the Polyline
diff --git a/MeshLib/Elements/Element.h b/MeshLib/Elements/Element.h
index 72308445b4f..383c3a6d26b 100644
--- a/MeshLib/Elements/Element.h
+++ b/MeshLib/Elements/Element.h
@@ -19,7 +19,7 @@
#include <limits>
#include <boost/optional.hpp>
-#include "Point.h"
+#include "GeoLib/Point.h"
#include "MeshLib/MeshEnums.h"
#include "MeshLib/Mesh.h"
@@ -155,8 +155,13 @@ public:
/// Returns true if these two indeces form an edge and false otherwise
virtual bool isEdge(unsigned i, unsigned j) const = 0;
- /// Returns true if the given point is not located outside of the element
- virtual bool isPntInElement(GeoLib::Point const& pnt) const = 0;
+ /**
+ * Checks if a point is inside the element.
+ * @param pnt a 3D GeoLib::Point object
+ * @param eps tolerance for numerical algorithm used or computing the property
+ * @return true if the point is not outside the element, false otherwise
+ */
+ virtual bool isPntInElement(GeoLib::Point const& pnt, double eps) const = 0;
/**
* Tests if the element is geometrically valid.
diff --git a/MeshLib/Elements/Face.h b/MeshLib/Elements/Face.h
index 6c767316171..d9a5ad0babf 100644
--- a/MeshLib/Elements/Face.h
+++ b/MeshLib/Elements/Face.h
@@ -59,14 +59,6 @@ public:
/// Destructor
virtual ~Face();
- /**
- * Check if the 3d GeoLib::Point is inside of the element.
- * @param pnt the 3d GeoLib::Point object
- * @param eps tolerance for numerical algorithm used or computing the property
- * @return true if the point is inside the element, false otherwise
- */
- virtual bool isPntInside(GeoLib::Point const& pnt, double eps = std::numeric_limits<double>::epsilon()) const = 0;
-
/**
* This method is pure virtual and is inherited from class @sa Element.
* It has to be implemented in the derived classes of class Face!
diff --git a/MeshLib/Elements/TemplateHex-impl.h b/MeshLib/Elements/TemplateHex-impl.h
index 1b863420a75..7a462b33579 100644
--- a/MeshLib/Elements/TemplateHex-impl.h
+++ b/MeshLib/Elements/TemplateHex-impl.h
@@ -140,14 +140,14 @@ bool TemplateHex<NNODES,CELLHEXTYPE>::isEdge(unsigned idx1, unsigned idx2) const
}
template <unsigned NNODES, CellType CELLHEXTYPE>
-bool TemplateHex<NNODES,CELLHEXTYPE>::isPntInElement(GeoLib::Point const& pnt) const
+bool TemplateHex<NNODES,CELLHEXTYPE>::isPntInElement(GeoLib::Point const& pnt, double eps) const
{
- return (GeoLib::isPointInTetrahedron(pnt, *_nodes[4], *_nodes[7], *_nodes[5], *_nodes[0]) ||
- GeoLib::isPointInTetrahedron(pnt, *_nodes[5], *_nodes[3], *_nodes[1], *_nodes[0]) ||
- GeoLib::isPointInTetrahedron(pnt, *_nodes[5], *_nodes[7], *_nodes[3], *_nodes[0]) ||
- GeoLib::isPointInTetrahedron(pnt, *_nodes[5], *_nodes[7], *_nodes[6], *_nodes[2]) ||
- GeoLib::isPointInTetrahedron(pnt, *_nodes[1], *_nodes[3], *_nodes[5], *_nodes[2]) ||
- GeoLib::isPointInTetrahedron(pnt, *_nodes[3], *_nodes[7], *_nodes[5], *_nodes[2]));
+ return (GeoLib::isPointInTetrahedron(pnt, *_nodes[4], *_nodes[7], *_nodes[5], *_nodes[0], eps) ||
+ GeoLib::isPointInTetrahedron(pnt, *_nodes[5], *_nodes[3], *_nodes[1], *_nodes[0], eps) ||
+ GeoLib::isPointInTetrahedron(pnt, *_nodes[5], *_nodes[7], *_nodes[3], *_nodes[0], eps) ||
+ GeoLib::isPointInTetrahedron(pnt, *_nodes[5], *_nodes[7], *_nodes[6], *_nodes[2], eps) ||
+ GeoLib::isPointInTetrahedron(pnt, *_nodes[1], *_nodes[3], *_nodes[5], *_nodes[2], eps) ||
+ GeoLib::isPointInTetrahedron(pnt, *_nodes[3], *_nodes[7], *_nodes[5], *_nodes[2], eps));
}
template <unsigned NNODES, CellType CELLHEXTYPE>
diff --git a/MeshLib/Elements/TemplateHex.h b/MeshLib/Elements/TemplateHex.h
index 33b7f6805d3..43231955e5c 100644
--- a/MeshLib/Elements/TemplateHex.h
+++ b/MeshLib/Elements/TemplateHex.h
@@ -104,8 +104,13 @@ public:
/// Returns true if these two indices form an edge and false otherwise
bool isEdge(unsigned i, unsigned j) const;
- /// Returns true if the given point is not located outside of the hexahedron
- bool isPntInElement(GeoLib::Point const& pnt) const;
+ /**
+ * Checks if a point is inside the element.
+ * @param pnt a 3D GeoLib::Point object
+ * @param eps tolerance for numerical algorithm used or computing the property
+ * @return true if the point is not outside the element, false otherwise
+ */
+ bool isPntInElement(GeoLib::Point const& pnt, double eps = std::numeric_limits<double>::epsilon()) const;
/**
* Tests if the element is geometrically valid.
diff --git a/MeshLib/Elements/TemplateLine-impl.h b/MeshLib/Elements/TemplateLine-impl.h
index b22ba9a7fa1..3c3a1499c66 100644
--- a/MeshLib/Elements/TemplateLine-impl.h
+++ b/MeshLib/Elements/TemplateLine-impl.h
@@ -56,12 +56,12 @@ TemplateLine<NNODES,CELLLINETYPE>::~TemplateLine()
{}
template <unsigned NNODES, CellType CELLLINETYPE>
-bool TemplateLine<NNODES,CELLLINETYPE>::isPntInElement(GeoLib::Point const& pnt) const
+bool TemplateLine<NNODES,CELLLINETYPE>::isPntInElement(GeoLib::Point const& pnt, double eps) const
{
double tmp;
double tmp_dst(0);
double const dist =MathLib::calcProjPntToLineAndDists(pnt.getCoords(), _nodes[0]->getCoords(), _nodes[1]->getCoords(), tmp, tmp_dst);
- return (dist < std::numeric_limits<double>::epsilon());
+ return (dist < eps);
}
template <unsigned NNODES, CellType CELLLINETYPE>
diff --git a/MeshLib/Elements/TemplateLine.h b/MeshLib/Elements/TemplateLine.h
index 3ae74c36b07..3a3d661baea 100644
--- a/MeshLib/Elements/TemplateLine.h
+++ b/MeshLib/Elements/TemplateLine.h
@@ -115,8 +115,13 @@ public:
return false;
}
- /// Returns true if pnt is located on the line segment and false otherwise
- bool isPntInElement(GeoLib::Point const& pnt) const;
+ /**
+ * Checks if a point is located on the line
+ * @param pnt a 3D GeoLib::Point object
+ * @param eps tolerance for numerical algorithm used or computing the property
+ * @return true if the point is located on the line, false otherwise
+ */
+ bool isPntInElement(GeoLib::Point const& pnt, double eps = std::numeric_limits<double>::epsilon()) const;
/**
* Tests if the element is geometrically valid.
diff --git a/MeshLib/Elements/TemplatePrism-impl.h b/MeshLib/Elements/TemplatePrism-impl.h
index 52cf5bb251a..b460682cd20 100644
--- a/MeshLib/Elements/TemplatePrism-impl.h
+++ b/MeshLib/Elements/TemplatePrism-impl.h
@@ -149,11 +149,11 @@ bool TemplatePrism<NNODES,CELLPRISMTYPE>::isEdge(unsigned idx1, unsigned idx2) c
}
template <unsigned NNODES, CellType CELLPRISMTYPE>
-bool TemplatePrism<NNODES,CELLPRISMTYPE>::isPntInElement(GeoLib::Point const& pnt) const
+bool TemplatePrism<NNODES,CELLPRISMTYPE>::isPntInElement(GeoLib::Point const& pnt, double eps) const
{
- return (GeoLib::isPointInTetrahedron(pnt, *_nodes[0], *_nodes[1], *_nodes[2], *_nodes[3]) ||
- GeoLib::isPointInTetrahedron(pnt, *_nodes[1], *_nodes[4], *_nodes[2], *_nodes[3]) ||
- GeoLib::isPointInTetrahedron(pnt, *_nodes[2], *_nodes[4], *_nodes[5], *_nodes[3]));
+ return (GeoLib::isPointInTetrahedron(pnt, *_nodes[0], *_nodes[1], *_nodes[2], *_nodes[3], eps) ||
+ GeoLib::isPointInTetrahedron(pnt, *_nodes[1], *_nodes[4], *_nodes[2], *_nodes[3], eps) ||
+ GeoLib::isPointInTetrahedron(pnt, *_nodes[2], *_nodes[4], *_nodes[5], *_nodes[3], eps));
}
template <unsigned NNODES, CellType CELLPRISMTYPE>
diff --git a/MeshLib/Elements/TemplatePrism.h b/MeshLib/Elements/TemplatePrism.h
index 4b32a2e359c..ea98f9ce71f 100644
--- a/MeshLib/Elements/TemplatePrism.h
+++ b/MeshLib/Elements/TemplatePrism.h
@@ -102,8 +102,13 @@ public:
/// Returns true if these two indeces form an edge and false otherwise
bool isEdge(unsigned i, unsigned j) const;
- /// Returns true if the given point is not located outside of the prism
- bool isPntInElement(GeoLib::Point const& pnt) const;
+ /**
+ * Checks if a point is inside the element.
+ * @param pnt a 3D GeoLib::Point object
+ * @param eps tolerance for numerical algorithm used or computing the property
+ * @return true if the point is not outside the element, false otherwise
+ */
+ bool isPntInElement(GeoLib::Point const& pnt, double eps = std::numeric_limits<double>::epsilon()) const;
/**
* Tests if the element is geometrically valid.
diff --git a/MeshLib/Elements/TemplatePyramid-impl.h b/MeshLib/Elements/TemplatePyramid-impl.h
index 823562a2995..a1a3f6ee541 100644
--- a/MeshLib/Elements/TemplatePyramid-impl.h
+++ b/MeshLib/Elements/TemplatePyramid-impl.h
@@ -151,10 +151,10 @@ bool TemplatePyramid<NNODES,CELLPYRAMIDTYPE>::isEdge(unsigned idx1, unsigned idx
}
template <unsigned NNODES, CellType CELLPYRAMIDTYPE>
-bool TemplatePyramid<NNODES,CELLPYRAMIDTYPE>::isPntInElement(GeoLib::Point const& pnt) const
+bool TemplatePyramid<NNODES,CELLPYRAMIDTYPE>::isPntInElement(GeoLib::Point const& pnt, double eps) const
{
- return (GeoLib::isPointInTetrahedron(pnt, *_nodes[0], *_nodes[1], *_nodes[2], *_nodes[4]) ||
- GeoLib::isPointInTetrahedron(pnt, *_nodes[0], *_nodes[2], *_nodes[3], *_nodes[4]));
+ return (GeoLib::isPointInTetrahedron(pnt, *_nodes[0], *_nodes[1], *_nodes[2], *_nodes[4], eps) ||
+ GeoLib::isPointInTetrahedron(pnt, *_nodes[0], *_nodes[2], *_nodes[3], *_nodes[4], eps));
}
template <unsigned NNODES, CellType CELLPYRAMIDTYPE>
diff --git a/MeshLib/Elements/TemplatePyramid.h b/MeshLib/Elements/TemplatePyramid.h
index 3c0a0d2f509..620075d88c4 100644
--- a/MeshLib/Elements/TemplatePyramid.h
+++ b/MeshLib/Elements/TemplatePyramid.h
@@ -100,8 +100,13 @@ public:
/// Returns true if these two indeces form an edge and false otherwise
bool isEdge(unsigned i, unsigned j) const;
- /// Returns true if the given point is not located outside of the pyramid
- bool isPntInElement(GeoLib::Point const& pnt) const;
+ /**
+ * Checks if a point is inside the element.
+ * @param pnt a 3D GeoLib::Point object
+ * @param eps tolerance for numerical algorithm used or computing the property
+ * @return true if the point is not outside the element, false otherwise
+ */
+ bool isPntInElement(GeoLib::Point const& pnt, double eps = std::numeric_limits<double>::epsilon()) const;
/**
* Tests if the element is geometrically valid.
diff --git a/MeshLib/Elements/TemplateQuad-impl.h b/MeshLib/Elements/TemplateQuad-impl.h
index 5c6e7bb4d4a..0c63f01fd35 100644
--- a/MeshLib/Elements/TemplateQuad-impl.h
+++ b/MeshLib/Elements/TemplateQuad-impl.h
@@ -106,10 +106,12 @@ bool TemplateQuad<NNODES,CELLQUADTYPE>::isEdge(unsigned idx1, unsigned idx2) con
}
template <unsigned NNODES, CellType CELLQUADTYPE>
-bool TemplateQuad<NNODES,CELLQUADTYPE>::isPntInside(GeoLib::Point const& pnt, double eps) const
+bool TemplateQuad<NNODES,CELLQUADTYPE>::isPntInElement(GeoLib::Point const& pnt, double eps) const
{
+ bool a (GeoLib::isPointInTriangle(pnt, *_nodes[0], *_nodes[1], *_nodes[2], eps));
+ bool b ( GeoLib::isPointInTriangle(pnt, *_nodes[0], *_nodes[2], *_nodes[3], eps));
return (GeoLib::isPointInTriangle(pnt, *_nodes[0], *_nodes[1], *_nodes[2], eps) ||
- GeoLib::isPointInTriangle(pnt, *_nodes[0], *_nodes[2], *_nodes[3], eps));
+ GeoLib::isPointInTriangle(pnt, *_nodes[0], *_nodes[2], *_nodes[3], eps));
}
template <unsigned NNODES, CellType CELLQUADTYPE>
diff --git a/MeshLib/Elements/TemplateQuad.h b/MeshLib/Elements/TemplateQuad.h
index 032e16eb70e..407c2b98056 100644
--- a/MeshLib/Elements/TemplateQuad.h
+++ b/MeshLib/Elements/TemplateQuad.h
@@ -89,16 +89,13 @@ public:
/// Returns true if these two indeces form an edge and false otherwise
bool isEdge(unsigned i, unsigned j) const;
- /// Returns true if the given point is not located outside of the quad
- bool isPntInElement(GeoLib::Point const& pnt) const { return isPntInside(pnt); }
-
/**
- * Check if the 3d GeoLib::Point is inside of the quad element.
- * @param pnt the 3d GeoLib::Point object
+ * Checks if a point is inside the element.
+ * @param pnt a 3D GeoLib::Point object
* @param eps tolerance for numerical algorithm used or computing the property
- * @return true if the point is inside the element, false otherwise
+ * @return true if the point is not outside the element, false otherwise
*/
- virtual bool isPntInside(GeoLib::Point const& pnt, double eps = std::numeric_limits<double>::epsilon()) const;
+ bool isPntInElement(GeoLib::Point const& pnt, double eps = std::numeric_limits<double>::epsilon()) const;
/**
* Tests if the element is geometrically valid.
diff --git a/MeshLib/Elements/TemplateTet-impl.h b/MeshLib/Elements/TemplateTet-impl.h
index c3c990a4409..5d70c1778a6 100644
--- a/MeshLib/Elements/TemplateTet-impl.h
+++ b/MeshLib/Elements/TemplateTet-impl.h
@@ -129,9 +129,9 @@ bool TemplateTet<NNODES,CELLTETTYPE>::isEdge(unsigned idx1, unsigned idx2) const
}
template <unsigned NNODES, CellType CELLTETTYPE>
-bool TemplateTet<NNODES,CELLTETTYPE>::isPntInElement(GeoLib::Point const& pnt) const
+bool TemplateTet<NNODES,CELLTETTYPE>::isPntInElement(GeoLib::Point const& pnt, double eps) const
{
- return GeoLib::isPointInTetrahedron(pnt, *_nodes[0], *_nodes[1], *_nodes[2], *_nodes[3]);
+ return GeoLib::isPointInTetrahedron(pnt, *_nodes[0], *_nodes[1], *_nodes[2], *_nodes[3], eps);
}
template <unsigned NNODES, CellType CELLTETTYPE>
diff --git a/MeshLib/Elements/TemplateTet.h b/MeshLib/Elements/TemplateTet.h
index 6b89f180506..7fee311b348 100644
--- a/MeshLib/Elements/TemplateTet.h
+++ b/MeshLib/Elements/TemplateTet.h
@@ -99,8 +99,13 @@ public:
/// Returns true if these two indeces form an edge and false otherwise
bool isEdge(unsigned i, unsigned j) const;
- /// Returns true if the given point is not located outside of the tetrahedron
- bool isPntInElement(GeoLib::Point const& pnt) const;
+ /**
+ * Checks if a point is inside the element.
+ * @param pnt a 3D GeoLib::Point object
+ * @param eps tolerance for numerical algorithm used or computing the property
+ * @return true if the point is not outside the element, false otherwise
+ */
+ bool isPntInElement(GeoLib::Point const& pnt, double eps = std::numeric_limits<double>::epsilon()) const;
/**
* Tests if the element is geometrically valid.
diff --git a/MeshLib/Elements/TemplateTri-impl.h b/MeshLib/Elements/TemplateTri-impl.h
index a4d2e0b735e..04f2f581e27 100644
--- a/MeshLib/Elements/TemplateTri-impl.h
+++ b/MeshLib/Elements/TemplateTri-impl.h
@@ -88,7 +88,7 @@ bool TemplateTri<NNODES,CELLTRITYPE>::isEdge(unsigned idx1, unsigned idx2) const
}
template <unsigned NNODES, CellType CELLTRITYPE>
-bool TemplateTri<NNODES,CELLTRITYPE>::isPntInside(GeoLib::Point const& pnt, double eps) const
+bool TemplateTri<NNODES,CELLTRITYPE>::isPntInElement(GeoLib::Point const& pnt, double eps) const
{
return GeoLib::isPointInTriangle(pnt, *_nodes[0], *_nodes[1], *_nodes[2], eps);
}
diff --git a/MeshLib/Elements/TemplateTri.h b/MeshLib/Elements/TemplateTri.h
index 8f340dd4255..0ddaf498348 100644
--- a/MeshLib/Elements/TemplateTri.h
+++ b/MeshLib/Elements/TemplateTri.h
@@ -93,16 +93,13 @@ public:
/// Returns true if these two indices form an edge and false otherwise
bool isEdge(unsigned idx1, unsigned idx2) const;
- /// Returns true if the given point is not located outside of the triangle
- bool isPntInElement(GeoLib::Point const& pnt) const { return isPntInside(pnt); }
-
/**
- * Check if the 3d GeoLib::Point is inside of the element.
- * @param pnt the 3d GeoLib::Point object
+ * Checks if a point is inside the element.
+ * @param pnt a 3D GeoLib::Point object
* @param eps tolerance for numerical algorithm used or computing the property
- * @return true if the point is inside the element, false otherwise
+ * @return true if the point is not outside the element, false otherwise
*/
- virtual bool isPntInside(GeoLib::Point const& pnt, double eps = std::numeric_limits<double>::epsilon()) const;
+ bool isPntInElement(GeoLib::Point const& pnt, double eps = std::numeric_limits<double>::epsilon()) const;
/**
* Tests if the element is geometrically valid
diff --git a/MeshLib/MeshEditing/Mesh2MeshPropertyInterpolation.cpp b/MeshLib/MeshEditing/Mesh2MeshPropertyInterpolation.cpp
index 281829ac3fc..f14311ad3f5 100644
--- a/MeshLib/MeshEditing/Mesh2MeshPropertyInterpolation.cpp
+++ b/MeshLib/MeshEditing/Mesh2MeshPropertyInterpolation.cpp
@@ -94,7 +94,7 @@ void Mesh2MeshPropertyInterpolation::interpolatePropertiesForMesh(Mesh *dest_mes
for (size_t j(0); j<n_nodes_in_vec; j++) {
MeshLib::Node const*const j_th_node((*i_th_vec)[j]);
if (elem_aabb.containsPoint(*j_th_node)) {
- if (dynamic_cast<MeshLib::Face*>(dest_elements[k])->isPntInside(* dynamic_cast<GeoLib::Point const*>(j_th_node), 30)) {
+ if (dest_elements[k]->isPntInElement(* dynamic_cast<GeoLib::Point const*>(j_th_node), 30)) {
dest_properties[k] += interpolated_src_node_properties[(*i_th_vec)[j]->getID()];
cnt++;
}
--
GitLab