summaryrefslogtreecommitdiff
path: root/ripser.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'ripser.cpp')
-rw-r--r--ripser.cpp6
1 files changed, 4 insertions, 2 deletions
diff --git a/ripser.cpp b/ripser.cpp
index 55e94b5..3f19858 100644
--- a/ripser.cpp
+++ b/ripser.cpp
@@ -577,6 +577,8 @@ public:
<< "\r" << std::flush;
#endif
+ reduction_matrix.append_column();
+
std::priority_queue<diameter_entry_t, std::vector<diameter_entry_t>,
greater_diameter_or_smaller_index<diameter_entry_t>>
working_reduction_column, working_coboundary;
@@ -618,8 +620,8 @@ public:
pivot_column_index.insert(
std::make_pair(get_entry(pivot), index_column_to_reduce));
- reduction_matrix.append_column();
-
+ const coefficient_t inverse =
+ multiplicative_inverse[get_coefficient(pivot)];
while (true) {
diameter_entry_t e = pop_pivot(working_reduction_column, modulus);
if (get_index(e) == -1) break;