From 6723f7687a27ac4e48a5d46bdac9cd4aa2478830 Mon Sep 17 00:00:00 2001 From: "jan.reininghaus" Date: Tue, 17 Dec 2013 10:17:07 +0000 Subject: new representation type heap_pivot_column git-svn-id: https://phat.googlecode.com/svn/trunk@147 8e3bb3c2-eed4-f18f-5264-0b6c94e6926d --- src/self_test.cpp | 15 +++++++++++++-- 1 file changed, 13 insertions(+), 2 deletions(-) (limited to 'src/self_test.cpp') diff --git a/src/self_test.cpp b/src/self_test.cpp index cfecd0e..21bb3bc 100644 --- a/src/self_test.cpp +++ b/src/self_test.cpp @@ -22,6 +22,7 @@ #include #include #include +#include #include #include @@ -36,6 +37,7 @@ int main( int argc, char** argv ) std::string test_data = argc > 1 ? argv[ 1 ] : "examples/torus.bin"; typedef phat::sparse_pivot_column Sparse; + typedef phat::heap_pivot_column Heap; typedef phat::full_pivot_column Full; typedef phat::bit_tree_pivot_column BitTree; typedef phat::vector_vector Vec_vec; @@ -57,6 +59,11 @@ int main( int argc, char** argv ) phat::boundary_matrix< Sparse > sparse_boundary_matrix = boundary_matrix; phat::compute_persistence_pairs< phat::chunk_reduction >( sparse_pairs, sparse_boundary_matrix ); + std::cout << "Running Chunk - Heap ..." << std::endl; + phat::persistence_pairs heap_pairs; + phat::boundary_matrix< Heap > heap_boundary_matrix = boundary_matrix; + phat::compute_persistence_pairs< phat::chunk_reduction >( heap_pairs, heap_boundary_matrix ); + std::cout << "Running Chunk - Full ..." << std::endl; phat::persistence_pairs full_pairs; phat::boundary_matrix< Full > full_boundary_matrix = boundary_matrix; @@ -82,8 +89,12 @@ int main( int argc, char** argv ) phat::boundary_matrix< Vec_list > vec_list_boundary_matrix = boundary_matrix; phat::compute_persistence_pairs< phat::chunk_reduction >( vec_list_pairs, vec_list_boundary_matrix ); - if( sparse_pairs != full_pairs ) { - std::cerr << "Error: sparse and full differ!" << std::endl; + if( sparse_pairs != heap_pairs ) { + std::cerr << "Error: sparse and heap differ!" << std::endl; + error = true; + } + if( heap_pairs != full_pairs ) { + std::cerr << "Error: heap and full differ!" << std::endl; error = true; } if( full_pairs != vec_vec_pairs ) { -- cgit v1.2.3