Remove automatic_dimension_set call from remove and prune functions

Make it public to be available
Modify tests accordingly


git-svn-id: svn+ssh://scm.gforge.inria.fr/svnroot/gudhi/branches/ST_automatic_dimension_set@2703 636b058d-ea47-450e-bf9e-a15bfbe3eedb

Former-commit-id: 2ff42d5327ae76b121390a913c4e548d6f59ea99

1 files changed, 13 insertions, 23 deletions

diff --git a/src/Simplex_tree/include/gudhi/Simplex_tree.h b/src/Simplex_tree/include/gudhi/Simplex_tree.h
index 5d1885b8..fd6cd72a 100644
--- a/src/Simplex_tree/include/gudhi/Simplex_tree.h
+++ b/src/Simplex_tree/include/gudhi/Simplex_tree.h
@@ -1161,13 +1161,11 @@ class Simplex_tree {
    * \post Some simplex tree functions require the filtration to be valid. `prune_above_filtration()`
    * function is not launching `initialize_filtration()` but returns the filtration modification information. If the
    * complex has changed , please call `initialize_filtration()` to recompute it.
+   * \post Be aware that `prune_above_filtration()` may change the simplex tree dimension (`automatic_dimension_set()`
+   * to be done).
    */
   bool prune_above_filtration(Filtration_value filtration) {
-    bool modified = rec_prune_above_filtration(root(), filtration);
-    if (modified) {
-      auto_dimension_set(dimension());
-    }
-    return modified;
+    return rec_prune_above_filtration(root(), filtration);
   }
 
  private:
@@ -1195,17 +1193,16 @@ class Simplex_tree {
     return modified;
   }
 
- private:
-  /** \brief Resets the Simplex_tree dimension.
-   * @param[in] old_dimension The former dimension value until the loop stopped when it is reached.
+ public:
+  /** \brief Deep search simplex tree dimension reset.
    * @return The dimension modification information.
-   * \pre Please check the simplex has not a too low dimension value.
-   * This cannot happen if set_dimension has not been performed.
+   * \pre Be sure the simplex tree has not a too low dimension value as the deep search stops when the former dimension
+   * has been reached (cf. `dimension()` and `set_dimension()` methods).
    */
-  bool auto_dimension_set(int old_dimension) {
+  bool automatic_dimension_set() {
     int new_dimension = -1;
     // Browse the tree from te left to the right as higher dimension cells are more likely on the left part of the tree
-    for (Simplex_handle sh : skeleton_simplex_range(old_dimension)) {
+    for (Simplex_handle sh : skeleton_simplex_range(dimension_)) {
 #ifdef DEBUG_TRACES
       for (auto vertex : simplex_vertex_range(sh)) {
         std::cout << " " << vertex;
@@ -1214,12 +1211,12 @@ class Simplex_tree {
 #endif  // DEBUG_TRACES
 
       int sh_dimension = dimension(sh);
-      if (sh_dimension >= old_dimension)
+      if (sh_dimension >= dimension_)
         // Stop browsing as soon as the dimension is reached, no need to go furter
         return false;
       new_dimension = std::max<int>(new_dimension, sh_dimension);
     }
-    set_dimension(new_dimension);
+    dimension_ = new_dimension;
     return true;
   }
 
@@ -1229,7 +1226,8 @@ class Simplex_tree {
    * @param[in] sh Simplex handle on the maximal simplex to remove.
    * \pre Please check the simplex has no coface before removing it.
    * \exception std::invalid_argument In debug mode, if sh has children.
-   * \post Be aware that removing is shifting data in a flat_map (initialize_filtration to be done).
+   * \post Be aware that removing is shifting data in a flat_map (`initialize_filtration()` to be done).
+   * \post Be aware that removing may change the simplex tree dimension (`automatic_dimension_set()` to be done).
    */
   void remove_maximal_simplex(Simplex_handle sh) {
     // Guarantee the simplex has no children
@@ -1244,17 +1242,9 @@ class Simplex_tree {
       // Special case when child is the root of the simplex tree, just remove it from members
       child->erase(sh);
     } else {
-      // Keep information before remove action
-      int sh_dim = dimension(sh);
-
       // Sibling is emptied : must be deleted, and its parent must point on his own Sibling
       child->oncles()->members().at(child->parent()).assign_children(child->oncles());
       delete child;
-
-      // No need to reset dimension in case maximal simplex is not the maximal dimension one
-      if (sh_dim >= dimension()) {
-        auto_dimension_set(sh_dim);
-      }
     }
   }