diff options
author | vrouvrea <vrouvrea@636b058d-ea47-450e-bf9e-a15bfbe3eedb> | 2017-01-24 10:28:32 +0000 |
---|---|---|
committer | vrouvrea <vrouvrea@636b058d-ea47-450e-bf9e-a15bfbe3eedb> | 2017-01-24 10:28:32 +0000 |
commit | 15a0852d4e182eb7340f82c68823106f69107e27 (patch) | |
tree | 6fd37c0390fdeb76390732ce9b33f72dea87389f /src/Bottleneck_distance | |
parent | f3ed36c07929334b6071de0ce5630da0b793f047 (diff) |
Last modification
git-svn-id: svn+ssh://scm.gforge.inria.fr/svnroot/gudhi/branches/bottleneck_misc_fixes@1991 636b058d-ea47-450e-bf9e-a15bfbe3eedb
Former-commit-id: 1e7e28db5a3f79ba5f65f4a3c025d23703638362
Diffstat (limited to 'src/Bottleneck_distance')
-rw-r--r-- | src/Bottleneck_distance/utils/alpha_rips_persistence_bottleneck_distance.cpp | 24 |
1 files changed, 15 insertions, 9 deletions
diff --git a/src/Bottleneck_distance/utils/alpha_rips_persistence_bottleneck_distance.cpp b/src/Bottleneck_distance/utils/alpha_rips_persistence_bottleneck_distance.cpp index 9043c947..6e9be97b 100644 --- a/src/Bottleneck_distance/utils/alpha_rips_persistence_bottleneck_distance.cpp +++ b/src/Bottleneck_distance/utils/alpha_rips_persistence_bottleneck_distance.cpp @@ -71,9 +71,9 @@ int main(int argc, char * argv[]) { Points_off_reader off_reader(off_file_points); - // ------------------------------------ + // -------------------------------------------- // Rips persistence - // ------------------------------------ + // -------------------------------------------- Rips_complex rips_complex(off_reader.get_point_cloud(), threshold, Euclidean_distance()); // Construct the Rips complex in a Simplex Tree @@ -90,14 +90,13 @@ int main(int argc, char * argv[]) { Persistent_cohomology rips_pcoh(rips_stree); // initializes the coefficient field for homology rips_pcoh.init_coefficients(p); - rips_pcoh.compute_persistent_cohomology(min_persistence); - rips_pcoh.output_diagram(); + // rips_pcoh.output_diagram(); - // ------------------------------------ + // -------------------------------------------- // Alpha persistence - // ------------------------------------ + // -------------------------------------------- Gudhi::alpha_complex::Alpha_complex<Kernel> alpha_complex(off_reader.get_point_cloud()); Simplex_tree alpha_stree; @@ -112,12 +111,15 @@ int main(int argc, char * argv[]) { Persistent_cohomology alpha_pcoh(alpha_stree); // initializes the coefficient field for homology alpha_pcoh.init_coefficients(p); - alpha_pcoh.compute_persistent_cohomology(min_persistence * min_persistence); - alpha_pcoh.output_diagram(); + // alpha_pcoh.output_diagram(); - for (int dim = 0; dim <= dim_max; dim ++) { + // -------------------------------------------- + // Bottleneck distance between both persistence + // -------------------------------------------- + double max_b_distance {}; + for (int dim = 0; dim < dim_max; dim ++) { std::vector< std::pair< Filtration_value , Filtration_value > > rips_intervals; std::vector< std::pair< Filtration_value , Filtration_value > > alpha_intervals; rips_intervals = rips_pcoh.intervals_in_dimension(dim); @@ -126,7 +128,11 @@ int main(int argc, char * argv[]) { double bottleneck_distance = Gudhi::persistence_diagram::bottleneck_distance(rips_intervals, alpha_intervals); std::cout << "In dimension " << dim << ", bottleneck distance = " << bottleneck_distance << std::endl; + if (bottleneck_distance > max_b_distance) + max_b_distance = bottleneck_distance; } + std::cout << "================================================================================" << std::endl; + std::cout << "Bottleneck distance is " << max_b_distance << std::endl; return 0; } |