diff --git a/Applications/FileIO/FEFLOW/FEFLOWMeshInterface.cpp b/Applications/FileIO/FEFLOW/FEFLOWMeshInterface.cpp
index fcbbd9f1d1cbe62029c29c064c596c09a6b36d41..22b47e28938137422ce4f7f721a0e612ede5329b 100644
--- a/Applications/FileIO/FEFLOW/FEFLOWMeshInterface.cpp
+++ b/Applications/FileIO/FEFLOW/FEFLOWMeshInterface.cpp
@@ -735,7 +735,7 @@ void FEFLOWMeshInterface::setMaterialIDs(
                 }
 
                 GeoLib::Polygon polygon(*poly, true);
-                if (polygon.isPntInPolygon(gpt[0], gpt[1], gpt[2]))
+                if (polygon.isPntInPolygon(gpt))
                 {
                     matId = j;
                     break;
diff --git a/Applications/Utils/MeshGeoTools/computeSurfaceNodeIDsInPolygonalRegion.cpp b/Applications/Utils/MeshGeoTools/computeSurfaceNodeIDsInPolygonalRegion.cpp
index a0e61a2d1c8ed504dc4ee1df9c6f18daf52e3c84..49ac89e93b28c39bc2741d94bf49d5c304b58f82 100644
--- a/Applications/Utils/MeshGeoTools/computeSurfaceNodeIDsInPolygonalRegion.cpp
+++ b/Applications/Utils/MeshGeoTools/computeSurfaceNodeIDsInPolygonalRegion.cpp
@@ -150,10 +150,10 @@ int main(int argc, char* argv[])
         std::vector<std::pair<std::size_t, double>> ids_and_areas;
         for (std::size_t k(0); k < all_sfc_nodes.size(); k++)
         {
-            MeshLib::Node const& pnt(*(all_sfc_nodes[k]));
-            if (polygon.isPntInPolygon(pnt[0], pnt[1], pnt[2]))
+            MeshLib::Node const& surface_node(*(all_sfc_nodes[k]));
+            if (polygon.isPntInPolygon(surface_node))
             {
-                ids_and_areas.emplace_back(pnt.getID(), areas[k]);
+                ids_and_areas.emplace_back(surface_node.getID(), areas[k]);
             }
         }
         if (ids_and_areas.empty())