message: Make it clearer that create_complex can/should be called without a second argument.
  Move the test restricting propagation to dim>1 earlier.
diff: @@ -228,11 +228,6 @@ class Alpha_complex {
  }
  public:
  - template <typename SimplicialComplexForAlpha>
  - bool create_complex(SimplicialComplexForAlpha& complex) {
  - typedef typename SimplicialComplexForAlpha::Filtration_value Filtration_value;
  - return create_complex(complex, std::numeric_limits<Filtration_value>::infinity());
  - }
  /** \brief Inserts all Delaunay triangulation into the simplicial complex.
  * It also computes the filtration values accordingly to the \ref createcomplexalgorithm
  @@ -240,7 +235,7 @@ class Alpha_complex {
  * \tparam SimplicialComplexForAlpha must meet `SimplicialComplexForAlpha` concept.
  *
  * @param[in] complex SimplicialComplexForAlpha to be created.
  - * @param[in] max_alpha_square maximum for alpha square value. Default value is +\f$\infty\f$.
  + * @param[in] max_alpha_square maximum for alpha square value. Default value is +\f$\infty\f$, and there is very little point using anything else since it does not save time.
  *
  * @return true if creation succeeds, false otherwise.
  *
  @@ -249,8 +244,8 @@ class Alpha_complex {
  *
  * Initialization can be launched once.
  */
  - template <typename SimplicialComplexForAlpha, typename Filtration_value>
  - bool create_complex(SimplicialComplexForAlpha& complex, Filtration_value max_alpha_square) {
  + template <typename SimplicialComplexForAlpha, typename Filtration_value = typename SimplicialComplexForAlpha::Filtration_value>
  + bool create_complex(SimplicialComplexForAlpha& complex, Filtration_value max_alpha_square = std::numeric_limits<Filtration_value>::infinity()) {
  // From SimplicialComplexForAlpha type required to insert into a simplicial complex (with or without subfaces).
  typedef typename SimplicialComplexForAlpha::Vertex_handle Vertex_handle;
  typedef typename SimplicialComplexForAlpha::Simplex_handle Simplex_handle;
  @@ -333,6 +328,8 @@ class Alpha_complex {
  std::cout << "filt(Sigma) is NaN : filt(Sigma) =" << complex.filtration(f_simplex) << std::endl;
  #endif // DEBUG_TRACES
  }
  + // No need to propagate further, unweighted points all have value 0
  + if (decr_dim > 1)
  propagate_alpha_filtration(complex, f_simplex, decr_dim);
  }
  }
  @@ -379,9 +376,6 @@ class Alpha_complex {
  #endif // DEBUG_TRACES
  // ### Else
  } else {
  - // No need to compute is_gabriel for dimension <= 2
  - // i.e. : Sigma = (3,1) => Tau = 1
  - if (decr_dim > 1) {
  // insert the Tau points in a vector for is_gabriel function
  Vector_of_CGAL_points pointVector;
  #ifdef DEBUG_TRACES
  @@ -422,7 +416,6 @@ class Alpha_complex {
  }
  }
  }
  - }
  };
  } // namespace alpha_complex
  