summaryrefslogtreecommitdiff
path: root/ripser.cpp
diff options
context:
space:
mode:
authorUlrich Bauer <ulrich.bauer@tum.de>2015-11-07 11:36:17 -0500
committerUlrich Bauer <ulrich.bauer@tum.de>2015-11-07 11:36:17 -0500
commit85c303ab2454ac17d74967fbfae3e9ee842165a2 (patch)
tree31874324c159cb9b8d911e1b457b8a19ead99913 /ripser.cpp
parentfc3ae90b8b0ed4de041867425edc16c5875f36ee (diff)
compare details
Diffstat (limited to 'ripser.cpp')
-rw-r--r--ripser.cpp14
1 files changed, 7 insertions, 7 deletions
diff --git a/ripser.cpp b/ripser.cpp
index 0f23a55..9f11c58 100644
--- a/ripser.cpp
+++ b/ripser.cpp
@@ -781,7 +781,7 @@ void compute_pairs(
entry_t pivot = make_entry(column_to_reduce, modulus - 1);
- std::cout << "reducing " << column_to_reduce << ": pivot " << std::flush;
+ std::cout << "reducing " << column_to_reduce << ":" << std::endl;
#ifdef ASSEMBLE_REDUCTION_MATRIX
reduction_matrix.append();
@@ -802,9 +802,9 @@ void compute_pairs(
const coefficient_t factor = modulus - get_coefficient(pivot);
-// std::priority_queue<entry_diameter_t, std::vector<entry_diameter_t>, decltype(comp) > eliminating_coboundary(comp);
+ std::priority_queue<entry_diameter_t, std::vector<entry_diameter_t>, decltype(comp) > eliminating_coboundary(comp);
-// std::cout << "w:" << get_column_vector(working_coboundary, modulus) << std::endl;
+ std::cout << "w:" << get_column_vector(working_coboundary, modulus) << std::endl;
#ifdef ASSEMBLE_REDUCTION_MATRIX
for (auto it = reduction_matrix.cbegin(j); it != reduction_matrix.cend(j); ++it)
@@ -844,7 +844,7 @@ void compute_pairs(
assert(coface_coefficient >= 0);
push_entry(working_coboundary, coface_index, coface_coefficient, coface_diameter);
-// push_entry(eliminating_coboundary, coface_index, coface_coefficient, coface_diameter);
+ push_entry(eliminating_coboundary, coface_index, coface_coefficient, coface_diameter);
}
}
}
@@ -854,12 +854,12 @@ void compute_pairs(
// std::cout << get_heap_vector(working_coboundary) << std::endl;
-// std::cout << "e:" << get_column_vector(eliminating_coboundary, modulus) << std::endl;
-// std::cout << "w:" << get_column_vector(working_coboundary, modulus) << std::endl << std::endl;
+ std::cout << "e:" << get_column_vector(eliminating_coboundary, modulus) << std::endl;
+ std::cout << "w:" << get_column_vector(working_coboundary, modulus) << std::endl << std::endl;
pivot = get_pivot(working_coboundary, modulus);
- std::cout << get_index(pivot) << " ";
+ std::cout << "pivot " << get_index(pivot) << std::endl;
if (get_index(pivot) != -1) {
auto pair = pivot_column_index.find(get_index(pivot));