From 8521d24c4584a6eb1149a604309ffa987862414f Mon Sep 17 00:00:00 2001
From: Ulrich Bauer <mail@ulrich-bauer.org>
Date: Mon, 28 Nov 2016 14:13:46 -0500
Subject: allow for negative edge filtration values

---
 ripser.cpp | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/ripser.cpp b/ripser.cpp
index 72e784c..efa672e 100644
--- a/ripser.cpp
+++ b/ripser.cpp
@@ -403,7 +403,7 @@ public:
 	}
 	
 	value_t compute_diameter(const index_t index, index_t dim) const {
-		value_t diam = 0;
+		value_t diam = -std::numeric_limits<value_t>::infinity();
 
 		vertices.clear();
 		get_simplex_vertices(index, dim, dist.size(), std::back_inserter(vertices));
@@ -682,7 +682,7 @@ public:
 
 				if (u != v) {
 #ifdef PRINT_PERSISTENCE_PAIRS
-					if (get_diameter(e) > 0) std::cout << " [0," << get_diameter(e) << ")" << std::endl;
+					if (get_diameter(e) != 0) std::cout << " [0," << get_diameter(e) << ")" << std::endl;
 #endif
 					dset.link(u, v);
 				} else
-- 
cgit v1.2.3