From 0191e039a08ccd68b9916fafb2da98260e5b4efb Mon Sep 17 00:00:00 2001
From: Marc Glisse <marc.glisse@inria.fr>
Date: Sat, 13 Apr 2019 00:08:23 +0200
Subject: Handle triangular matrix

---
 src/Rips_complex/include/gudhi/Sparse_rips_complex.h | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

(limited to 'src/Rips_complex/include')

diff --git a/src/Rips_complex/include/gudhi/Sparse_rips_complex.h b/src/Rips_complex/include/gudhi/Sparse_rips_complex.h
index 015c94d9..3aafc5ff 100644
--- a/src/Rips_complex/include/gudhi/Sparse_rips_complex.h
+++ b/src/Rips_complex/include/gudhi/Sparse_rips_complex.h
@@ -94,7 +94,7 @@ class Sparse_rips_complex {
   template <typename DistanceMatrix>
   Sparse_rips_complex(const DistanceMatrix& distance_matrix, double epsilon)
       : Sparse_rips_complex(boost::irange<Vertex_handle>(0, boost::size(distance_matrix)),
-                            [&](Vertex_handle i, Vertex_handle j) { return distance_matrix[j][i]; }, epsilon) {}
+                            [&](Vertex_handle i, Vertex_handle j) { return (i==j) ? 0 : (i<j) ? distance_matrix[j][i] : distance_matrix[i][j]; }, epsilon) {}
 
   /** \brief Fills the simplicial complex with the sparse Rips graph and
    * expands it with all the cliques, stopping at a given maximal dimension.
-- 
cgit v1.2.3