diff --git a/GeoLib/AnalyticalGeometry.cpp b/GeoLib/AnalyticalGeometry.cpp
index 00ae380c172347c42de1de72bfff56f81c608dcd..e4d8aae3bc84a5c268d032afc6e57c6275c9011d 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 3c11356b5e9a1b791dcf02f0794e93efe8c9f3db..a8a8cb204efaf7da2527fd755de985cbace1c71e 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 72308445b4fb6f9a92173e0d0db66451c3e053e2..383c3a6d26be9e1c2cf0bbd7855a8a02024c5cd6 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 6c767316171b038324f3d034ff2a6c119ed9576b..d9a5ad0babfb597da6c7960fb4adf0b36517229c 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 1b863420a75175e19ebf5263dc75c61bd700eb2c..7a462b3357980099c836350cebbd780b19a9decd 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 33b7f6805d36e6eb33d422be0ac8978f85339dc7..43231955e5ccd4d44aa4c7f1ecb47733cc6497f3 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 b22ba9a7fa12b4e638da2b1457aa995d63fd31fb..3c3a1499c666c1319e3de5ebf296941c2efad36b 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 3ae74c36b07ba6255a305ff0d7dbade197b55210..3a3d661baea87eedf4d177037a983a9b2bbc0ca6 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 52cf5bb251a24d39b08d695ad58bb05b079d2a36..b460682cd20d5014e2498d874771521e38b8046a 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 4b32a2e359cd3c68862baff380b92c0edf987acc..ea98f9ce71fb009c14a54cf665b834bdc5bfeedb 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 823562a2995ff41e9d2ebd939b9cf02e58092bb5..a1a3f6ee5413daf7e5dccb57056aad7dd6150109 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 3c0a0d2f509b6eb44516ac8bfcc66843664b8b2e..620075d88c404f532f2fc839c4965b50752da814 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 5c6e7bb4d4a829cb7b76284e9926fbe288dd6aff..0c63f01fd3573fba812ba8de1437a29067feb8ce 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 032e16eb70eef34c7bb3923768b26ffc32363d94..407c2b98056e693b2319d30436cd30b69e3d3cf4 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 c3c990a4409c3b69bb545fea0f604dd0a550cc91..5d70c1778a6e9756ce27911acfc46312051bf594 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 6b89f18050621f59e9460dfe705089790d021a45..7fee311b3484c277eaa9eb096630529e7cfb7d23 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 a4d2e0b735e29a4dc51fdb6cda2e83549ed7557a..04f2f581e27bd33c28c2bc6e62efeedce8a163db 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 8f340dd4255b8477719da1ad4c1db1638e802866..0ddaf4983489263af8efce6e282a867376b0413c 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 281829ac3fc8b0c19f58f7672eee8dcca305cda1..f14311ad3f5691c100a60a8e6a10b69d8d0f65f8 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++;
}