diff --git a/vcg/space/polygon3.h b/vcg/space/polygon3.h
index 031523df..069830c4 100644
--- a/vcg/space/polygon3.h
+++ b/vcg/space/polygon3.h
@@ -538,12 +538,12 @@ typename PolygonType::ScalarType PolyAspectRatio(const PolygonType &F,
 template<class PolygonType>
 typename PolygonType::ScalarType PolygonPointDistance(const PolygonType &F,
                                                       const vcg::Point3<typename PolygonType::ScalarType> &pos,
-                                                      vcg::Point3<typename PolygonType::ScalarType> &ClosestP)
+                                                      vcg::Point3<typename PolygonType::ScalarType> &ClosestP,
+                                                      typename PolygonType::ScalarType minD = std::numeric_limits<typename PolygonType::ScalarType>::max())
 {
     typedef typename PolygonType::ScalarType ScalarType;
     typedef typename PolygonType::CoordType CoordType;
 
-    ScalarType minD=std::numeric_limits<ScalarType>::max();
     CoordType bary=vcg::PolyBarycenter(F);
     for (size_t j=0;j<F.VN();j++)
     {