diff --git a/Geo/GFace.cpp b/Geo/GFace.cpp
index f4a558cf9fa04ba73cbfae323a1a860431ea1c2f..66c0d6dbe3427a9519eb5fd72484bed5c54a0118 100644
--- a/Geo/GFace.cpp
+++ b/Geo/GFace.cpp
@@ -522,6 +522,8 @@ void GFace::computeMeanPlane(const std::vector<MVertex*> &points)
 
 void GFace::computeMeanPlane(const std::vector<SPoint3> &points)
 {
+  if(points.empty()) return;
+
   // The concept of a mean plane computed in the sense of least
   // squares is fine for plane surfaces(!), but not really the best
   // one for non-plane surfaces. Indeed, imagine a quarter of a circle