diff --git a/Point_set_shape_detection_3/test/Point_set_shape_detection_3/test_scene.cpp b/Point_set_shape_detection_3/test/Point_set_shape_detection_3/test_scene.cpp
index 07fe4c9cceb..5282167816a 100644
--- a/Point_set_shape_detection_3/test/Point_set_shape_detection_3/test_scene.cpp
+++ b/Point_set_shape_detection_3/test/Point_set_shape_detection_3/test_scene.cpp
@@ -5,6 +5,7 @@
 #include <CGAL/Simple_cartesian.h>
 
 #include <CGAL/Shape_detection_3.h>
+#include <CGAL/Plane_regularization.h>
 #include <CGAL/Point_with_normal_3.h>
 #include <CGAL/property_map.h>
 
@@ -31,6 +32,8 @@ bool test_scene() {
   typedef CGAL::Shape_detection_3::Sphere<Traits>             Sphere;
   typedef CGAL::Shape_detection_3::Torus<Traits>              Torus;
 
+  typedef CGAL::Plane_regularization<Traits>                  Regularization;
+  
   Pwn_vector points;
 
   // Loads point set from a file. 
@@ -126,6 +129,10 @@ bool test_scene() {
     return false;
   }
 
+  // Test regularization
+  Regularization regularization (points, ransac);
+  regularization.run (50., 0.01);
+  
   Point_index_range pts = ransac.indices_of_unassigned_points();
 
   std::cout << " succeeded" << std::endl;