don't assemble columns in top dimension

1 files changed, 9 insertions, 8 deletions

diff --git a/ripser.cpp b/ripser.cpp
index e79791a..f099244 100644
--- a/ripser.cpp
+++ b/ripser.cpp
@@ -1101,14 +1101,15 @@ int main( int argc, char** argv ) {
             binomial_coeff, multiplicative_inverse
         );
         
-        assemble_columns_to_reduce(
-            columns_to_reduce,
-            pivot_column_index,
-            comp,
-            dim, n,
-            threshold,
-            binomial_coeff
-        );
+        if (dim < dim_max)
+            assemble_columns_to_reduce(
+                columns_to_reduce,
+                pivot_column_index,
+                comp,
+                dim, n,
+                threshold,
+                binomial_coeff
+            );
         
 //        if ( dim > 1 )
 //            diameters[dim] = std::vector<value_t>();