From ab2e89665e16be2dcd36e55df8543d4ff5974675 Mon Sep 17 00:00:00 2001
From: mcarrier <mcarrier@636b058d-ea47-450e-bf9e-a15bfbe3eedb>
Date: Mon, 30 Jan 2017 16:33:22 +0000
Subject: git-svn-id:
 svn+ssh://scm.gforge.inria.fr/svnroot/gudhi/branches/Nerve_GIC@2032
 636b058d-ea47-450e-bf9e-a15bfbe3eedb

Former-commit-id: e8901f1be678dbc71b8e6e7b5666cc475b573b1a
---
 src/Nerve_GIC/example/simple_GIC.cpp | 36 +++++++-----------------------------
 1 file changed, 7 insertions(+), 29 deletions(-)

(limited to 'src/Nerve_GIC/example/simple_GIC.cpp')

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;
   }
 
-- 
cgit v1.2.3