From 5694670b3e20f0cb935a751614ef12b6009a60c0 Mon Sep 17 00:00:00 2001
From: mathieu <mathieu.carriere3@gmail.com>
Date: Thu, 16 Jan 2020 15:58:15 -0500
Subject: fix to detect infinite persistence

---
 src/python/include/Persistent_cohomology_interface.h | 5 ++---
 1 file changed, 2 insertions(+), 3 deletions(-)

(limited to 'src')

diff --git a/src/python/include/Persistent_cohomology_interface.h b/src/python/include/Persistent_cohomology_interface.h
index 774eb56a..acc32b21 100644
--- a/src/python/include/Persistent_cohomology_interface.h
+++ b/src/python/include/Persistent_cohomology_interface.h
@@ -124,16 +124,15 @@ persistent_cohomology::Persistent_cohomology<FilteredComplex, persistent_cohomol
     std::vector<std::pair<int, std::pair<std::pair<double, int>, std::pair<double, int>>>> persistence;
     for (auto pair : persistent_pairs) {
 
-      int splx0, splx1;
-
       double f0 = stptr_->filtration(get<0>(pair));
       // Recursively get the top-dimensional cells / cofaces associated to the persistence generator
       std::vector<int> faces0; top_dimensional_cofaces(faces0, stptr_->key(get<0>(pair)));
       // Find the top-dimensional cell / coface with the same filtration value
       int cf; for (int i = 0; i < faces0.size(); i++){ if (stptr_->filtration(faces0[i]) == f0){cf = i; break;}}
       // Retrieve the index of the corresponding top-dimensional cell in the input data
-      splx0 = order[faces0[cf]];
+      int splx0 = order[faces0[cf]];
 
+      int splx1 = -1;
       if (isfinite(stptr_->filtration(get<1>(pair)))){
       double f1 = stptr_->filtration(get<1>(pair));
       // Recursively get the top-dimensional cells / cofaces associated to the persistence generator
-- 
cgit v1.2.3