diff options
| author | Marc Glisse <marc.glisse@inria.fr> | 2019-04-13 00:08:23 +0200 |
|---|---|---|
| committer | Marc Glisse <marc.glisse@inria.fr> | 2019-04-13 00:08:23 +0200 |
| commit | 0191e039a08ccd68b9916fafb2da98260e5b4efb (patch) | |
| tree | fde094e496e65102ed1280fbb74bdf04a266f395 /src | |
| parent | 9ceaa84d52d939a117d78b49fd19c8900387dadc (diff) | |
Handle triangular matrix
Diffstat (limited to 'src')
| -rw-r--r-- | src/Rips_complex/include/gudhi/Sparse_rips_complex.h | 2 |
1 files changed, 1 insertions, 1 deletions
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. |