diff --git a/Surface_modeling/include/CGAL/Deformation_Eigen_polar_closest_rotation_traits_3.h b/Surface_modeling/include/CGAL/Deformation_Eigen_polar_closest_rotation_traits_3.h
index 8d28dc5d279..66b250b3042 100644
--- a/Surface_modeling/include/CGAL/Deformation_Eigen_polar_closest_rotation_traits_3.h
+++ b/Surface_modeling/include/CGAL/Deformation_Eigen_polar_closest_rotation_traits_3.h
@@ -61,9 +61,9 @@ namespace CGAL {
       const Scalar th = std::sqrt(Eigen::NumTraits<Scalar>::dummy_precision());
 
       Eigen::SelfAdjointEigenSolver<Matrix> eig;
-      feclearexcept(FE_UNDERFLOW);
+      CGAL::feclearexcept(FE_UNDERFLOW);
       eig.computeDirect(A.transpose()*A);
-      if(fetestexcept(FE_UNDERFLOW) || eig.eigenvalues()(0)/eig.eigenvalues()(2)<th)
+      if(CGAL::fetestexcept(FE_UNDERFLOW) || eig.eigenvalues()(0)/eig.eigenvalues()(2)<th)
       { return false; }
 
       Vector S = eig.eigenvalues().cwiseSqrt();