diff --git a/Poisson_surface_reconstruction_3/include/CGAL/poisson_surface_reconstruction.h b/Poisson_surface_reconstruction_3/include/CGAL/poisson_surface_reconstruction.h
index 0972f5ca58f..1f43dd92d4f 100644
--- a/Poisson_surface_reconstruction_3/include/CGAL/poisson_surface_reconstruction.h
+++ b/Poisson_surface_reconstruction_3/include/CGAL/poisson_surface_reconstruction.h
@@ -35,9 +35,14 @@ namespace CGAL {
 
     Performs surface reconstruction as follows:
 
-    - compute the Poisson implicit function
-    - meshes the function with a user-defined precision
-    - outputs the result in a polygon mesh
+    - compute the Poisson implicit function, through a conjugate
+      gradient solver, represented as a piecewise linear function
+      stored on a 3D Delaunay mesh generated via Delaunay refinement
+    - meshes the function with a user-defined precision using another
+      round of Delaunay refinement: it contours the isosurface
+      corresponding to the isovalue of the median of the function
+      values at the input points
+    - outputs the result in a polygon mesh 
 
     This function relies mainly on the size parameter `spacing`. A
     reasonable solution is to use the average spacing of the input