No more use of threshold_ for prune_above_filtration because threshold may be not up to date

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

Former-commit-id: b17d8ef63bc173b4c4b6b0bde31f97d85672a754

1 files changed, 22 insertions, 26 deletions

diff --git a/src/Simplex_tree/include/gudhi/Simplex_tree.h b/src/Simplex_tree/include/gudhi/Simplex_tree.h
index 813e7924..92f4576d 100644
--- a/src/Simplex_tree/include/gudhi/Simplex_tree.h
+++ b/src/Simplex_tree/include/gudhi/Simplex_tree.h
@@ -1174,7 +1174,6 @@ class Simplex_tree {
  public:
   /** \brief Prune above filtration value given as parameter.
    * @param[in] filtration Maximum threshold value.
-   * \warning threshold_ is set from filtration given as parameter.
    * \warning The filtration must be valid. If the filtration has not been initialized yet, the method initializes it
    * (i.e. order the simplices). If the complex has changed since the last time the filtration was initialized, please
    * call `initialize_filtration()` to recompute it.
@@ -1182,33 +1181,30 @@ class Simplex_tree {
   void prune_above_filtration(Filtration_value filtration) {
     // No action if filtration is not stored
     if (Options::store_filtration) {
-      if (filtration < threshold_) {
-        threshold_ = filtration;
-        // Initialize filtration_vect_ if required
-        if (filtration_vect_.empty()) {
-          initialize_filtration();
-        }
+      // Initialize filtration_vect_ if required
+      if (filtration_vect_.empty()) {
+        initialize_filtration();
+      }
         
-        std::vector<std::vector<Vertex_handle>> simplex_list_to_removed;
-        // Loop in reverse mode until threshold is reached
-        // Do not erase while looping, because removing is shifting data in a flat_map
-        for (auto f_simplex = filtration_vect_.rbegin();
-             (f_simplex != filtration_vect_.rend()) && ((*f_simplex)->second.filtration() > threshold_);
-             f_simplex++) {
-          std::vector<Vertex_handle> simplex_to_remove;
-          for (auto vertex : simplex_vertex_range(*f_simplex))
-            simplex_to_remove.insert(simplex_to_remove.begin(), vertex);
-          simplex_list_to_removed.push_back(simplex_to_remove);
-        }
-        for (auto simplex_to_remove : simplex_list_to_removed) {
-          Simplex_handle sh = find_simplex(simplex_to_remove);
-          if (sh != null_simplex())
-            remove_maximal_simplex(sh);
-        }
-        // Re-initialize filtration_vect_ if dta were removed, because removing is shifting data in a flat_map
-        if (simplex_list_to_removed.size() > 0)
-          initialize_filtration();
+      std::vector<std::vector<Vertex_handle>> simplex_list_to_removed;
+      // Loop in reverse mode until threshold is reached
+      // Do not erase while looping, because removing is shifting data in a flat_map
+      for (auto f_simplex = filtration_vect_.rbegin();
+           (f_simplex != filtration_vect_.rend()) && ((*f_simplex)->second.filtration() > filtration);
+           f_simplex++) {
+        std::vector<Vertex_handle> simplex_to_remove;
+        for (auto vertex : simplex_vertex_range(*f_simplex))
+          simplex_to_remove.insert(simplex_to_remove.begin(), vertex);
+        simplex_list_to_removed.push_back(simplex_to_remove);
+      }
+      for (auto simplex_to_remove : simplex_list_to_removed) {
+        Simplex_handle sh = find_simplex(simplex_to_remove);
+        if (sh != null_simplex())
+          remove_maximal_simplex(sh);
       }
+      // Re-initialize filtration_vect_ if dta were removed, because removing is shifting data in a flat_map
+      if (simplex_list_to_removed.size() > 0)
+        initialize_filtration();
     }
   }