Ignore zero persistence in the input

1 files changed, 7 insertions, 1 deletions

diff --git a/geom_bottleneck/bottleneck/src/bottleneck.cpp b/geom_bottleneck/bottleneck/src/bottleneck.cpp
index 365995a..34f7ed4 100644
--- a/geom_bottleneck/bottleneck/src/bottleneck.cpp
+++ b/geom_bottleneck/bottleneck/src/bottleneck.cpp
@@ -538,7 +538,13 @@ bool readDiagramPointSet(const char* fname, std::vector<std::pair<double, double
 #endif
             return false;
         }
-        result.push_back(std::make_pair(x,y));
+        if ( x != y ) {
+            result.push_back(std::make_pair(x,y));
+        } else {
+#ifndef FOR_R_TDA
+            std::cerr << "Warning: in file " << fname << ", line number " << lineNumber << ", zero persistence point: \"" << line << "\"" << std::endl;
+#endif
+        }
     }
     f.close();
     return true;