diff --git a/Geo/GVertex.cpp b/Geo/GVertex.cpp
index 3bf0d3fa2944045360a06e2f7caf27053b10afef..4eb8f03ed355b7e8a11a346211908a50c348d8f6 100644
--- a/Geo/GVertex.cpp
+++ b/Geo/GVertex.cpp
@@ -73,3 +73,12 @@ MElement *GVertex::getMeshElement(unsigned int index)
     return points[index]; 
   return 0;
 }
+
+bool GVertex::isOnSeam(const GFace *gf) const
+{
+  std::list<GEdge*>::const_iterator eIter = l_edges.begin();
+  for (; eIter != l_edges.end(); eIter++) {
+    if ( (*eIter)->isSeam(gf) ) return true;
+  }
+  return false;
+}
diff --git a/Geo/GVertex.h b/Geo/GVertex.h
index 27db6204b6d77e2e06cc3810ea953fc0e582e3d3..65cd4ddf679a7ad3bac5982c8f1de49c9f6512d3 100644
--- a/Geo/GVertex.h
+++ b/Geo/GVertex.h
@@ -71,6 +71,9 @@ class GVertex : public GEntity
   // get the element at the given index
   MElement *getMeshElement(unsigned int index);
 
+  // return true if this vertex is on a seam of the given face
+  bool isOnSeam(const GFace *gf) const;
+
   std::vector<MPoint*> points;
 };