summaryrefslogtreecommitdiff
path: root/src/Nerve_GIC/example/simple_GIC.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'src/Nerve_GIC/example/simple_GIC.cpp')
-rw-r--r--src/Nerve_GIC/example/simple_GIC.cpp36
1 files changed, 7 insertions, 29 deletions
diff --git a/src/Nerve_GIC/example/simple_GIC.cpp b/src/Nerve_GIC/example/simple_GIC.cpp
index 33ebfcbb..d9297d04 100644
--- a/src/Nerve_GIC/example/simple_GIC.cpp
+++ b/src/Nerve_GIC/example/simple_GIC.cpp
@@ -1,18 +1,9 @@
-#include <gudhi/Rips_complex.h>
-// to construct Rips_complex from a OFF file of points
-#include <gudhi/Points_off_io.h>
-#include <gudhi/Simplex_tree.h>
-#include <gudhi/distance_functions.h>
#include <gudhi/GIC.h>
-#include <iostream>
-#include <string>
-#include <vector>
-
void usage(int nbArgs, char * const progName) {
std::cerr << "Error: Number of arguments (" << nbArgs << ") is not correct\n";
- std::cerr << "Usage: " << progName << " filename.off threshold dim_max [ouput_file.txt]\n";
- std::cerr << " i.e.: " << progName << " ../../data/points/alphacomplexdoc.off 60.0\n";
+ std::cerr << "Usage: " << progName << " filename.off threshold cover [ouput_file.txt]\n";
+ std::cerr << " i.e.: " << progName << " ../../data/points/test.off 1.5 test_cov \n";
exit(-1); // ----- >>
}
@@ -21,30 +12,20 @@ int main(int argc, char **argv) {
std::string off_file_name(argv[1]);
double threshold = atof(argv[2]);
- int dim_max = atoi(argv[3]);
+ std::string cover_file_name(argv[3]);
// Type definitions
- using Point = std::vector<float>;
- using Simplex_tree = Gudhi::Simplex_tree<>;
- using Filtration_value = Simplex_tree::Filtration_value;
- using Rips_complex = Gudhi::rips_complex::Rips_complex<Filtration_value>;
using Graph_t = boost::adjacency_list < boost::vecS, boost::vecS, boost::undirectedS,\
boost::property < vertex_filtration_t, Filtration_value >,\
boost::property < edge_filtration_t, Filtration_value > >;
- using Cover = std::map<int,int>;
// ----------------------------------------------------------------------------
// Init of a graph induced complex from an OFF file
// ----------------------------------------------------------------------------
- Gudhi::Points_off_reader<Point> off_reader(off_file_name);
- Rips_complex rips_complex_from_points(off_reader.get_point_cloud(), threshold, Euclidean_distance());
- Simplex_tree st;
- rips_complex_from_points.create_complex(st, 1);
- Cover C;
- Gudhi::graph_induced_complex::Graph_induced_complex GIC(st,C,dim_max);
- Simplex_tree stree;
- GIC.create_complex(stree);
+ Gudhi::graph_induced_complex::Graph_induced_complex GIC;
+ GIC.set_graph_simplex_tree(threshold, off_file_name); GIC.set_cover(cover_file_name); GIC.find_simplices();
+ Simplex_tree stree; GIC.create_complex(stree);
std::streambuf* streambufffer;
std::ofstream ouput_file_stream;
@@ -65,14 +46,11 @@ int main(int argc, char **argv) {
" - " << stree.num_simplices() << " simplices - " <<
stree.num_vertices() << " vertices." << std::endl;
- output_stream << "Iterator on graph induced complex simplices in the filtration order, with [filtration value]:" <<
- std::endl;
+ output_stream << "Iterator on graph induced complex simplices" << std::endl;
for (auto f_simplex : stree.filtration_simplex_range()) {
- output_stream << " ( ";
for (auto vertex : stree.simplex_vertex_range(f_simplex)) {
output_stream << vertex << " ";
}
- output_stream << ") -> " << "[" << stree.filtration(f_simplex) << "] ";
output_stream << std::endl;
}