1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
|
#ifndef CGAL_SKIN_SURFACE_WRITER_H
#define CGAL_SKIN_SURFACE_WRITER_H
#include <fstream>
#include <CGAL/IO/Polyhedron_iostream.h>
#include <CGAL/subdivide_skin_surface_mesh_3.h>
#include <CGAL/Skin_surface_refinement_policy_3.h>
template <class SkinSurface, class Polyhedron>
/// Write polyhedron with normals:
void write_polyhedron_with_normals(SkinSurface &skin,
Polyhedron &p,
std::ostream &out)
{
typedef typename Polyhedron::Vertex_iterator Vertex_iterator;
typedef typename Polyhedron::Facet_iterator Facet_iterator;
typedef typename Polyhedron::Halfedge_around_facet_circulator HFC;
typedef typename Polyhedron::Vertex_handle Vertex_handle;
typedef typename Polyhedron::Traits::Vector_3 Vector;
CGAL::Skin_surface_refinement_policy_3<SkinSurface, Polyhedron> policy(skin);
// Write header
out << "NOFF " << p.size_of_vertices ()
<< " " << p.size_of_facets()
<< " " << p.size_of_halfedges()
<< std::endl;
// Write vertices
for (Vertex_iterator vit = p.vertices_begin();
vit != p.vertices_end(); vit ++) {
Vector n = policy.normal(vit);
n = n/sqrt(n*n);
out << vit->point() << " " << n << std::endl;
}
// Write faces
CGAL::Inverse_index<Vertex_handle> index(p.vertices_begin(),
p.vertices_end());
for(Facet_iterator fi = p.facets_begin();
fi != p.facets_end(); ++fi) {
HFC hc = fi->facet_begin();
HFC hc_end = hc;
std::size_t n = circulator_size( hc);
out << n;
do {
Vertex_handle vh = (*hc).vertex();
out << " " << index[vh];
} while (++hc != hc_end);
out << "\n";
}
}
#endif // CGAL_SKIN_SURFACE_WRITER_H
|
