#include <gudhi/Points_off_io.h>

// For CGAL points type in dimension d
// cf. http://doc.cgal.org/latest/Kernel_d/classCGAL_1_1Point__d.html
#include <CGAL/Epick_d.h>

#include <iostream>
#include <string>
#include <vector>

typedef CGAL::Epick_d< CGAL::Dynamic_dimension_tag > Kernel;
typedef typename Kernel::Point_d Point_d;

void usage(int argc, char * const progName) {
  std::cerr << "Error: Number of arguments (" << argc << ") is not correct" << std::endl;
  std::cerr << "Usage: " << progName << " inputFile.off" << std::endl;
  exit(-1);
}

int main(int argc, char **argv) {
  if (argc != 2) usage(argc, (argv[0] - 1));

  std::string offInputFile(argv[1]);
  // Read the OFF file (input file name given as parameter) and triangulate points
  Gudhi::Points_off_reader<Point_d> off_reader(offInputFile);
  // Check the read operation was correct
  if (!off_reader.is_valid()) {
    std::cerr << "Unable to read file " << offInputFile << std::endl;
    exit(-1);
  }

  // Retrieve the triangulation
  std::vector<Point_d> point_cloud = off_reader.get_point_cloud();

  int n = 0;
  for (auto point : point_cloud) {
    std::cout << "Point[" << n << "] = ";
    for (int i = 0; i < point.dimension(); i++)
      std::cout << point[i] << " ";
    std::cout << "\n";
    ++n;
  }
  return 0;
}