diff --git a/Isosurfacing_3/include/CGAL/Isosurfacing_3/internal/dual_contouring_functors.h b/Isosurfacing_3/include/CGAL/Isosurfacing_3/internal/dual_contouring_functors.h
index 1f2fc7e612a..f0a94195c5a 100644
--- a/Isosurfacing_3/include/CGAL/Isosurfacing_3/internal/dual_contouring_functors.h
+++ b/Isosurfacing_3/include/CGAL/Isosurfacing_3/internal/dual_contouring_functors.h
@@ -180,9 +180,9 @@ bool cell_position_QEM(const typename Domain::Cell_descriptor& c,
   v_svd = x_hat + svd.solve(rhs - A * x_hat);
 
 
-  // bounding box
   if(constrain_to_cell)
   {
+    // @todo clamping back doesn't necessarily yield the optimal position within the cell
     v_svd[0] = std::clamp<FT>(v_svd[0], x_min, x_max);
     v_svd[1] = std::clamp<FT>(v_svd[1], y_min, y_max);
     v_svd[2] = std::clamp<FT>(v_svd[2], z_min, z_max);