diff --git a/src/Alpha_complex/test/Delaunay_complex_unit_test.h b/src/Alpha_complex/test/Delaunay_complex_unit_test.h
index e58cc1e4e..8870c6b97 100644
--- a/src/Alpha_complex/test/Delaunay_complex_unit_test.h
+++ b/src/Alpha_complex/test/Delaunay_complex_unit_test.h
@@ -24,7 +24,7 @@ using Simplex_handle = Simplex_tree::Simplex_handle;
 
 template<class CGAL_kernel>
 void compare_delaunay_complex_simplices() {
-  std::clog << "*****************************************************************************************************\n";
+  std::clog << "***************************************************************************************************\n";
   using Point = typename CGAL_kernel::Point_d;
   std::vector<Point> points;
   // 50 points on a 4-sphere
@@ -55,8 +55,14 @@ void compare_delaunay_complex_simplices() {
   std::clog << "Check simplices from alpha complex filtration values when output_squared_values is true\n";
   // Check all the simplices from alpha complex are in the Delaunay complex
   for (auto f_simplex : stree_from_alpha_complex.complex_simplex_range()) {
+    std::clog << "[";
+    for (auto vertex : stree_from_alpha_complex.simplex_vertex_range(f_simplex)) std::clog << vertex << " ";
+    std::clog << "] - ";
+
     Simplex_handle sh = stree_from_alpha_sqrt.find(stree_from_alpha_complex.simplex_vertex_range(f_simplex));
     BOOST_CHECK(sh != stree_from_alpha_sqrt.null_simplex());
+    std::clog << "alpha_sqrt = " << stree_from_alpha_sqrt.filtration(sh);
+    std::clog << " vs. sqrt(alpha) = " << std::sqrt(stree_from_alpha_complex.filtration(f_simplex)) << "\n";
     GUDHI_TEST_FLOAT_EQUALITY_CHECK(stree_from_alpha_sqrt.filtration(sh),
                                     std::sqrt(stree_from_alpha_complex.filtration(f_simplex)));
   }