diff --git a/Principal_component_analysis/test/Principal_component_analysis/fitting_3.C b/Principal_component_analysis/test/Principal_component_analysis/fitting_3.C
index 06a0b69c9e0..9c445c8f6ea 100644
--- a/Principal_component_analysis/test/Principal_component_analysis/fitting_3.C
+++ b/Principal_component_analysis/test/Principal_component_analysis/fitting_3.C
@@ -98,9 +98,9 @@ void test_3D_triangle_set(const unsigned int nb_triangles)
   FT quality;
   Plane plane;
   Point centroid;
-  quality = linear_least_squares_fitting_3(points.begin(),points.end(),plane);
-  quality = linear_least_squares_fitting_3(points.begin(),points.end(),plane,centroid);
-  quality = linear_least_squares_fitting_3(points.begin(),points.end(),plane,centroid,k);
+  quality = linear_least_squares_fitting_3(triangles.begin(),triangles.end(),plane);
+  quality = linear_least_squares_fitting_3(triangles.begin(),triangles.end(),plane,centroid);
+  quality = linear_least_squares_fitting_3(triangles.begin(),triangles.end(),plane,centroid,k);
   std::cout << "done (quality: " << quality << ")" << std::endl;
 
   Plane horizontal_plane(Point(0,0,0),Vector(0,0,1));