Fixing H2 sparse matrix threshold bug. Fixes #51

scikit-tda · Nov 13, 2018 · d6a68f4 · d6a68f4
1 parent dc9c50f
commit d6a68f4
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/ripser/ripser.cpp b/ripser/ripser.cpp
@@ -995,6 +995,8 @@ void ripser<sparse_distance_matrix>::assemble_columns_to_reduce(
 	}
 
 	simplices.swap(next_simplices);
+	std::sort(columns_to_reduce.begin(), columns_to_reduce.end(),
+				greater_diameter_or_smaller_index<diameter_index_t>());
 }
 
 ripserResults rips_dm(float* D, int N, int modulus, int dim_max, float threshold, int do_cocycles) {