11#ifndef SKELETON_BLOCKER_SKELETON_BLOCKER_OFF_IO_H_
12#define SKELETON_BLOCKER_SKELETON_BLOCKER_OFF_IO_H_
14#include <gudhi/Off_reader.h>
22namespace skeleton_blocker {
27template<
typename Complex>
31 typedef typename Complex::Point Point;
33 const bool load_only_points_;
38 load_only_points_(load_only_points) { }
40 void init(
int dim,
int num_vertices,
int num_faces,
int num_edges) {
45 void point(
const std::vector<double>& point) {
46 complex_.
add_vertex(Point(point.begin(), point.end()));
49 void maximal_face(
const std::vector<int>& face) {
50 if (!load_only_points_) {
51 for (
size_t i = 0; i < face.size(); ++i)
52 for (
size_t j = i + 1; j < face.size(); ++j) {
64template<
typename Complex>
68 typedef typename Complex::Simplex Simplex;
69 typedef typename Complex::Point Point;
71 const bool load_only_points_;
72 std::vector<Point> points_;
73 std::vector<Simplex> maximal_faces_;
78 load_only_points_(load_only_points) { }
80 void init(
int dim,
int num_vertices,
int num_faces,
int num_edges) {
81 maximal_faces_.reserve(num_faces);
82 points_.reserve(num_vertices);
85 void point(
const std::vector<double>& point) {
86 points_.emplace_back(point.begin(), point.end());
89 void maximal_face(
const std::vector<int>& face) {
90 if (!load_only_points_) {
94 maximal_faces_.emplace_back(s);
99 complex_ = make_complex_from_top_faces<Complex>(maximal_faces_.begin(), maximal_faces_.end(),
100 points_.begin(), points_.end());
107template<
typename Complex>
116 bool read_only_points =
false,
bool is_flag =
false) : valid_(false) {
117 std::ifstream stream(name_file);
118 if (stream.is_open()) {
119 if (is_flag || read_only_points) {
122 valid_ = off_reader.
read(off_visitor);
126 valid_ = off_reader.
read(off_visitor);
142template<
typename Complex>
143class Skeleton_blocker_off_writer {
150 Skeleton_blocker_off_writer(
const std::string & name_file,
const Complex& save_complex) {
153 std::ofstream stream(name_file);
154 if (stream.is_open()) {
157 stream << save_complex.num_vertices() <<
" " << num_triangles <<
" 0 \n";
161 std::map<Vertex_handle, size_t> vertex_num;
162 size_t current_vertex = 0;
165 vertex_num[v] = current_vertex++;
166 const auto& pt(save_complex.
point(v));
175 stream << vertex_num[x] <<
" ";
180 std::cerr <<
"could not open file " << name_file << std::endl;
187namespace skbl = skeleton_blocker;
OFF file reader top class visitor.
Definition: Off_reader.h:29
bool read(OffVisitor &off_visitor)
Read an OFF file and calls the following methods :
Definition: Off_reader.h:51
Edge_handle add_edge_without_blockers(Vertex_handle a, Vertex_handle b)
Adds an edge between vertices a and b without blockers.
Definition: Skeleton_blocker_complex.h:566
Complex_triangle_around_vertex_range triangle_range(Vertex_handle v) const
Range over triangles around a vertex of the simplicial complex. Methods .begin() and ....
Definition: Skeleton_blocker_complex.h:1358
Complex_vertex_range vertex_range() const
Returns a Complex_vertex_range over all vertices of the complex.
Definition: Skeleton_blocker_complex.h:1282
Class that represents a geometric complex that can be simplified. The class allows access to points o...
Definition: Skeleton_blocker_geometric_complex.h:29
Vertex_handle add_vertex()
Add a vertex to the complex with a default constructed associated point.
Definition: Skeleton_blocker_geometric_complex.h:85
const Point & point(Vertex_handle v) const
Returns the Point associated to the vertex v.
Definition: Skeleton_blocker_geometric_complex.h:101
Off reader visitor that can be passed to Off_reader to read a Skeleton_blocker_complex.
Definition: Skeleton_blocker_off_io.h:28
Class that allows to load a Skeleton_blocker_complex from an off file.
Definition: Skeleton_blocker_off_io.h:108
Skeleton_blocker_off_reader(const std::string &name_file, Complex &read_complex, bool read_only_points=false, bool is_flag=false)
Definition: Skeleton_blocker_off_io.h:115
bool is_valid() const
Definition: Skeleton_blocker_off_io.h:134
Off reader visitor that can be passed to Off_reader to read a Skeleton_blocker_complex.
Definition: Skeleton_blocker_off_io.h:65
Gudhi namespace.
Definition: SimplicialComplexForAlpha.h:14
Definition: SkeletonBlockerDS.h:56
Handle type for the vertices of a cell complex.
Definition: VertexHandle.h:15