diff --git a/Principal_component_analysis/include/CGAL/linear_least_squares_fitting_3.h b/Principal_component_analysis/include/CGAL/linear_least_squares_fitting_3.h
index 60d136dbfa4..b607b83cdf5 100644
--- a/Principal_component_analysis/include/CGAL/linear_least_squares_fitting_3.h
+++ b/Principal_component_analysis/include/CGAL/linear_least_squares_fitting_3.h
@@ -32,12 +32,7 @@
 #include <CGAL/linear_least_squares_fitting_tetrahedra_3.h>
 #include <CGAL/linear_least_squares_fitting_spheres_3.h>
 
-#ifdef CGAL_EIGEN3_ENABLED
-#include <CGAL/Eigen_diagonalize_traits.h>
-#else
-#include <CGAL/Diagonalize_traits.h>
-#endif
-
+#include <CGAL/Default_diagonalize_traits.h>
 
 #include <CGAL/Dimension.h>
 
@@ -85,11 +80,7 @@ linear_least_squares_fitting_3(InputIterator first,
   typedef typename std::iterator_traits<InputIterator>::value_type Value_type;
   typedef typename Kernel_traits<Value_type>::Kernel Kernel;
   return CGAL::linear_least_squares_fitting_3(first,beyond,object,centroid,tag,Kernel(),
-#ifdef CGAL_EIGEN3_ENABLED
-					      Eigen_diagonalize_traits<typename Kernel::FT, 3>());
-#else
-  					      Diagonalize_traits<typename Kernel::FT, 3>());
-#endif
+					      Default_diagonalize_traits<typename Kernel::FT, 3>());
 
 }
 
@@ -109,11 +100,7 @@ linear_least_squares_fitting_3(InputIterator first,
   typedef typename Kernel_traits<Value_type>::Kernel Kernel;
   typename Kernel::Point_3 centroid; // not used by caller
   return CGAL::linear_least_squares_fitting_3(first,beyond,object,centroid,tag,Kernel(),
-#ifdef CGAL_EIGEN3_ENABLED
-					      Eigen_diagonalize_traits<typename Kernel::FT, 3>());
-#else
-  					      Diagonalize_traits<typename Kernel::FT, 3>());
-#endif
+					      Default_diagonalize_traits<typename Kernel::FT, 3>());
 
 }