summaryrefslogtreecommitdiff
path: root/src/Bottleneck_distance/example/bottleneck_read_file_example.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'src/Bottleneck_distance/example/bottleneck_read_file_example.cpp')
-rw-r--r--src/Bottleneck_distance/example/bottleneck_read_file_example.cpp71
1 files changed, 32 insertions, 39 deletions
diff --git a/src/Bottleneck_distance/example/bottleneck_read_file_example.cpp b/src/Bottleneck_distance/example/bottleneck_read_file_example.cpp
index ceedccc5..4c74b66e 100644
--- a/src/Bottleneck_distance/example/bottleneck_read_file_example.cpp
+++ b/src/Bottleneck_distance/example/bottleneck_read_file_example.cpp
@@ -4,7 +4,7 @@
*
* Authors: Francois Godi, small modifications by Pawel Dlotko
*
- * Copyright (C) 2015 INRIA (France)
+ * Copyright (C) 2015 INRIA
*
* This program is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
@@ -28,49 +28,42 @@
#include <sstream>
#include <string>
-std::vector< std::pair<double, double> > read_diagram_from_file( const char* filename )
-{
- std::ifstream in;
- in.open( filename );
- std::vector< std::pair<double, double> > result;
- if ( !in.is_open() )
- {
- std::cerr << "File : " << filename << " do not exist. The program will now terminate \n";
- throw "File do not exist \n";
- }
+std::vector< std::pair<double, double> > read_diagram_from_file(const char* filename) {
+ std::ifstream in;
+ in.open(filename);
+ std::vector< std::pair<double, double> > result;
+ if (!in.is_open()) {
+ std::cerr << "File : " << filename << " do not exist. The program will now terminate \n";
+ throw "File do not exist \n";
+ }
- std::string line;
- while (!in.eof())
- {
- getline(in,line);
- if ( line.length() != 0 )
- {
- std::stringstream lineSS;
- lineSS << line;
- double beginn, endd;
- lineSS >> beginn;
- lineSS >> endd;
- result.push_back( std::make_pair( beginn , endd ) );
- }
+ std::string line;
+ while (!in.eof()) {
+ getline(in, line);
+ if (line.length() != 0) {
+ std::stringstream lineSS;
+ lineSS << line;
+ double beginn, endd;
+ lineSS >> beginn;
+ lineSS >> endd;
+ result.push_back(std::make_pair(beginn, endd));
}
- in.close();
- return result;
-} //read_diagram_from_file
+ }
+ in.close();
+ return result;
+} // read_diagram_from_file
-int main( int argc , char** argv )
-{
- if ( argc < 3 )
- {
- std::cout << "To run this program please provide as an input two files with persistence diagrams. Each file " <<
- "should contain a birth-death pair per line. Third, optional parameter is an error bound on a bottleneck" <<
- " distance (set by default to zero). The program will now terminate \n";
+int main(int argc, char** argv) {
+ if (argc < 3) {
+ std::cout << "To run this program please provide as an input two files with persistence diagrams. Each file " <<
+ "should contain a birth-death pair per line. Third, optional parameter is an error bound on a bottleneck" <<
+ " distance (set by default to zero). The program will now terminate \n";
}
- std::vector< std::pair< double , double > > diag1 = read_diagram_from_file( argv[1] );
- std::vector< std::pair< double , double > > diag2 = read_diagram_from_file( argv[2] );
+ std::vector< std::pair< double, double > > diag1 = read_diagram_from_file(argv[1]);
+ std::vector< std::pair< double, double > > diag2 = read_diagram_from_file(argv[2]);
double tolerance = 0.;
- if ( argc == 4 )
- {
- tolerance = atof( argv[3] );
+ if (argc == 4) {
+ tolerance = atof(argv[3]);
}
double b = Gudhi::persistence_diagram::bottleneck_distance(diag1, diag2, tolerance);
std::cout << "The distance between the diagrams is : " << b << ". The tolerance is : " << tolerance << std::endl;