diff --git a/GeoLib/PointVec.cpp b/GeoLib/PointVec.cpp
index dab41d47257a30ac8089e80339a1fdb13ab14219..ff6514d04d0797ab82679a3273859cf136fe99a5 100644
--- a/GeoLib/PointVec.cpp
+++ b/GeoLib/PointVec.cpp
@@ -207,14 +207,4 @@ void PointVec::correctNameIDMapping()
 	}
 }
 
-std::vector<GeoLib::Point*>* PointVec::getSubset(const std::vector<std::size_t> &subset) const
-{
-	std::vector<GeoLib::Point*> *new_points (new std::vector<GeoLib::Point*>(subset.size()));
-
-	const std::size_t nPoints(subset.size());
-	for (std::size_t i = 0; i < nPoints; i++)
-		(*new_points)[i] = new GeoLib::Point(*(*this->_data_vec)[subset[i]]);
-
-	return new_points;
-}
 } // end namespace
diff --git a/GeoLib/PointVec.h b/GeoLib/PointVec.h
index d768d2528507da77d51bcfcd8618acf1ada19b3b..3fbed3af5329984727b70f1eaf08ec3a32694179 100644
--- a/GeoLib/PointVec.h
+++ b/GeoLib/PointVec.h
@@ -100,9 +100,6 @@ public:
 
 	const GeoLib::AABB<GeoLib::Point>& getAABB () const;
 
-	/// Returns a subset of this point vector containing only the points specified in "subset"
-	std::vector<GeoLib::Point*>* getSubset(const std::vector<std::size_t> &subset) const;
-
 private:
 	/**
 	 * Removes points out of the given point set that have the (nearly) same coordinates, i.e.